MagicalRsq-X: A cross-cohort transferable genotype imputation quality metric.

Since genotype imputation was introduced, researchers have been relying on the estimated imputation quality from imputation software to perform post-imputation quality control (QC). However, this quality estimate (denoted as Rsq) performs less well for lower-frequency variants. We recently published MagicalRsq, a machine-learning-based imputation quality calibration, which leverages additional typed markers from the same cohort and outperforms Rsq as a QC metric. In this work, we extended the original MagicalRsq to allow cross-cohort model training and named the new model MagicalRsq-X. We removed the cohort-specific estimated minor allele frequency and included linkage disequilibrium scores and recombination rates as additional features. Leveraging whole-genome sequencing data from TOPMed, specifically participants in the BioMe, JHS, WHI, and MESA studies, we performed comprehensive cross-cohort evaluations for predominantly European and African ancestral individuals based on their inferred global ancestry with the 1000 Genomes and Human Genome Diversity Project data as reference. Our results suggest MagicalRsq-X outperforms Rsq in almost every setting, with 7.3%-14.4% improvement in squared Pearson correlation with true R2, corresponding to 85-218 K variant gains. We further developed a metric to quantify the genetic distances of a target cohort relative to a reference cohort and showed that such metric largely explained the performance of MagicalRsq-X models. Finally, we found MagicalRsq-X saved up to 53 known genome-wide significant variants in one of the largest blood cell trait GWASs that would be missed using the original Rsq for QC. In conclusion, MagicalRsq-X shows superiority for post-imputation QC and benefits genetic studies by distinguishing well and poorly imputed lower-frequency variants.

Jasmonate activates a CsMPK6-CsMYC2 module that regulates the expression of ß-citraurin biosynthetic genes and fruit coloration in orange (Citrus sinensis).

Carotenoids are natural pigments that influence the color of citrus fruit. The red-colored carotenoid ß-citraurin is responsible for the peel color in "Newhall" orange (Citrus sinensis). Although jasmonates are known to regulate the biosynthesis and accumulation of carotenoids, their effects on ß-citraurin biosynthesis in citrus fruit remain unclear. Here, we determined that treatment with methyl jasmonate (MeJA) significantly promotes fruit coloration and ß-citraurin production in "Newhall" orange. A MeJA treatment induced the expression of CsMYC2, which encodes a transcription factor that serves as a master regulator of jasmonate responses. CsMYC2 bound the promoter of the gene that encodes carotenoid cleavage dioxygenase 4b (CsCCD4b), the key gene for ß-citraurin biosynthesis, and the promoters of genes that encode phytoene synthase (CsPSY), lycopene ß-cyclase (CsLCYb), and ß-carotene hydroxylase (CsBCH) and induced their expression. In addition, CsMYC2 promoted CsMPK6 expression. Notably, we found that CsMPK6 interacted with CsMYC2 and that this interaction decreased the stability and DNA-binding activity of CsMYC2. Thus, we conclude that negative feedback regulation attenuates JA signaling during the jasmonate-induced coloration of citrus fruit. Together, our findings indicate that jasmonates induce ß-citraurin biosynthesis in citrus by activating a CsMPK6-CsMYC2 cascade, thereby affecting fruit coloration.

S-nitrosylation of AMPKγ impairs coronary collateral circulation and disrupts VSMC reprogramming.

Collateral circulation is essential for blood resupply to the ischemic heart, which is dictated by the contractile phenotypic restoration of vascular smooth muscle cells (VSMC). Here we investigate whether S-nitrosylation of AMP-activated protein kinase (AMPK), a key regulator of the VSMC phenotype, impairs collateral circulation. In rats with collateral growth and development, nitroglycerin decreases coronary collateral blood flow (CCBF), inhibits vascular contractile phenotypic restoration, and increases myocardial infarct size, accompanied by reduced AMPK activity in the collateral zone. Nitric oxide (NO) S-nitrosylates human recombinant AMPKγ1 at cysteine 131 and decreases AMP sensitivity of AMPK. In VSMCs, exogenous expression of S-nitrosylation-resistant AMPKγ1 or deficient NO synthase (iNOS) prevents the disruption of VSMC reprogramming. Finally, hyperhomocysteinemia or hyperglycemia increases AMPKγ1 S-nitrosylation, prevents vascular contractile phenotypic restoration, reduces CCBF, and increases the infarct size of the heart in Apoe-/- mice, all of which is rescued in Apoe-/-/iNOSsm-/- mice or Apoe-/- mice with enforced expression of the AMPKγ1-C130A mutant following RI/MI. We conclude that nitrosative stress disrupts coronary collateral circulation during hyperhomocysteinemia or hyperglycemia through AMPK S-nitrosylation.

The BTB-BACK-TAZ domain protein MdBT2 reduces drought resistance by weakening the positive regulatory effect of MdHDZ27 on apple drought tolerance via ubiquitination.

Drought stress is one of the dominating challenges to the growth and productivity in crop plants. Elucidating the molecular mechanisms of plants responses to drought stress is fundamental to improve fruit quality. However, such molecular mechanisms are poorly understood in apple (Malus domestica Borkh.). In this study, we explored that the BTB-BACK-TAZ protein, MdBT2, negatively modulates the drought tolerance of apple plantlets. Moreover, we identified a novel Homeodomain-leucine zipper (HD-Zip) transcription factor, MdHDZ27, using a yeast two-hybrid (Y2H) screen with MdBT2 as the bait. Overexpression of MdHDZ27 in apple plantlets, calli, and tomato plantlets enhanced their drought tolerance by promoting the expression of drought tolerance-related genes [responsive to dehydration 29A (MdRD29A) and MdRD29B]. Biochemical analyses demonstrated that MdHDZ27 directly binds to and activates the promoters of MdRD29A and MdRD29B. Furthermore, in vitro and in vivo assays indicate that MdBT2 interacts with and ubiquitinates MdHDZ27, via the ubiquitin/26S proteasome pathway. This ubiquitination results in the degradation of MdHDZ27 and weakens the transcriptional activation of MdHDZ27 on MdRD29A and MdRD29B. Finally, a series of transgenic analyses in apple plantlets further clarified the role of the relationship between MdBT2 and MdHDZ27, as well as the effect of their interaction on drought resistance in apple plantlets. Collectively, our findings reveal a novel mechanism by which the MdBT2-MdHDZ27 regulatory module controls drought tolerance, which is of great significance for enhancing the drought resistance of apple and other plants.

Pleiotropic modifiers of age-related diabetes and neonatal intestinal obstruction in cystic fibrosis.

Individuals with cystic fibrosis (CF) develop complications of the gastrointestinal tract influenced by genetic variants outside of CFTR. Cystic fibrosis-related diabetes (CFRD) is a distinct form of diabetes with a variable age of onset that occurs frequently in individuals with CF, while meconium ileus (MI) is a severe neonatal intestinal obstruction affecting ∼20% of newborns with CF. CFRD and MI are slightly correlated traits with previous evidence of overlap in their genetic architectures. To better understand the genetic commonality between CFRD and MI, we used whole-genome-sequencing data from the CF Genome Project to perform genome-wide association. These analyses revealed variants at 11 loci (6 not previously identified) that associated with MI and at 12 loci (5 not previously identified) that associated with CFRD. Of these, variants at SLC26A9, CEBPB, and PRSS1 associated with both traits; variants at SLC26A9 and CEBPB increased risk for both traits, while variants at PRSS1, the higher-risk alleles for CFRD, conferred lower risk for MI. Furthermore, common and rare variants within the SLC26A9 locus associated with MI only or CFRD only. As expected, different loci modify risk of CFRD and MI; however, a subset exhibit pleiotropic effects indicating etiologic and mechanistic overlap between these two otherwise distinct complications of CF.

MagicalRsq: Machine-learning-based genotype imputation quality calibration.

Whole-genome sequencing (WGS) is the gold standard for fully characterizing genetic variation but is still prohibitively expensive for large samples. To reduce costs, many studies sequence only a subset of individuals or genomic regions, and genotype imputation is used to infer genotypes for the remaining individuals or regions without sequencing data. However, not all variants can be well imputed, and the current state-of-the-art imputation quality metric, denoted as standard Rsq, is poorly calibrated for lower-frequency variants. Here, we propose MagicalRsq, a machine-learning-based method that integrates variant-level imputation and population genetics statistics, to provide a better calibrated imputation quality metric. Leveraging WGS data from the Cystic Fibrosis Genome Project (CFGP), and whole-exome sequence data from UK BioBank (UKB), we performed comprehensive experiments to evaluate the performance of MagicalRsq compared to standard Rsq for partially sequenced studies. We found that MagicalRsq aligns better with true R2 than standard Rsq in almost every situation evaluated, for both European and African ancestry samples. For example, when applying models trained from 1,992 CFGP sequenced samples to an independent 3,103 samples with no sequencing but TOPMed imputation from array genotypes, MagicalRsq, compared to standard Rsq, achieved net gains of 1.4 million rare, 117k low-frequency, and 18k common variants, where net gains were gained numbers of correctly distinguished variants by MagicalRsq over standard Rsq. MagicalRsq can serve as an improved post-imputation quality metric and will benefit downstream analysis by better distinguishing well-imputed variants from those poorly imputed. MagicalRsq is freely available on GitHub.

TOP-LD: A tool to explore linkage disequilibrium with TOPMed whole-genome sequence data.

Current publicly available tools that allow rapid exploration of linkage disequilibrium (LD) between markers (e.g., HaploReg and LDlink) are based on whole-genome sequence (WGS) data from 2,504 individuals in the 1000 Genomes Project. Here, we present TOP-LD, an online tool to explore LD inferred with high-coverage (∼30×) WGS data from 15,578 individuals in the NHLBI Trans-Omics for Precision Medicine (TOPMed) program. TOP-LD provides a significant upgrade compared to current LD tools, as the TOPMed WGS data provide a more comprehensive representation of genetic variation than the 1000 Genomes data, particularly for rare variants and in the specific populations that we analyzed. For example, TOP-LD encompasses LD information for 150.3, 62.2, and 36.7 million variants for European, African, and East Asian ancestral samples, respectively, offering 2.6- to 9.1-fold increase in variant coverage compared to HaploReg 4.0 or LDlink. In addition, TOP-LD includes tens of thousands of structural variants (SVs). We demonstrate the value of TOP-LD in fine-mapping at the GGT1 locus associated with gamma glutamyltransferase in the African ancestry participants in UK Biobank. Beyond fine-mapping, TOP-LD can facilitate a wide range of applications that are based on summary statistics and estimates of LD. TOP-LD is freely available online.

Admix-kit: an integrated toolkit and pipeline for genetic analyses of admixed populations.

SUMMARY: Admixed populations, with their unique and diverse genetic backgrounds, are often underrepresented in genetic studies. This oversight not only limits our understanding but also exacerbates existing health disparities. One major barrier has been the lack of efficient tools tailored for the special challenges of genetic studies of admixed populations. Here, we present admix-kit, an integrated toolkit and pipeline for genetic analyses of admixed populations. Admix-kit implements a suite of methods to facilitate genotype and phenotype simulation, association testing, genetic architecture inference, and polygenic scoring in admixed populations. AVAILABILITY AND IMPLEMENTATION: Admix-kit package is open-source and available at https://github.com/KangchengHou/admix-kit. Additionally, users can use the pipeline designed for admixed genotype simulation available at https://github.com/UW-GAC/admix-kit_workflow.

Super interactive promoters provide insight into cell type-specific regulatory networks in blood lineage cell types.

Existing studies of chromatin conformation have primarily focused on potential enhancers interacting with gene promoters. By contrast, the interactivity of promoters per se, while equally critical to understanding transcriptional control, has been largely unexplored, particularly in a cell type-specific manner for blood lineage cell types. In this study, we leverage promoter capture Hi-C data across a compendium of blood lineage cell types to identify and characterize cell type-specific super-interactive promoters (SIPs). Notably, promoter-interacting regions (PIRs) of SIPs are more likely to overlap with cell type-specific ATAC-seq peaks and GWAS variants for relevant blood cell traits than PIRs of non-SIPs. Moreover, PIRs of cell-type-specific SIPs show enriched heritability of relevant blood cell trait (s), and are more enriched with GWAS variants associated with blood cell traits compared to PIRs of non-SIPs. Further, SIP genes tend to express at a higher level in the corresponding cell type. Importantly, SIP subnetworks incorporating cell-type-specific SIPs and ATAC-seq peaks help interpret GWAS variants. Examples include GWAS variants associated with platelet count near the megakaryocyte SIP gene EPHB3 and variants associated lymphocyte count near the native CD4 T-Cell SIP gene ETS1. Interestingly, around 25.7% ~ 39.6% blood cell traits GWAS variants residing in SIP PIR regions disrupt transcription factor binding motifs. Importantly, our analysis shows the potential of using promoter-centric analyses of chromatin spatial organization data to identify biologically important genes and their regulatory regions.

Direct Hot-Electron Transfer at the Au Nanoparticle/Monolayer Transition-Metal Dichalcogenide Interface Observed with Ultrahigh Spatiotemporal Resolution.

Plasmon-induced hot-electron transfer at the metallic nanoparticle/semiconductor interface is the basis of plasmon-enhanced photocatalysis and energy harvesting. However, limited by the nanoscale size of hot spots and femtosecond time scale of hot-electron transfer, direct observation is still challenging. Herein, by using spatiotemporal-resolved photoemission electron microscopy with a two-color pump-probe beamline, we directly observed such a process with a concise system, the Au nanoparticle/monolayer transition-metal dichalcogenide (TMD) interface. The ultrafast hot-electron transfer from Au nanoparticles to monolayer TMDs and the plasmon-enhanced transfer process were directly measured and verified through an in situ comparison with the Au film/TMD interface and free TMDs. The lifetime at the Au nanoparticle/MoSe2 interface decreased from 410 to 42 fs, while the photoemission intensities exhibited a 27-fold increase compared to free MoSe2. We also measured the evolution of hot electrons in the energy distributions, indicating the hot-electron injection and decay happened in an ultrafast time scale of ∼50 fs without observable electron cooling.

Transcriptome-wide association study in UK Biobank Europeans identifies associations with blood cell traits.

Previous genome-wide association studies (GWAS) of hematological traits have identified over 10 000 distinct trait-specific risk loci. However, at these loci, the underlying causal mechanisms remain incompletely characterized. To elucidate novel biology and better understand causal mechanisms at known loci, we performed a transcriptome-wide association study (TWAS) of 29 hematological traits in 399 835 UK Biobank (UKB) participants of European ancestry using gene expression prediction models trained from whole blood RNA-seq data in 922 individuals. We discovered 557 gene-trait associations for hematological traits distinct from previously reported GWAS variants in European populations. Among the 557 associations, 301 were available for replication in a cohort of 141 286 participants of European ancestry from the Million Veteran Program. Of these 301 associations, 108 replicated at a strict Bonferroni adjusted threshold ($\alpha$= 0.05/301). Using our TWAS results, we systematically assigned 4261 out of 16 900 previously identified hematological trait GWAS variants to putative target genes. Compared to coloc, our TWAS results show reduced specificity and increased sensitivity in external datasets to assign variants to target genes.

Whole genome sequence analysis of platelet traits in the NHLBI Trans-Omics for Precision Medicine (TOPMed) initiative.

Platelets play a key role in thrombosis and hemostasis. Platelet count (PLT) and mean platelet volume (MPV) are highly heritable quantitative traits, with hundreds of genetic signals previously identified, mostly in European ancestry populations. We here utilize whole genome sequencing (WGS) from NHLBI's Trans-Omics for Precision Medicine initiative (TOPMed) in a large multi-ethnic sample to further explore common and rare variation contributing to PLT (n = 61 200) and MPV (n = 23 485). We identified and replicated secondary signals at MPL (rs532784633) and PECAM1 (rs73345162), both more common in African ancestry populations. We also observed rare variation in Mendelian platelet-related disorder genes influencing variation in platelet traits in TOPMed cohorts (not enriched for blood disorders). For example, association of GP9 with lower PLT and higher MPV was partly driven by a pathogenic Bernard-Soulier syndrome variant (rs5030764, p.Asn61Ser), and the signals at TUBB1 and CD36 were partly driven by loss of function variants not annotated as pathogenic in ClinVar (rs199948010 and rs571975065). However, residual signal remained for these gene-based signals after adjusting for lead variants, suggesting that additional variants in Mendelian genes with impacts in general population cohorts remain to be identified. Gene-based signals were also identified at several genome-wide association study identified loci for genes not annotated for Mendelian platelet disorders (PTPRH, TET2, CHEK2), with somatic variation driving the result at TET2. These results highlight the value of WGS in populations of diverse genetic ancestry to identify novel regulatory and coding signals, even for well-studied traits like platelet traits.

Whole-genome sequencing association analysis of quantitative red blood cell phenotypes: The NHLBI TOPMed program.

Whole-genome sequencing (WGS), a powerful tool for detecting novel coding and non-coding disease-causing variants, has largely been applied to clinical diagnosis of inherited disorders. Here we leveraged WGS data in up to 62,653 ethnically diverse participants from the NHLBI Trans-Omics for Precision Medicine (TOPMed) program and assessed statistical association of variants with seven red blood cell (RBC) quantitative traits. We discovered 14 single variant-RBC trait associations at 12 genomic loci, which have not been reported previously. Several of the RBC trait-variant associations (RPN1, ELL2, MIDN, HBB, HBA1, PIEZO1, and G6PD) were replicated in independent GWAS datasets imputed to the TOPMed reference panel. Most of these discovered variants are rare/low frequency, and several are observed disproportionately among non-European Ancestry (African, Hispanic/Latino, or East Asian) populations. We identified a 3 bp indel p.Lys2169del (g.88717175_88717177TCT[4]) (common only in the Ashkenazi Jewish population) of PIEZO1, a gene responsible for the Mendelian red cell disorder hereditary xerocytosis (MIM: 194380), associated with higher mean corpuscular hemoglobin concentration (MCHC). In stepwise conditional analysis and in gene-based rare variant aggregated association analysis, we identified several of the variants in HBB, HBA1, TMPRSS6, and G6PD that represent the carrier state for known coding, promoter, or splice site loss-of-function variants that cause inherited RBC disorders. Finally, we applied base and nuclease editing to demonstrate that the sentinel variant rs112097551 (nearest gene RPN1) acts through a cis-regulatory element that exerts long-range control of the gene RUVBL1 which is essential for hematopoiesis. Together, these results demonstrate the utility of WGS in ethnically diverse population-based samples and gene editing for expanding knowledge of the genetic architecture of quantitative hematologic traits and suggest a continuum between complex trait and Mendelian red cell disorders.

Whole-genome sequencing in diverse subjects identifies genetic correlates of leukocyte traits: The NHLBI TOPMed program.

Many common and rare variants associated with hematologic traits have been discovered through imputation on large-scale reference panels. However, the majority of genome-wide association studies (GWASs) have been conducted in Europeans, and determining causal variants has proved challenging. We performed a GWAS of total leukocyte, neutrophil, lymphocyte, monocyte, eosinophil, and basophil counts generated from 109,563,748 variants in the autosomes and the X chromosome in the Trans-Omics for Precision Medicine (TOPMed) program, which included data from 61,802 individuals of diverse ancestry. We discovered and replicated 7 leukocyte trait associations, including (1) the association between a chromosome X, pseudo-autosomal region (PAR), noncoding variant located between cytokine receptor genes (CSF2RA and CLRF2) and lower eosinophil count; and (2) associations between single variants found predominantly among African Americans at the S1PR3 (9q22.1) and HBB (11p15.4) loci and monocyte and lymphocyte counts, respectively. We further provide evidence indicating that the newly discovered eosinophil-lowering chromosome X PAR variant might be associated with reduced susceptibility to common allergic diseases such as atopic dermatitis and asthma. Additionally, we found a burden of very rare FLT3 (13q12.2) variants associated with monocyte counts. Together, these results emphasize the utility of whole-genome sequencing in diverse samples in identifying associations missed by European-ancestry-driven GWASs.

Comparative transcriptome analysis reveals genes involved in trichome development and metabolism in tobacco.

BACKGROUND: The glandular trichomes of tobacco (Nicotiana tabacum) can efficiently produce secondary metabolites. They act as natural bioreactors, and their natural products function to protect plants against insect-pests and pathogens and are also components of industrial chemicals. To clarify the molecular mechanisms of tobacco glandular trichome development and secondary metabolic regulation, glandular trichomes and glandless trichomes, as well as other different developmental tissues, were used for RNA sequencing and analysis. RESULTS: By comparing glandless and glandular trichomes with other tissues, we obtained differentially expressed genes. They were obviously enriched in KEGG pathways, such as cutin, suberine, and wax biosynthesis, flavonoid and isoflavonoid biosynthesis, terpenoid biosynthesis, and plant-pathogen interaction. In particular, the expression levels of genes related to the terpenoid, flavonoid, and wax biosynthesis pathway mainly showed down-regulation in glandless trichomes, implying that they lack the capability to synthesize certain exudate compounds. Among the differentially expressed genes, 234 transcription factors were found, including AP2-ERFs, MYBs, bHLHs, WRKYs, Homeoboxes (HD-ZIP), and C2H2-ZFs. These transcription factor and genes that highly expressed in trichomes or specially expressed in GT or GLT. Following the overexpression of R2R3-MYB transcription factor Nitab4.5_0011760g0030.1 in tobacco, an increase in the number of branched glandular trichomes was observed. CONCLUSIONS: Our data provide comprehensive gene expression information at the transcriptional level and an understanding of the regulatory pathways involved in glandular trichome development and secondary metabolism.

SnapHiC: a computational pipeline to identify chromatin loops from single-cell Hi-C data.

Single-cell Hi-C (scHi-C) analysis has been increasingly used to map chromatin architecture in diverse tissue contexts, but computational tools to define chromatin loops at high resolution from scHi-C data are still lacking. Here, we describe Single-Nucleus Analysis Pipeline for Hi-C (SnapHiC), a method that can identify chromatin loops at high resolution and accuracy from scHi-C data. Using scHi-C data from 742 mouse embryonic stem cells, we benchmark SnapHiC against a number of computational tools developed for mapping chromatin loops and interactions from bulk Hi-C. We further demonstrate its use by analyzing single-nucleus methyl-3C-seq data from 2,869 human prefrontal cortical cells, which uncovers cell type-specific chromatin loops and predicts putative target genes for noncoding sequence variants associated with neuropsychiatric disorders. Our results indicate that SnapHiC could facilitate the analysis of cell type-specific chromatin architecture and gene regulatory programs in complex tissues.

MdPRX34L, a class III peroxidase gene, activates the immune response in apple to the fungal pathogen Botryosphaeria dothidea.

MAIN CONCLUSION: MdPRX34L enhanced resistance to Botryosphaeria dothidea by increasing salicylic acid (SA) and abscisic acid (ABA) content as well as the expression of related defense genes. The class III peroxidase (PRX) multigene family is involved in complex biological processes. However, the molecular mechanism of PRXs in the pathogen defense of plants against Botryosphaeria dothidea (B. dothidea) remains unclear. Here, we cloned the PRX gene MdPRX34L, which was identified as a positive regulator of the defense response to B. dothidea, from the apple cultivar 'Royal Gala.' Overexpression of MdPRX34L in apple calli decreased sensitivity to salicylic acid (SA) and abscisic acid（ABA). Subsequently, overexpression of MdPRX34L in apple calli increased resistance to B. dothidea infection. In addition, SA contents and the expression levels of genes related to SA synthesis and signaling in apple calli overexpressing MdPRX34L were higher than those in the control after inoculation, suggesting that MdPRX34L enhances resistance to B. dothidea via the SA pathway. Interestingly, infections in apple calli by B. dothidea caused an increase in endogenous levels of ABA followed by induction of ABA-related genes expression. These findings suggest a potential mechanism by which MdPRX34L enhances plant-pathogen defense against B. dothidea by regulating the SA and ABA pathways.

The transcriptional regulatory module CsHB5-CsbZIP44 positively regulates abscisic acid-mediated carotenoid biosynthesis in citrus (Citrus spp.).

Carotenoids contribute to fruit coloration and are valuable sources of provitamin A in the human diet. Abscisic acid (ABA) plays an essential role in fruit coloration during citrus fruit ripening, but little is known about the underlying mechanisms. Here, we identified a novel bZIP transcription activator called CsbZIP44, which serves as a central regulator of ABA-mediated citrus carotenoid biosynthesis. CsbZIP44 directly binds to the promoters of four carotenoid metabolism-related genes (CsDXR, CsGGPPs, CsBCH1 and CsNCED2) and activates their expression. Furthermore, our research indicates that CsHB5, a positive regulator of ABA and carotenoid-driven processes, activates the expression of CsbZIP44 by binding to its promoter. Additionally, CsHB5 interacts with CsbZIP44 to form a transcriptional regulatory module CsHB5-CsbZIP44, which is responsive to ABA induction and promotes carotenoid accumulation in citrus. Interestingly, we also discover a positive feedback regulation loop between the ABA signal and carotenoid biosynthesis mediated by the CsHB5-CsbZIP44 transcriptional regulatory module. Our findings show that CsHB5-CsbZIP44 precisely modulates ABA signal-mediated carotenoid metabolism, providing an effective strategy for quality improvement of citrus fruit and other crops.

eSCAN: scan regulatory regions for aggregate association testing using whole-genome sequencing data.

Multiple statistical methods for aggregate association testing have been developed for whole-genome sequencing (WGS) data. Many aggregate variants in a given genomic window and ignore existing knowledge to define test regions, resulting in many identified regions not clearly linked to genes, and thus, limiting biological understanding. Functional information from new technologies (such as Hi-C and its derivatives), which can help link enhancers to their effector genes, can be leveraged to predefine variant sets for aggregate testing in WGS data. Here, we propose the eSCAN (scan the enhancers) method for genome-wide assessment of enhancer regions in sequencing studies, combining the advantages of dynamic window selection in SCANG (SCAN the Genome), a previously developed method, with the advantages of incorporating putative regulatory regions from annotation. eSCAN, by searching in putative enhancers, increases statistical power and aids mechanistic interpretation, as demonstrated by extensive simulation studies. We also apply eSCAN for blood cell traits using NHLBI Trans-Omics for Precision Medicine WGS data. Results from real data analysis show that eSCAN is able to capture more significant signals, and these signals are of shorter length (indicating higher resolution fine-mapping capability) and drive association of larger regions detected by other methods.

Phytochrome B mediates dim-light-reduced insect resistance by promoting the ethylene pathway in rice.

Increasing planting density is one of the most effective ways to improve crop yield. However, one major factor that limits crop planting density is the weakened immunity of plants to pathogens and insects caused by dim light (DL) under shade conditions. The molecular mechanism underlying how DL compromises plant immunity remains unclear. Here, we report that DL reduces rice (Oryza sativa) resistance against brown planthopper (BPH; Nilaparvata lugens) by elevating ethylene (ET) biosynthesis and signaling in a Phytochrome B (OsPHYB)-dependent manner. The DL-reduced BPH resistance is relieved in osphyB mutants, but aggravated in OsPHYB overexpressing plants. Further, we found that DL reduces the nuclear accumulation of OsphyB, thus alleviating Phytochrome Interacting Factor Like14 (OsPIL14) degradation, consequently leading to the up-regulation of 1-Aminocyclopropane-1-Carboxylate Oxidase1 (OsACO1) and an increase in ET levels. In addition, we found that nuclear OsphyB stabilizes Ethylene Insensitive Like2 (OsEIL2) by competitively interacting with EIN3 Binding F-Box Protein (OsEBF1) to enhance ET signaling in rice, which contrasts with previous findings that phyB blocks ET signaling by facilitating Ethylene Insensitive3 (EIN3) degradation in other plant species. Thus, enhanced ET biosynthesis and signaling reduces BPH resistance under DL conditions. Our findings provide insights into the molecular mechanism of the light-regulated ET pathway and host-insect interactions and potential strategies for sustainable insect management.