High-throughput characterization of functional variants highlights heterogeneity and polygenicity underlying lung cancer susceptibility.

Genome-wide association studies (GWASs) have identified numerous lung cancer risk-associated loci. However, decoding molecular mechanisms of these associations is challenging since most of these genetic variants are non-protein-coding with unknown function. Here, we implemented massively parallel reporter assays (MPRAs) to simultaneously measure the allelic transcriptional activity of risk-associated variants. We tested 2,245 variants at 42 loci from 3 recent GWASs in East Asian and European populations in the context of two major lung cancer histological types and exposure to benzo(a)pyrene. This MPRA approach identified one or more variants (median 11 variants) with significant effects on transcriptional activity at 88% of GWAS loci. Multimodal integration of lung-specific epigenomic data demonstrated that 63% of the loci harbored multiple potentially functional variants in linkage disequilibrium. While 22% of the significant variants showed allelic effects in both A549 (adenocarcinoma) and H520 (squamous cell carcinoma) cell lines, a subset of the functional variants displayed a significant cell-type interaction. Transcription factor analyses nominated potential regulators of the functional variants, including those with cell-type-specific expression and those predicted to bind multiple potentially functional variants across the GWAS loci. Linking functional variants to target genes based on four complementary approaches identified candidate susceptibility genes, including those affecting lung cancer cell growth. CRISPR interference of the top functional variant at 20q13.33 validated variant-to-gene connections, including RTEL1, SOX18, and ARFRP1. Our data provide a comprehensive functional analysis of lung cancer GWAS loci and help elucidate the molecular basis of heterogeneity and polygenicity underlying lung cancer susceptibility.

Genetic susceptibility to chronic diseases leads to heart failure among Europeans: the influence of leukocyte telomere length.

BACKGROUND: Genetic susceptibility to various chronic diseases has been shown to influence heart failure (HF) risk. However, the underlying biological pathways, particularly the role of leukocyte telomere length (LTL), are largely unknown. We investigated the impact of genetic susceptibility to chronic diseases and various traits on HF risk, and whether LTL mediates or modifies the pathways. METHODS: We conducted prospective cohort analyses on 404 883 European participants from the UK Biobank, including 9989 incident HF cases. Multivariable Cox regression was used to estimate associations between HF risk and 24 polygenic risk scores (PRSs) for various diseases or traits previously generated using a Bayesian approach. We assessed multiplicative interactions between the PRSs and LTL previously measured in the UK Biobank using quantitative PCR. Causal mediation analyses were conducted to estimate the proportion of the total effect of PRSs acting indirectly through LTL, an integrative marker of biological aging. RESULTS: We identified 9 PRSs associated with HF risk, including those for various cardiovascular diseases or traits, rheumatoid arthritis (P = 1.3E-04), and asthma (P = 1.8E-08). Additionally, longer LTL was strongly associated with decreased HF risk (P-trend = 1.7E-08). Notably, LTL strengthened the asthma-HF relationship significantly (P-interaction = 2.8E-03). However, LTL mediated only 1.13% (P < 0.001) of the total effect of the asthma PRS on HF risk. CONCLUSIONS: Our findings shed light onto the shared genetic susceptibility between HF risk, asthma, rheumatoid arthritis, and other traits. Longer LTL strengthened the genetic effect of asthma in the pathway to HF. These results support consideration of LTL and PRSs in HF risk prediction.

Estimating the overall fraction of phenotypic variance attributed to high-dimensional predictors measured with error.

In prospective genomic studies (e.g., DNA methylation, metagenomics, and transcriptomics), it is crucial to estimate the overall fraction of phenotypic variance (OFPV) attributed to the high-dimensional genomic variables, a concept similar to heritability analyses in genome-wide association studies (GWAS). Unlike genetic variants in GWAS, these genomic variables are typically measured with error due to technical limitation and temporal instability. While the existing methods developed for GWAS can be used, ignoring measurement error may severely underestimate OFPV and mislead the design of future studies. Assuming that measurement error variances are distributed similarly between causal and noncausal variables, we show that the asymptotic attenuation factor equals to the average intraclass correlation coefficients of all genomic variables, which can be estimated based on a pilot study with repeated measurements. We illustrate the method by estimating the contribution of microbiome taxa to body mass index and multiple allergy traits in the American Gut Project. Finally, we show that measurement error does not cause meaningful bias when estimating the correlation of effect sizes for two traits.

Low-intensity ultrasound ameliorates brain organoid integration and rescues microcephaly deficits.

Human brain organoids represent a remarkable platform for modeling neurological disorders and a promising brain repair approach. However, the effects of physical stimulation on their development and integration remain unclear. Here, we report that low-intensity ultrasound significantly increases neural progenitor cell proliferation and neuronal maturation in cortical organoids. Histological assays and single-cell gene expression analyses reveal that low-intensity ultrasound improves the neural development in cortical organoids. Following organoid grafts transplantation into the injured somatosensory cortices of adult mice, longitudinal electrophysiological recordings and histological assays reveal that ultrasound-treated organoid grafts undergo advanced maturation. They also exhibit enhanced pain-related gamma-band activity and more disseminated projections into the host brain than the untreated groups. Finally, low-intensity ultrasound ameliorates neuropathological deficits in a microcephaly brain organoid model. Hence, low-intensity ultrasound stimulation advances the development and integration of brain organoids, providing a strategy for treating neurodevelopmental disorders and repairing cortical damage.

Association of oral microbiome with oral human papillomavirus infection: a population study of the National Health and Nutrition Examination Survey, 2009-2012.

BACKGROUND: Oral human papillomavirus(HPV) infection and the oral microbiome are associated with oropharyngeal cancer. However, population-based data on the association of oral microbiome with oral HPV infection are limited. METHOD: We performed a cross-sectional analysis of 5,496 participants aged 20-59 in National Health and Nutrition Examination Surveys(NHANES):2009-2012. The association between either oral microbiome alpha diversity or beta diversity and oral HPV infection was assessed using multivariable logistic regression or principal coordinate analyses(PCoA) and multivariate analysis of variance(PERMANOVA). RESULTS: For alpha diversity, we found a lower number of observed Amplicon sequence variants(ASVs) (adjusted odds ratio[aOR] = 0.996; 95%CI = 0.992-0.999) and reduced Faith's Phylogenetic Diversity(aOR = 0.95; 95%CI = 0.90-0.99) associated with high-risk oral HPV infection in the overall population. This trend was observed in males for both high-risk and any oral HPV infection. Beta diversity showed differentiation of oral microbiome community by high-risk oral HPV infection as measured by Bray-Curtis dissimilarity (R2 = 0.054%; P = .029) and unweighted UniFrac distance (R2 = 0.046%; P = .045) among the overall population, and associations were driven by males. CONCLUSIONS: Both oral microbiome alpha diversity(within-sample richness and phylogenetic diversity) and beta diversity(heterogeneous dispersion of oral microbiome community) are associated with HPV infection. Longitudinal studies are needed to characterize the role of the microbiome in the natural history of oral HPV infection.

Rare germline deleterious variants increase susceptibility for lung cancer.

Although multiple common susceptibility loci for lung cancer (LC) have been identified by genome-wide association studies, they can explain only a small portion of heritability. The etiological contribution of rare deleterious variants (RDVs) to LC risk is not fully characterized and may account for part of the missing heritability. Here, we sequenced the whole exomes of 2777 participants from the Environment and Genetics in Lung cancer Etiology study, a homogenous population including 1461 LC cases and 1316 controls. In single-variant analyses, we identified a new RDV, rs77187983 [EHBP1, odds ratio (OR) = 3.13, 95% confidence interval (CI) = 1.34-7.30, P = 0.008] and replicated two previously reported RDVs, rs11571833 (BRCA2, OR = 2.18; 95% CI = 1.25-3.81, P = 0.006) and rs752672077 (MPZL2, OR = 3.70, 95% CI = 1.04-13.15, P = 0.044). In gene-based analyses, we confirmed BRCA2 (P = 0.007) and ATM (P = 0.014) associations with LC risk and identified TRIB3 (P = 0.009), involved in maintaining genome stability and DNA repair, as a new candidate susceptibility gene. Furthermore, cases were enriched with RDVs in homologous recombination repair [carrier frequency (CF) = 22.9% versus 19.5%, P = 0.017] and Fanconi anemia (CF = 12.5% versus 10.2%, P = 0.036) pathways. Our results were not significant after multiple testing corrections but were enriched in cases versus controls from large scale public biobank resources, including The Cancer Genome Atlas, FinnGen and UK Biobank. Our study identifies novel candidate genes and highlights the importance of RDVs in DNA repair-related genes for LC susceptibility. These findings improve our understanding of LC heritability and may contribute to the development of risk stratification and prevention strategies.

Cell-type-specific meQTLs extend melanoma GWAS annotation beyond eQTLs and inform melanocyte gene-regulatory mechanisms.

Although expression quantitative trait loci (eQTLs) have been powerful in identifying susceptibility genes from genome-wide association study (GWAS) findings, most trait-associated loci are not explained by eQTLs alone. Alternative QTLs, including DNA methylation QTLs (meQTLs), are emerging, but cell-type-specific meQTLs using cells of disease origin have been lacking. Here, we established an meQTL dataset by using primary melanocytes from 106 individuals and identified 1,497,502 significant cis-meQTLs. Multi-QTL colocalization with meQTLs, eQTLs, and mRNA splice-junction QTLs from the same individuals together with imputed methylome-wide and transcriptome-wide association studies identified candidate susceptibility genes at 63% of melanoma GWAS loci. Among the three molecular QTLs, meQTLs were the single largest contributor. To compare melanocyte meQTLs with those from malignant melanomas, we performed meQTL analysis on skin cutaneous melanomas from The Cancer Genome Atlas (n = 444). A substantial proportion of meQTL probes (45.9%) in primary melanocytes is preserved in melanomas, while a smaller fraction of eQTL genes is preserved (12.7%). Integration of melanocyte multi-QTLs and melanoma meQTLs identified candidate susceptibility genes at 72% of melanoma GWAS loci. Beyond GWAS annotation, meQTL-eQTL colocalization in melanocytes suggested that 841 unique genes potentially share a causal variant with a nearby methylation probe in melanocytes. Finally, melanocyte trans-meQTLs identified a hotspot for rs12203592, a cis-eQTL of a transcription factor, IRF4, with 131 candidate target CpGs. Motif enrichment and IRF4 ChIP-seq analysis demonstrated that these target CpGs are enriched in IRF4 binding sites, suggesting an IRF4-mediated regulatory network. Our study highlights the utility of cell-type-specific meQTLs.

A 584 bp deletion in CTRB2 inhibits chymotrypsin B2 activity and secretion and confers risk of pancreatic cancer.

Genome-wide association studies (GWASs) have discovered 20 risk loci in the human genome where germline variants associate with risk of pancreatic ductal adenocarcinoma (PDAC) in populations of European ancestry. Here, we fine-mapped one such locus on chr16q23.1 (rs72802365, p = 2.51 × 10-17, OR = 1.36, 95% CI = 1.31-1.40) and identified colocalization (PP = 0.87) with aberrant exon 5-7 CTRB2 splicing in pancreatic tissues (pGTEx = 1.40 × 10-69, ßGTEx = 1.99; pLTG = 1.02 × 10-30, ßLTG = 1.99). Imputation of a 584 bp structural variant overlapping exon 6 of CTRB2 into the GWAS datasets resulted in a highly significant association with pancreatic cancer risk (p = 2.83 × 10-16, OR = 1.36, 95% CI = 1.31-1.42), indicating that it may underlie this signal. Exon skipping attributable to the deletion (risk) allele introduces a premature stop codon in exon 7 of CTRB2, yielding a truncated chymotrypsinogen B2 protein that lacks chymotrypsin activity, is poorly secreted, and accumulates intracellularly in the endoplasmic reticulum (ER). We propose that intracellular accumulation of a nonfunctional chymotrypsinogen B2 protein leads to ER stress and pancreatic inflammation, which may explain the increased pancreatic cancer risk in carriers of CTRB2 exon 6 deletion alleles.

Phenotypic and genetically predicted leucocyte telomere length and lung cancer risk in the prospective UK Biobank.

We investigated phenotypic leucocyte telomere length (LTL), genetically predicted LTL (gTL), and lung cancer risk among 371 890 participants, including 2829 incident cases, from the UK Biobank. Using multivariable Cox regression, we found dose-response relationships between longer phenotypic LTL (p-trendcontinuous=2.6×10-5), longer gTL predicted using a polygenic score with 130 genetic instruments (p-trendcontinuous=4.2×10-10), and overall lung cancer risk, particularly for adenocarcinoma. The associations were prominent among never smokers. Mendelian Randomization analyses supported causal associations between longer telomere length and lung cancer (HRper 1 SD gTL=1.87, 95% CI: 1.49 to 2.36, p=4.0×10-7), particularly adenocarcinoma (HRper 1 SD gTL=2.45, 95%CI: 1.69 to 3.57, p=6.5×10-6).

Epigenome-wide association study of lung cancer among never smokers in two prospective cohorts in Shanghai, China.

BACKGROUND: The aetiology of lung cancer among individuals who never smoked remains elusive, despite 15% of lung cancer cases in men and 53% in women worldwide being unrelated to smoking. Epigenetic alterations, particularly DNA methylation (DNAm) changes, have emerged as potential drivers. Yet, few prospective epigenome-wide association studies (EWAS), primarily focusing on peripheral blood DNAm with limited representation of never smokers, have been conducted. METHODS: We conducted a nested case-control study of 80 never-smoking incident lung cancer cases and 83 never-smoking controls within the Shanghai Women's Health Study and Shanghai Men's Health Study. DNAm was measured in prediagnostic oral rinse samples using Illumina MethylationEPIC array. Initially, we conducted an EWAS to identify differentially methylated positions (DMPs) associated with lung cancer in the discovery sample of 101 subjects. The top 50 DMPs were further evaluated in a replication sample of 62 subjects, and results were pooled using fixed-effect meta-analysis. RESULTS: Our study identified three DMPs significantly associated with lung cancer at the epigenome-wide significance level of p<8.22×10-8. These DMPs were identified as cg09198866 (MYH9; TXN2), cg01411366 (SLC9A10) and cg12787323. Furthermore, examination of the top 1000 DMPs indicated significant enrichment in epithelial regulatory regions and their involvement in small GTPase-mediated signal transduction pathways. Additionally, GrimAge acceleration was identified as a risk factor for lung cancer (OR=1.19 per year; 95% CI 1.06 to 1.34). CONCLUSIONS: While replication in a larger sample size is necessary, our findings suggest that DNAm patterns in prediagnostic oral rinse samples could provide novel insights into the underlying mechanisms of lung cancer in never smokers.

A promoter polymorphism defines distinct roles in anther development for Col-0 and Ler-0 alleles of Arabidopsis ACYL-COA BINDING PROTEIN3.

Acyl-CoA-Binding Proteins (ACBPs) bind acyl-CoA esters and function in lipid metabolism. Although acbp3-1, the ACBP3 mutant in Arabidopsis thaliana ecotype Col-0, displays normal floral development, the acbp3-2 mutant from ecotype Ler-0 characterized herein exhibits defective adaxial anther lobes and improper sporocyte formation. To understand these differences and identify the role of ERECTA in ACBP3 function, the acbp3 mutants and acbp3-erecta (er) lines were analyzed by microscopy for anther morphology and high-performance liquid chromatography for lipid composition. Defects in Landsberg anther development were related to the ERECTA-mediated pathway because the progenies of acbp3-2 × La-0 and acbp3-1 × er-1 in Col-0 showed normal anthers, contrasting to that of acbp3-2 in Ler-0. Polymorphism in the regulatory region of ACBP3 enabled its function in anther development in Ler-0 but not Col-0 which harbored an AT-repeat insertion. ACBP3 expression and anther development in acbp3-2 were restored using ACBP3pro (Ler)::ACBP3 not ACBP3pro (Col)::ACBP3. SPOROCYTELESS (SPL), a sporocyte formation regulator activated ACBP3 transcription in Ler-0 but not Col-0. For anther development, the ERECTA-related role of ACBP3 is required in Ler-0, but not Col-0. The disrupted promoter regulatory region for SPL binding in Col-0 eliminates the role of ACBP3 in anther development.

Boosting Triticeae crop grain yield by manipulating molecular modules to regulate inflorescence architecture: insights and knowledge from other cereal crops.

One of the challenges for global food security is to reliably and sustainably improve the grain yield of cereal crops. One solution is to modify the architecture of the grain-bearing inflorescence to optimize for grain number and size. Cereal inflorescences are complex structures, with determinacy, branching patterns, and spikelet/floret growth patterns that vary by species. Recent decades have witnessed rapid advancements in our understanding of the genetic regulation of inflorescence architecture in rice, maize, wheat, and barley. Here, we summarize current knowledge on key genetic factors underlying the different inflorescence morphologies of these crops and model plants (Arabidopsis and tomato), focusing particularly on the regulation of inflorescence meristem determinacy and spikelet meristem identity and determinacy. We also discuss strategies to identify and utilize these superior alleles to optimize inflorescence architecture and, ultimately, improve crop grain yield.

Real-time, random-access organ screening for carbapenem-resistant organisms (CRO) reduces CRO-associated, donor-derived infection mortality in lung transplant recipients.

PURPOSE: Donor-derived infection (DDI) has become an important factor affecting the prognosis of lung transplantation patients. The risks versus benefits of using donor organs infected with multidrug-resistant organisms (MDRO), especially carbapenem-resistant organisms (CRO), are frequently debated. Traditional microbial culture and antimicrobial susceptibility testing at present fail to meet the needs of quick CRO determination for donor lungs before acquisition. In this study, we explored a novel screening method by using Xpert® Carba-R assay for CRO in donor lungs in a real-time manner to reduce CRO-associated DDI mortality. METHODS: This study was registered on chictr.org.cn (ChiCTR2100053687) on November 2021. In the Xpert Carba-R screening group, donor lungs were screened for CRO infection by the Xpert Carba-R test on bronchoalveolar fluid (BALF) before acquisition. If the result was negative, donor lung acquisition and subsequent lung transplantation were performed. In the thirty-five potential donors, nine (25.71%) with positive Xpert Carba-R results in BALF were declined for lung transplantation. Twenty-six recipients and the matching CRO-negative donor lungs (74.29%) were included in the Xpert Carba-R screening group. In the control group, nineteen recipients underwent lung transplants without Xpert Carba-R screening. The incidence and mortality of CRO-associated DDI were collected and contrasted between the two groups. RESULTS: Multivariate analysis showed that CRO-related death due to DDI within 60 days was significantly lower in the Xpert Carba-R screening group than that in the control group (OR = 0.05, 95% CI 0.003-0.74, p = 0.03). CONCLUSION: Real-time CRO screening of donor lungs before transplantation at the point of care by the Xpert Carba-R helps clinicians formulate lung transplantation strategies quickly and reduces the risk of subsequent CRO infection improving the prognosis of lung transplantation.

Genomic variation in 3,010 diverse accessions of Asian cultivated rice.

Here we analyse genetic variation, population structure and diversity among 3,010 diverse Asian cultivated rice (Oryza sativa L.) genomes from the 3,000 Rice Genomes Project. Our results are consistent with the five major groups previously recognized, but also suggest several unreported subpopulations that correlate with geographic location. We identified 29 million single nucleotide polymorphisms, 2.4 million small indels and over 90,000 structural variations that contribute to within- and between-population variation. Using pan-genome analyses, we identified more than 10,000 novel full-length protein-coding genes and a high number of presence-absence variations. The complex patterns of introgression observed in domestication genes are consistent with multiple independent rice domestication events. The public availability of data from the 3,000 Rice Genomes Project provides a resource for rice genomics research and breeding.

Effect of black cloth ointment on hypertrophic scar formation: An investigation using integrated network pharmacology and animal assay.

BACKGROUND: Hypertrophic scars (HS) are a common disfiguring condition in daily clinical encounters which brings a lot of anxieties and concerns to patients, but the treatment options of HS are limited. Black cloth ointment (BCO), as a cosmetic ointment applicable to facial scars, has shown promising therapeutic effects for facial scarring. However, the molecular mechanisms underlying its therapeutic effects remain unclear. MATERIAL AND METHODS: Network pharmacology was first applied to analyze the major active components of BCO and the related signaling pathways. Subsequently, rabbit ear scar model was successfully established to determine the pharmacological effects of BCO and its active component ß-elemene on HS. Finally, the molecular mechanism of BCO and ß-elemene was analyzed by Western blot. RESULTS: Through the network pharmacology, it showed that ß-elemene was the main active ingredient of BCO, and it could significantly improve the pathological structure of HS and reduce collagen deposition. BCO and ß-elemene could increase the expression of ER stress-related markers and promote the increase of apoptotic proteins in the Western blot experiment and induce the apoptosis of myofibroblasts. CONCLUSIONS: Our findings indicate that the material basis for the scar-improving effects of the BCO is ß-elemene, and cellular apoptosis is the key mechanism through which the BCO and ß-elemene exert their effects.

Adaptive Weighted Data Fusion for Line Structured Light and Photometric Stereo Measurement System.

Line structured light (LSL) measurement systems can obtain high accuracy profiles, but the overall clarity relies greatly on the sampling interval of the scanning process. Photometric stereo (PS), on the other hand, is sensitive to tiny features but has poor geometrical accuracy. Cooperative measurement with these two methods is an effective way to ensure precision and clarity results. In this paper, an LSL-PS cooperative measurement system is brought out. The calibration methods used in the LSL and PS measurement system are given. Then, a data fusion algorithm with adaptive weights is proposed, where an error function that contains the 3D point cloud matching error and normal vector error is established. The weights, which are based on the angles of adjacent normal vectors, are also added to the error function. Afterward, the fusion results can be obtained by solving linear equations. From the experimental results, it can be seen that the proposed method has the advantages of both the LSL and PS methods. The 3D reconstruction results have the merits of high accuracy and high clarity.

Prospects for rice in 2050.

A key to achieve the goals put forward in the UN's 2030 Agenda for Sustainable Development, it will need transformative change to our agrifood systems. We must mount to the global challenge to achieve food security in a sustainable manner in the context of climate change, population growth, urbanization, and depletion of natural resources. Rice is one of the major staple cereal crops that has contributed, is contributing, and will still contribute to the global food security. To date, rice yield has held pace with increasing demands, due to advances in both fundamental and biological studies, as well as genomic and molecular breeding practices. However, future rice production depends largely on the planting of resilient cultivars that can acclimate and adapt to changing environmental conditions. This Special Issue highlight with reviews and original research articles the exciting and growing field of rice-environment interactions that could benefit future rice breeding. We also outline open questions and propose future directions of 2050 rice research, calling for more attentions to develop environment-resilient rice especially hybrid rice, upland rice and perennial rice.

Combined genome-wide association study and epistasis analysis reveal multifaceted genetic architectures of plant height in Asian cultivated rice.

Plant height (PH) in rice (Oryza sativa) is an important trait for its adaptation and agricultural performance. Discovery of the semi-dwarf1 (SD1) mutation initiated the Green Revolution, boosting rice yield and fitness, but the underlying genetic regulation of PH in rice remains largely unknown. Here, we performed genome-wide association study (GWAS) and identified 12 non-repetitive QTL/genes regulating PH variation in 619 Asian cultivated rice accessions. One of these was an SD1 structural variant, not normally detected in standard GWAS analyses. Given the strong effect of SD1 on PH, we also divided 619 accessions into subgroups harbouring distinct SD1 haplotypes, and found a further 85 QTL/genes for PH, revealing genetic heterogeneity that may be missed by analysing a broad, diverse population. Moreover, we uncovered two epistatic interaction networks of PH-associated QTL/genes in the japonica (Geng)-dominant SD1NIP subgroup. In one of them, the hub QTL/gene qphSN1.4/GAMYB interacted with qphSN3.1/OsINO80, qphSN3.4/HD16/EL1, qphSN6.2/LOC_Os06g11130, and qphSN10.2/MADS56. Sequence variations in GAMYB and MADS56 were associated with their expression levels and PH variations, and MADS56 was shown to physically interact with MADS57 to coregulate expression of gibberellin (GA) metabolic genes OsGA2ox3 and Elongated Uppermost Internode1 (EUI1). Our study uncovered the multifaceted genetic architectures of rice PH, and provided novel and abundant genetic resources for breeding semi-dwarf rice and new candidates for further mechanistic studies on regulation of PH in rice.

Genetic architecture of seed glycerolipids in Asian cultivated rice.

Glycerolipids are essential for rice development and grain quality but its genetic regulation remains unknown. Here we report its genetic base using metabolite-based genome-wide association study and metabolite-based quantitative traits locus (QTL) analyses based on lipidomic profiles of seeds from 587 Asian cultivated rice accessions and 103 chromosomal segment substitution lines, respectively. We found that two genes encoding phosphatidylcholine (PC):diacylglycerol cholinephosphotransferase (OsLP1) and granule-bound starch synthase I (Waxy) contribute to variations in saturated triacylglycerol (TAG) and lyso-PC contents, respectively. We demonstrated that allelic variation in OsLP1 sequence between indica and japonica results in different enzymatic preference for substrate PC-16:0/16:0 and different saturated TAG levels. Further evidence demonstrated that OsLP1 also affects heading date, and that co-selection of OsLP1 and a flooding-tolerant QTL in Aus results in the abundance of saturated TAGs associated with flooding tolerance. Moreover, we revealed that the sequence polymorphisms in Waxy has pleiotropic effects on lyso-PC and amylose content. We proposed that rice seed glycerolipids have been unintentionally shaped during natural and artificial selection for adaptive or import seed quality traits. Collectively, our findings provide valuable genetic resources for rice improvement and evolutionary insights into seed glycerolipid variations in rice.

SUITOR: Selecting the number of mutational signatures through cross-validation.

For de novo mutational signature analysis, the critical first step is to decide how many signatures should be expected in a cancer genomics study. An incorrect number could mislead downstream analyses. Here we present SUITOR (Selecting the nUmber of mutatIonal signaTures thrOugh cRoss-validation), an unsupervised cross-validation method that requires little assumptions and no numerical approximations to select the optimal number of signatures without overfitting the data. In vitro studies and in silico simulations demonstrated that SUITOR can correctly identify signatures, some of which were missed by other widely used methods. Applied to 2,540 whole-genome sequenced tumors across 22 cancer types, SUITOR selected signatures with the smallest prediction errors and almost all signatures of breast cancer selected by SUITOR were validated in an independent breast cancer study. SUITOR is a powerful tool to select the optimal number of mutational signatures, facilitating downstream analyses with etiological or therapeutic importance.