A non-canonical role for a small nucleolar RNA in ribosome biogenesis and senescence.

Cellular senescence is an irreversible state of cell-cycle arrest induced by various stresses, including aberrant oncogene activation, telomere shortening, and DNA damage. Through a genome-wide screen, we discovered a conserved small nucleolar RNA (snoRNA), SNORA13, that is required for multiple forms of senescence in human cells and mice. Although SNORA13 guides the pseudouridylation of a conserved nucleotide in the ribosomal decoding center, loss of this snoRNA minimally impacts translation. Instead, we found that SNORA13 negatively regulates ribosome biogenesis. Senescence-inducing stress perturbs ribosome biogenesis, resulting in the accumulation of free ribosomal proteins (RPs) that trigger p53 activation. SNORA13 interacts directly with RPL23, decreasing its incorporation into maturing 60S subunits and, consequently, increasing the pool of free RPs, thereby promoting p53-mediated senescence. Thus, SNORA13 regulates ribosome biogenesis and the p53 pathway through a non-canonical mechanism distinct from its role in guiding RNA modification. These findings expand our understanding of snoRNA functions and their roles in cellular signaling.

The enormous repetitive Antarctic krill genome reveals environmental adaptations and population insights.

Antarctic krill (Euphausia superba) is Earth's most abundant wild animal, and its enormous biomass is vital to the Southern Ocean ecosystem. Here, we report a 48.01-Gb chromosome-level Antarctic krill genome, whose large genome size appears to have resulted from inter-genic transposable element expansions. Our assembly reveals the molecular architecture of the Antarctic krill circadian clock and uncovers expanded gene families associated with molting and energy metabolism, providing insights into adaptations to the cold and highly seasonal Antarctic environment. Population-level genome re-sequencing from four geographical sites around the Antarctic continent reveals no clear population structure but highlights natural selection associated with environmental variables. An apparent drastic reduction in krill population size 10 mya and a subsequent rebound 100 thousand years ago coincides with climate change events. Our findings uncover the genomic basis of Antarctic krill adaptations to the Southern Ocean and provide valuable resources for future Antarctic research.

African lungfish genome sheds light on the vertebrate water-to-land transition.

Lungfishes are the closest extant relatives of tetrapods and preserve ancestral traits linked with the water-to-land transition. However, their huge genome sizes have hindered understanding of this key transition in evolution. Here, we report a 40-Gb chromosome-level assembly of the African lungfish (Protopterus annectens) genome, which is the largest genome assembly ever reported and has a contig and chromosome N50 of 1.60 Mb and 2.81 Gb, respectively. The large size of the lungfish genome is due mainly to retrotransposons. Genes with ultra-long length show similar expression levels to other genes, indicating that lungfishes have evolved high transcription efficacy to keep gene expression balanced. Together with transcriptome and experimental data, we identified potential genes and regulatory elements related to such terrestrial adaptation traits as pulmonary surfactant, anxiolytic ability, pentadactyl limbs, and pharyngeal remodeling. Our results provide insights and key resources for understanding the evolutionary pathway leading from fishes to humans.

Target-directed microRNA degradation regulates developmental microRNA expression and embryonic growth in mammals.

MicroRNAs (miRNAs) are post-transcriptional regulators of gene expression that play critical roles in development and disease. Target-directed miRNA degradation (TDMD), a pathway in which miRNAs that bind to specialized targets with extensive complementarity are rapidly decayed, has emerged as a potent mechanism of controlling miRNA levels. Nevertheless, the biological role and scope of miRNA regulation by TDMD in mammals remains poorly understood. To address these questions, we generated mice with constitutive or conditional deletion of Zswim8, which encodes an essential TDMD factor. Loss of Zswim8 resulted in developmental defects in the heart and lungs, growth restriction, and perinatal lethality. Small RNA sequencing of embryonic tissues revealed widespread miRNA regulation by TDMD and greatly expanded the known catalog of miRNAs regulated by this pathway. These experiments also uncovered novel features of TDMD-regulated miRNAs, including their enrichment in cotranscribed clusters and examples in which TDMD underlies "arm switching," a phenomenon wherein the dominant strand of a miRNA precursor changes in different tissues or conditions. Importantly, deletion of two miRNAs, miR-322 and miR-503, rescued growth of Zswim8-null embryos, directly implicating the TDMD pathway as a regulator of mammalian body size. These data illuminate the broad landscape and developmental role of TDMD in mammals.

RBM33 directs the nuclear export of transcripts containing GC-rich elements.

Although splicing is a major driver of RNA nuclear export, many intronless RNAs are efficiently exported to the cytoplasm through poorly characterized mechanisms. For example, GC-rich sequences promote nuclear export in a splicing-independent manner, but how GC content is recognized and coupled to nuclear export is unknown. Here, we developed a genome-wide screening strategy to investigate the mechanism of export of NORAD, an intronless cytoplasmic long noncoding RNA (lncRNA). This screen revealed an RNA binding protein, RBM33, that directs the nuclear export of NORAD and numerous other transcripts. RBM33 directly binds substrate transcripts and recruits components of the TREX-NXF1/NXT1 RNA export pathway. Interestingly, high GC content emerged as the feature that specifies RBM33-dependent nuclear export. Accordingly, RBM33 directly binds GC-rich elements in target transcripts. These results provide a broadly applicable strategy for the genetic dissection of nuclear export mechanisms and reveal a long-sought nuclear export pathway for transcripts with GC-rich sequences.

Algorithm for optimized mRNA design improves stability and immunogenicity.

Messenger RNA (mRNA) vaccines are being used to combat the spread of COVID-19 (refs. 1-3), but they still exhibit critical limitations caused by mRNA instability and degradation, which are major obstacles for the storage, distribution and efficacy of the vaccine products4. Increasing secondary structure lengthens mRNA half-life, which, together with optimal codons, improves protein expression5. Therefore, a principled mRNA design algorithm must optimize both structural stability and codon usage. However, owing to synonymous codons, the mRNA design space is prohibitively large-for example, there are around 2.4 × 10632 candidate mRNA sequences for the SARS-CoV-2 spike protein. This poses insurmountable computational challenges. Here we provide a simple and unexpected solution using the classical concept of lattice parsing in computational linguistics, where finding the optimal mRNA sequence is analogous to identifying the most likely sentence among similar-sounding alternatives6. Our algorithm LinearDesign finds an optimal mRNA design for the spike protein in just 11 minutes, and can concurrently optimize stability and codon usage. LinearDesign substantially improves mRNA half-life and protein expression, and profoundly increases antibody titre by up to 128 times in mice compared to the codon-optimization benchmark on mRNA vaccines for COVID-19 and varicella-zoster virus. This result reveals the great potential of principled mRNA design and enables the exploration of previously unreachable but highly stable and efficient designs. Our work is a timely tool for vaccines and other mRNA-based medicines encoding therapeutic proteins such as monoclonal antibodies and anti-cancer drugs7,8.

Fluctuation of lysosomal protein degradation in neural stem cells of the postnatal mouse brain.

Lysosomes are intracellular organelles responsible for degrading diverse macromolecules delivered from several pathways, including the endo-lysosomal and autophagic pathways. Recent reports have suggested that lysosomes are essential for regulating neural stem cells in developing, adult and aged brains. However, the activity of these lysosomes has yet to be monitored in these brain tissues. Here, we report the development of a new probe to measure lysosomal protein degradation in brain tissue by immunostaining. Our results indicate that lysosomal protein degradation fluctuates in neural stem cells of the hippocampal dentate gyrus, depending on age and brain disorders. Neural stem cells increase their lysosomal activity during hippocampal development in the dentate gyrus, but aging and aging-related disease reduce lysosomal activity. In addition, physical exercise increases lysosomal activity in neural stem cells and astrocytes in the dentate gyrus. We therefore propose that three different stages of lysosomal activity exist: the state of increase during development, the stable state during adulthood and the state of reduction due to damage caused by either age or disease.

Carbon starvation raises capacities in bacterial antibiotic resistance and viral auxiliary carbon metabolism in soils.

Organic carbon availability in soil is crucial for shaping microbial communities, yet, uncertainties persist concerning microbial adaptations to carbon levels and the ensuing ecological and evolutionary consequences. We investigated organic carbon metabolism, antibiotic resistance, and virus-host interactions in soils subjected to 40 y of chemical and organic fertilization that led to contrasting carbon availability: carbon-poor and carbon-rich soils, respectively. Carbon-poor soils drove the enrichment of putative genes involved in organic matter decomposition and exhibited specialization in utilizing complex organic compounds, reflecting scramble competition. This specialization confers a competitive advantage of microbial communities in carbon-poor soils but reduces their buffering capacity in terms of organic carbon metabolisms, making them more vulnerable to environmental fluctuations. Additionally, in carbon-poor soils, viral auxiliary metabolic genes linked to organic carbon metabolism increased host competitiveness and environmental adaptability through a strategy akin to "piggyback the winner." Furthermore, putative antibiotic resistance genes, particularly in low-abundance drug categories, were enriched in carbon-poor soils as an evolutionary consequence of chemical warfare (i.e., interference competition). This raises concerns about the potential dissemination of antibiotic resistance from conventional agriculture that relies on chemical-only fertilization. Consequently, carbon starvation resulting from long-term chemical-only fertilization increases microbial adaptations to competition, underscoring the importance of implementing sustainable agricultural practices to mitigate the emergence and spread of antimicrobial resistance and to increase soil carbon storage.

miR-26 suppresses adipocyte progenitor differentiation and fat production by targeting Fbxl19.

Fat storage in adult mammals is a highly regulated process that involves the mobilization of adipocyte progenitor cells (APCs) that differentiate to produce new adipocytes. Here we report a role for the broadly conserved miR-26 family of microRNAs (miR-26a-1, miR-26a-2, and miR-26b) as major regulators of APC differentiation and adipose tissue mass. Deletion of all miR-26-encoding loci in mice resulted in a dramatic expansion of adipose tissue in adult animals fed normal chow. Conversely, transgenic overexpression of miR-26a protected mice from high-fat diet-induced obesity. These effects were attributable to a cell-autonomous function of miR-26 as a potent inhibitor of APC differentiation. miR-26 blocks adipogenesis, at least in part, by repressing expression of Fbxl19, a conserved miR-26 target without a previously known role in adipocyte biology that encodes a component of SCF-type E3 ubiquitin ligase complexes. These findings have therefore revealed a novel pathway that plays a critical role in regulating adipose tissue formation in vivo and suggest new potential therapeutic targets for obesity and related disorders.

Zn2+-dependent association of cysteine-rich protein with virion orchestrates morphogenesis of rod-shaped viruses.

The majority of rod-shaped and some filamentous plant viruses encode a cysteine-rich protein (CRP) that functions in viral virulence; however, the roles of these CRPs in viral infection remain largely unknown. Here, we used barley stripe mosaic virus (BSMV) as a model to investigate the essential role of its CRP in virus morphogenesis. The CRP protein γb directly interacts with BSMV coat protein (CP), the mutations either on the His-85 site in γb predicted to generate a potential CCCH motif or on the His-13 site in CP exposed to the surface of the virions abolish the zinc-binding activity and their interaction. Immunogold-labeling assays show that γb binds to the surface of rod-shaped BSMV virions in a Zn2+-dependent manner, which enhances the RNA binding activity of CP and facilitates virion assembly and stability, suggesting that the Zn2+-dependent physical association of γb with the virion is crucial for BSMV morphogenesis. Intriguingly, the tightly binding of diverse CRPs to their rod-shaped virions is a general feature employed by the members in the families Virgaviridae (excluding the genus Tobamovirus) and Benyviridae. Together, these results reveal a hitherto unknown role of CRPs in the assembly and stability of virus particles, and expand our understanding of the molecular mechanism underlying virus morphogenesis.

Natural variation in the prolyl 4-hydroxylase gene PtoP4H9 contributes to perennial stem growth in Populus.

Perennial trees must maintain stem growth throughout their entire lifespan to progressively increase in size as they age. The overarching question of the molecular mechanisms that govern stem perennial growth in trees remains largely unanswered. Here we deciphered the genetic architecture that underlies perennial growth trajectories using genome-wide association studies (GWAS) for measures of growth traits across years in a natural population of Populus tomentosa. By analyzing the stem growth trajectory, we identified PtoP4H9, encoding prolyl 4-hydroxylase 9, which is responsible for the natural variation in the growth rate of diameter at breast height (DBH) across years. Quantifying the dynamic genetic contribution of PtoP4H9 loci to stem growth showed that PtoP4H9 played a pivotal role in stem growth regulation. Spatiotemporal expression analysis showed that PtoP4H9 was highly expressed in cambium tissues of poplars of various ages. Overexpression and knockdown of PtoP4H9 revealed that it altered cell expansion to regulate cell wall modification and mechanical characteristics, thereby promoting stem growth in Populus. We showed that natural variation in PtoP4H9 occurred in a BASIC PENTACYSTEINE transcription factor PtoBPC1-binding promoter element controlling PtoP4H9 expression. The geographic distribution of PtoP4H9 allelic variation was consistent with the modes of selection among populations. Altogether, our study provides important genetic insights into dynamic stem growth in Populus, and we confirmed PtoP4H9 as a potential useful marker for breeding or genetic engineering of poplars.

In situ formation of cocatalytic sites boosts single-atom catalysts for nitrogen oxide reduction.

Nitrogen oxide (NOx) pollution presents a severe threat to the environment and human health. Catalytic reduction of NOx with H2 using single-atom catalysts poses considerable potential in the remediation of air pollution; however, the unfavorable process of H2 dissociation limits its practical application. Herein, we report that the in situ formation of PtTi cocatalytic sites (which are stabilized by Pt-Ti bonds) over Pt1/TiO2 significantly increases NOx conversion by reducing the energy barrier of H2 activation. We demonstrate that two H atoms of H2 molecule are absorbed by adjacent Pt atoms in Pt-O and Pt-Ti, respectively, which can promote the cleave of H-H bonds. Besides, PtTi sites facilitate the adsorption of NO molecules and further lower the activation barrier of the whole de-NOx reaction. Extending the concept to Pt1/Nb2O5 and Pd1/TiO2 systems also sees enhanced catalytic activities, demonstrating that engineering the cocatalytic sites can be a general strategy for the design of high-efficiency catalysts that can benefit environmental sustainability.

The MAPK-Alfin-like 7 module negatively regulates ROS scavenging genes to promote NLR-mediated immunity.

Nucleotide-binding leucine-rich repeat (NLR) receptor-mediated immunity includes rapid production of reactive oxygen species (ROS) and transcriptional reprogramming, which is controlled by transcription factors (TFs). Although some TFs have been reported to participate in NLR-mediated immune response, most TFs are transcriptional activators, and whether and how transcriptional repressors regulate NLR-mediated plant defenses remains largely unknown. Here, we show that the Alfin-like 7 (AL7) interacts with N NLR and functions as a transcriptional repressor. Knockdown and knockout of AL7 compromise N NLR-mediated resistance against tobacco mosaic virus, whereas AL7 overexpression enhances defense, indicating a positive regulatory role for AL7 in immunity. AL7 binds to the promoters of ROS scavenging genes to inhibit their transcription during immune responses. Mitogen-activated protein kinases (MAPKs), salicylic acid-induced protein kinase (SIPK), and wound-induced protein kinase (WIPK) directly interact with and phosphorylate AL7, which impairs the AL7-N interaction and enhances its DNA binding activity, which promotes ROS accumulation and enables immune activation. In addition to N, AL7 is also required for the function of other Toll interleukin 1 receptor/nucleotide-binding/leucine-rich repeats (TNLs) including Roq1 and RRS1-R/RPS4. Our findings reveal a hitherto unknown MAPK-AL7 module that negatively regulates ROS scavenging genes to promote NLR-mediated immunity.

Loss of Dis3l2 partially phenocopies Perlman syndrome in mice and results in up-regulation of Igf2 in nephron progenitor cells.

Loss of function of the DIS3L2 exoribonuclease is associated with Wilms tumor and the Perlman congenital overgrowth syndrome. LIN28, a Wilms tumor oncoprotein, triggers the DIS3L2-mediated degradation of the precursor of let-7, a microRNA that inhibits Wilms tumor development. These observations have led to speculation that DIS3L2-mediated tumor suppression is attributable to let-7 regulation. Here we examine new DIS3L2-deficient cell lines and mouse models, demonstrating that DIS3L2 loss has no effect on mature let-7 levels. Rather, analysis of Dis3l2-null nephron progenitor cells, a potential cell of origin of Wilms tumors, reveals up-regulation of Igf2, a growth-promoting gene strongly associated with Wilms tumorigenesis. These findings nominate a new potential mechanism underlying the pathology associated with DIS3L2 deficiency.

An Immunocompetent Mongolian Gerbil Model for Hepatitis E Virus Genotype 1 Infection.

BACKGROUND & AIMS: Hepatitis E virus (HEV), primarily genotype 1 (HEV-1), causes approximately 20.1 million infections, 44,000 deaths, and 3000 stillbirths annually. Current evidence indicates that HEV-1 is only transmitted in humans. Here, we evaluated whether Mongolian gerbils can serve as animal models for HEV-1 infection. METHODS: Mongolian gerbils were used for HEV-1 and hepatitis E virus genotype 3 infection experiments. HEV infection parameters, including detection of HEV RNA and HEV antigen, liver function assessment, and histopathology, were evaluated. RESULTS: We adapted a clinical isolate of HEV-1 for Mongolian gerbils by serial passaging in feces of aged male gerbils. The gerbil-adapted strain obtained at passage 3 induced a robust, acute HEV infection, characterized by stable fecal virus shedding, elevated liver enzymes, histopathologic changes in the liver, and seroconversion to anti-HEV. An infectious complementary DNA clone of the adapted virus was generated. HEV-1-infected pregnant gerbils showed a high rate of maternal mortality and vertical transmission. HEV RNA or antigens were detected in the liver, kidney, intestine, placenta, testis, and fetus liver. Liver and placental transcriptomic analyses indicated activation of host immunity. Tacrolimus prolonged HEV-1 infection, whereas ribavirin cleared infection. The protective efficacy of a licensed HEV vaccine was validated using this model. CONCLUSIONS: HEV-1 efficiently infected Mongolian gerbils. This HEV-1 infection model will be valuable for investigating hepatitis E immunopathogenesis and evaluating vaccines and antivirals against HEV.

Optimization of immunosuppression strategies for the establishment of chronic hepatitis E virus infection in rabbits.

Chronic hepatitis E mostly occurs in organ transplant recipients and can lead to rapid liver fibrosis and cirrhosis. Previous studies found that the development of chronic hepatitis E virus (HEV) infection is linked to the type of immunosuppressant used. Animal models are crucial for the study of pathogenesis of chronic hepatitis E. We previously established a stable chronic HEV infection rabbit model using cyclosporine A (CsA), a calcineurin inhibitor (CNI)-based immunosuppressant. However, the immunosuppression strategy and timing may be optimized, and how different types of immunosuppressants affect the establishment of chronic HEV infection in this model is still unknown. Here, we showed that chronic HEV infection can be established in 100% of rabbits when CsA treatment was started at HEV challenge or even 4 weeks after. Tacrolimus or prednisolone treatment alone also contributed to chronic HEV infection, resulting in 100% and 77.8% chronicity rates, respectively, while mycophenolate mofetil (MMF) only led to a 28.6% chronicity rate. Chronic HEV infection was accompanied with a persistent activation of innate immune response evidenced by transcriptome analysis. The suppressed adaptive immune response evidenced by low expression of genes related to cytotoxicity (like perforin and FasL) and low anti-HEV seroconversion rates may play important roles in causing chronic HEV infection. By analyzing HEV antigen concentrations with different infection outcomes, we also found that HEV antigen levels could indicate chronic HEV infection development. This study optimized the immunosuppression strategies for establishing chronic HEV infection in rabbits and highlighted the potential association between the development of chronic HEV infection and immunosuppressants.IMPORTANCEOrgan transplant recipients are at high risk of chronic hepatitis E and generally receive a CNI-based immunosuppression regimen containing CNI (tacrolimus or CsA), MMF, and/or corticosteroids. Previously, we established stable chronic HEV infection in a rabbit model by using CsA before HEV challenge. In this study, we further optimized the immunosuppression strategies for establishing chronic HEV infection in rabbits. Chronic HEV infection can also be established when CsA treatment was started at the same time or even 4 weeks after HEV challenge, clearly indicating the risk of progression to chronic infection under these circumstances and the necessity of HEV screening for both the recipient and the donor preoperatively. CsA, tacrolimus, or prednisolone instead of MMF significantly contributed to chronic HEV infection. HEV antigen in acute infection phase indicates the development of chronic infection. Our results have important implications for understanding the potential association between chronic HEV infection and immunosuppressants.

Transcription factors PuPRE6/PuMYB12 and histone deacetylase PuHDAC9-like regulate sucrose levels in pear.

The improvement of fruit quality, in particular sugar content, has been a major goal of plant breeding programmes for many years. Here, 2 varieties of the Ussurian pear (Pyrus ussuriensis), Nanguo, and its high-sucrose accumulation bud sport, Nanhong, were used to study the molecular mechanisms regulating sucrose transport in fruits. Comparative transcriptome analysis showed that in Nanhong fruit, an MYB transcription factor, PuMYB12, and a sucrose transporter protein, PuSUT4-like, were expressed at higher levels, while a paclobutrazol resistance transcription factor, PuPRE6, and a histone deacetylase (HDAC), PuHDAC9-like, were expressed at lower levels in Nanguo fruit. PuSUT4-like silencing and overexpression experiments in Nanguo pear showed that PuSUT4-like is essential for sucrose transportation. PuPRE6 and PuMYB12 act as antagonistic complexes to regulate PuSUT4-like transcription and sucrose accumulation. The histone deacetylation levels of the PuMYB12 and PuSUT4-like promoters were higher in Nanguo fruit than in Nanhong fruit, and Y1H assays showed that HDAC PuHDAC9-like bound directly to the promoters of PuMYB12 and PuSUT4-like. Our results uncovered transcription regulation and epigenetic mechanisms underlying sucrose accumulation in pears.

The extracellular cyclophilin A-integrin ß2 complex as a therapeutic target of viral pneumonia.

The acute respiratory virus infection can induce uncontrolled inflammatory responses, such as cytokine storm and viral pneumonia, which are the major causes of death in clinical cases. Cyclophilin A (CypA) is mainly distributed in the cytoplasm of resting cells and released into the extracellular space in response to inflammatory stimuli. Extracellular CypA (eCypA) is upregulated and promotes inflammatory response in severe COVID-19 patients. However, how eCypA promotes virus-induced inflammatory response remains elusive. Here, we observe that eCypA is induced by influenza A and B viruses and SARS-CoV-2 in cells, mice, or patients. Anti-CypA mAb reduces pro-inflammatory cytokines production, leukocytes infiltration, and lung injury in virus-infected mice. Mechanistically, eCypA binding to integrin ß2 triggers integrin activation, thereby facilitating leukocyte trafficking and cytokines production via the focal adhesion kinase (FAK)/GTPase and FAK/ERK/P65 pathways, respectively. These functions are suppressed by the anti-CypA mAb that specifically blocks eCypA-integrin ß2 interaction. Overall, our findings reveal that eCypA-integrin ß2 signaling mediates virus-induced inflammatory response, indicating that eCypA is a potential target for antibody therapy against viral pneumonia.

LazySampling and LinearSampling: fast stochastic sampling of RNA secondary structure with applications to SARS-CoV-2.

Many RNAs fold into multiple structures at equilibrium, and there is a need to sample these structures according to their probabilities in the ensemble. The conventional sampling algorithm suffers from two limitations: (i) the sampling phase is slow due to many repeated calculations; and (ii) the end-to-end runtime scales cubically with the sequence length. These issues make it difficult to be applied to long RNAs, such as the full genomes of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). To address these problems, we devise a new sampling algorithm, LazySampling, which eliminates redundant work via on-demand caching. Based on LazySampling, we further derive LinearSampling, an end-to-end linear time sampling algorithm. Benchmarking on nine diverse RNA families, the sampled structures from LinearSampling correlate better with the well-established secondary structures than Vienna RNAsubopt and RNAplfold. More importantly, LinearSampling is orders of magnitude faster than standard tools, being 428× faster (72 s versus 8.6 h) than RNAsubopt on the full genome of SARS-CoV-2 (29 903 nt). The resulting sample landscape correlates well with the experimentally guided secondary structure models, and is closer to the alternative conformations revealed by experimentally driven analysis. Finally, LinearSampling finds 23 regions of 15 nt with high accessibilities in the SARS-CoV-2 genome, which are potential targets for COVID-19 diagnostics and therapeutics.

LinearCoFold and LinearCoPartition: linear-time algorithms for secondary structure prediction of interacting RNA molecules.

Many RNAs function through RNA-RNA interactions. Fast and reliable RNA structure prediction with consideration of RNA-RNA interaction is useful, however, existing tools are either too simplistic or too slow. To address this issue, we present LinearCoFold, which approximates the complete minimum free energy structure of two strands in linear time, and LinearCoPartition, which approximates the cofolding partition function and base pairing probabilities in linear time. LinearCoFold and LinearCoPartition are orders of magnitude faster than RNAcofold. For example, on a sequence pair with combined length of 26,190 nt, LinearCoFold is 86.8× faster than RNAcofold MFE mode, and LinearCoPartition is 642.3× faster than RNAcofold partition function mode. Surprisingly, LinearCoFold and LinearCoPartition's predictions have higher PPV and sensitivity of intermolecular base pairs. Furthermore, we apply LinearCoFold to predict the RNA-RNA interaction between SARS-CoV-2 genomic RNA (gRNA) and human U4 small nuclear RNA (snRNA), which has been experimentally studied, and observe that LinearCoFold's prediction correlates better with the wet lab results than RNAcofold's.