Origin and evolution of the octoploid strawberry genome.

Cultivated strawberry emerged from the hybridization of two wild octoploid species, both descendants from the merger of four diploid progenitor species into a single nucleus more than 1 million years ago. Here we report a near-complete chromosome-scale assembly for cultivated octoploid strawberry (Fragaria × ananassa) and uncovered the origin and evolutionary processes that shaped this complex allopolyploid. We identified the extant relatives of each diploid progenitor species and provide support for the North American origin of octoploid strawberry. We examined the dynamics among the four subgenomes in octoploid strawberry and uncovered the presence of a single dominant subgenome with significantly greater gene content, gene expression abundance, and biased exchanges between homoeologous chromosomes, as compared with the other subgenomes. Pathway analysis showed that certain metabolomic and disease-resistance traits are largely controlled by the dominant subgenome. These findings and the reference genome should serve as a powerful platform for future evolutionary studies and enable molecular breeding in strawberry.

Haplotype-phased genome and evolution of phytonutrient pathways of tetraploid blueberry.

BACKGROUND: Highbush blueberry (Vaccinium corymbosum) has long been consumed for its unique flavor and composition of health-promoting phytonutrients. However, breeding efforts to improve fruit quality in blueberry have been greatly hampered by the lack of adequate genomic resources and a limited understanding of the underlying genetics encoding key traits. The genome of highbush blueberry has been particularly challenging to assemble due, in large part, to its polyploid nature and genome size. FINDINGS: Here, we present a chromosome-scale and haplotype-phased genome assembly of the cultivar "Draper," which has the highest antioxidant levels among a diversity panel of 71 cultivars and 13 wild Vaccinium species. We leveraged this genome, combined with gene expression and metabolite data measured across fruit development, to identify candidate genes involved in the biosynthesis of important phytonutrients among other metabolites associated with superior fruit quality. Genome-wide analyses revealed that both polyploidy and tandem gene duplications modified various pathways involved in the biosynthesis of key phytonutrients. Furthermore, gene expression analyses hint at the presence of a spatial-temporal specific dominantly expressed subgenome including during fruit development. CONCLUSIONS: These findings and the reference genome will serve as a valuable resource to guide future genome-enabled breeding of important agronomic traits in highbush blueberry.

Autogenous cross-regulation of Quaking mRNA processing and translation balances Quaking functions in splicing and translation.

Quaking protein isoforms arise from a single Quaking gene and bind the same RNA motif to regulate splicing, translation, decay, and localization of a large set of RNAs. However, the mechanisms by which Quaking expression is controlled to ensure that appropriate amounts of each isoform are available for such disparate gene expression processes are unknown. Here we explore how levels of two isoforms, nuclear Quaking-5 (Qk5) and cytoplasmic Qk6, are regulated in mouse myoblasts. We found that Qk5 and Qk6 proteins have distinct functions in splicing and translation, respectively, enforced through differential subcellular localization. We show that Qk5 and Qk6 regulate distinct target mRNAs in the cell and act in distinct ways on their own and each other's transcripts to create a network of autoregulatory and cross-regulatory feedback controls. Morpholino-mediated inhibition of Qk translation confirms that Qk5 controls Qk RNA levels by promoting accumulation and alternative splicing of Qk RNA, whereas Qk6 promotes its own translation while repressing Qk5. This Qk isoform cross-regulatory network responds to additional cell type and developmental controls to generate a spectrum of Qk5/Qk6 ratios, where they likely contribute to the wide range of functions of Quaking in development and cancer.

Distinct and shared functions of ALS-associated proteins TDP-43, FUS and TAF15 revealed by multisystem analyses.

The RNA-binding protein (RBP) TAF15 is implicated in amyotrophic lateral sclerosis (ALS). To compare TAF15 function to that of two ALS-associated RBPs, FUS and TDP-43, we integrate CLIP-seq and RNA Bind-N-Seq technologies, and show that TAF15 binds to ∼4,900 RNAs enriched for GGUA motifs in adult mouse brains. TAF15 and FUS exhibit similar binding patterns in introns, are enriched in 3' untranslated regions and alter genes distinct from TDP-43. However, unlike FUS and TDP-43, TAF15 has a minimal role in alternative splicing. In human neural progenitors, TAF15 and FUS affect turnover of their RNA targets. In human stem cell-derived motor neurons, the RNA profile associated with concomitant loss of both TAF15 and FUS resembles that observed in the presence of the ALS-associated mutation FUS R521G, but contrasts with late-stage sporadic ALS patients. Taken together, our findings reveal convergent and divergent roles for FUS, TAF15 and TDP-43 in RNA metabolism.

Quaking promotes monocyte differentiation into pro-atherogenic macrophages by controlling pre-mRNA splicing and gene expression.

A hallmark of inflammatory diseases is the excessive recruitment and influx of monocytes to sites of tissue damage and their ensuing differentiation into macrophages. Numerous stimuli are known to induce transcriptional changes associated with macrophage phenotype, but posttranscriptional control of human macrophage differentiation is less well understood. Here we show that expression levels of the RNA-binding protein Quaking (QKI) are low in monocytes and early human atherosclerotic lesions, but are abundant in macrophages of advanced plaques. Depletion of QKI protein impairs monocyte adhesion, migration, differentiation into macrophages and foam cell formation in vitro and in vivo. RNA-seq and microarray analysis of human monocyte and macrophage transcriptomes, including those of a unique QKI haploinsufficient patient, reveal striking changes in QKI-dependent messenger RNA levels and splicing of RNA transcripts. The biological importance of these transcripts and requirement for QKI during differentiation illustrates a central role for QKI in posttranscriptionally guiding macrophage identity and function.

Quaking and PTB control overlapping splicing regulatory networks during muscle cell differentiation.

Alternative splicing contributes to muscle development, but a complete set of muscle-splicing factors and their combinatorial interactions are unknown. Previous work identified ACUAA ("STAR" motif) as an enriched intron sequence near muscle-specific alternative exons such as Capzb exon 9. Mass spectrometry of myoblast proteins selected by the Capzb exon 9 intron via RNA affinity chromatography identifies Quaking (QK), a protein known to regulate mRNA function through ACUAA motifs in 3' UTRs. We find that QK promotes inclusion of Capzb exon 9 in opposition to repression by polypyrimidine tract-binding protein (PTB). QK depletion alters inclusion of 406 cassette exons whose adjacent intron sequences are also enriched in ACUAA motifs. During differentiation of myoblasts to myotubes, QK levels increase two- to threefold, suggesting a mechanism for QK-responsive exon regulation. Combined analysis of the PTB- and QK-splicing regulatory networks during myogenesis suggests that 39% of regulated exons are under the control of one or both of these splicing factors. This work provides the first evidence that QK is a global regulator of splicing during muscle development in vertebrates and shows how overlapping splicing regulatory networks contribute to gene expression programs during differentiation.

Divergent roles of ALS-linked proteins FUS/TLS and TDP-43 intersect in processing long pre-mRNAs.

FUS/TLS (fused in sarcoma/translocated in liposarcoma) and TDP-43 are integrally involved in amyotrophic lateral sclerosis (ALS) and frontotemporal dementia. We found that FUS/TLS binds to RNAs from >5,500 genes in mouse and human brain, primarily through a GUGGU-binding motif. We identified a sawtooth-like binding pattern, consistent with co-transcriptional deposition of FUS/TLS. Depletion of FUS/TLS from the adult nervous system altered the levels or splicing of >950 mRNAs, most of which are distinct from RNAs dependent on TDP-43. Abundance of only 45 RNAs was reduced after depletion of either TDP-43 or FUS/TLS from mouse brain, but among these were mRNAs that were transcribed from genes with exceptionally long introns and that encode proteins that are essential for neuronal integrity. Expression levels of a subset of these were lowered after TDP-43 or FUS/TLS depletion in stem cell-derived human neurons and in TDP-43 aggregate-containing motor neurons in sporadic ALS, supporting a common loss-of-function pathway as one component underlying motor neuron death from misregulation of TDP-43 or FUS/TLS.

Integrative genome-wide analysis reveals cooperative regulation of alternative splicing by hnRNP proteins.

Understanding how RNA binding proteins control the splicing code is fundamental to human biology and disease. Here, we present a comprehensive study to elucidate how heterogeneous nuclear ribonucleoparticle (hnRNP) proteins, among the most abundant RNA binding proteins, coordinate to regulate alternative pre-mRNA splicing (AS) in human cells. Using splicing-sensitive microarrays, crosslinking and immunoprecipitation coupled with high-throughput sequencing (CLIP-seq), and cDNA sequencing, we find that more than half of all AS events are regulated by multiple hnRNP proteins and that some combinations of hnRNP proteins exhibit significant synergy, whereas others act antagonistically. Our analyses reveal position-dependent RNA splicing maps, in vivo consensus binding sites, a surprising level of cross- and autoregulation among hnRNP proteins, and the coordinated regulation by hnRNP proteins of dozens of other RNA binding proteins and genes associated with cancer. Our findings define an unprecedented degree of complexity and compensatory relationships among hnRNP proteins and their splicing targets that likely confer robustness to cells.

A post-transcriptional regulatory switch in polypyrimidine tract-binding proteins reprograms alternative splicing in developing neurons.

Many metazoan gene transcripts exhibit neuron-specific splicing patterns, but the developmental control of these splicing events is poorly understood. We show that the splicing of a large group of exons is reprogrammed during neuronal development by a switch in expression between two highly similar polypyrimidine tract-binding proteins, PTB and nPTB (neural PTB). PTB is a well-studied regulator of alternative splicing, but nPTB is a closely related paralog whose functional relationship to PTB is unknown. In the brain, nPTB protein is specifically expressed in post-mitotic neurons, whereas PTB is restricted to neuronal precursor cells (NPC), glia, and other nonneuronal cells. Interestingly, nPTB mRNA transcripts are found in NPCs and other nonneuronal cells, but in these cells nPTB protein expression is repressed. This repression is due in part to PTB-induced alternative splicing of nPTB mRNA, leading to nonsense-mediated decay (NMD). However, we find that even properly spliced mRNA fails to express nPTB protein when PTB is present, indicating contributions from additional post-transcriptional mechanisms. The PTB-controlled repression of nPTB results in a mutually exclusive pattern of expression in the brain, where the loss of PTB in maturing neurons allows the synthesis of nPTB in these cells. To examine the consequences of this switch, we used splicing-sensitive microarrays to identify different sets of exons regulated by PTB, nPTB, or both proteins. During neuronal differentiation, the splicing of these exon sets is altered as predicted from the observed changes in PTB and nPTB expression. These data show that the post-transcriptional switch from PTB to nPTB controls a widespread alternative splicing program during neuronal development.

Cell type and culture condition-dependent alternative splicing in human breast cancer cells revealed by splicing-sensitive microarrays.

Growing evidence indicates that alternative or aberrant pre-mRNA splicing takes place during the development, progression, and metastasis of breast cancer. However, which splicing changes that might contribute directly to tumorigenesis or cancer progression remain to be elucidated. We used splicing-sensitive microarrays to detect differences in alternative splicing between two breast cancer cell lines, MCF7 (estrogen receptor positive) and MDA-MB-231 (estrogen receptor negative), as well as cultured human mammary epithelial cells. Several splicing alterations in genes, including CD44, FAS, RBM9, hnRNPA/B, APLP2, and MYL6, were detected by the microarray and verified by reverse transcription-PCR. We also compared splicing in these breast cancer cells cultured in either two-dimensional flat dishes or in three-dimensional Matrigel conditions. Only a subset of the splicing differences that distinguish MCF7 cells from MDA-MB-231 cells under two-dimensional culture condition is retained under three-dimensional conditions, suggesting that alternative splicing events are influenced by the geometry of the culture conditions of these cells. Further characterization of splicing patterns of several genes in MCF7 cells grown in Matrigel and in xenograft in nude mice shows that splicing is similar under both conditions. Thus, our oligonucleotide microarray can effectively detect changes in alternative splicing in different cells or in the same cells grown in different environments. Our findings also illustrate the potential for understanding gene expression with resolution of alternative splicing in the study of breast cancer.

Detection and measurement of alternative splicing using splicing-sensitive microarrays.

Splicing and alternative splicing are major processes in the interpretation and expression of genetic information for metazoan organisms. The study of splicing is moving from focused attention on the regulatory mechanisms of a selected set of paradigmatic alternative splicing events to questions of global integration of splicing regulation with genome and cell function. For this reason, parallel methods for detecting and measuring alternative splicing are necessary. We have adapted the splicing-sensitive oligonucleotide microarrays used to estimate splicing efficiency in yeast to the study of alternative splicing in vertebrate cells and tissues. We use gene models incorporating knowledge about splicing to design oligonucleotides specific for discriminating alternatively spliced mRNAs from each other. Here we present the main strategies for design, application, and analysis of spotted oligonucleotide arrays for detection and measurement of alternative splicing. We demonstrate these strategies using a two-intron yeast gene that has been altered to produce different amounts of alternatively spliced RNAs, as well as by profiling alternative splicing in NCI 60 cancer cell lines.

Exploring functional relationships between components of the gene expression machinery.

Eukaryotic gene expression requires the coordinated activity of many macromolecular machines including transcription factors and RNA polymerase, the spliceosome, mRNA export factors, the nuclear pore, the ribosome and decay machineries. Yeast carrying mutations in genes encoding components of these machineries were examined using microarrays to measure changes in both pre-mRNA and mRNA levels. We used these measurements as a quantitative phenotype to ask how steps in the gene expression pathway are functionally connected. A multiclass support vector machine was trained to recognize the gene expression phenotypes caused by these mutations. In several cases, unexpected phenotype assignments by the computer revealed functional roles for specific factors at multiple steps in the gene expression pathway. The ability to resolve gene expression pathway phenotypes provides insight into how the major machineries of gene expression communicate with each other.