Search

Background Recent studies using single molecule, real-time (SMRT) sequencing revealed a substantial population of analyzable long cell-free DNA (cfDNA) in plasma. Potential clinical utilities of such long cfDNA in pregnancy and cancer have been demonstrated. However, the performance of different long-read sequencing platforms for the analysis of long cfDNA remains unknown. Methods Size biases of SMRT sequencing by Pacific Biosciences (PacBio) and nanopore sequencing by Oxford Nanopore Technologies (ONT) were evaluated using artificial mixtures of sonicated human and mouse DNA of different sizes. cfDNA from plasma samples of pregnant women at different trimesters, hepatitis B carriers, and patients with hepatocellular carcinoma were sequenced with the 2 platforms. Results Both platforms showed biases to sequence longer (1500 bp vs 200 bp) DNA fragments, with PacBio showing a stronger bias (5-fold overrepresentation of long fragments vs 2-fold in ONT). Percentages of cfDNA fragments 500 bp were around 6-fold higher in PacBio compared with ONT. End motif profiles of cfDNA from PacBio and ONT were similar, yet exhibited platform-dependent patterns. Tissue-of-origin analysis based on single-molecule methylation patterns showed comparable performance on both platforms. Conclusions SMRT sequencing generated data with higher percentages of long cfDNA compared with nanopore sequencing. Yet, a higher number of long cfDNA fragments eligible for the tissue-of-origin analysis could be obtained from nanopore sequencing due to its much higher throughput. When analyzing the size and end motif of cfDNA, one should be aware of the analytical characteristics and possible biases of the sequencing platforms being used.

Obesity represents a major health challenge because it substantially increases the risk of metabolic diseases. Capsaicin, the major active ingredient of Capsicum spp., has been reported to possess anti-obesity activity. Hereon, the effect of capsaicin on glucose uptake and consumption in hepatocytes was extensively studied. Capsaicin was shown to accelerate the glucose uptake/consumption and the ATP production of hepatocytes. The elevation of intracellular Ca2+ was thought to be a potential mechanism. By transcriptome analysis, 78, 146 and 507 differentially expressed genes (DEGs) were identified between capsaicin and the control group for 4 h, 12 h and 24 h treatments. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that most of the DEGs were involved in canonical pathways, like MAPK and PI3K-AKT signaling pathways. Clustering analysis showed that many DEGs were associated with glucose and amino acid metabolism. The variation trend in genes related to glucose and amino acid metabolism (like CTH, VEGFA, PCK2 and IGFBP3) in the quantitative PCR (q-PCR) assay was consistent with the transcriptome data. These results demonstrated that capsaicin efficiently accelerated the glucose uptake and consumption of hepatocytes.

Background Jagged ends of plasma DNA are a recently recognized class of fragmentomic markers for cell-free DNA, reflecting the activity of nucleases. A number of recent studies have also highlighted the importance of jagged ends in the context of pregnancy and oncology. However, knowledge regarding the generation of jagged ends is incomplete. Methods Jaggedness of plasma DNA was analyzed based on Jag-seq, which utilized the differential methylation signals introduced by the DNA end-repair process. We investigated the jagged ends in plasma DNA using mouse models by deleting the deoxyribonuclease 1 (Dnase1), DNA fragmentation factor subunit beta (Dffb), or deoxyribonuclease 1 like 3 (Dnase1l3) gene. Results Aberrations in the profile of plasma DNA jagged ends correlated with the type of nuclease that had been genetically deleted, depending on nucleosomal structures. The deletion of Dnase1l3 led to a significant reduction of jaggedness for those plasma DNA molecules involving more than 1 nucleosome (e.g., size ranges 240-290 bp, 330-380 bp, and 420-470 bp). However, less significant effects of Dnase1 and Dffb deletions were observed regarding different sizes of DNA fragments. Interestingly, the aberration in plasma DNA jagged ends related to multinucleosomes was observed in human subjects with familial systemic lupus erythematosus with Dnase1l3 deficiency and human subjects with sporadic systemic lupus erythematosus. Conclusions Detailed understanding of the relationship between nuclease and plasma DNA jaggedness has opened up avenues for biomarker development.

Cell-free DNA (cfDNA) fragmentation is nonrandom, at least partially mediated by various DNA nucleases, forming characteristic cfDNA end motifs. However, there is a paucity of tools for deciphering the relative contributions of cfDNA cleavage patterns related to underlying fragmentation factors. In this study, through non-negative matrix factorization algorithm, we used 256 5′ 4-mer end motifs to identify distinct types of cfDNA cleavage patterns, referred to as “founder” end-motif profiles (F-profiles). F-profiles were associated with different DNA nucleases based on whether such patterns were disrupted in nuclease-knockout mouse models. Contributions of individual F-profiles in a cfDNA sample could be determined by deconvolutional analysis. We analyzed 93 murine cfDNA samples of different nuclease-deficient mice and identified six types of F-profiles. F-profiles I, II, and III were linked to deoxyribonuclease 1 like 3 (DNASE1L3), deoxyribonuclease 1 (DNASE1), and DNA fragmentation factor subunit beta (DFFB), respectively. We revealed that 42.9% of plasma cfDNA molecules were attributed to DNASE1L3-mediated fragmentation, whereas 43.4% of urinary cfDNA molecules involved DNASE1-mediated fragmentation. We further demonstrated that the relative contributions of F-profiles were useful to inform pathological states, such as autoimmune disorders and cancer. Among the six F-profiles, the use of F-profile I could inform the human patients with systemic lupus erythematosus. F-profile VI could be used to detect individuals with hepatocellular carcinoma, with an area under the receiver operating characteristic curve of 0.97. F-profile VI was more prominent in patients with nasopharyngeal carcinoma undergoing chemoradiotherapy. We proposed that this profile might be related to oxidative stress.

Cell-free DNA (cfDNA) fragmentation patterns contain important molecular information linked to tissues of origin. We explored the possibility of using fragmentation patterns to predict cytosine-phosphate-guanine (CpG) methylation of cfDNA, obviating the use of bisulfite treatment and associated risks of DNA degradation. This study investigated the cfDNA cleavage profile surrounding a CpG (i.e., within an 11-nucleotide [nt] window) to analyze cfDNA methylation. The cfDNA cleavage proportion across positions within the window appeared nonrandom and exhibited correlation with methylation status. The mean cleavage proportion was ∼twofold higher at the cytosine of methylated CpGs than unmethylated ones in healthy controls. In contrast, the mean cleavage proportion rapidly decreased at the 1-nt position immediately preceding methylated CpGs. Such differential cleavages resulted in a characteristic change in relative presentations of CGN and NCG motifs at 5′ ends, where N represented any nucleotide. CGN/NCG motif ratios were correlated with methylation levels at tissue-specific methylated CpGs (e.g., placenta or liver) (Pearson’s absolute r > 0.86). cfDNA cleavage profiles were thus informative for cfDNA methylation and tissue-of-origin analyses. Using CG-containing end motifs, we achieved an area under a receiver operating characteristic curve (AUC) of 0.98 in differentiating patients with and without hepatocellular carcinoma and enhanced the positive predictive value of nasopharyngeal carcinoma screening (from 19.6 to 26.8%). Furthermore, we elucidated the feasibility of using cfDNA cleavage patterns to deduce CpG methylation at single CpG resolution using a deep learning algorithm and achieved an AUC of 0.93. FRAGmentomics-based Methylation Analysis (FRAGMA) presents many possibilities for noninvasive prenatal, cancer, and organ transplantation assessment.

BackgroundDouble-stranded DNA in plasma is known to carry single-stranded ends, called jagged ends. Plasma DNA jagged ends are biomarkers for pathophysiologic states such as pregnancy and cancer. It remains unknown whether urinary cell-free DNA (cfDNA) molecules have jagged ends.MethodsJagged ends of cfDNA were detected by incorporating unmethylated cytosines during a DNA end-repair process, followed by bisulfite sequencing. Incorporation of unmethylated cytosines during the repair of the jagged ends lowered the apparent methylation levels measured by bisulfite sequencing and were used to calculate a jagged end index. This approach is called jagged end analysis by sequencing.ResultsThe jagged end index of urinary cfDNA was higher than that of plasma DNA. The jagged end index profile of plasma DNA displayed several strongly oscillating major peaks at intervals of approximately 165 bp (i.e., nucleosome size) and weakly oscillating minor peaks with periodicities of approximately 10 bp. In contrast, the urinary DNA jagged end index profile showed weakly oscillating major peaks but strongly oscillating minor peaks. The jagged end index was generally higher in nucleosomal linker DNA regions. Patients with bladder cancer (n = 46) had lower jagged end indexed of urinary DNA than participants without bladder cancer (n = 39). The area under the curve for differentiating between patients with and without bladder cancer was 0.83.ConclusionsJagged ends represent a property of urinary cfDNA. The generation of jagged ends might be related to nucleosomal structures, with enrichment in linker DNA regions. Jagged ends of urinary DNA could potentially serve as a new biomarker for bladder cancer detection.

Plasma DNA fragmentomics is an emerging area of research covering plasma DNA sizes, end points, and nucleosome footprints. In the present study, we found a significant increase in the diversity of plasma DNA end motifs in patients with hepatocellular carcinoma (HCC). Compared with patients without HCC, patients with HCC showed a preferential pattern of 4-mer end motifs. In particular, the abundance of plasma DNA motif CCCA was much lower in patients with HCC than in subjects without HCC. The aberrant end motifs were also observed in patients with other cancer types, including colorectal cancer, lung cancer, nasopharyngeal carcinoma, and head and neck squamous cell carcinoma. We further observed that the profile of plasma DNA end motifs originating from the same organ, such as the liver, placenta, and hematopoietic cells, generally clustered together. The profile of end motifs may therefore serve as a class of biomarkers for liquid biopsy in oncology, noninvasive prenatal testing, and transplantation monitoring. Significance: Plasma DNA molecules originating from the liver, HCC and other cancers, placenta, and hematopoietic cells each harbor a set of characteristic plasma DNA end motifs. Such markers carry tissue-of-origin information and represent a new class of biomarkers in the nascent field of fragmentomics. This article is highlighted in the In This Issue feature, p. 627

Background Analysis of circulating tumor DNA has become increasingly important as a tool for cancer care. However, the focus of previous studies has been on short fragments of DNA. Also, bisulfite sequencing, a conventional approach for methylation analysis, causes DNA degradation, which is not ideal for the assessment of long DNA properties and methylation patterns. This study attempted to overcome such obstacles by single-molecule sequencing. Methods Single-molecule real-time (SMRT) sequencing was used to sequence plasma DNA. We performed fragment size and direct methylation analysis for each molecule. A methylation score concerning single-molecule methylation patterns was used for cancer detection. Results A substantial proportion of plasma DNA was longer than 1 kb with a median of 16% in hepatocellular carcinoma (HCC) patients, hepatitis B virus carriers, and healthy individuals. The longest plasma DNA molecule in the HCC patients was 39.8 kb. Tumoral cell-free DNA (cfDNA) was generally shorter than nontumoral cfDNA. The longest tumoral cfDNA was 13.6 kb. Tumoral cfDNA had lower methylation levels compared with nontumoral cfDNA (median: 59.3% vs 76.9%). We developed and analyzed a metric reflecting single-molecule methylation patterns associated with cancer, named the HCC methylation score. HCC patients displayed significantly higher HCC methylation scores than those without HCC. Interestingly, compared to using short cfDNA (area under the receiver operating characteristic [ROC] curve, AUC: 0.75), the use of long cfDNA molecules greatly enhanced the discriminatory power (AUC: 0.91). Conclusions A previously unidentified long cfDNA population was revealed in cancer patients. The presence and direct methylation analysis of these molecules open new possibilities for cancer liquid biopsy.

Significance We revealed a large population of long cell-free DNA molecules (up to 23,635 bp in length) in maternal plasma and developed an approach which leveraged the abundance of CpG sites on long molecules to deduce the tissue of origin of individual plasma DNA molecules based on single-molecule methylation analysis. We illustrated how such an approach may be utilized to achieve noninvasive prenatal testing of monogenic diseases. We also revealed a reduction in amounts of such long cell-free DNA molecules and a different end motif profile in maternal plasma DNA from pregnancies with preeclampsia. Hence, long cell-free DNA molecules represent a valuable resource of biomarker development for pregnancy-associated disorders., In the field of circulating cell-free DNA, most of the studies have focused on short DNA molecules (e.g.

Significance Single molecule real-time (SMRT) sequencing theoretically offers the opportunity to directly assess certain base modifications of native DNA molecules without any prior chemical/enzymatic conversions and PCR amplification, using kinetic signals of a DNA polymerase. However, the kinetic signal changes caused by 5mC modification are extremely subtle. Hence, the robust genome-wide measurement of 5mC modification has not been achieved. We enhanced 5mC detection using SMRT sequencing by holistically analyzing kinetic signals of a DNA polymerase and sequence context for every base within a measurement window. We employed a convolutional neural network to train a methylation classification model, leading to genome-wide 5mC detection. The sensitivity and specificity reached 90% and 94%, with a 99% correlation of overall methylation level with bisulfite sequencing., 5-Methylcytosine (5mC) is an important type of epigenetic modification. Bisulfite sequencing (BS-seq) has limitations, such as severe DNA degradation. Using single molecule real-time sequencing, we developed a methodology to directly examine 5mC. This approach holistically examined kinetic signals of a DNA polymerase (including interpulse duration and pulse width) and sequence context for every nucleotide within a measurement window, termed the holistic kinetic (HK) model. The measurement window of each analyzed double-stranded DNA molecule comprised 21 nucleotides with a cytosine in a CpG site in the center. We used amplified DNA (unmethylated) and M.SssI-treated DNA (methylated) (M.SssI being a CpG methyltransferase) to train a convolutional neural network. The area under the curve for differentiating methylation states using such samples was up to 0.97. The sensitivity and specificity for genome-wide 5mC detection at single-base resolution reached 90% and 94%, respectively. The HK model was then tested on human–mouse hybrid fragments in which each member of the hybrid had a different methylation status. The model was also tested on human genomic DNA molecules extracted from various biological samples, such as buffy coat, placental, and tumoral tissues. The overall methylation levels deduced by the HK model were well correlated with those by BS-seq (r = 0.99; P < 0.0001) and allowed the measurement of allele-specific methylation patterns in imprinted genes. Taken together, this methodology has provided a system for simultaneous genome-wide genetic and epigenetic analyses.

Maize was domesticated from lowland teosinte (Zea mays ssp. parviglumis), but the contribution of highland teosinte (Zea mays ssp. mexicana, hereafter mexicana) to modern maize is not clear. Here, two genomes for Mo17 (a modern maize inbred) and mexicana are assembled using a meta-assembly strategy after sequencing of 10 lines derived from a maize-teosinte cross. Comparative analyses reveal a high level of diversity between Mo17, B73, and mexicana, including three Mb-size structural rearrangements. The maize spontaneous mutation rate is estimated to be 2.17 × 10−8 ~3.87 × 10−8 per site per generation with a nonrandom distribution across the genome. A higher deleterious mutation rate is observed in the pericentromeric regions, and might be caused by differences in recombination frequency. Over 10% of the maize genome shows evidence of introgression from the mexicana genome, suggesting that mexicana contributed to maize adaptation and improvement. Our data offer a rich resource for constructing the pan-genome of Zea mays and genetic improvement of modern maize varieties., Maize was domesticated from wild lowland progenitors that co-existed with upland subspecies in Southwestern Mexico. Here Yang et al. use a meta-assembly approach to assemble an upland mexicana genome and find evidence of introgression suggesting it contributed to modern maize adaptation

International audience; Plasma DNA fragmentomics is an emerging area in cell-free DNA diagnostics and research. In murine models, it has been shown that the extracellular DNase, DNASE1L3, plays a role in the fragmentation of plasma DNA. In humans, DNASE1L3 deficiency causes familial monogenic systemic lupus erythematosus with childhood onset and anti-dsDNA reactivity. In this study, we found that human patients with DNASE1L3 disease-associated gene variations showed aberrations in size and a reduction of a ''CC'' end motif of plasma DNA. Furthermore, we demonstrated that DNA from DNASE1L3-digested cell nuclei showed a median length of 153 bp with CC motif frequencies resembling plasma DNA from healthy individuals. Adeno-associated virus-based transduction of Dnase1l3 into Dnase1l3-deficient mice restored the end motif profiles to those seen in the plasma DNA of wild-type mice. Our findings demonstrate that DNASE1L3 is an important player in the fragmentation of plasma DNA, which appears to act in a cell-extrinsic manner to regulate plasma DNA size and motif frequency.

BACKGROUND Cellular mitochondrial DNA (mtDNA) is organized as circular, covalently closed and double-stranded DNA. Studies have demonstrated the presence of short mtDNA fragments in plasma. It is not known whether circular mtDNA might concurrently exist with linear mtDNA in plasma. METHODS We elucidated the topology of plasma mtDNA using restriction enzyme BfaI cleavage signatures on mtDNA fragment ends to differentiate linear and circular mtDNA. mtDNA fragments with both ends carrying BfaI cleavage signatures were defined as circular-derived mtDNA, whereas those with no cleavage signature or with 1 cleavage signature were defined as linear-derived mtDNA. An independent assay using exonuclease V to remove linear DNA followed by restriction enzyme MspI digestion was used for confirming the conclusions based on BfaI cleavage analysis. We analyzed the presence of BfaI cleavage signatures on plasma DNA ends in nonhematopoietically and hematopoietically derived DNA molecules by sequencing plasma DNA of patients with liver transplantation and bone marrow transplantation. RESULTS Both linear and circular mtDNA coexisted in plasma. In patients with liver transplantation, donor-derived (i.e., liver) mtDNA molecules were mainly linear (median fraction, 91%; range, 75%–97%), whereas recipient-derived (i.e., hematopoietic) mtDNA molecules were mainly circular (median fraction, 88%; range, 77%–93%). The proportion of linear mtDNA was well correlated with liver DNA contribution in the plasma DNA pool (r = 0.83; P value = 0.0008). Consistent data were obtained from a bone marrow transplantation recipient in whom the donor-derived (i.e., hematopoietic) mtDNA molecules were predominantly circular. CONCLUSIONS Linear and circular mtDNA molecules coexist in plasma and may have different tissue origins.

Significance Circulating DNA in plasma has many diagnostic applications, including noninvasive prenatal testing and cancer liquid biopsy. Plasma DNA consists of short fragments of DNA. However, there is little information about mechanisms that are involved in the fragmentation of plasma DNA. We showed that mice in which Dnase1l3 had been deleted showed aberrations in the fragmentation of plasma DNA. We also observed a change in the ranked frequencies of end motifs of plasma DNA caused by the Dnase1l3 deletion. In Dnase1l3−/− mice pregnant with Dnase1l3+/− fetuses, we observed a partial reversal of the plasma DNA aberrations. This study has thus linked the fields of nuclease biology and circulating nucleic acids and has opened up avenues for future research., Circulating DNA in plasma consists of short DNA fragments. The biological processes generating such fragments are not well understood. DNASE1L3 is a secreted DNASE1-like nuclease capable of digesting DNA in chromatin, and its absence causes anti-DNA responses and autoimmunity in humans and mice. We found that the deletion of Dnase1l3 in mice resulted in aberrations in the fragmentation of plasma DNA. Such aberrations included an increase in short DNA molecules below 120 bp, which was positively correlated with anti-DNA antibody levels. We also observed an increase in long, multinucleosomal DNA molecules and decreased frequencies of the most common end motifs found in plasma DNA. These aberrations were independent of anti-DNA response, suggesting that they represented a primary effect of DNASE1L3 loss. Pregnant Dnase1l3−/− mice carrying Dnase1l3+/− fetuses showed a partial restoration of normal frequencies of plasma DNA end motifs, suggesting that DNASE1L3 from Dnase1l3-proficient fetuses could enter maternal systemic circulation and affect both fetal and maternal DNA fragmentation in a systemic as well as local manner. However, the observed shortening of circulating fetal DNA relative to maternal DNA was not affected by the deletion of Dnase1l3. Collectively, our findings demonstrate that DNASE1L3 plays a role in circulating plasma DNA homeostasis by enhancing fragmentation and influencing end-motif frequencies. These results support a distinct role of DNASE1L3 as a regulator of the physical form and availability of cell-free DNA and may have important implications for the mechanism whereby this enzyme prevents autoimmunity.

Significance We identified differentiating molecular characteristics of plasma EBV DNA between nasopharyngeal carcinoma (NPC) patients and non-NPC subjects. Sequencing-based analysis revealed higher amounts of plasma EBV DNA and generally longer fragment lengths of plasma viral molecules in NPC patients than in non-NPC subjects. Based on these findings, we have developed a highly accurate blood-based test for screening of NPC. Such an approach is shown to enhance the positive predictive value and demonstrate a superior performance for NPC screening. It also obviates the need of a follow-up blood sample and therefore allows single time-point testing. We believe that this more clinically practical protocol would facilitate NPC screening on a population scale., Circulating tumor-derived DNA testing for cancer screening has recently been demonstrated in a prospective study on identification of nasopharyngeal carcinoma (NPC) among 20,174 asymptomatic individuals. Plasma EBV DNA, a marker for NPC, was detected using real-time PCR. While plasma EBV DNA was persistently detectable in 97.1% of the NPCs identified, ∼5% of the general population had transiently detectable plasma EBV DNA. We hypothesized that EBV DNA in plasma of subjects with or without NPC may have different molecular characteristics. We performed target-capture sequencing of plasma EBV DNA and identified differences in the abundance and size profiles of EBV DNA molecules within plasma of NPC and non-NPC subjects. NPC patients had significantly higher amounts of plasma EBV DNA, which showed longer fragment lengths. Cutoff values were established from an exploratory dataset and tested in a validation sample set. Adopting an algorithm that required a sample to concurrently pass cutoffs for EBV DNA counting and size measurements, NPCs were detected at a positive predictive value (PPV) of 19.6%. This represented superior performance compared with the PPV of 11.0% in the prospective screening study, which required participants with an initially detectable plasma EBV DNA result to be retested within 4 weeks. The observed differences in the molecular nature of EBV DNA molecules in plasma of subjects with or without NPC were successfully translated into a sequencing-based test that had a high PPV for NPC screening and achievable through single time-point testing.

An ideal plant height is essential for crop yield. Some Polycomb group (PcG) genes mutants exhibit a dwarf phenotype in rice. To determine how PcG genes regulate plant height, we investigated the phenotypes of the emf2b mutant and OsEMF2b, OsFIE2 and OsCLF RNAi transgenic plants; they all exhibited dwarf phenotype. Further analyses indicated that stem elongation at jointing stage was seriously inhibited in emf2b and RNAi transgenic plants. Reduced cell expansion and cell division of stem caused significant reduction of stem length during mature period of OsEMF2b, OsFIE2 and OsCLF RNAi transgenic plants. Transcription analysis revealed that cell division, cell expansion and plant hormones related genes differentially expressed between emf2b and WT. In addition, PcG genes mutants weakened GA signal and GA concentration and leaded to suppresseion of plant height. Analysis of differentially expressed genes revealed that 109 up-regulated and 19 down-regulated genes were identified in both emf2b and fie2. H3K27me3-modified sites were observed in 95 of the 109 up-regulated genes, and some of them were up-regulated in OsFIE2, OsCLF and OsEMF2b RNAi transgenic plants, and their H3K27me3 levels were reduced in emf2b. Moreover, OsEMF2b interacted with OsCLF. Therefore, we speculated that these PcG genes, OsFIE2, OsCLF and OsEMF2b, may work as a PRC2 to regulate rice height.

We isolated Pseudomonas putida (P. putida) strain 1A00316 from Antarctica. This bacterium has a high efficiency against Meloidogyne incognita (M. incognita) in vitro and under greenhouse conditions. The complete genome of P. putida 1A00316 was sequenced using PacBio single molecule real-time (SMRT) technology. A comparative genomic analysis of 16 Pseudomonas strains revealed that although P. putida 1A00316 belonged to P. putida, it was phenotypically more similar to nematicidal Pseudomonas fluorescens (P. fluorescens) strains. We characterized the diversity and specificity of nematicidal factors in P. putida 1A00316 with comparative genomics and functional analysis and found that P. putida 1A00316 has diverse nematicidal factors including protein alkaline metalloproteinase AprA and two secondary metabolites, hydrogen cyanide and cyclo-(l-isoleucyl-l-proline). We show for the first time that cyclo-(l-isoleucyl-l-proline) exhibit nematicidal activity in P. putida. Interestingly, our study had not detected common nematicidal factors such as 2,4-diacetylphloroglucinol (2,4-DAPG) and pyrrolnitrin in P. putida 1A00316. The results of the present study reveal the diversity and specificity of nematicidal factors in P. putida strain 1A00316.

Long intergenic non-coding RNAs (lincRNAs) may play widespread roles in gene regulation and other biological processes, however, a systematic examination of the functions of lincRNAs in the biological responses of rice to phosphate (Pi) starvation has not been performed. Here, we used a computational method to predict the functions of lincRNAs in Pi-starved rice. Overall, 3,170 lincRNA loci were identified using RNA sequencing data from the roots and shoots of control and Pi-starved rice. A competing endogenous RNA (ceRNA) network was constructed for each tissue by considering the competing relationships between lincRNAs and genes, and the correlations between the expression levels of RNAs in ceRNA pairs. Enrichment analyses showed that most of the communities in the networks were related to the biological processes of Pi starvation. The lincRNAs in the two tissues were individually functionally annotated based on the ceRNA networks, and the differentially expressed lincRNAs were biologically meaningful. For example, XLOC_026030 was upregulated from 3 days after Pi starvation, and its functional annotation was ‘cellular response to Pi starvation’. In conclusion, we systematically annotated lincRNAs in rice and identified those involved in the biological response to Pi starvation.