A designer synthetic chromosome fragment functions in moss.

Rapid advances in DNA synthesis techniques have enabled the assembly and engineering of viral and microbial genomes, presenting new opportunities for synthetic genomics in multicellular eukaryotic organisms. These organisms, characterized by larger genomes, abundant transposons and extensive epigenetic regulation, pose unique challenges. Here we report the in vivo assembly of chromosomal fragments in the moss Physcomitrium patens, producing phenotypically virtually wild-type lines in which one-third of the coding region of a chromosomal arm is replaced by redesigned, chemically synthesized fragments. By eliminating 55.8% of a 155 kb endogenous chromosomal region, we substantially simplified the genome without discernible phenotypic effects, implying that many transposable elements may minimally impact growth. We also introduced other sequence modifications, such as PCRTag incorporation, gene locus swapping and stop codon substitution. Despite these substantial changes, the complex epigenetic landscape was normally established, albeit with some three-dimensional conformation alterations. The synthesis of a partial multicellular eukaryotic chromosome arm lays the foundation for the synthetic moss genome project (SynMoss) and paves the way for genome synthesis in multicellular organisms.

Decoding the Heterogeneity and Specialized Function of Translation Machinery Through Ribosome Profiling in Yeast Mutants of Initiation Factors.

The nuanced heterogeneity and specialized functions of translation machinery are increasingly recognized as crucial for precise translational regulation. Here, high-throughput ribosomal profiling (ribo-seq) is used to analyze the specialized roles of eukaryotic initiation factors (eIFs) in the budding yeast. By examining changes in ribosomal distribution across the genome resulting from knockouts of eIF4A, eIF4B, eIF4G1, CAF20, or EAP1, or knockdowns of eIF1, eIF1A, eIF4E, or PAB1, two distinct initiation-factor groups, the "looping" and "scanning" groups are discerned, based on similarities in the ribosomal landscapes their perturbation induced. The study delves into the cis-regulatory sequence features of genes influenced predominantly by each group, revealing that genes more dependent on the looping-group factors generally have shorter transcripts and poly(A) tails. In contrast, genes more dependent on the scanning-group factors often possess upstream open reading frames and exhibit a higher GC content in their 5' untranslated regions. From the ribosomal RNA fragments identified in the ribo-seq data, ribosomal heterogeneity associated with perturbation of specific initiation factors is further identified, suggesting their potential roles in regulating ribosomal components. Collectively, the study illuminates the complexity of translational regulation driven by heterogeneity and specialized functions of translation machinery, presenting potential approaches for targeted gene translation manipulation.

A translational regulator MHZ9 modulates ethylene signaling in rice.

Ethylene plays essential roles in rice growth, development and stress adaptation. Translational control of ethylene signaling remains unclear in rice. Here, through analysis of an ethylene-response mutant mhz9, we identified a glycine-tyrosine-phenylalanine (GYF) domain protein MHZ9, which positively regulates ethylene signaling at translational level in rice. MHZ9 is localized in RNA processing bodies. The C-terminal domain of MHZ9 interacts with OsEIN2, a central regulator of rice ethylene signaling, and the N-terminal domain directly binds to the OsEBF1/2 mRNAs for translational inhibition, allowing accumulation of transcription factor OsEIL1 to activate the downstream signaling. RNA-IP seq and CLIP-seq analyses reveal that MHZ9 associates with hundreds of RNAs. Ribo-seq analysis indicates that MHZ9 is required for the regulation of ~ 90% of genes translationally affected by ethylene. Our study identifies a translational regulator MHZ9, which mediates translational regulation of genes in response to ethylene, facilitating stress adaptation and trait improvement in rice.

A high-resolution transcriptomic atlas depicting nitrogen fixation and nodule development in soybean.

Although root nodules are essential for biological nitrogen fixation in legumes, the cell types and molecular regulatory mechanisms contributing to nodule development and nitrogen fixation in determinate nodule legumes, such as soybean (Glycine max), remain incompletely understood. Here, we generated a single-nucleus resolution transcriptomic atlas of soybean roots and nodules at 14 days post inoculation (dpi) and annotated 17 major cell types, including six that are specific to nodules. We identified the specific cell types responsible for each step in the ureides synthesis pathway, which enables spatial compartmentalization of biochemical reactions during soybean nitrogen fixation. By utilizing RNA velocity analysis, we reconstructed the differentiation dynamics of soybean nodules, which differs from those of indeterminate nodules in Medicago truncatula. Moreover, we identified several putative regulators of soybean nodulation and two of these genes, GmbHLH93 and GmSCL1, were as-yet uncharacterized in soybean. Overexpression of each gene in soybean hairy root systems validated their respective roles in nodulation. Notably, enrichment for cytokinin-related genes in soybean nodules led to identification of the cytokinin receptor, GmCRE1, as a prominent component of the nodulation pathway. GmCRE1 knockout in soybean resulted in a striking nodule phenotype with decreased nitrogen fixation zone and depletion of leghemoglobins, accompanied by downregulation of nodule-specific gene expression, as well as almost complete abrogation of biological nitrogen fixation. In summary, this study provides a comprehensive perspective of the cellular landscape during soybean nodulation, shedding light on the underlying metabolic and developmental mechanisms of soybean nodule formation.

Distance-dependent inhibition of translation initiation by downstream out-of-frame AUGs is consistent with a Brownian ratchet process of ribosome scanning.

BACKGROUND: Eukaryotic ribosomes are widely presumed to scan mRNA for the AUG codon to initiate translation in a strictly 5'-3' movement (i.e., strictly unidirectional scanning model), so that ribosomes initiate translation exclusively at the 5' proximal AUG codon (i.e., the first-AUG rule). RESULTS: We generate 13,437 yeast variants, each with an ATG triplet placed downstream (dATGs) of the annotated ATG (aATG) codon of a green fluorescent protein. We find that out-of-frame dATGs can inhibit translation at the aATG, but with diminishing strength over increasing distance between aATG and dATG, undetectable beyond ~17 nt. This phenomenon is best explained by a Brownian ratchet mechanism of ribosome scanning, in which the ribosome uses small-amplitude 5'-3' and 3'-5' oscillations with a net 5'-3' movement to scan the AUG codon, thereby leading to competition for translation initiation between aAUG and a proximal dAUG. This scanning model further predicts that the inhibitory effect induced by an out-of-frame upstream AUG triplet (uAUG) will diminish as uAUG approaches aAUG, which is indeed observed among the 15,586 uATG variants generated in this study. Computational simulations suggest that each triplet is scanned back and forth approximately ten times until the ribosome eventually migrates to downstream regions. Moreover, this scanning process could constrain the evolution of sequences downstream of the aATG to minimize proximal out-of-frame dATG triplets in yeast and humans. CONCLUSIONS: Collectively, our findings uncover the basic process by which eukaryotic ribosomes scan for initiation codons, and how this process could shape eukaryotic genome evolution.

Evidence for a mouse origin of the SARS-CoV-2 Omicron variant.

The rapid accumulation of mutations in the SARS-CoV-2 Omicron variant that enabled its outbreak raises questions as to whether its proximal origin occurred in humans or another mammalian host. Here, we identified 45 point mutations that Omicron acquired since divergence from the B.1.1 lineage. We found that the Omicron spike protein sequence was subjected to stronger positive selection than that of any reported SARS-CoV-2 variants known to evolve persistently in human hosts, suggesting a possibility of host-jumping. The molecular spectrum of mutations (i.e., the relative frequency of the 12 types of base substitutions) acquired by the progenitor of Omicron was significantly different from the spectrum for viruses that evolved in human patients but resembled the spectra associated with virus evolution in a mouse cellular environment. Furthermore, mutations in the Omicron spike protein significantly overlapped with SARS-CoV-2 mutations known to promote adaptation to mouse hosts, particularly through enhanced spike protein binding affinity for the mouse cell entry receptor. Collectively, our results suggest that the progenitor of Omicron jumped from humans to mice, rapidly accumulated mutations conducive to infecting that host, then jumped back into humans, indicating an inter-species evolutionary trajectory for the Omicron outbreak.

Host-specific asymmetric accumulation of mutation types reveals that the origin of SARS-CoV-2 is consistent with a natural process.

The capacity of RNA viruses to adapt to new hosts and rapidly escape the host immune system is largely attributable to de novo genetic diversity that emerges through mutations in RNA. Although the molecular spectrum of de novo mutations-the relative rates at which various base substitutions occur-are widely recognized as informative toward understanding the evolution of a viral genome, little attention has been paid to the possibility of using molecular spectra to infer the host origins of a virus. Here, we characterize the molecular spectrum of de novo mutations for SARS-CoV-2 from transcriptomic data obtained from virus-infected cell lines, enabled by the use of sporadic junctions formed during discontinuous transcription as molecular barcodes. We find that de novo mutations are generated in a replication-independent manner, typically on the genomic strand, and highly dependent on mutagenic mechanisms specific to the host cellular environment. De novo mutations will then strongly influence the types of base substitutions accumulated during SARS-CoV-2 evolution, in an asymmetric manner favoring specific mutation types. Consequently, similarities between the mutation spectra of SARS-CoV-2 and the bat coronavirus RaTG13, which have accumulated since their divergence strongly suggest that SARS-CoV-2 evolved in a host cellular environment highly similar to that of bats before its zoonotic transfer into humans. Collectively, our findings provide data-driven support for the natural origin of SARS-CoV-2.

Single-cell transcriptome atlas of the leaf and root of rice seedlings.

As a multicellular organism, rice flourishes relying on gene expression diversity among cells of various functions. However, cellular-resolution transcriptome features are yet to be fully recognized, let alone cell-specific transcriptional responses to environmental stimuli. In this study, we apply single-cell RNA sequencing to both shoot and root of rice seedlings growing in Kimura B nutrient solution or exposed to various abiotic stresses and characterize transcriptomes for a total of 237,431 individual cells. We identify 15 and 9 cell types in the leaf and root, respectively, and observe that common transcriptome features are often shared between leaves and roots in the same tissue layer, except for endodermis or epidermis. Abiotic stress stimuli alter gene expression largely in a cell type-specific manner, but for a given cell type, different stresses often trigger transcriptional regulation of roughly the same set of genes. Besides, we detect proportional changes in cell populations in response to abiotic stress and investigate the underlying molecular mechanisms through single-cell reconstruction of the developmental trajectory. Collectively, our study represents a benchmark-setting data resource of single-cell transcriptome atlas for rice seedlings and an illustration of exploiting such resources to drive discoveries in plant biology.

The memory of neuronal mitochondrial stress is inherited transgenerationally via elevated mitochondrial DNA levels.

The memory of stresses experienced by parents can be passed on to descendants as a forecast of the challenges to come. Here, we discovered that the neuronal mitochondrial perturbation-induced systemic mitochondrial unfolded protein response (UPRmt) in Caenorhabditis elegans can be transmitted to offspring over multiple generations. The transgenerational activation of UPRmt is mediated by maternal inheritance of elevated levels of mitochondrial DNA (mtDNA), which causes the proteostasis stress within mitochondria. Furthermore, results from intercrossing studies using wild C. elegans strains further support that maternal inheritance of higher levels of mtDNA can induce the UPRmt in descendants. The mitokine Wnt signalling pathway is required for the transmission of elevated mtDNA levels across generations, thereby conferring lifespan extension and stress resistance to offspring. Collectively, our results reveal that the nervous system can transmit stress signals across generations by increasing mtDNA in the germline, enabling descendants to better cope with anticipated challenges.

Integrated global analysis reveals a vitamin E-vitamin K1 sub-network, downstream of COLD1, underlying rice chilling tolerance divergence.

Rice, a staple food with tropical/subtropical origination, is susceptible to cold stress, one of the major constraints on its yield and distribution. Asian cultivated rice consists of two subspecies with diverged chilling tolerance to adapt to different environments. The mechanism underlying this divergence remains obscure with a few known factors, including membrane protein CHILLING-TOLERANCE DIVERGENCE 1 (COLD1). Here, we reveal a vitamin E-vitamin K1 sub-network responsible for chilling tolerance divergence through global analyses. Rice genome regions responsible for tolerance divergence are identified with chromosome segment substitution lines (CSSLs). Comparative transcriptomic and metabolomic analysis of chilling-tolerant CSSL4-1 and parent lines uncovered a vitamin E-vitamin K1 sub-network in chloroplast with tocopherol (vitamin E) mediating chloroplast-to-nucleus signaling. COLD1, located in the substitution segment in CSSL4-1, is confirmed as its upstream regulator by transgenic material analysis. Our work uncovers a pathway downstream of COLD1, through which rice modulates chilling tolerance for thermal adaptation, with potential utility in crop improvement.

The Ghd7 transcription factor represses ARE1 expression to enhance nitrogen utilization and grain yield in rice.

The genetic improvement of nitrogen use efficiency (NUE) of crops is vital for grain productivity and sustainable agriculture. However, the regulatory mechanism of NUE remains largely elusive. Here, we report that the rice Grain number, plant height, and heading date7 (Ghd7) gene genetically acts upstream of ABC1 REPRESSOR1 (ARE1), a negative regulator of NUE, to positively regulate nitrogen utilization. As a transcriptional repressor, Ghd7 directly binds to two Evening Element-like motifs in the promoter and intron 1 of ARE1, likely in a cooperative manner, to repress its expression. Ghd7 and ARE1 display diurnal expression patterns in an inverse oscillation manner, mirroring a regulatory scheme based on these two loci. Analysis of a panel of 2656 rice varieties suggests that the elite alleles of Ghd7 and ARE1 have undergone diversifying selection during breeding. Moreover, the allelic distribution of Ghd7 and ARE1 is associated with the soil nitrogen deposition rate in East Asia and South Asia. Remarkably, the combination of the Ghd7 and ARE1 elite alleles substantially improves NUE and yield performance under nitrogen-limiting conditions. Collectively, these results define a Ghd7-ARE1-based regulatory mechanism of nitrogen utilization, providing useful targets for genetic improvement of rice NUE.

The Missing Expression Level-Evolutionary Rate Anticorrelation in Viruses Does Not Support Protein Function as a Main Constraint on Sequence Evolution.

One of the central goals in molecular evolutionary biology is to determine the sources of variation in the rate of sequence evolution among proteins. Gene expression level is widely accepted as the primary determinant of protein evolutionary rate, because it scales with the extent of selective constraints imposed on a protein, leading to the well-known negative correlation between expression level and protein evolutionary rate (the E-R anticorrelation). Selective constraints have been hypothesized to entail the maintenance of protein function, the avoidance of cytotoxicity caused by protein misfolding or nonspecific protein-protein interactions, or both. However, empirical tests evaluating the relative importance of these hypotheses remain scarce, likely due to the nontrivial difficulties in distinguishing the effect of a deleterious mutation on a protein's function versus its cytotoxicity. We realized that examining the sequence evolution of viral proteins could overcome this hurdle. It is because purifying selection against mutations in a viral protein that result in cytotoxicity per se is likely relaxed, whereas purifying selection against mutations that impair viral protein function persists. Multiple analyses of SARS-CoV-2 and nine other virus species revealed a complete absence of any E-R anticorrelation. As a control, the E-R anticorrelation does exist in human endogenous retroviruses where purifying selection against cytotoxicity is present. Taken together, these observations do not support the maintenance of protein function as the main constraint on protein sequence evolution in cellular organisms.

Disome-seq reveals widespread ribosome collisions that promote cotranslational protein folding.

BACKGROUND: The folding of proteins is challenging in the highly crowded and sticky environment of a cell. Regulation of translation elongation may play a crucial role in ensuring the correct folding of proteins. Much of our knowledge regarding translation elongation comes from the sequencing of mRNA fragments protected by single ribosomes by ribo-seq. However, larger protected mRNA fragments have been observed, suggesting the existence of an alternative and previously hidden layer of regulation. RESULTS: In this study, we performed disome-seq to sequence mRNA fragments protected by two stacked ribosomes, a product of translational pauses during which the 5'-elongating ribosome collides with the 3'-paused one. We detected widespread ribosome collisions that are related to slow ribosome release when stop codons are at the A-site, slow peptide bond formation from proline, glycine, asparagine, and cysteine when they are at the P-site, and slow leaving of polylysine from the exit tunnel of ribosomes. The structure of disomes obtained by cryo-electron microscopy suggests a different conformation from the substrate of the ribosome-associated protein quality control pathway. Collisions occurred more frequently in the gap regions between α-helices, where a translational pause can prevent the folding interference from the downstream peptides. Paused or collided ribosomes are associated with specific chaperones, which can aid in the cotranslational folding of the nascent peptides. CONCLUSIONS: Therefore, cells use regulated ribosome collisions to ensure protein homeostasis.

Synchronized replication of genes encoding the same protein complex in fast-proliferating cells.

DNA replication perturbs the dosage balance among genes; at mid-S phase, early-replicating genes have doubled their copies while late-replicating ones have not. Dosage imbalance among genes, especially within members of a protein complex, is toxic to cells. However, the molecular mechanisms that cells use to deal with such imbalance remain not fully understood. Here, we validate at the genomic scale that the dosage between early- and late-replicating genes is imbalanced in HeLa cells. We propose the synchronized replication hypothesis that genes sensitive to stoichiometric relationships will be replicated simultaneously to maintain stoichiometry. In support of this hypothesis, we observe that genes encoding the same protein complex have similar replication timing but mainly in fast-proliferating cells such as embryonic stem cells and cancer cells. We find that the synchronized replication observed in cancer cells, but not in slow-proliferating differentiated cells, is due to convergent evolution during tumorigenesis that restores synchronized replication timing within protein complexes. Taken together, our study reveals that the demand for dosage balance during S phase plays an important role in the optimization of the replication-timing program; this selection is relaxed during differentiation as the cell cycle prolongs and is restored during tumorigenesis as the cell cycle shortens.

Impact of poly(A)-tail G-content on Arabidopsis PAB binding and their role in enhancing translational efficiency.

BACKGROUND: Polyadenylation plays a key role in producing mature mRNAs in eukaryotes. It is widely believed that the poly(A)-binding proteins (PABs) uniformly bind to poly(A)-tailed mRNAs, regulating their stability and translational efficiency. RESULTS: We observe that the homozygous triple mutant of broadly expressed Arabidopsis thaliana PABs, AtPAB2, AtPAB4, and AtPAB8, is embryonic lethal. To understand the molecular basis, we characterize the RNA-binding landscape of these PABs. The AtPAB-binding efficiency varies over one order of magnitude among genes. To identify the sequences accounting for the variation, we perform poly(A)-seq that directly sequences the full-length poly(A) tails. More than 10% of poly(A) tails contain at least one guanosine (G); among them, the G-content varies from 0.8 to 28%. These guanosines frequently divide poly(A) tails into interspersed A-tracts and therefore cause the variation in the AtPAB-binding efficiency among genes. Ribo-seq and genome-wide RNA stability assays show that AtPAB-binding efficiency of a gene is positively correlated with translational efficiency rather than mRNA stability. Consistently, genes with stronger AtPAB binding exhibit a greater reduction in translational efficiency when AtPAB is depleted. CONCLUSIONS: Our study provides a new mechanism that translational efficiency of a gene can be regulated through the G-content-dependent PAB binding, paving the way for a better understanding of poly(A) tail-associated regulation of gene expression.

Overdosage of Balanced Protein Complexes Reduces Proliferation Rate in Aneuploid Cells.

Cells with complex aneuploidies display a wide range of phenotypic abnormalities. However, the molecular basis for this has been mainly studied in trisomic (2n + 1) and disomic (n + 1) cells. To determine how karyotype affects proliferation in cells with complex aneuploidies, we generated 92 2n + x yeast strains in which each diploid cell has between 3 and 12 extra chromosomes. Genome-wide and, for individual protein complexes, proliferation defects are caused by the presence of protein complexes in which all subunits are balanced at the 3-copy level. Proteomics revealed that over 50% of 3-copy members of imbalanced complexes were expressed at only 2n protein levels, whereas members of complexes in which all subunits are stoichiometrically balanced at 3 copies per cell had 3n protein levels. We validated this finding using orthogonal datasets from yeast and from human cancers. Taken together, our study provides an explanation of how aneuploidy affects phenotype.

The lysine catabolite saccharopine impairs development by disrupting mitochondrial homeostasis.

Amino acid catabolism is frequently executed in mitochondria; however, it is largely unknown how aberrant amino acid metabolism affects mitochondria. Here we report the requirement for mitochondrial saccharopine degradation in mitochondrial homeostasis and animal development. In Caenorhbditis elegans, mutations in the saccharopine dehydrogenase (SDH) domain of the bi-functional enzyme α-aminoadipic semialdehyde synthase AASS-1 greatly elevate the lysine catabolic intermediate saccharopine, which causes mitochondrial damage by disrupting mitochondrial dynamics, leading to reduced adult animal growth. In mice, failure of mitochondrial saccharopine oxidation causes lethal mitochondrial damage in the liver, leading to postnatal developmental retardation and death. Importantly, genetic inactivation of genes that raise the mitochondrial saccharopine precursors lysine and α-ketoglutarate strongly suppresses SDH mutation-induced saccharopine accumulation and mitochondrial abnormalities in C. elegans Thus, adequate saccharopine catabolism is essential for mitochondrial homeostasis. Our study provides mechanistic and therapeutic insights for understanding and treating hyperlysinemia II (saccharopinuria), an aminoacidopathy with severe developmental defects.