The Genetic Landscape of Familial Pulmonary Fibrosis.

Rationale and Objectives: Up to 20% of idiopathic interstitial lung disease is familial, referred to as familial pulmonary fibrosis (FPF). An integrated analysis of FPF genetic risk was performed by comprehensively evaluating for genetic rare variants (RVs) in a large cohort of FPF kindreds. Methods: Whole-exome sequencing and/or candidate gene sequencing from affected individuals in 569 FPF kindreds was performed, followed by cosegregation analysis in large kindreds, gene burden analysis, gene-based risk scoring, cell-type enrichment analysis, and coexpression network construction. Measurements and Main Results: It was found that 14.9-23.4% of genetic risk in kindreds could be explained by RVs in genes previously linked to FPF, predominantly telomere-related genes. New candidate genes were identified in a small number of families-including SYDE1, SERPINB8, GPR87, and NETO1-and tools were developed for evaluation and prioritization of RV-containing genes across kindreds. Several pathways were enriched for RV-containing genes in FPF, including focal adhesion and mitochondrial complex I assembly. By combining single-cell transcriptomics with prioritized candidate genes, expression of RV-containing genes was discovered to be enriched in smooth muscle cells, type II alveolar epithelial cells, and endothelial cells. Conclusions: In the most comprehensive FPF genetic study to date, the prevalence of RVs in known FPF-related genes was defined, and new candidate genes and pathways relevant to FPF were identified. However, new RV-containing genes shared across multiple kindreds were not identified, thereby suggesting that heterogeneous genetic variants involving a variety of genes and pathways mediate genetic risk in most FPF kindreds.

Bioinformatics resources facilitate understanding and harnessing clinical research of SARS-CoV-2.

The coronavirus disease 2019 (COVID-19) pandemic, caused by the coronavirus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has created an unprecedented threat to public health. The pandemic has been sweeping the globe, impacting more than 200 countries, with more outbreaks still lurking on the horizon. At the time of the writing, no approved drugs or vaccines are available to treat COVID-19 patients, prompting an urgent need to decipher mechanisms underlying the pathogenesis and develop curative treatments. To fight COVID-19, researchers around the world have provided specific tools and molecular information for SARS-CoV-2. These pieces of information can be integrated to aid computational investigations and facilitate clinical research. This paper reviews current knowledge, the current status of drug development and various resources for key steps toward effective treatment of COVID-19, including the phylogenetic characteristics, genomic conservation and interaction data. The final goal of this paper is to provide information that may be utilized in bioinformatics approaches and aid target prioritization and drug repurposing. Several SARS-CoV-2-related tools/databases were reviewed, and a web-portal named OverCOVID (http://bis.zju.edu.cn/overcovid/) is constructed to provide a detailed interpretation of SARS-CoV-2 basics and share a collection of resources that may contribute to therapeutic advances. These information could improve researchers' understanding of SARS-CoV-2 and help to accelerate the development of new antiviral treatments.

Proteogenomic Characterization of the Pathogenic Fungus Aspergillus flavus Reveals Novel Genes Involved in Aflatoxin Production.

Aspergillus flavus (A. flavus), a pathogenic fungus, can produce carcinogenic and toxic aflatoxins that are a serious agricultural and medical threat worldwide. Attempts to decipher the aflatoxin biosynthetic pathway have been hampered by the lack of a high-quality genome annotation for A. flavus. To address this gap, we performed a comprehensive proteogenomic analysis using high-accuracy mass spectrometry data for this pathogen. The resulting high-quality data set confirmed the translation of 8724 previously predicted genes and identified 732 novel proteins, 269 splice variants, 447 single amino acid variants, 188 revised genes. A subset of novel proteins was experimentally validated by RT-PCR and synthetic peptides. Further functional annotation suggested that a number of the identified novel proteins may play roles in aflatoxin biosynthesis and stress responses in A. flavus. This comprehensive strategy also identified a wide range of posttranslational modifications (PTMs), including 3461 modification sites from 1765 proteins. Functional analysis suggested the involvement of these modified proteins in the regulation of cellular metabolic and aflatoxin biosynthetic pathways. Together, we provided a high-quality annotation of A. flavus genome and revealed novel insights into the mechanisms of aflatoxin production and pathogenicity in this pathogen.

OsbHLH98 regulates leaf angle in rice through transcriptional repression of OsBUL1.

Leaf angle is an important agronomic trait in cereals that helps determine plant yield by affecting planting density. However, the regulation mechanism of leaf angle remained elusive. Here, we show that OsbHLH98, a rice bHLH transcription factor, negatively regulates leaf angle. osbhlh98 mutant leaves formed a larger leaf angle, whereas transgenic plants overexpressing OsbHLH98 exhibited a slight reduction in leaf angle. We determined that the changes in leaf angle resulted from increased number and size of parenchyma cells on the adaxial side of the lamina joint in osbhlh98 mutants. Experiments using reporter constructs showed that OsbHLH98 is expressed on the adaxial side of lamina joints, consistent with its proposed function in regulating leaf angle. Furthermore, we established by chromatin immunoprecipitation and CUT&RUN that OsBUL1 is a direct downstream target of OsbHLH98. Transactivation assays and reverse transcription quantitative polymerase chain reaction (RT-qPCR) analysis indicated that OsbHLH98 represses OsBUL1 transcription. Our results demonstrate that OsbHLH98 negatively regulates leaf angle by counteracting brassinosteroid-induced cell elongation via the repression of OsBUL1 transcription. The characterization of OsbHLH98 and its role in determining leaf angle will lay the foundation to develop the ideal plant architecture for adaptation to high planting density.

Antioxidant-related catalase CTA1 regulates development, aflatoxin biosynthesis, and virulence in pathogenic fungus Aspergillus flavus.

Reactive oxygen species (ROS) induce the synthesis of a myriad of secondary metabolites, including aflatoxins. It raises significant concern as it is a potent environmental contaminant. In Aspergillus flavus., antioxidant enzymes link ROS stress response with coordinated gene regulation of aflatoxin biosynthesis. In this study, we characterized the function of a core component of the antioxidant enzyme catalase (CTA1) of A. flavus. Firstly, we verified the presence of cta1 corresponding protein (CTA1) by Western blot analysis and mass-spectrometry based analysis. Then, the functional study revealed that the growth, sporulation and sclerotia formation significantly increased, while aflatoxins production and virulence were decreased in the cta1 deletion mutant as compared with the WT and complementary strains. Furthermore, the absence of the cta1 gene resulted in a significant rise in the intracellular ROS level, which in turn added to the oxidative stress level of cells. A further quantitative proteomics investigation hinted that in vivo, CTA1 might maintain the ROS level to facilitate the aflatoxin synthesis. All in all, the pleiotropic phenotype of A. flavus CTA1 deletion mutant revealed that the antioxidant system plays a crucial role in fungal development, aflatoxins biosynthesis and virulence.

Toward a next-generation atlas of RNA secondary structure.

RNA structure plays a crucial role in gene maturation, regulation and function. Determining the form and frequency of RNA folds is essential for a better understanding of how RNA exerts its functions. Low-throughput studies have focused on RNA primary sequences and expression levels, but with an emphasis on relatively small numbers of transcripts. However, with the recent advent of high-throughput technologies, it is realistic to begin analyzing RNA secondary structures on a genome-wide scale. Here, we review genome-wide RNA secondary structure profiles as well as advances in computational structure predictions. We further discuss the novel characteristics of RNA secondary structure across messenger RNAs. Probing RNA secondary structure by high-throughput sequencing will enable us to build atlases of RNA secondary structures, an important step in helping us to understand the versatility of RNA functions in diverse cellular processes.

Genome-wide characterization and evolution analysis of miniature inverted-repeat transposable elements (MITEs) in moso bamboo (Phyllostachys heterocycla).

MAIN CONCLUSION: Moso bamboo MITEs were genome-wide identified first time, and data shows that MITEs contribute to the genomic diversity and differentiation of bamboo. Miniature inverted-repeat transposable elements (MITEs) are widespread in animals and plants. There are a large number of transposable elements in moso bamboo (Phyllostachys heterocycla var. pubescens) genome, but the genome-wide information of moso bamboo MITEs is not known yet. Here we identified 362 MITE families with a total of 489,592 MITE-related sequences, accounting for 4.74 % of the moso bamboo genome. The 362 MITE families are clustered into six known and one unknown super-families. Our analysis indicated that moso bamboo MITEs preferred to reside in or near the genes that might be involved in regulation of host gene expression. Of the seven super-families, three might undergo major expansion event twice, respectively, during 8-11 million years ago (mya) ago and 22-28 mya ago; two might experience a long expansion period from 6 to 13 mya. Almost 1/3 small RNAs might be derived from the MITE sequences. Some MITE families generate small RNAs mainly from the terminals, while others predominantly from the central region. Given the high copy number of MITEs, many siRNAs and miRNAs derived from MITE sequences and the preferential insertion of MITE into gene regions, MITEs may contribute to the genomic diversity and differentiation of bamboo.

Identification and Analysis of Medicago truncatula Auxin Transporter Gene Families Uncover their Roles in Responses to Sinorhizobium meliloti Infection.

Auxin transport plays a pivotal role in the interaction between legume species and nitrogen-fixing bacteria to form symbioses. Auxin influx carriers auxin resistant 1/like aux 1 (AUX/LAX), efflux carriers pin-formed (PIN) and efflux/conditional P-glycoprotein (PGP/ABCB) are three major protein families participating in auxin polar transport. We used the latest Medicago truncatula genome sequence to characterize and analyze the M. truncatula LAX (MtLAX), M. truncatula PIN (MtPIN) and M. truncatula ABCB (MtABCB) families. Transient expression experiments indicated that three representative auxin transporters (MtLAX3, MtPIN7 and MtABCB1) showed cell plasma membrane localizations. The expression of most MtLAX, MtPIN and MtABCB genes was up-regulated in the roots and was down-regulated in the shoots by Sinorhizobium meliloti infection in the wild type (WT). However, the expression of these genes was down-regulated in both the roots and shoots of an infection-resistant mutant, dmi3. The different expression patterns between the WT and the mutant roots indicated that auxin relocation may be involved in rhizobial infection responses. Furthermore, IAA contents were significantly up-regulated in the shoots and down-regulated in the roots after Sinorhizobium meliloti infection in the WT. Inoculation of roots with rhizobia may reduce the auxin loading from shoots to roots by inhibiting the expression of most auxin transporter genes. However, the rate of change of gene expression and IAA contents in the dmi3 mutant were obviously lower than in the WT. The identification and expression analysis of auxin transporter genes helps us to understand the roles of auxin in the regulation of nodule formation in M. truncatula.

sRNA profiling in Aspergillus flavus reveals differentially expressed miRNA-like RNAs response to water activity and temperature.

Small non-coding RNA (sRNA) in various organisms remains a mysterious subject. Although microRNAs (miRNAs) have been intensively investigated in plants and animals, the study of miRNAs in fungi has been limited. Only microRNA-like RNAs (milRNAs) have been reported in several filamentous fungi. In this study, Illumina deep sequencing was performed to characterize the sRNA in Aspergillus flavus and to evaluate their responses to water activity and temperature. Global expression analysis showed an extensively differential expression of sRNA loci in A. flavus under different temperature or water activities. In addition, a total of 135 milRNAs were identified in A. flavus. The milRNA profiles obtained in deep sequencing were further validated by RT-qPCR assay. The presence and differential expression of milRNAs under different temperature or water activities in A. flavus imply that milRNAs might play important roles in the mycotoxin biosynthesis and mycelium growth in fungi A. flavus.

Pan-cancer Analysis Reveals Cancer-dependent Expression of SOX17 and Associated Clinical Outcomes.

BACKGROUND/AIM: SRY-box containing gene 17 (SOX17) plays a pivotal role in cancer onset and progression and is considered a potential target for cancer diagnosis and treatment. However, the expression pattern of SOX17 in cancer and its clinical relevance remains unknown. Here, we explored the relationship between the expression of SOX17 and drug response by examining SOX17 expression patterns across multiple cancer types. MATERIALS AND METHODS: Single-cell and bulk RNA-seq analyses were used to explore the expression profile of SOX17. Analysis results were verified with qPCR and immunohistochemistry. Survival, drug response, and co-expression analyses were performed to illustrate its correlation with clinical outcomes. RESULTS: The results revealed that abnormal expression of SOX17 is highly heterogenous across multiple cancer types, indicating that SOX17 manifests as a cancer type-dependent feature. Furthermore, the expression pattern of SOX17 is also associated with cancer prognosis in certain cancer types. Strong SOX17 expression correlates with the potency of small molecule drugs that affect PI3K/mTOR signaling. FGF18, a gene highly relevant to SOX17, is involved in the PI3K-AKT signaling pathway. Single-cell RNA-seq analysis demonstrated that SOX17 is mainly expressed in endothelial cells and barely expressed in other cells but spreads to other cell types during the development of ovarian cancer. CONCLUSION: Our study revealed the expression pattern of SOX17 in pan-cancer through bulk and single-cell RNA-seq analyses and determined that SOX17 is related to the diagnosis, staging, and prognosis of some tumors. These findings have clinical implications and may help identify mechanistic pathways amenable to pharmacological interventions.

Molecular cloning and characterization of OsCHR4, a rice chromatin-remodeling factor required for early chloroplast development in adaxial mesophyll.

Mi-2 protein, the central component of the NuRD nucleosome remodeling and histone deacetylase complex, plays a role in transcriptional repression in animals. Mi-2-like genes have been reported in Arabidopsis, though their function in monocots remains largely unknown. In the present study, a rice Mi-2-like gene, OsCHR4 (Oryza sativa Chromatin Remodeling 4, LOC_Os07g03450), was cloned from a rice mutant with adaxial albino leaves. The Oschr4 mutant exhibited defective chloroplasts in adaxial mesophyll, but not in abaxial mesophyll. Ultrastructural observations indicated that proplastid growth and/or thylakoid membrane formation in adaxial mesophyll cells was blocked in the Oschr4 mutant. Subcellular localization revealed that OsCHR4::GFP fusion protein was targeted to the nuclei. OsCHR4 was mainly expressed in the root meristem, flower, vascular bundle, and mesophyll cells by promoter::GUS analysis in transgenic rice. The transcripts of some nuclear- and plastid-encoded genes required for early chloroplast development and photosynthesis were decreased in the adaxial albino mesophyll of the Oschr4 mutant. These observations provide evidence that OsCHR4, the rice Mi-2-like protein, plays an important role in early chloroplast development in adaxial mesophyll cells. The results increase our understanding of the molecular mechanism underlying tissue-specific chloroplast development in plants.

MyBioNet: interactively visualize, edit and merge biological networks on the web.

SUMMARY: MyBioNet is a web-based application for biological network analysis, which provides user-friendly web interfaces to visualize, edit and merge biological networks. In addition, MyBioNet integrated KEGG metabolic network data from 1366 organisms and allows users to search and navigate interesting networks. AVAILABILITY AND IMPLEMENTATION: All KEGG metabolic network data are organized and stored in the MySQL database. MyBioNet is implemented in Flex/Actionscript and PHP languages and deployed on an Apache web server. MyBioNet is accessible through all the Flash-embedded browsers at http://bis.zju.edu.cn/mybionet/. CONTACT: mchen@zju.edu.cn.

OsPHF1 regulates the plasma membrane localization of low- and high-affinity inorganic phosphate transporters and determines inorganic phosphate uptake and translocation in rice.

PHOSPHATE TRANSPORTER TRAFFIC FACILITATOR1 (PHF1) is known to regulate the plasma membrane localization of PHT1;1, a high-affinity inorganic phosphate (Pi) transporter in Arabidopsis (Arabidopsis thaliana). OsPHF1, a rice (Oryza sativa) gene homologous to AtPHF1, was isolated and found to regulate the localization of both low- and high-affinity Pi transporters to the plasma membrane. Three OsPHF1 allelic mutants carrying one-point mutations at the fifth WD-repeat motif and two at the transmembrane helix, respectively, showed arsenate resistance and severely reduced Pi accumulation. The data indicate that mutation of OsPHF1 results in the endoplasmic reticulum retention of the low-affinity Pi transporter OsPT2 and high-affinity Pi transporter OsPT8. Mutation of OsPHF1 also reduced Pi accumulation in plants exhibiting excessive shoot Pi accumulation due to the overexpression of OsPHR2. However, the transcript level of OsPHF1 itself is not controlled by OsPHR2. Overexpression of OsPHF1 increased Pi accumulation in both roots and shoots in a solution culture with Pi-supplied condition. These results indicate that the role of OsPHF1 is unique in the localization of both low- and high-affinity Pi transporters on the plasma membrane in rice and determines Pi uptake and translocation in rice. The similar function of PHF1 required to facilitate PHT1 transit through the endoplasmic reticulum between Arabidopsis and rice provides an example of expectations from what one would deduce from sequence comparisons to extend knowledge from Arabidopsis to crops.

Origin and evolutionary analysis of the plant-specific TIFY transcription factor family.

A substantial number of transcription factor families have been investigated from all kingdoms of life, but a particular class of plant-specific TIFY transcription factors, characterized by a highly conserved TIFY domain, lacks a systemic analysis of its origin and evolutionary relationships among different plant species. After exhaustive genome-wide searches against 14 genomes, TIFY transcription factors were identified and classified into four subfamilies TIFY, PPD, JAZ and ZML according to their different domain architectures. Results show that the TIFY domain of the ZML subfamily possesses a core "TLS[F/Y]XG" motif rather than the "TIFYXG" motif that is dominant in the other three subfamilies. A comprehensive survey of the TIFY family allowed us to discover a new group within the JAZ subfamily and to identify several novel conserved motifs via phylogenetic analysis. Evolutional analysis indicates that whole genome duplication and tandem duplication contributed to the expansion of the TIFY family in plants.

Single-Cell RNA-Sequencing Atlas Reveals the Tumor Microenvironment of Metastatic High-Grade Serous Ovarian Carcinoma.

Ovarian cancer is the most common and lethal gynecological tumor in women worldwide. High-grade serous ovarian carcinoma (HGSOC) is one of the histological subtypes of epithelial ovarian cancer, accounting for 70%. It often occurs at later stages associated with a more fatal prognosis than endometrioid carcinomas (EC), another subtype of epithelial ovarian cancer. However, the molecular mechanism and biology underlying the metastatic HGSOC (HG_M) immunophenotype remain poorly elusive. Here, we performed single-cell RNA sequencing analyses of primary HGSOC (HG_P) samples, metastatic HGSOC (HG_M) samples, and endometrioid carcinomas (EC) samples. We found that ERBB2 and HOXB-AS3 genes were more amplified in metastasis tumors than in primary tumors. Notably, high-grade serous ovarian cancer metastases are accompanied by dysregulation of multiple pathways. Malignant cells with features of epithelial-mesenchymal transition (EMT) affiliated with poor overall survival were identified. In addition, cancer-associated fibroblasts with EMT-program were enriched in HG_M, participating in angiogenesis and immune regulation, such as IL6/STAT3 pathway activity. Compared with ECs, HGSOCs exhibited higher T cell infiltration. PRDM1 regulators may be involved in T cell exhaustion in ovarian cancer. The CX3CR1_macro subpopulation may play a role in promoting tumor progression in ovarian cancer with high expression of BAG3, IL1B, and VEGFA. The new targets we discovered in this study will be useful in the future, providing guidance on the treatment of ovarian cancer.

Small RNAs in angiosperms: sequence characteristics, distribution and generation.

High-throughput sequencing (HTS) has opened up a new era for small RNA (sRNA) exploration. Using HTS data for a global survey of sRNAs in 26 angiosperms, elevated GC contents were detected in the monocots, whereas the 5(')-terminal compositions were quite uniform among the angiosperms. Chromosome-wide distribution patterns of sRNAs were investigated by using scrolling-window analysis. We performed de novo natural antisense transcript (NAT) prediction, and found that the overlapping regions of trans-NATs, but not cis-NATs, were hotspots for sRNA generation. One cis-NAT generates phased natural antisense short interfering RNAs (nat-siRNAs) specifically from flowers in Arabidopsis, while one in rice produces phased nat-siRNAs from grains, suggesting their organ-specific regulatory roles.

A novel phosphoinositide kinase Fab1 regulates biosynthesis of pathogenic aflatoxin in Aspergillus flavus.

Aspergillus flavus (A. flavus) is one of the most important model environmental fungi which can produce a potent toxin and carcinogen known as aflatoxin. Aflatoxin contamination causes massive agricultural economic loss and a critical human health issue each year. Although a functional vacuole has been highlighted for its fundamental importance in fungal virulence, the molecular mechanisms of the vacuole in regulating the virulence of A. flavus remain largely unknown. Here, we identified a novel vacuole-related protein in A. flavus, the ortholog of phosphatidylinositol-3-phosphate-5-kinase (Fab1) in Saccharomyces cerevisiae. This kinase was located at the vacuolar membrane, and loss of fab1 function was found to affect the growth, conidia and sclerotial development, cellular acidification and metal ion homeostasis, aflatoxin production and pathogenicity of A. flavus. Further functional analysis revealed that Fab1 was required to maintain the vacuole size and cell morphology. Additional quantitative proteomic analysis suggested that Fab1 was likely to play an important role in maintaining vacuolar/cellular homeostasis, with vacuolar dysregulation upon fab1 deletion leading to impaired aflatoxin synthesis in this fungus. Together, these results provide insight into the molecular mechanisms by which this pathogen produces aflatoxin and mediates its pathogenicity, and may facilitate dissection of the vacuole-mediated regulatory network in A. flavus.

Auxin-related gene families in abiotic stress response in Sorghum bicolor.

Sorghum, a C4 model plant, has been studied to develop an understanding of the molecular mechanism of resistance to stress. The auxin-response genes, auxin/indole-3-acetic acid (Aux/IAA), auxin-response factor (ARF), Gretchen Hagen3 (GH3), small auxin-up RNAs, and lateral organ boundaries (LBD), are involved in growth/development and stress/defense responses in Arabidopsis and rice, but they have not been studied in sorghum. In the present paper, the chromosome distribution, gene duplication, promoters, intron/exon, and phylogenic relationships of Aux/IAA, ARF, GH3, and LBD genes in sorghum are presented. Furthermore, real-time PCR analysis demonstrated these genes are differently expressed in leaf/root of sorghum and indicated the expression profile of these gene families under IAA, brassinosteroid (BR), salt, and drought treatments. The SbGH3 and SbLBD genes, expressed in low level under natural condition, were highly induced by salt and drought stress consistent with their products being involved in both abiotic stresses. Three genes, SbIAA1, SbGH3-13, and SbLBD32, were highly induced under all the four treatments, IAA, BR, salt, and drought. The analysis provided new evidence for role of auxin in stress response, implied there are cross talk between auxin, BR and abiotic stress signaling pathways.

Functional analysis of the structural domain of ARF proteins in rice (Oryza sativa L.).

Auxin response factors (ARFs) are key regulators of plant growth and development. Through interaction with auxin/indole acetic acid (Aux/IAA) proteins, they influence the expression of auxin response genes. An ARF gene family has been predicted in rice, but the functions of the individual structural domains of the OsARFs remain obscure. Bioinformatics was used to analyse the position of the DNA-binding domain (DBD), middle region (MR), and C-terminal dimerization domain (CTD) of OsARFs, and experimentally confirmed the presence of a classical monopartite nuclear localization signal (NLS) in the DBD. The DBD was shown to contribute to nuclear localization of OsARF proteins in addition to its known DNA-binding function. Interactions between 14 integrated OsARFs and 15 OsIAA proteins were tested using yeast two-hybrid assays. It was found that eight OsARF activators interacted with the 15 OsIAA proteins, while six OsARF repressors did not. The interactions between the MR+CTD or CTD of 10 OsARFs and 15 OsIAA proteins were also tested and the results were consistent with those of each intact OsARF, although some slight differences in interaction intensity were observed by α-galactosidase quantitative assays. The truncated CTD of OsARF11 did not interact with any OsIAA, implying that the CTD is required for ARF-IAA dimerization, and that the MR influences the interaction intensity in yeast. A subset of the interactions in yeast were also observed in tobacco plants using firefly luciferase complementation imaging assays, indicating that these interactions are specific in plants, and might have a special role in the auxin signalling response. This study provides new insight into the structure of OsARF proteins and ARF-Aux/IAA interactions.

Identification and characterization of pineapple leaf lncRNAs in crassulacean acid metabolism (CAM) photosynthesis pathway.

Long noncoding RNAs (lncRNAs) have been identified in many mammals and plants and are known to play crucial roles in multiple biological processes. Pineapple is an important tropical fruit and a good model for studying the plant evolutionary adaptation to the dry environment and the crassulacean acid metabolism (CAM) photosynthesis strategy; however, the lncRNAs involved in CAM pathway remain poorly characterized. Here, we analyzed the available RNA-seq data sets derived from 26 pineapple leaf samples at 13 time points and identified 2,888 leaf lncRNAs, including 2,046 long intergenic noncoding RNAs (lincRNAs) and 842 long noncoding natural antisense transcripts (lncNATs). Pineapple leaf lncRNAs are expressed in a highly tissue-specific manner. Co-expression analysis of leaf lncRNA and mRNA revealed that leaf lncRNAs are preferentially associated with photosynthesis genes. We further identified leaf lncRNAs that potentially function as competing endogenous RNAs (ceRNAs) of two CAM photosynthesis pathway genes, PPCK and PEPC, and revealed their diurnal expression pattern in leaves. Moreover, we found that 48% of lncRNAs exhibit diurnal expression patterns in leaves, suggesting their important roles in CAM. This study conducted a comprehensive genome-wide analysis of leaf lncRNAs and identified their role in gene expression regulation of the CAM photosynthesis pathway in pineapple.