Evolutionarily informed deep learning methods for predicting relative transcript abundance from DNA sequence

Deep learning methodologies have revolutionized prediction in many fields and show potential to do the same in molecular biology and genetics. However, applying these methods in their current forms ignores evolutionary dependencies within biological systems and can result in false positives and spurious conclusions. We developed two approaches that account for evolutionary relatedness in machine learning models: (i) gene-family–guided splitting and (ii) ortholog contrasts. The first approach accounts for evolution by constraining model training and testing sets to include different gene families. The second approach uses evolutionarily informed comparisons between orthologous genes to both control for and leverage evolutionary divergence during the training process. The two approaches were explored and validated within the context of mRNA expression level prediction and have the area under the ROC curve (auROC) values ranging from 0.75 to 0.94. Model weight inspections showed biologically interpretable patterns, resulting in the hypothesis that the 3′ UTR is more important for fine-tuning mRNA abundance levels while the 5′ UTR is more important for large-scale changes

Debunking the myth of the endogenous antiangiogenic Vegfaxxxb transcripts

In this opinion article we critically assess evidence for the existence of a family of antiangiogenic vascular endothelial growth factor (Vegfaxxxb) transcripts, arising from the use of a phylogenetically conserved alternative distal splice site within exon 8 of the VEGFA gene. We explain that prior evidence for Vegfaxxxb transcripts in tissues rests heavily upon flawed RT-PCR methodologies, with the extensive use of 5′-tailing in primer design being the main issue. Furthermore, our analysis of large RNA-seq data sets (human and ovine) fails to identify a single Vegfaxxxb transcript. Therefore, we challenge the very existence of Vegfaxxxb transcripts, which further questions the physiological relevance of studies based on the use of ‘anti-VEGFAxxxb’ antibodies. Our analysis has implications for the proposed therapeutic use of isoform-specific anti-VEGFA strategies for treating cancer and retinopathies

A noncanonical auxin-sensing mechanism is required for organ morphogenesis in Arabidopsis

Tissue patterning in multicellular organisms is the output of precise spatio–temporal regulation of gene expression coupled with changes in hormone dynamics. In plants, the hormone auxin regulates growth and development at every stage of a plant's life cycle. Auxin signaling occurs through binding of the auxin molecule to a TIR1/AFB F-box ubiquitin ligase, allowing interaction with Aux/IAA transcriptional repressor proteins. These are subsequently ubiquitinated and degraded via the 26S proteasome, leading to derepression of auxin response factors (ARFs). How auxin is able to elicit such a diverse range of developmental responses through a single signaling module has not yet been resolved. Here we present an alternative auxin-sensing mechanism in which the ARF ARF3/ETTIN controls gene expression through interactions with process-specific transcription factors. This noncanonical hormone-sensing mechanism exhibits strong preference for the naturally occurring auxin indole 3-acetic acid (IAA) and is important for coordinating growth and patterning in diverse developmental contexts such as gynoecium morphogenesis, lateral root emergence, ovule development, and primary branch formation. Disrupting this IAA-sensing ability induces morphological aberrations with consequences for plant fitness. Therefore, our findings introduce a novel transcription factor-based mechanism of hormone perception in plants. Note that there is a CORRIGENDUM to this article: http://eprints.whiterose.ac.uk/132306/ http://genesdev.cshlp.org/content/31/17/1821.ful

Deleterious Mutation Burden and Its Association with Complex Traits in Sorghum (Sorghum bicolor)

Sorghum (Sorghum bicolor L.) is a major food cereal for millions of people worldwide. The sorghum genome, like other species, accumulates deleterious mutations, likely impacting its fitness. The lack of recombination, drift, and the coupling with favorable loci impede the removal of deleterious mutations from the genome by selection. To study how deleterious variants impact phenotypes, we identified putative deleterious mutations among ∼5.5 M segregating variants of 229 diverse biomass sorghum lines. We provide the whole-genome estimate of the deleterious burden in sorghum, showing that ∼33% of nonsynonymous substitutions are putatively deleterious. The pattern of mutation burden varies appreciably among racial groups. Across racial groups, the mutation burden correlated negatively with biomass, plant height, specific leaf area (SLA), and tissue starch content (TSC), suggesting that deleterious burden decreases trait fitness. Putatively deleterious variants explain roughly one-half of the genetic variance. However, there is only moderate improvement in total heritable variance explained for biomass (7.6%) and plant height (average of 3.1% across all stages). There is no advantage in total heritable variance for SLA and TSC. The contribution of putatively deleterious variants to phenotypic diversity therefore appears to be dependent on the genetic architecture of traits. Overall, these results suggest that incorporating putatively deleterious variants into genomic models slightly improves prediction accuracy because of extensive linkage. Knowledge of deleterious variants could be leveraged for sorghum breeding through either genome editing and/or conventional breeding that focuses on the selection of progeny with fewer deleterious alleles

MiR-662 is associated with metastatic relapse in early-stage breast cancer and promotes metastasis by stimulating cancer cell stemness

Background Breast cancer (BC) metastasis, which often occurs in bone, contributes substantially to mortality. MicroRNAs play a fundamental role in BC metastasis, although microRNA-regulated mechanisms driving metastasis progression remain poorly understood. Methods MiRome analysis in serum from BC patients was performed by TaqMan™ low-density array. MiR-662 was overexpressed following MIMIC-transfection or lentivirus transduction. Animal models were used to investigate the role of miR-662 in BC (bone) metastasis. The effect of miR-662-overexpressing BC cell conditioned medium on osteoclastogenesis was investigated. ALDEFLUOR assays were performed to study BC stemness. RNA-sequencing transcriptomic analysis of miR-662-overexpressing BC cells was performed to evaluate gene expression changes. Results High levels of hsa-miR-662 (miR-662) in serum from BC patients, at baseline (time of surgery), were associated with future recurrence in bone. At an early-stage of the metastatic disease, miR-662 could mask the presence of BC metastases in bone by inhibiting the differentiation of bone-resorbing osteoclasts. Nonetheless, metastatic miR-662-overexpressing BC cells then progressed as overt osteolytic metastases thanks to increased stem cell-like traits. Conclusions MiR-662 is involved in BC metastasis progression, suggesting it may be used as a prognostic marker to identify BC patients at high risk of metastasis

Machine learning-enabled phenotyping for GWAS and TWAS of WUE traits in 869 field-grown sorghum accessions

Sorghum (Sorghum bicolor) is a model C4 crop made experimentally tractable by extensive genomic and genetic resources. Biomass sorghum is studied as a feedstock for biofuel and forage. Mechanistic modeling suggests that reducing stomatal conductance (gs) could improve sorghum intrinsic water use efficiency (iWUE) and biomass production. Phenotyping to discover genotype-to-phenotype associations remains a bottleneck in understanding the mechanistic basis for natural variation in gs and iWUE. This study addressed multiple methodological limitations. Optical tomography and a machine learning tool were combined to measure stomatal density (SD). This was combined with rapid measurements of leaf photosynthetic gas exchange and specific leaf area (SLA). These traits were the subject of genome-wide association study and transcriptome-wide association study across 869 field-grown biomass sorghum accessions. The ratio of intracellular to ambient CO2 was genetically correlated with SD, SLA, gs, and biomass production. Plasticity in SD and SLA was interrelated with each other and with productivity across wet and dry growing seasons. Moderate-to-high heritability of traits studied across the large mapping population validated associations between DNA sequence variation or RNA transcript abundance and trait variation. A total of 394 unique genes underpinning variation in WUE-related traits are described with higher confidence because they were identified in multiple independent tests. This list was enriched in genes whose Arabidopsis (Arabidopsis thaliana) putative orthologs have functions related to stomatal or leaf development and leaf gas exchange, as well as genes with nonsynonymous/missense variants. These advances in methodology and knowledge will facilitate improving C4 crop WUE

Transcriptional responses of winter barley to cold indicate nucleosome remodelling as a specific feature of crown tissues

We report a series of microarray-based comparisons of gene expression in the leaf and crown of the winter barley cultivar Luxor, following the exposure of young plants to various periods of low (above and below zero) temperatures. A transcriptomic analysis identified genes which were either expressed in both the leaf and crown, or specifically in one or the other. Among the former were genes responsible for calcium and abscisic acid signalling, polyamine synthesis, late embryogenesis abundant proteins and dehydrins. In the crown, the key organ for cereal overwintering, cold treatment induced transient changes in the transcription of nucleosome assembly genes, and especially H2A and HTA11, which have been implicated in cold sensing in Arabidopsis thaliana. In the leaf, various heat-shock proteins were induced. Differences in expression pattern between the crown and leaf were frequent for genes involved in certain pathways responsible for osmolyte production (sucrose and starch, raffinose, γ-aminobutyric acid metabolism), sugar signalling (trehalose metabolism) and secondary metabolism (lignin synthesis). The action of proteins with antifreeze activity, which were markedly induced during hardening, was demonstrated by a depression in the ice nucleation temperature