Deep whole-genome analysis of 494 hepatocellular carcinomas.

Over half of hepatocellular carcinoma (HCC) cases diagnosed worldwide are in China1-3. However, whole-genome analysis of hepatitis B virus (HBV)-associated HCC in Chinese individuals is limited4-8, with current analyses of HCC mainly from non-HBV-enriched populations9,10. Here we initiated the Chinese Liver Cancer Atlas (CLCA) project and performed deep whole-genome sequencing (average depth, 120×) of 494 HCC tumours. We identified 6 coding and 28 non-coding previously undescribed driver candidates. Five previously undescribed mutational signatures were found, including aristolochic-acid-associated indel and doublet base signatures, and a single-base-substitution signature that we termed SBS_H8. Pentanucleotide context analysis and experimental validation confirmed that SBS_H8 was distinct to the aristolochic-acid-associated SBS22. Notably, HBV integrations could take the form of extrachromosomal circular DNA, resulting in elevated copy numbers and gene expression. Our high-depth data also enabled us to characterize subclonal clustered alterations, including chromothripsis, chromoplexy and kataegis, suggesting that these catastrophic events could also occur in late stages of hepatocarcinogenesis. Pathway analysis of all classes of alterations further linked non-coding mutations to dysregulation of liver metabolism. Finally, we performed in vitro and in vivo assays to show that fibrinogen alpha chain (FGA), determined as both a candidate coding and non-coding driver, regulates HCC progression and metastasis. Our CLCA study depicts a detailed genomic landscape and evolutionary history of HCC in Chinese individuals, providing important clinical implications.

LncRNA CCTT-mediated RNA-DNA and RNA-protein interactions facilitate the recruitment of CENP-C to centromeric DNA during kinetochore assembly.

Kinetochore assembly on centromeres is central for chromosome segregation, and defects in this process cause mitotic errors and aneuploidy. Besides the well-established protein network, emerging evidence suggests the involvement of regulatory RNA in kinetochore assembly; however, it has remained elusive about the identity of such RNA, let alone its mechanism of action in this critical process. Here, we report CCTT, a previously uncharacterized long non-coding RNA (lncRNA) transcribed from the arm of human chromosome 17, which plays a vital role in kinetochore assembly. We show that CCTT highly localizes to all centromeres via the formation of RNA-DNA triplex and specifically interacts with CENP-C to help engage this blueprint protein in centromeres, and consequently, CCTT loss triggers extensive mitotic errors and aneuploidy. These findings uncover a non-centromere-derived lncRNA that recruits CENP-C to centromeres and shed critical lights on the function of centromeric DNA sequences as anchor points for kinetochore assembly.

Extending the spectrum of fully integrated photonics to submicrometre wavelengths.

Integrated photonics has profoundly affected a wide range of technologies underpinning modern society1-4. The ability to fabricate a complete optical system on a chip offers unrivalled scalability, weight, cost and power efficiency5,6. Over the last decade, the progression from pure III-V materials platforms to silicon photonics has significantly broadened the scope of integrated photonics, by combining integrated lasers with the high-volume, advanced fabrication capabilities of the commercial electronics industry7,8. Yet, despite remarkable manufacturing advantages, reliance on silicon-based waveguides currently limits the spectral window available to photonic integrated circuits (PICs). Here, we present a new generation of integrated photonics by directly uniting III-V materials with silicon nitride waveguides on Si wafers. Using this technology, we present a fully integrated PIC at photon energies greater than the bandgap of silicon, demonstrating essential photonic building blocks, including lasers, amplifiers, photodetectors, modulators and passives, all operating at submicrometre wavelengths. Using this platform, we achieve unprecedented coherence and tunability in an integrated laser at short wavelength. Furthermore, by making use of this higher photon energy, we demonstrate superb high-temperature performance and kHz-level fundamental linewidths at elevated temperatures. Given the many potential applications at short wavelengths, the success of this integration strategy unlocks a broad range of new integrated photonics applications.

Construction of the single-diamond-structured titania scaffold-Recreation of the holy grail photonic structure.

As one of the most stunning biological nanostructures, the single-diamond (SD) surface discovered in beetles and weevils exoskeletons possesses the widest complete photonic bandgap known to date and is renowned as the "holy grail" of photonic materials. However, the synthesis of SD is difficult due to its thermodynamical instability compared to the energetically favoured bicontinuous double diamond and other easily formed lattices; thus, the artificial fabrication of SD has long been a formidable challenge. Herein, we report a bottom-up approach to fabricate SD titania networks via a one-pot cooperative assembly scenario employing the diblock copolymer poly(ethylene oxide)-block-polystyrene as a soft template and titanium diisopropoxide bis(acetylacetonate) as an inorganic precursor in a mixed solvent, in which the SD scaffold was obtained by kinetically controlled nucleation and growth in the skeletal channels of the diamond minimal surface formed by the polymer matrix. Electron crystallography investigations revealed the formation of tetrahedrally connected SD frameworks with the space group Fd [Formula: see text] m in a polycrystalline anatase form. A photonic bandgap calculation showed that the resulting SD structure has a wide and complete bandgap. This work solves the complex synthetic enigmas and offers a frontier in hyperbolic surfaces, biorelevant materials, next-generation optical devices, etc.

Chromatin remodeling analysis reveals the RdDM pathway responds to low-phosphorus stress in maize.

Chromatin in eukaryotes folds into a complex three-dimensional (3D) structure that is essential for controlling gene expression and cellular function and is dynamically regulated in biological processes. Studies on plant phosphorus signaling have concentrated on single genes and gene interactions. It is critical to expand the existing signaling pathway in terms of its 3D structure. In this study, low-Pi treatment led to greater chromatin volume. Furthermore, low-Pi stress increased the insulation score and the number of TAD-like domains, but the effects on the A/B compartment were not obvious. The methylation levels of target sites (hereafter as RdDM levels) peaked at specific TAD-like boundaries, whereas RdDM peak levels at conserved TAD-like boundaries shifted and decreased sharply. The distribution pattern of RdDM sites originating from the Helitron transposons matched that of genome-wide RdDM sites near TAD-like boundaries. RdDM pathway genes were upregulated in the middle or early stages and downregulated in the later stages under low-Pi conditions. The RdDM pathway mutant ddm1a showed increased tolerance to low-Pi stress, with shortened and thickened roots contributing to higher Pi uptake from the shallow soil layer. ChIP-seq results revealed that ZmDDM1A could bind to Pi- and root development-related genes. Strong associations were found between interacting genes in significantly different chromatin-interaction regions and root traits. These findings not only expand the mechanisms by which plants respond to low-Pi stress through the RdDM pathway but also offer a crucial framework for the analysis of biological issues using 3D genomics.

Fracture-driven power amplification in a hydrogel launcher.

Robotic tasks that require robust propulsion abilities such as jumping, ejecting or catapulting require power-amplification strategies where kinetic energy is generated from pre-stored energy. Here we report an engineered accumulated strain energy-fracture power-amplification method that is inspired by the pressurized fluidic squirting mechanism of Ecballium elaterium (squirting cucumber plants). We realize a light-driven hydrogel launcher that harnesses fast liquid vapourization triggered by the photothermal response of an embedded graphene suspension. This vapourization leads to appreciable elastic energy storage within the surrounding hydrogel network, followed by rapid elastic energy release within 0.3 ms. These soft hydrogel robots achieve controlled launching at high velocity with a predictable trajectory. The accumulated strain energy-fracture method was used to create an artificial squirting cucumber that disperses artificial seeds over metres, which can further achieve smart seeding through an integrated radio-frequency identification chip. This power-amplification strategy provides a basis for propulsive motion to advance the capabilities of miniaturized soft robotic systems.

Differentiation of Type 17 Mucosal-Associated Invariant T Cells in Circulation Contributes to the Severity of Sepsis.

Mucosal-associated invariant T (MAIT) cells are essential in defending against infection. Sepsis is a systemic inflammatory response to infection and a leading cause of death. The relationship between the overall competency of the host immune response and disease severity is not fully elucidated. This study identified a higher proportion of circulating MAIT17 with expression of IL-17A and retinoic acid receptor-related orphan receptor γt in patients with sepsis. The proportion of MAIT17 was correlated with the severity of sepsis. Single-cell RNA-sequencing analysis revealed an enhanced expression of lactate dehydrogenase A (LDHA) in MAIT17 in patients with sepsis. Cell-culture experiments demonstrated that phosphoinositide 3-kinase-LDHA signaling was required for retinoic acid receptor-related orphan receptor γt expression in MAIT17. Finally, the elevated levels of plasma IL-18 promoted the differentiation of circulating MAIT17 cells in sepsis. In summary, this study reveals a new role of circulating MAIT17 in promoting sepsis severity and suggests the phosphoinositide 3-kinase-LDHA signaling as a driving force in MAIT17 responses.

hsa_circ_0007919 promotes pancreatic cancer metastasis by modulating Sp1-mediated THBS1 transcription.

CircRNAs are abnormally expressed in various cancers and play an important role in the occurrence and development of cancers. However, their biological functions and the underlying molecular mechanisms in pancreatic cancer (PC) metastasis are incompletely understood. Differentially expressed circRNAs were identified by second-generation transcriptome sequencing in three pairs of PC tissues and adjacent tissues. The expression and prognostic significance of hsa_circ_0007919 were evaluated by qRT-PCR and Kaplan-Meier survival curves. Gain- and loss-of-function assays were conducted to detect the role of hsa_circ_0007919 in PC metastasis in vitro. A lung metastasis model and IHC experiments were conducted to confirm the effects of hsa_circ_0007919 on tumor metastasis in vivo. Mechanistically, RNA immunoprecipitation and chromatin immunoprecipitation assays were conducted to explore the interplay among hsa_circ_0007919, Sp1, and the THBS1 promoter. hsa_circ_0007919 was significantly upregulated in PC tissues and cells and was correlated with lymph node metastasis, TNM stage, and poor prognosis. Knockdown of hsa_circ_0007919 significantly suppressed the migration and invasion of PC cells in vitro and inhibited tumor metastasis in vivo. However, overexpression of hsa_circ_0007919 exerted the opposite effects. Mechanistically, hsa_circ_0007919 could recruit the transcription factor Sp1 to inhibit THBS1 transcription, thereby facilitating PC metastasis. hsa_circ_0007919 can promote the metastasis of PC by inhibiting THBS1 expression. hsa_circ_0007919 may be a potential therapeutic target in PC.

UBA5 inhibition restricts lung adenocarcinoma via blocking macrophage M2 polarization and cisplatin resistance.

UBA5, a ubiquitin-like activated enzyme involved in ufmylation and sumoylation, presents a viable target for pancreatic and breast cancer treatments, yet its role in lung adenocarcinoma (LUAD) remains underexplored. This study reveals UBA5's tumor-promoting effect in LUAD, as evidenced by its upregulation in patients and positive correlation with TNM stages. Elevated UBA5 levels predict poor outcomes for these patients. Pharmacological inhibition of UBA5 using DKM 2-93 significantly curtails the growth of A549, H1299, and cisplatin-resistant A549 (A549/DDP) LUAD cells in vitro. Additionally, UBA5 knockdown via shRNA lentivirus suppresses tumor growth both in vitro and in vivo. High UBA5 expression adversely alters the tumor immune microenvironment, affecting immunostimulators, MHC molecules, chemokines, receptors, and immune cell infiltration. Notably, UBA5 expression correlates positively with M2 macrophage infiltration, the predominant immune cells in LUAD. Co-culture experiments further demonstrate that UBA5 knockdown directly inhibits M2 macrophage polarization and lactate production in LUAD. Moreover, in vivo studies show reduced M2 macrophage infiltration following UBA5 knockdown. UBA5 expression is also associated with increased tumor heterogeneity, including tumor mutational burden, microsatellite instability, neoantigen presence, and homologous recombination deficiency. Experiments indicate that UBA5 overexpression promotes cisplatin resistance in vitro, whereas UBA5 inhibition enhances cisplatin sensitivity in both in vitro and in vivo settings. Overall, these findings suggest that targeting UBA5 inhibits LUAD by impeding cancer cell proliferation, M2 macrophage polarization, and cisplatin resistance.

Recruitment of transcription factor ETS1 to activated accessible regions promotes the transcriptional program of cilia genes.

Defects in cilia genes, which are critical for cilia formation and function, can cause complicated ciliopathy syndromes involving multiple organs and tissues; however, the underlying regulatory mechanisms of the networks of cilia genes in ciliopathies remain enigmatic. Herein, we have uncovered the genome-wide redistribution of accessible chromatin regions and extensive alterations of expression of cilia genes during Ellis-van Creveld syndrome (EVC) ciliopathy pathogenesis. Mechanistically, the distinct EVC ciliopathy-activated accessible regions (CAAs) are shown to positively regulate robust changes in flanking cilia genes, which are a key requirement for cilia transcription in response to developmental signals. Moreover, a single transcription factor, ETS1, can be recruited to CAAs, leading to prominent chromatin accessibility reconstruction in EVC ciliopathy patients. In zebrafish, the collapse of CAAs driven by ets1 suppression subsequently causes defective cilia proteins, resulting in body curvature and pericardial oedema. Our results depict a dynamic landscape of chromatin accessibility in EVC ciliopathy patients, and uncover an insightful role for ETS1 in controlling the global transcriptional program of cilia genes by reprogramming the widespread chromatin state.

Metal Clusters Confined in Chiral Zeolitic Imidazolate Framework for Circularly Polarized-Luminescence Inks.

Chiroptical activities arising in nanoclusters (NCs) are emerging as one of the most dynamic areas of modern science. However, devising an overarching strategy that is capable of concurrently enhancing the photoluminescence (PL) and circularly polarized luminescence (CPL) of metal NCs remains a formidable challenge. Herein, gold and silver nanoclusters (AuNCs, AgNCs) are endowed with CPL, for the first time, through a universal host-guest approach─centered around perturbing a chiral microenvironment within chiral hosts, simultaneously enhancing emissions. Remarkably, the photoluminescence quantum yield (PLQY) of AuNCs has undergone an increase of over 200 times upon confinement, escalating from 0.05% to 12%, and demonstrates a CPL response. Moreover, a three-dimensional (3D) model termed "NCs@CMOF" featuring CPL activity is created using metal cluster-based assembly inks through the process of 3D printing. This work introduces a potentially straightforward and versatile approach for achieving both PL enhancement and CPL activities in metal clusters.

Inhibition of UBA52 induces autophagy via EMC6 to suppress hepatocellular carcinoma tumorigenesis and progression.

Ubiquitin A-52 residue ribosomal protein fusion product 1 (UBA52) has a role in the occurrence and development of tumours. However, the mechanism by which UBA52 regulates hepatocellular carcinoma (HCC) tumorigenesis and progression remains poorly understood. By using the Cell Counting Kit (CCK-8), colony formation, wound healing and Transwell assays, we assessed the effects of UBA52 knockdown and overexpression on the proliferation and migration of HCC cells in vitro. By establishing subcutaneous and metastatic tumour models in nude mice, we evaluated the effects of UBA52 on HCC cell proliferation and migration in vivo. Through bioinformatic analysis of data from the Gene Expression Profiling Interactive Analysis (GEPIA) and The Cancer Genome Atlas (TCGA) databases, we discovered that UBA52 is associated with autophagy. In addition, we discovered that HCC tissues with high UBA52 expression had a poor prognosis in patients. Moreover, knockdown of UBA52 reduced HCC cell growth and metastasis both in vitro and in vivo. Mechanistically, knockdown of UBA52 induced autophagy through EMC6 in HCC cells. These findings suggest that UBA52 promoted the proliferation and migration of HCC cells through autophagy regulation via EMC6 and imply that UBA52 may be a viable novel treatment target for HCC patients.

A heterogeneous tumor immune microenvironment of uncommon epidermal growth factor receptor mutant non-small cell lung cancer.

Common epidermal growth factor receptor (EGFR) mutations are usually not considered for immunotherapy in non-small cell lung cancer (NSCLC) due to poor efficacy. However, whether uncommon EGFR mutations are suitable for immunotherapy has not been thoroughly studied. Thus, we explored the tumor immune microenvironment (TME) features in uncommon EGFR mutant NSCLC. In this study, a total of 41 patients with EGFR mutations were included, the majority (85.4%) of whom were stage I. Among them, 22 patients harbored common mutations, while 19 patients presented with uncommon mutations. Compared with common mutations, uncommon mutations exhibited more infiltrating T cells and fewer M2 macrophages, upregulated expression of antigen processing and a presentation pathway. Unsupervised clustering based on the mIF profile identified two classes with heterogeneous TME in uncommon mutations. Class 1 featured the absence of PD-1+ cytotoxic T cell infiltration, and class 2 displayed a hotter TME because of the downregulated expression of hypoxia (p < 0.001), oxidative phosphorylation (p = 0.009), and transforming growth factor beta signaling (p = 0.01) pathways as well as increased expression of CTLA4 (p = 0.001) and PDCD1 (p = 0.004). The association of CTLA4 and PDCD1 with TME profiles was validated in a TCGA lung adenocarcinoma cohort with uncommon EGFR mutations. Our study reveals the distinct and heterogeneous TME features in uncommon EGFR mutant NSCLC.

A Pt1Pd Single-Atom Alloy Nanozyme with Boosted Enzyme-Like Activity for Efficient Photo-Mediated Tumor Therapy.

Single-atom nanozymes (SAzymes) are emerging natural enzyme mimics and have attracted much attention in the biomedical field. SAzymes with Metal─Nx sites designed on carbon matrixes are currently the mainstream in research. It is of great significance to further expand the types of SAzymes to enrich the nanozyme library. Single-atom alloys (SAAs) are a material in which single-atom metal sites are dispersed onto another active metal matrix, and currently, there is limited research on their enzyme-like catalytic performance. In this work, a biodegradable Pt1Pd SAA is fabricated via a simple galvanic replacement strategy, and for the first time reveals its intrinsic enzyme-like catalytic performance including catalase-, oxidase-, and peroxidase-like activities, as well as its photodynamic effect. Experimental characterizations demonstrate that the introduction of single-atom Pt sites contributes to enhancing the affinity of Pt1Pd single-atom alloy nanozyme (SAAzyme) toward substrates, thus exhibiting boosted catalytic efficiency. In vitro and in vivo experiments demonstrate that Pt1Pd SAAzyme exhibits a photo-controlled therapeutic effect, with a tumor inhibition rate of up to 100%. This work provides vital guidance for opening the research direction of SAAs in enzyme-like catalysis.

Multiple strategies, including 6mA methylation, affecting plant alternative splicing in allopolyploid peanut.

Alternative splicing (AS), an important post-transcriptional regulation mechanism in eukaryotes, can significantly increase transcript diversity and contribute to gene expression regulation and many other complicated developmental processes. While plant gene AS events are well described, few studies have investigated the comprehensive regulation machinery of plant AS. Here, we use multi-omics to analyse peanut AS events. Using long-read isoform sequencing, 146 464 full-length non-chimeric transcripts were obtained, resulting in annotation corrections for 1782 genes and the identification of 4653 new loci. Using Iso-Seq RNA sequences, 271 776 unique splice junctions were identified, 82.49% of which were supported by transcriptome data. We characterized 50 977 polyadenylation sites for 23 262 genes, 12 369 of which had alternative polyadenylation sites. AS allows differential regulation of the same gene by miRNAs at the isoform level coupled with polyadenylation. In addition, we identified many long non-coding RNAs and fusion transcripts. There is a suppressed effect of 6mA on AS and gene expression. By analysis of chromatin structures, the genes located in the boundaries of topologically associated domains, proximal chromosomal telomere regions, inter- or intra-chromosomal loops were found to have more unique splice isoforms, higher expression, lower 6mA and more transposable elements (TEs) in their gene bodies than the other genes, indicating that chromatin interaction, 6mA and TEs play important roles in AS and gene expression. These results greatly refine the peanut genome annotation and contribute to the study of gene expression and regulation in peanuts. This work also showed AS is associated with multiple strategies for gene regulation.

Longitudinal Fluctuations in Treatment Response After OnabotulinumToxinA and Sacral Neuromodulation for Refractory Urgency Incontinence.

PURPOSE: We compared fluctuations in treatment response after onabotulinumtoxinA and sacral neuromodulation for urgency incontinence using Markov models. MATERIALS AND METHODS: We fit data from a randomized trial to Markov models to compare transitions of success/failure over 6 months between 200 U onabotulinumtoxinA and sacral neuromodulation. Objective failure was <50% reduction in urgency incontinence episodes from baseline; subjective failure "strongly disagree" to "neutral" to the Patient Global Symptom Control questionnaire. RESULTS: Of the 357 participants (median baseline daily urgency incontinence episodes 4.7 [IQR 3.7-6.0]) 61% vs 51% and 3.2% vs 6.1% reported persistent states of objective success and failure over 6 months after onabotulinumtoxinA vs sacral neuromodulation. Participants receiving onabotulinumtoxinA vs sacral neuromodulation had lower 30-day transition probabilities from objective and subjective success to failure (10% vs 14%, ratio 0.75 [95% CI 0.55-0.95]; 14% vs 21%, ratio 0.70 [95% CI 0.51-0.89]). The 30-day transition probability from objective and subjective failure to success did not differ between onabotulinumtoxinA and sacral neuromodulation (40% vs 36%, ratio 1.11 [95% CI 0.73-1.50]; 18% vs 17%, ratio 1.14 [95% CI 0.65-1.64]). CONCLUSIONS: Over 6 months after treatment, 2 in 5 women's symptoms fluctuate. Within these initial 6 months, women receiving onabotulinumtoxinA transitioned from success to failure over 30 days less often than sacral neuromodulation. For both treatments, there was an almost 20%-40% probability over 30 days that women returned to subjective and objective success after failure. Markov models add important information to longitudinal models on how symptoms fluctuate after urgency incontinence treatment.

Lithium chloride promotes mesenchymal-epithelial transition in murine cutaneous wound healing via inhibiting CXCL9 and IGF2.

Cutaneous wound healing is a challenge in plastic and reconstructive surgery. In theory, cells undergoing mesenchymal transition will achieve re-epithelialization through mesenchymal-epithelial transition at the end of wound healing. But in fact, some pathological stimuli will inhibit this biological process and result in scar formation. If mesenchymal-epithelial transition can be activated at the corresponding stage, the ideal wound healing may be accomplished. Two in vivo skin defect mouse models and dermal-derived mesenchymal cells were used to evaluate the effect of lithium chloride in wound healing. The mesenchymal-epithelial transition was detected by immunohistochemistry staining. In vivo, differentially expressed genes were analysed by transcriptome analyses and the subsequent testing was carried out. We found that lithium chloride could promote murine cutaneous wound healing and facilitate mesenchymal-epithelial transition in vivo and in vitro. In lithium chloride group, scar area was smaller and the collagen fibres are also orderly arranged. The genes related to mesenchyme were downregulated and epithelial mark genes were activated after intervention. Moreover, transcriptome analyses suggested that this effect might be related to the inhibition of CXCL9 and IGF2, subsequent assays demonstrated it. Lithium chloride can promote mesenchymal-epithelial transition via downregulating CXCL9 and IGF2 in murine cutaneous wound healing, the expression of IGF2 is regulated by ß-catenin. It may be a potential promising therapeutic drug for alleviating postoperative scar and promoting re-epithelialization in future.

Designing future peanut: the power of genomics-assisted breeding.

KEY MESSAGE: Integrating GAB methods with high-throughput phenotyping, genome editing, and speed breeding hold great potential in designing future smart peanut cultivars to meet market and food supply demands. Cultivated peanut (Arachis hypogaea L.), a legume crop greatly valued for its nourishing food, cooking oil, and fodder, is extensively grown worldwide. Despite decades of classical breeding efforts, the actual on-farm yield of peanut remains below its potential productivity due to the complicated interplay of genotype, environment, and management factors, as well as their intricate interactions. Integrating modern genomics tools into crop breeding is necessary to fast-track breeding efficiency and rapid progress. When combined with speed breeding methods, this integration can substantially accelerate the breeding process, leading to faster access of improved varieties to farmers. Availability of high-quality reference genomes for wild diploid progenitors and cultivated peanuts has accelerated the process of gene/quantitative locus discovery, developing markers and genotyping assays as well as a few molecular breeding products with improved resistance and oil quality. The use of new breeding tools, e.g., genomic selection, haplotype-based breeding, speed breeding, high-throughput phenotyping, and genome editing, is probable to boost genetic gains in peanut. Moreover, renewed attention to efficient selection and exploitation of targeted genetic resources is also needed to design high-quality and high-yielding peanut cultivars with main adaptation attributes. In this context, the combination of genomics-assisted breeding (GAB), genome editing, and speed breeding hold great potential in designing future improved peanut cultivars to meet market and food supply demands.

Genome-wide association studies dissect low-phosphorus stress response genes underling field and seedling traits in maize.

Phosphorus (P) is an essential element for plant growth, and its deficiency can cause decreased crop yield. This study systematically evaluated the low-phosphate (Pi) response traits in a large population at maturity and seedling stages, and explored candidate genes and their interrelationships with specific traits. The results revealed a greater sensitivity of seedling maize to low-Pi stress compared to that at maturity stage. The phenotypic response patterns to low-Pi stress at different stages were independent. Chlorophyll content was found to be a potential indicator for screening low-Pi-tolerant materials in the field. A total of 2900 and 1446 significantly associated genes at the maturity and seedling stages were identified, respectively. Among these genes, 972 were uniquely associated with maturity traits, while 330 were specifically detected at the seedling stage under low-Pi stress. Moreover, 768 and 733 genes were specifically associated with index values (low-Pi trait/normal-Pi trait) at maturity and seedling stage, respectively. Genetic network diagrams showed that the low-Pi response gene Zm00001d022226 was specifically associated with multiple primary P-related traits under low-Pi conditions. A total of 963 out of 2966 genes specifically associated with traits under low-Pi conditions or index values were found to be induced by low-Pi stress. Notably, ZmSPX4.1 and ZmSPX2 were sharply up-regulated in response to low-Pi stress across different lines or tissues. These findings advance our understanding of maize's response to low-Pi stress at different developmental stages, shedding light on the genes and pathways implicated in this response.

Mining for QTL controlling maize low-phosphorus response genes combined with deep resequencing of RIL parental genomes and in silico GWAS analysis.

KEY MESSAGE: Extensive and comprehensive phenotypic data from a maize RIL population under both low- and normal-Pi treatments were used to conduct QTL mapping. Additionally, we integrated parental resequencing data from the RIL population, GWAS results, and transcriptome data to identify candidate genes associated with low-Pi stress in maize. Phosphorus (Pi) is one of the essential nutrients that greatly affect the maize yield. However, the genes underlying the QTL controlling maize low-Pi response remain largely unknown. In this study, a total of 38 traits at both seedling and maturity stages were evaluated under low- and normal-Pi conditions using a RIL population constructed from X178 (tolerant) and 9782 (sensitive), and most traits varied significantly between low- and normal-Pi treatments. Twenty-nine QTLs specific to low-Pi conditions were identified after excluding those with common intervals under both low- and normal-Pi conditions. Furthermore, 45 additional QTLs were identified based on the index value ((Trait_under_LowPi-Trait_under_NormalPi)/Trait_under_NormalPi) of each trait. These 74 QTLs collectively were classified as Pi-dependent QTLs. Additionally, 39 Pi-dependent QTLs were clustered in nine HotspotQTLs. The Pi-dependent QTL interval contained 19,613 unique genes, 6,999 of which exhibited sequence differences with non-synonymous mutation sites between X178 and 9782. Combined with in silico GWAS results, 277 consistent candidate genes were identified, with 124 genes located within the HotspotQTL intervals. The transcriptome analysis revealed that 21 genes, including the Pi transporter ZmPT7 and the strigolactones pathway-related gene ZmPDR1, exhibited consistent low-Pi stress response patterns across various maize inbred lines or tissues. It is noteworthy that ZmPDR1 in maize roots can be sharply up-regulated by low-Pi stress, suggesting its potential importance as a candidate gene for responding to low-Pi stress through the strigolactones pathway.