Time‐specific and pleiotropic quantitative trait loci coordinately modulate stem growth in Populus

In perennial woody plants, the coordinated increase of stem height and diameter during juvenile growth improves competitiveness (i.e. access to light); however, the factors underlying variation in stem growth remain unknown in trees. Here, we used linkage‐linkage disequilibrium (linkage‐LD) mapping to decipher the genetic architecture underlying three growth traits during juvenile stem growth. We used two Populus populations: a linkage mapping population comprising a full‐sib family of 1,200 progeny and an association mapping panel comprising 435 unrelated individuals from nearly the entire natural range of Populus tomentosa. We mapped 311 quantitative trait loci (QTL) for three growth traits at 12 timepoints to 42 regions in 17 linkage groups. Of these, 28 regions encompassing 233 QTL were annotated as 27 segmental homology regions (SHRs). Using SNPs identified by whole‐genome re‐sequencing of the 435‐member association mapping panel, we identified significant SNPs (P ≤ 9.4 × 10^?7) within 27 SHRs that affect stem growth at nine timepoints with diverse additive and dominance patterns, and these SNPs exhibited complex allelic epistasis over the juvenile growth period. Nineteen genes linked to potential causative alleles that have time‐specific or pleiotropic effects, and mostly overlapped with significant signatures of selection within SHRs between climatic regions represented by the association mapping panel. Five genes with potential time‐specific effects showed species‐specific temporal expression profiles during the juvenile stages of stem growth in five representative Populus species. Our observations revealed the importance of considering temporal genetic basis of complex traits, which will facilitate the molecular design of tree ideotypes.     

Sl‐lncRNA15492 interacts with Sl‐miR482a and affects Solanum lycopersicum immunity against Phytophthora infestans

Long non‐coding RNAs (lncRNAs) are involved in the resistance of plants to infection by pathogens via interactions with microRNAs (miRNAs). Long non‐coding RNAs are cleaved by miRNAs to produce phased small interfering RNAs (phasiRNAs), which, as competing endogenous RNAs (ceRNAs), function as decoys for mature miRNAs, thus inhibiting their expression, and contain pre‐miRNA sequences to produce mature miRNAs. However, whether lncRNAs and miRNAs mediate other molecular mechanisms during plant resistance to pathogens is unknown. In this study, as a positive regulator, Sl‐lncRNA15492 from tomato (Solanum lycopersicum Zaofen No. 2) plants affected tomato resistance to Phytophthora infestans. Gain‐ and loss‐of‐function experiments and RNA ligase‐mediated 5′‐amplification of cDNA ends (RLM‐5′ RACE) also revealed that Sl‐miR482a was negatively involved in tomato resistance by targeting Sl‐NBS‐LRR genes and that silencing of Sl‐NBS‐LRR1 decreased tomato resistance. Sl‐lncRNA15492 inhibited the expression of mature Sl‐miR482a, whose precursor was located within the antisense sequence of Sl‐lncRNA15492. Further degradome analysis and additional RLM‐5′ RACE experiments verified that mature Sl‐miR482a could also cleave Sl‐lncRNA15492. These results provide a mechanism by which lncRNAs might inhibit precursor miRNA expression through antisense strands of lncRNAs, and demonstrate that Sl‐lncRNA15492 and Sl‐miR482a mutually inhibit the maintenance of Sl‐NBS‐LRR1 homeostasis during tomato resistance to P. infestans.     

Long noncoding RNAs involve in resistance to Verticillium dahliae,a fungal disease in cotton

Long noncoding RNAs (lncRNAs) have several known functions in plant development, but their possible roles in responding to plant disease remain largely unresolved. In this study, we described a comprehensive disease‐responding lncRNA profiles in defence against a cotton fungal disease Verticillium dahliae. We further revealed the conserved and specific characters of disease‐responding process between two cotton species. Conservatively for two cotton species, we found the expression dominance of induced lncRNAs in the Dt subgenome, indicating a biased induction pattern in the co‐existing subgenomes of allotetraploid cotton. Comparative analysis of lncRNA expression and their proposed functions in resistant Gossypium barbadense cv. ‘7124’ versus susceptible Gossypium hirsutum cv. ‘YZ1’ revealed their distinct disease response mechanisms. Species‐specific (LS) lncRNAs containing more SNPs displayed a fiercer inducing level postinfection than the species‐conserved (core) lncRNAs. Gene Ontology enrichment of LS lncRNAs and core lncRNAs indicates distinct roles in the process of biotic stimulus. Further functional analysis showed that two core lncRNAs, GhlncNAT‐ANX2‐ and GhlncNAT‐RLP7‐silenced seedlings, displayed an enhanced resistance towards V. dahliae and Botrytis cinerea, possibly associated with the increased expression of LOX1 and LOX2. This study represents the first characterization of lncRNAs involved in resistance to fungal disease and provides new clues to elucidate cotton disease response mechanism.     

Roles of small RNAs in soybean defense against Phytophthora sojae infection

The genus Phytophthora consists of many notorious pathogens of crops and forestry trees. At present, battling Phytophthora diseases is challenging due to a lack of understanding of their pathogenesis. We investigated the role of small RNAs in regulating soybean defense in response to infection by Phytophthora sojae, the second most destructive pathogen of soybean. Small RNAs, including microRNAs (miRNAs) and small interfering RNAs (siRNAs), are universal regulators that repress target gene expression in eukaryotes. We identified known and novel small RNAs that differentially accumulated during P. sojae infection in soybean roots. Among them, miR393 and miR166 were induced by heat‐inactivated P. sojae hyphae, indicating that they may be involved in soybean basal defense. Indeed, knocking down the level of mature miR393 led to enhanced susceptibility of soybean to P. sojae; furthermore, the expression of isoflavonoid biosynthetic genes was drastically reduced in miR393 knockdown roots. These data suggest that miR393 promotes soybean defense against P. sojae. In addition to miRNAs, P. sojae infection also resulted in increased accumulation of phased siRNAs (phasiRNAs) that are predominantly generated from canonical resistance genes encoding nucleotide binding‐leucine rich repeat proteins and genes encoding pentatricopeptide repeat‐containing proteins. This work identifies specific miRNAs and phasiRNAs that regulate defense‐associated genes in soybean during Phytophthora infection.     

LncKdm2b controls self‐renewal of embryonic stem cells via activating expression of transcription factor Zbtb3

Divergent long noncoding RNAs (lncRNAs) represent a major lncRNA biotype in mouse and human genomes. The biological and molecular functions of the divergent lncRNAs remain largely unknown. Here, we show that lncKdm2b, a divergent lncRNA for Kdm2b gene, is conserved among five mammalian species and highly expressed in embryonic stem cells (ESCs) and early embryos. LncKdm2b knockout impairs ESC self‐renewal and causes early embryonic lethality. LncKdm2b can activate Zbtb3 by promoting the assembly and ATPase activity of Snf2‐related CREBBP activator protein (SRCAP) complex in trans. Zbtb3 potentiates the ESC self‐renewal in a Nanog‐dependent manner. Finally, Zbtb3 deficiency impairs the ESC self‐renewal and early embryonic development. Therefore, our findings reveal that lncRNAs may represent an additional layer of the regulation of ESC self‐renewal and early embryogenesis.     

Deregulated lncRNA expression profile in the mouse lung adenocarcinomas with KRAS‐G12D mutation and P53 knockout

Recent studies have demonstrated that aberrant long non‐coding RNAs (lncRNAs) expression are suggested to be closely associated with multiple human diseases, lung cancer included. However, the roles of lncRNAs in lung cancer are not well understood. In this study, we used microarrays to investigate the aberrantly expressed lncRNAs in the mouse lung adenocarcinoma with P53 knockout and the Kras^G12D mutation. Results revealed that 6424 lncRNAs were differentially expressed (≥ 2‐fold change, P < .05). Two hundred and ten lncRNAs showed more than 8‐fold change and conserved across human and were further analysed in the primary mouse lung adenocarcinoma KP cells, which were isolated from the p53 knockout and the Kras^G12D mutation mice. Among all the 210 lncRNAs, 11 lncRNAs' expression was regulated by P53, 33 lncRNAs by KRAS and 13 lncRNAs by hypoxia in the primary KP cells, respectively. NONMMUT015812, which was remarkably up‐regulated in the mouse lung adenocarcinoma and negatively regulated by the P53 re‐expression, was detected to analyse its cellular function. Results showed that knockdown of NONMMUT015812 by shRNAs decreased proliferation and migration abilities of KP cells. Among those aberrantly expressed lncRNAs in the mouse lung adenocarcinoma, NONMMUT015812 was a potential oncogene.     

Comprehensive analysis of a ceRNA network reveals potential prognostic cytoplasmic lncRNAs involved in HCC progression

The aberrant expression of long noncoding RNAs (lncRNAs) has drawn increasing attention in the field of hepatocellular carcinoma (HCC) biology. In the present study, we obtained the expression profiles of lncRNAs, microRNAs (miRNAs), and messenger RNAs (mRNAs) in 371 HCC tissues and 50 normal tissues from The Cancer Genome Atlas (TCGA) and identified hepatocarcinogenesis-specific differentially expressed genes (DEGs, log fold change ≥ 2, FDR < 0.01), including 753 lncRNAs, 97 miRNAs, and 1,535 mRNAs. Because the specific functions of lncRNAs are closely related to their intracellular localizations and because the cytoplasm is the main location for competitive endogenous RNA (ceRNA) action, we analyzed not only the interactions among these DEGs but also the distributions of lncRNAs (cytoplasmic, nuclear or both). Then, an HCC-associated deregulated ceRNA network consisting of 37 lncRNAs, 10 miRNAs, and 26 mRNAs was constructed after excluding those lncRNAs located only in the nucleus. Survival analysis of this network demonstrated that 15 lncRNAs, 3 miRNAs, and 16 mRNAs were significantly correlated with the overall survival of HCC patients (p < 0.01). Through multivariate Cox regression and lasso analysis, a risk score system based on 13 lncRNAs was constructed, which showed good discrimination and predictive ability for HCC patient survival time. This ceRNA network-construction approach, based on lncRNA distribution, not only narrowed the scope of target lncRNAs but also provided specific candidate molecular biomarkers for evaluating the prognosis of HCC, which will help expand our understanding of the ceRNA mechanisms involved in the early development of HCC.     

Three vital RNA functions and interactions in the process of silk gland apoptosis in silkworm Bombyx mori

The Silkworm Bombyx mori is an important insect in terms of economics and a model organism with a complete metamorphosis. The economic importance of silkworms is dependent on the functions of the silkgland, a specialized organ that synthesizes silk proteins. The silk gland undergoes massive degeneration during the larval to pupal stage, which involves in cell apoptosis. In this paper, high throughput sequencing was used to detect the expression of messenger RNA (mRNA), long noncoding RNA (lncRNA), and microRNA (miRNA) from silk glands of Day 3 in the fifth instar larvae (L5D3) and the spinning 36h (sp36h). We analyzed the Gene Ontology (GO) functions of target genes of the differentially expressed lncRNAs and miRNAs. We investigated the regulations of mRNA, lncRNA, and miRNA on silk gland apoptosis in L5D3 and sp36h. In total, 10,947 lncRNAs were detected in the silk gland and the index number TCONS‐00021360 lncRNA may be involved in the process of apoptosis. In addition, 344 miRNAs targeted 285 mRNAs were related to the death process under GO entry. The results indicated that miRNAs play an important role in the molecular regulation of the silk gland apoptosis compared with that of lncRNAs. Finally, we screened 746 lncRNAs and 20 miRNAs that might interact with BmDredd, and drew an interaction network among them.