Long noncoding RNAs (lncRNAs) have emerged as playing crucial roles in abiotic stress responsive regulation, however, the mechanism of lncRNAs underlying drought-tolerance remains largely unknown in cassava, an important tropical and sub-tropical root crop of remarkable drought tolerance. In this study, a total of 833 high-confidence lncRNAs, including 652 intergenic and 181 anti-sense lncRNAs, were identified in cassava leaves and root using strand-specific RNA-seq technology, of which 124 were drought-responsive. Trans-regulatory co-expression network revealed that lncRNAs exhibited tissue-specific expression patterns and they preferred to function differently in distinct tissues: e.g., cell-related metabolism, cell wall, and RNA regulation of transcription in folded leaf (FL); degradation of major carbohydrate (CHO) metabolism, calvin cycle and light reaction, light signaling, and tetrapyrrole synthesis in full expanded leaf (FEL); synthesis of major CHO metabolism, nitrogen-metabolism, photosynthesis, and redox in bottom leaf (BL); and hormone metabolism, secondary metabolism, calcium signaling, and abiotic stress in root (RT). In addition, 27 lncRNA-mRNA pairs referred to cis-acting regulation were identified, and these lncRNAs regulated the expression of their neighboring genes mainly through hormone metabolism, RNA regulation of transcription, and signaling of receptor kinase. Besides, 11 lncRNAs were identified acting as putative target mimics of known miRNAs in cassava. Finally, five drought-responsive lncRNAs and 13 co-expressed genes involved in trans-acting, cis-acting, or target mimic regulation were selected and confirmed by qRT-PCR. These findings provide a comprehensive view of cassava lncRNAs in response to drought stress, which will enable in-depth functional analysis in the future.

Coffea arabica L. (C. arabica) is an economically important agricultural crop and the most popular beverage worldwide. To analyze genetic diversity and provide genetic resources for the selection and breeding of superior varieties of C. arabica, 61 cultivated Arabica coffee accessions were analyzed in the study, including 12 resequencing accessions from previous research and 49 accessions that were resequenced in this study. Single nucleotide polymorphisms (SNPs) and insertion–deletions (InDels) were statistically analyzed. Based on SNP variations, a genetic structure analysis, phylogenetic tree construction, and principal component analysis were performed for the 61 coffee accessions. The results showed that a total of 805.46 Gb of raw whole-genome resequencing data was obtained from the 61 coffee accessions, with 781.29 Gb of high-quality sequencing data after filtering. In total, 7,013,820 SNP sites and 1,074,329 InDel sites were detected. The average sequencing depth ranged from 6.69× to 19.35×, and the coverage ranged from 85.49% to 96.43%. The population genetic structure and phylogenetic analysis of the 61 coffee accessions revealed four lineages, suggesting that they had at least four ancestral genetic components. Catimor exhibited the highest genetic diversity, while Geisha had the lowest genetic diversity. The selective sweep analysis indicated that among the selected genes in Catimor, disease-resistance genes were significantly more numerous than in other coffee varieties. The genome resequencing data and genetic markers identified from the 61 cultivated Arabica coffee materials provided insights into the genetic variation in Arabica coffee germplasm and facilitated extensive genetic research.

The MADS-box gene family constitutes a vital group of transcription factors that play significant roles in regulating plant growth, development, and signal transduction processes. However, research on the MADS-box genes in cassava (Manihot esculenta) has been relatively limited. To gain deeper insights into the functions of the MADS-box genes in cassava development, in this study, we undertook a comprehensive genome-wide identification of the MADS-box gene family in cassava. We identified a total of 86 MADS-box genes with complete domains in the cassava genome, designated as MeMADS01 to MeMADS86. Through bioinformatic analyses, we investigated the basic physicochemical properties, conserved motifs, chromosomal locations, and phylogenetic relationships of the cassava MADS-box genes. The MADS-box gene family of cassava exhibited conservation, as well as species-specific characteristics, with intron loss being a predominant mode of evolution for the MADS-box genes. Expression pattern variations in the MeMADS genes across different tissues offer insights into their potential functions. Time-ordered gene co-expression network (TO-GCN), transcriptome data, and RT-qPCR analysis suggested the responsiveness of the MADS-box genes to drought stress. Meanwhile, MeMADS12 might be involved in regulating flowering under drought conditions via an ABA (abscisic acid)-dependent pathway. These findings provide valuable resources for a deeper understanding of the biological roles of the MADS-box genes in cassava.