Search

U12-type or minor introns are found in most multicellular eukaryotes and constitute ∼0.5% of all introns in species with a minor spliceosome. Although the biological significance for the evolutionary conservation of U12-type introns is debated, mutations disrupting U12 splicing cause developmental defects in both plants and animals. In human hematopoietic stem cells, U12 splicing defects disrupt proper differentiation of myeloid lineages and are associated with myelodysplastic syndrome, predisposing individuals to acute myeloid leukemia. Mutants in the maize ortholog of RNA binding motif protein 48 (RBM48) have aberrant U12-type intron splicing. Human RBM48 was recently purified biochemically as part of the minor spliceosome and shown to recognize the 5′ end of the U6atac snRNA. In this report, we use CRISPR/Cas9-mediated ablation of RBM48 in human K-562 cells to show the genetic function of RBM48. RNA-seq analysis comparing wild-type and mutant K-562 genotypes found that 48% of minor intron-containing genes have significant U12-type intron retention in RBM48 mutants. Comparing these results to maize rbm48 mutants defined a subset of minor intron-containing genes disrupted in both species. Mutations in the majority of these orthologous minor intron-containing genes have been reported to cause developmental defects in both plants and animals. Our results provide genetic evidence that the primary defect of human RBM48 mutants is aberrant U12-type intron splicing, while a comparison of human and maize RNA-seq data identifies candidate genes likely to mediate mutant phenotypes of U12-type splicing defects.

Plants grown in spaceflight experience novel environmental signals, including those associated with microgravity and ionizing radiation. Spaceflight triggers a response involving transcriptional re-programming and altered cell morphology, though many aspects of this response remain uncharacterized. We analyzed the spaceflight-induced transcriptome with a focus on genes that undergo alternative splicing to examine differential splicing associated with spaceflight—an unstudied characteristic of the molecular response to spaceflight exposure. RNA sequence data obtained during the APEX03 spaceflight experiment that was collected from two Arabidopsis thaliana ecotypes at two seedling stages grown onboard the International Space Station, or as ground controls at Kennedy Space Center, were re-examined to detect alternative splicing differences induced by spaceflight. Presence/absence variation analysis was used to identify putative expression-level differences in alternatively spliced isoforms between spaceflight and ground controls and was followed by analysis of significant differential alternative splicing. This study provides the first evidence of a role for alternative splicing in the molecular processes of physiological adaptation to the spaceflight environment.

The spiny mouse (Acomys species) has emerged as an exciting research organism due to its remarkable ability to undergo scarless regeneration of skin wounds and ear punches. Excitingly, Acomys species demonstrate scar-free healing in a wide-range of tissues beyond the skin. In this perspective article, we discuss published findings from a variety of tissues to highlight how this emerging research organism could shed light on numerous clinically relevant human diseases. We also discuss the challenges of working with this emerging research organism and suggest strategies for future Acomys-inspired research.

Jasmonate is an essential phytohormone regulating plant growth, development, and defense. Alternative splicing (AS) in jasmonate ZIM‐domain (JAZ) repressors is well‐characterized and plays an important role in jasmonate signaling regulation. However, it is unknown whether other genes in the jasmonate signaling pathway are regulated by AS. We explore the potential for AS regulation in three Arabidopsis genotypes (WT, jaz2, jaz7) in response to methyl jasmonate (MeJA) treatment with respect to: (a) differential AS, (b) differential miRNA targeted AS, and (c) AS isoforms with novel functions. AS events identified from transcriptomic data were validated with proteomic data. Protein interaction networks identified two genes, SKIP and ALY4 whose products have both DNA‐ and RNA‐binding affinities, as potential key regulators mediating jasmonate signaling and AS regulation. We observed cases where AS alone, or AS and transcriptional regulation together, can influence gene expression in response to MeJA. Twenty‐one genes contain predicted miRNA target sites subjected to AS, which implies that AS is coupled to miRNA regulation. We identified 30 cases where alternatively spliced isoforms may have novel functions. For example, AS of bHLH160 generates an isoform without a basic domain, which may convert it from an activator to a repressor. Our study identified potential key regulators in AS regulation of jasmonate signaling pathway. These findings highlight the importance of AS regulation in the jasmonate signaling pathway, both alone and in collaboration with other regulators. Significance statement By exploring alternative splicing, we demonstrate its regulation in the jasmonate signaling pathway alone or in collaboration with other posttranscriptional regulations such as nonsense and microRNA‐mediated decay. A signal transduction network model for alternative splicing in jasmonate signaling pathway was generated, contributing to our understanding for this important, prevalent, but relatively unexplored regulatory mechanism in plants.

U12-type or minor introns are found in most multicellular eukaryotes and constitute ∼0.5% of all introns in species with a minor spliceosome. Although the biological significance for evolutionary conservation of U12-type introns is debated, mutations disrupting U12 splicing cause developmental defects in both plants and animals. In human hematopoietic stem cells, U12 splicing defects disrupt proper differentiation of myeloid lineages and are associated with myelodysplastic syndrome (MDS), predisposing individuals to acute myeloid leukemia. Mutants in the maize ortholog of RNA Binding Motif Protein48 (RBM48) have aberrant U12-type intron splicing. Human RBM48 was recently purified biochemically as part of the minor spliceosome and shown to recognize the 5’ end of the U6atac snRNA. In this report, we use CRISPR/Cas9-mediated ablation of RBM48 in human K-562 cells to show the genetic function of RBM48. RNA-seq analysis comparing wild-type and mutant K-562 genotypes found that 48% of minor intron containing genes (MIGs) have significant U12-type intron retention in RBM48 mutants. Comparing these results to maize rbm48 mutants defined a subset of MIGs disrupted in both species. Mutations in the majority of these orthologous MIGs have been reported to cause developmental defects in both plants and animals. Our results provide genetic evidence that the primary defect of human RBM48 mutants is aberrant U12-type intron splicing, while a comparison of human and maize RNA-seq data identifies candidate genes likely to mediate mutant phenotypes of U12-type splicing defects.

Polyploidy has played a pivotal and recurring role in angiosperm evolution. Allotetraploids arise from hybridization between species and possess duplicated gene copies (homeologs) that serve redundant roles immediately after polyploidization. Although polyploidization is a major contributor to plant evolution, it remains poorly understood. We describe an analytical approach for assessing homeolog-specific expression that begins with de novo assembly of parental transcriptomes and effectively (i) reduces redundancy in de novo assemblies, (ii) identifies putative orthologs, (iii) isolates common regions between orthologs, and (iv) assesses homeolog-specific expression using a robust Bayesian Poisson-Gamma model to account for sequence bias when mapping polyploid reads back to parental references. Using this novel methodology, we examine differential homeolog contributions to the transcriptome in the recently formed allopolyploids Tragopogon mirus and T. miscellus (Compositae). Notably, we assess a larger Tragopogon gene set than previous studies of this system. Using carefully identified orthologous regions and filtering biased orthologs, we find in both allopolyploids largely balanced expression with no strong parental bias. These new methods can be used to examine homeolog expression in any tetrapolyploid system without requiring a reference genome.

International audience; Compared to traditionally used plastid or ribosomal markers, nuclear gene markers provide many advantages for molecular systematics of plants, and high-throughput sequencing technologies are making large nuclear datasets available at an unprecedented rate. We used targeted amplicon sequencing of 44 nuclear genes to construct a time-calibrated phylogeny of genus Metrosideros (Myrtaceae), evaluate recent systematic revisions, and assess whether phylogenetic signal within the Hawaiian Archipelago is correlated with island biogeography or morphological diversification. We generated a final dataset of 40 nuclear genes for 187 specimens, used multiple search heuristics and species-tree analysis to estimate a phylogeny, and incorporated new fossils for the genus to estimate divergence times across the dataset. All analyses supported the monophyly of Metrosideros, including Carpolepis and Tepualia. Hawaiian Metrosideros were monophyletic and dated to 3.1 MYA using new fossils for the genus, which is intermediate to previous estimates based on nuclear ribosomal/chloroplast loci and calibrated with island ages. Within the Hawaiian Metrosideros clade, we observed short branch lengths and unresolved relationships, and phylogenetic patterns were not concordant with biogeographic hypotheses of island progression, or the delineation of taxa or morphotypes. Average nucleotide diversity was relatively consistent across the Hawaiian Islands with the exception of slightly lower diversity on Kaua'i. These results provide a data-rich estimate of the timing of a single introduction of Metrosideros to Hawai'i and highlight the need for molecular markers with higher evolutionary rates for resolution of relationships within this recent radiation.

Despite its economic importance as a bioenergy crop and key role in riparian ecosystems, little is known about genetic diversity and adaptation of the eastern cottonwood (Populus deltoides). Here, we report the first population genomics study for this species, conducted on a sample of 425 unrelated individuals collected in 13 states of the southeastern United States. The trees were genotyped by targeted resequencing of 18,153 genes and 23,835 intergenic regions, followed by the identification of single nucleotide polymorphisms (SNPs). This natural P. deltoides population showed low levels of subpopulation differentiation (F ST = 0.022–0.106), high genetic diversity (θW = 0.00100, π = 0.00170), a large effective population size (N e ≈ 32,900), and low to moderate levels of linkage disequilibrium. Additionally, genomewide scans for selection (Tajima's D), subpopulation differentiation (XTX), and environmental association analyses with eleven climate variables carried out with two different methods (LFMM and BAYENV2) identified genes putatively involved in local adaptation. Interestingly, many of these genes were also identified as adaptation candidates in another poplar species, Populus trichocarpa, indicating possible convergent evolution. This study constitutes the first assessment of genetic diversity and local adaptation in P. deltoides throughout the southern part of its range, information we expect to be of use to guide management and breeding strategies for this species in future, especially in the face of climate change.

Alternative splicing (AS) is a major source of transcript and proteome diversity, but examining AS in species without well-annotated reference genomes remains difficult. Research on both human and mouse has demonstrated the advantages of using Iso-Seq™ data for isoform-level transcriptome analysis, including the study of AS and gene fusion. We applied Iso-Seq™ to investigate AS in Amborella trichopoda, a phylogenetically pivotal species that is sister to all other living angiosperms. Our data show that, compared with RNA-Seq data, the Iso-Seq™ platform provides better recovery on large transcripts, new gene locus identification and gene model correction. Reference-based AS detection with Iso-Seq™ data identifies AS within a higher fraction of multi-exonic genes than observed for published RNA-Seq analysis (45.8% vs. 37.5%). These data demonstrate that the Iso-Seq™ approach is useful for detecting AS events. Using the Iso-Seq-defined transcript collection in Amborella as a reference, we further describe a pipeline for detection of AS isoforms from PacBio Iso-Seq™ without using a reference sequence (de novo). Results using this pipeline show a 66%-76% overall success rate in identifying AS events. This de novoAS detection pipeline provides a method to accurately characterize and identify bona fide alternatively spliced transcripts in any nonmodel system that lacks a reference genome sequence. Hence, our pipeline has huge potential applications and benefits to the broader biology community.

Significance The avocado is a nutritious, economically important fruit species that occupies an unresolved position near the earliest evolutionary branchings of flowering plants. Our nuclear genome sequences of Mexican and Hass variety avocados inform ancient evolutionary relationships and genome doublings and the admixed nature of Hass and provide a look at how pathogen interactions have shaped the avocado’s more recent genomic evolutionary history., The avocado, Persea americana, is a fruit crop of immense importance to Mexican agriculture with an increasing demand worldwide. Avocado lies in the anciently diverged magnoliid clade of angiosperms, which has a controversial phylogenetic position relative to eudicots and monocots. We sequenced the nuclear genomes of the Mexican avocado race, P. americana var. drymifolia, and the most commercially popular hybrid cultivar, Hass, and anchored the latter to chromosomes using a genetic map. Resequencing of Guatemalan and West Indian varieties revealed that ∼39% of the Hass genome represents Guatemalan source regions introgressed into a Mexican race background. Some introgressed blocks are extremely large, consistent with the recent origin of the cultivar. The avocado lineage experienced 2 lineage-specific polyploidy events during its evolutionary history. Although gene-tree/species-tree phylogenomic results are inconclusive, syntenic ortholog distances to other species place avocado as sister to the enormous monocot and eudicot lineages combined. Duplicate genes descending from polyploidy augmented the transcription factor diversity of avocado, while tandem duplicates enhanced the secondary metabolism of the species. Phenylpropanoid biosynthesis, known to be elicited by Colletotrichum (anthracnose) pathogen infection in avocado, is one enriched function among tandems. Furthermore, transcriptome data show that tandem duplicates are significantly up- and down-regulated in response to anthracnose infection, whereas polyploid duplicates are not, supporting the general view that collections of tandem duplicates contribute evolutionarily recent “tuning knobs” in the genome adaptive landscapes of given species.

Background Natural rubber, an indispensable commodity used in approximately 40,000 products, is fundamental to the tire industry. The rubber tree species Hevea brasiliensis (Willd. ex Adr. de Juss.) Muell-Arg., which is native the Amazon rainforest, is the major producer of latex worldwide. Rubber tree breeding is time consuming, expensive and requires large field areas. Thus, genetic studies could optimize field evaluations, thereby reducing the time and area required for these experiments. In this work, transcriptome sequencing was used to identify a full set of transcripts and to evaluate the gene expression involved in the different cold-response strategies of the RRIM600 (cold-resistant) and GT1 (cold-tolerant) genotypes. Results We built a comprehensive transcriptome using multiple database sources, which resulted in 104,738 transcripts clustered in 49,304 genes. The RNA-seq data from the leaf tissues sampled at four different times for each genotype were used to perform a gene-level expression analysis. Differentially expressed genes (DEGs) were identified through pairwise comparisons between the two genotypes for each time series of cold treatments. DEG annotation revealed that RRIM600 and GT1 exhibit different chilling tolerance strategies. To cope with cold stress, the RRIM600 clone upregulates genes promoting stomata closure, photosynthesis inhibition and a more efficient reactive oxygen species (ROS) scavenging system. The transcriptome was also searched for putative molecular markers (single nucleotide polymorphisms (SNPs) and microsatellites) in each genotype. and a total of 27,111 microsatellites and 202,949 (GT1) and 156,395 (RRIM600) SNPs were identified in GT1 and RRIM600. Furthermore, a search for alternative splicing (AS) events identified a total of 20,279 events. Conclusions The elucidation of genes involved in different chilling tolerance strategies associated with molecular markers and information regarding AS events provides a powerful tool for further genetic and genomic analyses of rubber tree breeding. Electronic supplementary material The online version of this article (10.1186/s12864-019-5852-5) contains supplementary material, which is available to authorized users.

The African spiny mouse, Acomys spp., is capable of scar-free dermal wound healing. Here, we have performed a comprehensive analysis of gene expression throughout wound healing following full-thickness excisional dermal wounds in both Acomys cahirinus and Mus musculus. Additionally, we provide an annotated, de novo transcriptome assembly of A. cahirinus skin and skin wounds. Using a novel computational comparative RNA-Seq approach along with pathway and co-expression analyses, we identify enrichment of regeneration associated genes as well as upregulation of genes directly related to muscle development or function. Our RT-qPCR data reveals induction of the myogenic regulatory factors, as well as upregulation of embryonic myosin, starting between days 14 and 18 post-wounding in A. cahirinus. In contrast, the myogenic regulatory factors remain downregulated, embryonic myosin is only modestly upregulated, and no new muscle fibers of the panniculus carnosus are generated in M. musculus wounds. Additionally, we show that Col6a1, a key component of the satellite cell niche, is upregulated in A. cahirinus compared to M. musculus. Our data also demonstrate that the macrophage profile and inflammatory response is different between species, with A. cahirinus expressing significantly higher levels of Il10. We also demonstrate differential expression of the upstream regulators Wnt7a, Wnt2 and Wnt6 during wound healing. Our analyses demonstrate that A. cahirinus is capable of de novo skeletal muscle regeneration of the panniculus carnosus following removal of the extracellular matrix. We believe this study represents the first detailed analysis of de novo skeletal muscle regeneration observed in an adult mammal.

Natural rubber is an indispensable commodity used in approximately 40,000 products and is fundamental to the tire industry. Among the species that produce latex, the rubber tree [Hevea brasiliensis (Willd. ex Adr. de Juss.) Muell-Arg.], a species native to the Amazon rainforest, is the major producer of latex used worldwide. The Amazon Basin presents optimal conditions for rubber tree growth, but the occurrence of South American leaf blight, which is caused by the fungus Microcyclus ulei (P. Henn) v. Arx, limits rubber tree production. Currently, rubber tree plantations are located in scape regions that exhibit suboptimal conditions such as high winds and cold temperatures. Rubber tree breeding programs aim to identify clones that are adapted to these stress conditions. However, rubber tree breeding is time-consuming, taking more than 20 years to develop a new variety. It is also expensive and requires large field areas. Thus, genetic studies could optimize field evaluations, thereby reducing the time and area required for these experiments. Transcriptome sequencing using next-generation sequencing (RNA-seq) is a powerful tool to identify a full set of transcripts and for evaluating gene expression in model and non-model species. In this study, we constructed a comprehensive transcriptome to evaluate the cold response strategies of the RRIM600 (cold-resistant) and GT1 (cold-tolerant) genotypes. Furthermore, we identified putative microsatellite (SSR) and single-nucleotide polymorphism (SNP) markers. Alternative splicing, which is an important mechanism for plant adaptation under abiotic stress, was further identified, providing an important database for further studies of cold tolerance.

We report here the genome sequences of three newly isolated phages that infect Mycobacterium smegmatis mc 2 155. Phages Findley, Hurricane, and TBond007 were discovered in geographically distinct locations and are related to cluster K mycobacteriophages, with Findley being similar to subcluster K2 phages and Hurricane and TBond007 being similar to subcluster K3 phages.

The maize W22 inbred has served as a platform for maize genetics since the mid twentieth century. To streamline maize genome analyses, we have sequenced and de novo assembled a W22 reference genome using short-read sequencing technologies. We show that significant structural heterogeneity exists in comparison to the B73 reference genome at multiple scales, from transposon composition and copy number variation to single-nucleotide polymorphisms. The generation of this reference genome enables accurate placement of thousands of Mutator (Mu) and Dissociation (Ds) transposable element insertions for reverse and forward genetics studies. Annotation of the genome has been achieved using RNA-seq analysis, differential nuclease sensitivity profiling and bisulfite sequencing to map open reading frames, open chromatin sites and DNA methylation profiles, respectively. Collectively, the resources developed here integrate W22 as a community reference genome for functional genomics and provide a foundation for the maize pan-genome.

One difficulty when identifying alternative splicing (AS) events in plants is distinguishing functional AS from splicing noise. One way to add confidence to the validity of a splice isoform is to observe that it is conserved across evolutionarily related species. We use a high throughput method to identify junction-based conserved AS events from RNA-Seq data across nine plant species, including five grass monocots (maize, sorghum, rice, Brachpodium, and foxtail millet), plus two nongrass monocots (banana and African oil palm), the eudicot Arabidopsis, and the basal angiosperm Amborella. In total, 9804 AS events were found to be conserved between two or more species studied. In grasses containing large regions of conserved synteny, the frequency of conserved AS events is twice that observed for genes outside of conserved synteny blocks. In plant-specific RS and RS2Z subfamilies of the serine/arginine (SR) splice-factor proteins, we observe both conservation and divergence of AS events after the whole genome duplication in maize. In addition, plant-specific RS and RS2Z splice-factor subfamilies are highly connected with R2R3-MYB in STRING functional protein association networks built using genes exhibiting conserved AS. Furthermore, we discovered that functional protein association networks constructed around genes harboring conserved AS events are enriched for phosphatases, kinases, and ubiquitylation genes, which suggests that AS may participate in regulating signaling pathways. These data lay the foundation for identifying and studying conserved AS events in the monocots, particularly across grass species, and this conserved AS resource identifies an additional layer between genotype to phenotype that may impact future crop improvement efforts.

RNA splicing of U12-type introns functions in human cell differentiation, but it is not known whether this class of introns has a similar role in plants. The maize ROUGH ENDOSPERM3 (RGH3) protein is orthologous to the human splicing factor, ZRSR2. ZRSR2 mutations are associated with myelodysplastic syndrome (MDS) and cause U12 splicing defects. Maize rgh3 mutants have aberrant endosperm cell differentiation and proliferation. We found that most U12-type introns are retained or misspliced in rgh3. Genes affected in rgh3 and ZRSR2 mutants identify cell cycle and protein glycosylation as common pathways disrupted. Transcripts with retained U12-type introns can be found in polysomes, suggesting that splicing efficiency can alter protein isoforms. The rgh3 mutant protein disrupts colocalization with a known ZRSR2-interacting protein, U2AF2. These results indicate conserved function for RGH3/ZRSR2 in U12 splicing and a deeply conserved role for the minor spliceosome to promote cell differentiation from stem cells to terminal fates.

Next-generation sequencing (NGS) has revolutionized molecular systematics as well as population and conservation genetics. It is now possible to obtain enormous amounts of gene sequence data from any species in a short time at low cost. More technological advances are on the horizon, ensuring that this trend will continue throughout the coming decade. These rapid advances provide unprecedented opportunities in systematics; they also pose new challenges, requiring that the next generation of systematists be well-versed in new skill sets (e.g., bioinformatics). As examples of the potential of NGS, it is now possible to develop genetic resources for any plant system that poses intriguing evolutionary questions. During the next decade many new "evolutionary model systems" will become available as systematists rapidly develop the necessary genetic/genomic frameworks for many previously unstudied plants. Phylogenetic reconstruction will be conducted at an unprecedented pace at both deep and fine scales with datasets of numerous taxa and genes—this includes rapid progress on assembling a more com- prehensive Tree of Life for green plants. For example, complete plastid genome sequencing is now routinely facilitating analyses of hundreds of taxa at deep levels, as well as enabling complete plastid genome phylogeographic analyses at the population level. Gene capture methods hold enormous promise for the rapid and inexpensive analyses of complete plastid genomes, as well as studies of hundreds of selected (targeted) nuclear loci. NGS has also had a big impact on population genetics, initially by dramatically simplifying microsatellite marker development, but more recently by opening new possibilities through vari- ous genotyping-by-sequencing (GBS) approaches that have great potential to expand on the types of questions that can be addressed at the population level. Transcriptome sequencing has enabled the construction of large datasets of nuclear genes while also providing a wealth of plastid and mitochondrial genes. NGS has also facilitated probe development for studies of chromosomes using FISH (fluorescence in situ hybridization). NGS is also making the rapid sequencing of complete nuclear genomes routine, thus transforming our field and opening up new avenues of systematic endeavor in comparative genomics. However, even as sequencing costs drop and technological advances make complete nuclear genome sequencing more com - monplace, genome assembly will remain a major challenge.

Strawberry (Fragaria L.) genotypes bear remarkable phenotypic similarity, even across ploidy levels. Additionally, breeding programs seek to introgress alleles from wild germplasm, so objective molecular description of genetic variation has great value. In this report, a high-throughput, robust protocol for generating highly-informative simple sequence repeat (SSR) patterns is presented to address these issues. The methods are comparable to SSR use in DNA typing in humans and are based on identification of high-number repeats composed of tetra- through nona-nucleotide repeat units found in the Fragaria vesca genome sequence. Individual SSR-containing regions were examined for variability over a range of 219 strawberry genotypes. A single-fluorophore secondary labeling strategy was devised that allows simultaneous amplification of eight SSR regions in a single PCR reaction. The approach yields reproducible, highly-variable, complex patterns (Shannon-Weaver Index 7.09–13.88). The technique may be applied to detect closely-related individuals across ploidy levels, including full sibling progeny in an inter-related octoploid pedigree. Genetic diversity among various cultivars and progenitor wild species in the United States Department of Agriculture-Agricultural Research Service Fragaria Supercore collection was also evaluated. The results build on known relationships, and also raise questions about accepted relationships between several genotypes.

Summary New hybrid species might be expected to show patterns of gene expression intermediate to those shown by parental species [1, 2]. "Transcriptomic shock" may also occur, in which gene expression is disrupted; this may be further modified by whole genome duplication (causing allopolyploidy) [3–16]. "Shock" can include instantaneous partitioning of gene expression between parental copies of genes among tissues [16–19]. These effects have not previously been studied at a population level in a natural allopolyploid plant species. Here, we survey tissue-specific expression of 144 duplicated gene pairs derived from different parental species (homeologs) in two natural populations of 40-generation-old allotetraploid Tragopogon miscellus (Asteraceae) plants. We compare these results with patterns of allelic expression in both in vitro "hybrids" and hand-crossed F 1 hybrids between the parental diploids T. dubius and T. pratensis , and with patterns of homeolog expression in synthetic (S 1 ) allotetraploids. Partitioning of expression was frequent in natural allopolyploids, but F 1 hybrids and S 1 allopolyploids showed less partitioning of expression than the natural allopolyploids and the in vitro "hybrids" of diploid parents. Our results suggest that regulation of gene expression is relaxed in a concerted manner upon hybridization, and new patterns of partitioned expression subsequently emerge over the generations following allopolyploidization.

Strawberry (Fragaria spp.) is a valuable fruit crop as well as an outstanding system for studying functional genomics in plants. The goal of this study was to substantially increase and analyze the available expressed sequence information in the genus by examining the transcriptome of the cultivated strawberry (Fragaria × ananassa Duchesne). To maximize transcript diversity and discovery, plants representing an octoploid strawberry cultivar were subjected to a broad range of treatments. Plant materials were pooled by tissue type. cDNA pools were sequenced by the Roche-454 GS-FLX system and assembled into over 32,000 contigs. Predictions of cellular localization and function were made by associating assembled contigs to annotated homologs, and the tissue pool tags provided a means to assess the overall expression pattern for any given transcript. Contigs comprised of reads originating from only one organ type and those present equally in all plant organs were both identifi ed. Bacterial and fungal sequences found in the strawberry samples provide a metagenomic survey of the microbial community of a greenhouse strawberry plant. This study utilized an innovative assembly strategy on pooled tissues, thus providing a foundation for developing tissue-specifi c tools, an opportunity to identify alleles for marker-assisted selection, a reference of strawberry gene annotations, and a basis for comparative transcriptomics between cultivated strawberry, its diploid ancestors, and the wider Rosaceae family. C ULTIVATED STRAWBERRY is a valuable crop, prized for its fl avorful and nutritious fruits. Strawberry is an important crop for many regional economies worldwide and off ers high levels of healthful compounds. Species within the genus Fragaria occur naturally throughout the northern hemisphere and represent many levels of ploidy, including wild diploids and the wild and cultivated octoploids (reviewed in Folta and Davis, 2006; Hummer and Hancock, 2009). Recent attention has been brought to the diploid strawberry species Fragaria vesca Coville because of its utility as a highly tractable experimental system within the economically important Rosaceae family (Folta and Davis, 2006). Strawberry is emerging as a powerful functional genomics system. Fragaria vesca’s haploid genome size of ~200 Mb is among the smallest of the rosaceous crop plants (Shulaev et al., 2008), and its complete genome sequence will be available in 2010, providing an invaluable resource with reference to its octoploid, cultivated relatives. Most importantly, several genotypes of both diploid and octoploid strawberry are readily transformable, making it an outstanding system for functional studies. Protocols exist for diploid and octoploid lines that allow effi cient regeneration of transgenic materials in a period of weeks to months (Folta and Dhingra, 2006; Mezzetti and

Uncovering conserved alternative splicing (AS) events can identify AS events that perform important functions. This is especially useful for identifying premature stop codon containing (PTC) AS isoforms that may regulate protein expression by being targets for nonsense mediated decay. This report discusses the identification of a PTC containing splice isoform of the TFIIIA gene that is highly conserved in land plants. TFIIIA is essential for RNA Polymerase III-based transcription of 5S rRNA in eukaryotes. Two independent groups have determined that the PTC containing alternative exon is ultraconserved and is coupled with nonsense-mediated mRNA decay. The alternative exon appears to have been derived by the exonization of 5S ribosomal RNA (5S rRNA) within the gene of its own transcription regulator, TFIIIA. This provides the first evidence of ancient exaptation of 5S rRNA in plants, suggesting a novel gene regulation model mediated by the AS of an anciently exonized non-coding element.

Rice transcription factor RF2a binds to the BoxII cis element of the promoter of rice tungro bacilliform virus and activates promoter expression. The acidic acid-rich domain of RF2a is a transcription activator and has been partially characterized (Dai et al., 2003). The RF2a acidic domain (A; amino acids 49-116) was fused with the synthetic zinc finger ZF-TF 2C7 and was co-introduced with a reporter gene into transgenic Arabidopsis plants. Expression of the reporter gene was increased up to seven times by the effector. In transient assays in tobacco BY-2 protoplasts, we identified a subdomain comprising amino acids 56-84 (A5) that was equally as effective as an activator as the entire acidic domain. A chemically inducible system was used to show determined that A and A5 domains are equally as effective in transcription activation as the well-characterized VP16 activation domain. Bioinformatics analyses revealed that the A5 domain is present only in b-ZIP transcription factors. In dicots, the A domain contains an insertion of four amino acids that is not present in monocot proteins. The A5 domain, and similar domains in other b-ZIP transcription factors, is predicted to form an anti-parallel beta sheet structure.

Sequence capture technologies, pioneered in mammalian genomes, enable the resequencing of targeted genomic regions. Most capture protocols require blocking DNA, the production of which in large quantities can prove challenging. A blocker-free, two-stage capture protocol was developed using NimbleGen arrays. The first capture depletes the library of repetitive sequences, while the second enriches for target loci. This strategy was used to resequence non-repetitive portions of an approximately 2.2 Mb chromosomal interval and a set of 43 genes dispersed in the 2.3 Gb maize genome. This approach achieved approximately 1800-3000-fold enrichment and 80-98% coverage of targeted bases. More than 2500 SNPs were identified in target genes. Low rates of false-positive SNP predictions were obtained, even in the presence of captured paralogous sequences. Importantly, it was possible to recover novel sequences from non-reference alleles. The ability to design novel repeat-subtraction and target capture arrays makes this technology accessible in any species.

Advances in next-generation sequencing technology have facilitated the discovery of single nucleotide polymorphisms (SNPs). Sequenom-based SNP-typing assays were developed for 1359 maize SNPs identified via comparative next-generation transcriptomic sequencing. Approximately 75% of these SNPs were successfully converted into genetic markers that can be scored reliably and used to generate a SNP-based genetic map by genotyping recombinant inbred lines from the intermated B73 × Mo17 population. The quantitative nature of Sequenom-based SNP assays led to the development of a time- and cost-efficient strategy to genetically map mutants via quantitative bulked segregant analysis. This strategy was used to rapidly map the loci associated with several dozen recessive mutants. Because a mutant can be mapped using as few as eight multiplexed sets of SNP assays on a bulk of as few as 20 mutant F2 individuals, this strategy is expected to be widely adopted for mapping in many species.

Recent genome-wide estimates of the alternative splicing (AS) rate of plant genes indicate that AS is more prevalent than originally expected (Reddy 2007), although substantially fewer plant genes undergo AS than has been observed for human genes. Approximately 20% of multiexon genes in several plant species, including Arabidopsis, rice, moss, and maize, are predicted to be alternatively spliced based on EST/cDNA evidence (Campbell et al. 2006; Wang and Brendel 2006; Barbazuk et al. 2008). The functional significance of almost all known and predicted plant AS events has yet to be determined; however, the results of a few functional analyses conducted demonstrate roles for AS in important plant processes such as photosynthesis, defense response, flowering, and cereal grain quality (Reddy 2007), as well as some metabolic pathways (Gorlach et al. 1995) and mRNA processing (Golovkin and Reddy 1996; Kalyna et al. 2006). Identification of conserved AS events among multiple evolutionarily distant species will prioritize AS events for functional characterization and help identify relevant cis- and trans-regulatory factors (Xing and Lee 2006; Reddy 2007; Barbazuk et al. 2008). AS events that occur within the translated regions of mRNAs may contribute to the diversity of the proteome (Stamm et al. 2005). Some splice isoforms containing a premature stop codon (PTC) are often not translated. Instead, they are targeted for nonsense-mediated mRNA decay (NMD), which is an RNA surveillance system that recognizes PTC-containing mRNAs and targets them for degradation (Maquat 2004). Lewis et al. (2003) proposed that coupled AS and NMD, also termed RUST (regulated unproductive splicing and translation), may function to regulate protein expression by generating NMD-targeted isoforms. Recently, ultraconserved elements in mammals associated with cassette exons in some splicing activator SR (serine-arginine rich) proteins have been discovered, and transcript isoforms that include these conserved sequences contain PTCs and are subjected to NMD. The conserved nature of these elements in SR proteins suggests that their unproductive splicing is functionally important for the autoregulation or homeostatic control of these splicing regulators (Lareau et al. 2007; Ni et al. 2007). Conserved AS–NMD events have also been reported within SR proteins in plants (Kalyna et al. 2006). A genome-wide comparative analysis for AS events conserved between Arabidopsis and rice (Oryza sativa sp. Nipponbare), which diverged 140–150 Mya (Chaw et al. 2004), identified a conserved cassette exon splicing event associated with highly conserved sequences within the TFIIIA gene (transcription factor for polymerase III A). In eukaryotes, 5S ribosomal RNA (5S rRNA) is transcribed by RNA polymerase III (Pol III) and is required for the assembly of the large ribosomal subunit (Szymanski et al. 2003). TFIIIA is essential for the Pol III-based transcription of 5S rRNA. TFIIIA was the first zinc finger (ZF) protein characterized and the first protein found to bind both 5S rDNA and 5S rRNA (Engelke et al. 1980; Pelham and Brown 1980). The dual functions of TFIIIA suggest a negative feedback autoregulation model for 5S rRNA transcription (Cassiday and Maher III 2002). A few cis- and trans-factors have been implicated in the developmental regulation of the TFIIIA in Xenopus (Pfaff and Taylor 1992, 1998; Griffin et al. 2003; Penberthy et al. 2003); however, the transcriptional regulation of TFIIIA is not well understood. Post-transcriptional regulation of TFIIIA has never been reported in any species. The current analysis provides evidence of post-transcriptional autoregulation of TFIIIA in plants mediated by coupling a cassette exon AS event in TFIIIA pre-mRNA with NMD. This cassette exon is highly conserved in land plants and exhibits striking sequence and predicted secondary structure similarity to 5S rRNA. Similar to exonization of Alu elements in mammalian genomes (Sorek et al. 2002), the cassette exon may have derived from a 5S rDNA integration into the TFIIIA DNA and acquired a role in regulating TFIIIA homeostasis.

Alternative splicing (AS) creates multiple mRNA transcripts from a single gene. While AS is known to contribute to gene regulation and proteome diversity in animals, the study of its importance in plants is in its early stages. However, recently available plant genome and transcript sequence data sets are enabling a global analysis of AS in many plant species. Results of genome analysis have revealed differences between animals and plants in the frequency of alternative splicing. The proportion of plant genes that have one or more alternative transcript isoforms is ∼20%, indicating that AS in plants is not rare, although this rate is approximately one-third of that observed in human. The majority of plant AS events have not been functionally characterized, but evidence suggests that AS participates in important plant functions, including stress response, and may impact domestication and trait selection. The increasing availability of plant genome sequence data will enable larger comparative analyses that will identify functionally important plant AS events based on their evolutionary conservation, determine the influence of genome duplication on the evolution of AS, and discover plant-specific cis-elements that regulate AS. This review summarizes recent analyses of AS in plants, discusses the importance of further analysis, and suggests directions for future efforts.

The contribution of epigenetic alterations to natural variation for gene transcription levels remains unclear. In this study, we investigated the functional targets of the maize chromomethylase ZMET2 in multiple inbred lines to determine whether epigenetic changes conditioned by this chromomethylase are conserved or variable within the species. Gene expression microarrays were hybridized with RNA samples from the inbred lines B73 and Mo17 and from near-isogenic derivatives containing the loss-of-function allele zmet2-m1. A set of 126 genes that displayed statistically significant differential expression in zmet2 mutants relative to wild-type plants in at least one of the two genetic backgrounds was identified. Analysis of the transcript levels in both wild-type and mutant individuals revealed that only 10% of these genes were affected in zmet2 mutants in both B73 and Mo17 genetic backgrounds. Over 80% of the genes with expression patterns affected by zmet2 mutations display variation for gene expression between wild-type B73 and Mo17 plants. Further analysis was performed for 7 genes that were transcriptionally silent in wild-type B73, but expressed in B73 zmet2-m1, wild-type Mo17, and Mo17 zmet2-m1 lines. Mapping experiments confirmed that the expression differences in wild-type B73 relative to Mo17 inbreds for these genes were caused by cis-acting regulatory variation. Methylation-sensitive PCR and bisulfite sequencing demonstrated that for 5 of these genes the CpNpG methylation in the wild-type B73 genetic background was substantially decreased in the B73 zmet2-m1 mutant and in wild-type Mo17. A survey of eight maize inbreds reveals that each of these 5 genes exhibit transcriptionally silent and methylated states in some inbred lines and unmethylated, expressed states in other inbreds, providing evidence for natural variation in epigenetic states for some maize genes.

A massively parallel pyro-sequencing technology commercialized by 454 Life Sciences Corporation was used to sequence the transcriptomes of shoot apical meristems isolated from two inbred lines of maize using laser capture microdissection (LCM). A computational pipeline that uses the POLYBAYES polymorphism detection system was adapted for 454 ESTs and used to detect SNPs (single nucleotide polymorphisms) between the two inbred lines. Putative SNPs were computationally identified using 260 000 and 280 000 454 ESTs from the B73 and Mo17 inbred lines, respectively. Over 36 000 putative SNPs were detected within 9980 unique B73 genomic anchor sequences (MAGIs). Stringent post-processing reduced this number to > 7000 putative SNPs. Over 85% (94/110) of a sample of these putative SNPs were successfully validated by Sanger sequencing. Based on this validation rate, this pilot experiment conservatively identified > 4900 valid SNPs within > 2400 maize genes. These results demonstrate that 454-based transcriptome sequencing is an excellent method for the high-throughput acquisition of gene-associated SNPs.

Alternative splicing (AS) plays important roles in many plant functions, but its conservation across the plant kingdom is not known. We describe a methodology to identify AS events and identify conserved AS events across large phylogenetic distances using RNA-Seq datasets. We applied this methodology to transcriptome data from nine angiosperms including Amborella, the single sister species to all other extant flowering plants. AS events within 40-70% of the expressed multi-exonic genes per species were found, 27,120 of which are conserved among two or more of the taxa studied. While many events are species specific, many others are shared across long evolutionary distances suggesting they have functional significance. Conservation of AS event data provides an estimate of the number of ancestral AS events present at each node of the tree representing the nine species studied. Furthermore, the presence or absence of AS isoforms between species with different whole genome duplication (WGD) histories provides the opportunity to examine the impact of WDG on AS potential. Examining AS in gene families identifies those with high rates of AS, and conservation can distinguish ancient events vs. recent or species specific adaptations. The MADS-box and SR protein families are found to represent families with low and high occurrences of AS, respectively, yet their AS events were likely present in the MRCA of angiosperms.

454 DNA sequencing technology achieves significant throughput relative to traditional approaches. More than 261,000 ESTs were generated by 454 Life Sciences from cDNA isolated using laser capture microdissection (LCM) from the developmentally important shoot apical meristem (SAM) of maize (Zea mays L.). This single sequencing run annotated >25,000 maize genomic sequences and also captured ∼400 expressed transcripts for which homologous sequences have not yet been identified in other species. Approximately 70% of the ESTs generated in this study had not been captured during a previous EST project conducted using a cDNA library constructed from hand-dissected apex tissue that is highly enriched for SAMs. In addition, at least 30% of the 454-ESTs do not align to any of the ∼648,000 extant maize ESTs using conservative alignment criteria. These results indicate that the combination of LCM and the deep sequencing possible with 454 technology enriches for SAM transcripts not present in current EST collections. RT-PCR was used to validate the expression of 27 genes whose expression had been detected in the SAM via LCM-454 technology, but that lacked orthologs in GenBank. Significantly, transcripts from ∼74% (20/27) of these validated SAM-expressed “orphans” were not detected in meristem-rich immature ears. We conclude that the coupling of LCM and 454 sequencing technologies facilitates the discovery of rare, possibly cell-type-specific transcripts. [The sequence data from this study have been submitted to GenBank under accession nos. DW724699–DW985434.]

Maize (Zea mays) possesses a large, highly repetitive genome, and subsequently a number of reduced-representation sequencing approaches have been used to try and enrich for gene space while eluding difficulties associated with repetitive DNA. This article documents the ability of publicly available maize expressed sequence tag and Genome Survey Sequences (GSSs; many of which were isolated through the use of reduced representation techniques) to recognize and provide coverage of 78 maize full-length cDNAs (FLCs). All 78 FLCs in the dataset were identified by at least three GSSs, indicating that the majority of maize genes have been identified by at least one currently available GSS. Both methyl-filtration and high-Cot enrichment methods provided a 7- to 8-fold increase in gene discovery rates as compared to random sequencing. The available maize GSSs aligned to 75% of the FLC nucleotides used to perform searches, while the expressed sequence tag sequences aligned to 73% of the nucleotides. Our data suggest that at least approximately 95% of maize genes have been tagged by at least one GSS. While the GSSs are very effective for gene identification, relatively few (18%) of the FLCs are completely represented by GSSs. Analysis of the overlap of coverage and bias due to position within a gene suggest that RescueMu, methyl-filtration, and high-Cot methods are at least partially nonredundant.

Premise of the study: Primers were developed to amplify low-copy nuclear genes in Hawaiian Metrosideros (Myrtaceae). Methods and Results: Data from a pooled 454 Titanium run of the partial transcriptomes of four Metrosideros taxa were used to identify the loci of interest. Ten exon-primed intron-crossing (EPIC) markers were amplified and sequenced directly with success in Metrosideros, as well as in a representative selection of Myrtaceae, including Syzygium, Psidium, and Melaleuca for most of the markers. The loci amplified ranged between 500 and 1100 bp, and up to 117 polymorphic sites were observed within an individual gene alignment. Two introns contained microsatellites in some of the species. Conclusions: These novel primer pairs should be useful for phylogenetic analysis and population genetics of a broad range of Myrtaceae, particularly the diverse fleshy-fruited tribes Syzygieae and Myrteae.

Contemporary methods to assay gene expression depend on a stable set of reference transcripts for accurate quantitation. A lack of well-tested reference genes slows progress in characterizing gene expression in high-value specialty crops. In this study, a set of strawberry (Fragaria spp.) constitutively expressed reference genes has been identified by merging digital gene expression data with expression profiling. Constitutive reference candidates were validated using quantitative PCR and hybridization. Several transcripts have been identified that show improved stability across tissues relative to traditional reference transcripts. Results are similar between commercial octoploid strawberry and the diploid model. Our findings also show that while some never-before-used references are appropriate for most applications, even the most stable reference transcripts require careful assessment across the diverse tissues and fruit developmental states before being adopted as controls., Facultad de Ciencias Exactas, Instituto de Fisiología Vegetal

Associations between arbuscular mycorrhizal (AM) fungi and plants are an ancient and widespread plant microbe symbioses. Most land plants can associate with this specialized group of soil fungi (in the Glomeromycota), which enhance plant nutrient uptake in return for C derived from plant photosynthesis. Elucidating the mechanisms involved in the symbiosis between obligate symbionts such as AM fungi and plant roots is challenging because AM fungal transcripts in roots are in low abundance and reference genomes for the fungi have not been available. A deep sequencing metatranscriptomics approach was applied to a wild-type tomato and a tomato mutant (Solanum lycopersicum L. cultivar RioGrande 76R) incapable of supporting a functional AM symbiosis, revealing novel AM fungal and microbial transcripts expressed in colonized roots. We confirm transcripts known to be mycorrhiza associated and report the discovery of more than 500 AM fungal and novel plant transcripts associated with mycorrhizal tomato roots including putative Zn, Fe, aquaporin, and carbohydrate transporters as well as mycorrhizal-associated alternative gene splicing. This analysis provides a fundamental step toward identifying the molecular mechanisms of mineral and carbohydrate exchange during the symbiosis. The utility of this metatranscriptomic approach to explore an obligate biotrophic interaction is illustrated, especially as it relates to agriculturally relevant biological processes. A rbuscular mycorrhizal (AM) fungi are important root symbionts that associate with the majority of land plants including most agricultural species (Smith and Read, 2008). They are obligate mutualistic biotrophs that provide an additional (fungal) pathway of mineral nutrient (mainly inorganic P, N, S, and Zn) uptake from the soil (Allen and Shachar-Hill, 2009; Govindarajulu et al., 2005; Javot et al., 2007), enhance drought tolerance (Aroca et al., 2008), and increased pathogen protection (Liu et al., 2007). In return for soil-derived nutrients, the plant supplies C to the fungus in the form of photosynthesis-derived sugars (Pfeffer et al., 1999). Establishment