VSSI-GGD: A Variation Sparse EEG Source Imaging Approach Based on Generalized Gaussian Distribution.

Electroencephalographic (EEG) source imaging (ESI) is a powerful method for studying brain functions and surgical resection of epileptic foci. However, accurately estimating the location and extent of brain sources remains challenging due to noise and background interference in EEG signals. To reconstruct extended brain sources, we propose a new ESI method called Variation Sparse Source Imaging based on Generalized Gaussian Distribution (VSSI-GGD). VSSI-GGD uses the generalized Gaussian prior as a sparse constraint on the spatial variation domain and embeds it into the Bayesian framework for source estimation. Using a variational technique, we approximate the intractable true posterior with a Gaussian density. Through convex analysis, the Bayesian inference problem is transformed entirely into a series of regularized L2p -norm ( ) optimization problems, which are efficiently solved with the ADMM algorithm. Imaging results of numerical simulations and human experimental dataset analysis reveal the superior performance of VSSI-GGD, which provides higher spatial resolution with clear boundaries compared to benchmark algorithms. VSSI-GGD can potentially serve as an effective and robust spatiotemporal EEG source imaging method. The source code of VSSI-GGD is available at https://github.com/Mashirops/VSSI-GGD.git.

The Idesia polycarpa genome provides insights into its evolution and oil biosynthesis.

The deciduous tree Idesia polycarpa can provide premium edible oil with high polyunsaturated fatty acid contents. Here, we generate its high-quality reference genome, which is ∼1.21 Gb, comprising 21 pseudochromosomes and 42,086 protein-coding genes. Phylogenetic and genomic synteny analyses show that it diverged with Populus trichocarpa about 16.28 million years ago. Notably, most fatty acid biosynthesis genes are not only increased in number in its genome but are also highly expressed in the fruits. Moreover, we identify, through genome-wide association analysis and RNA sequencing, the I. polycarpa SUGAR TRANSPORTER 5 (IpSTP5) gene as a positive regulator of high oil accumulation in the fruits. Silencing of IpSTP5 by virus-induced gene silencing causes a significant reduction of oil content in the fruits, suggesting it has the potential to be used as a molecular marker to breed the high-oil-content cultivars. Our results collectively lay the foundation for breeding the elite cultivars of I. polycarpa.

TNF and IFNγ-induced cell death requires IRF1 and ELAVL1 to promote CASP8 expression.

TNFα and IFNγ (TNF/IFNγ) synergistically induce caspase-8 activation and cancer cell death. However, the mechanism of IFNγ in promoting TNF-initiated caspase-8 activation in cancer cells is poorly understood. Here, we found that in addition to CASP8, CYLD is transcriptionally upregulated by IFNγ-induced transcription factor IRF1. IRF1-mediated CASP8 and CYLD upregulation additively mediates TNF/IFNγ-induced cancer cell death. Clinically, the expression levels of TNF, IFNγ, CYLD, and CASP8 in melanoma tumors are increased in patients responsive to immune checkpoint blockade (ICB) therapy after anti-PD-1 treatment. Accordingly, our genetic screen revealed that ELAVL1 (HuR) is required for TNF/IFNγ-induced caspase-8 activation. Mechanistically, ELAVL1 binds CASP8 mRNA and extends its stability to sustain caspase-8 expression both in IFNγ-stimulated and in basal conditions. Consequently, ELAVL1 determines death receptors-initiated caspase-8-dependent cell death triggered from stimuli including TNF and TRAIL by regulating basal/stimulated caspase-8 levels. As caspase-8 is a master regulator in cell death and inflammation, these results provide valuable clues for tumor immunotherapy and inflammatory diseases.

The chromosome-level assembly of the wild diploid alfalfa genome provides insights into the full landscape of genomic variations between cultivated and wild alfalfa.

Alfalfa (Medicago sativa L.) is one of the most important forage legumes in the world, including autotetraploid (M. sativa ssp. sativa) and diploid alfalfa (M. sativa ssp. caerulea, progenitor of autotetraploid alfalfa). Here, we reported a high-quality genome of ZW0012 (diploid alfalfa, 769 Mb, contig N50 = 5.5 Mb), which was grouped into the Northern group in population structure analysis, suggesting that our genome assembly filled a major gap among the members of M. sativa complex. During polyploidization, large phenotypic differences occurred between diploids and tetraploids, and the genetic information underlying its massive phenotypic variations remains largely unexplored. Extensive structural variations (SVs) were identified between ZW0012 and XinJiangDaYe (an autotetraploid alfalfa with released genome). We identified 71 ZW0012-specific PAV genes and 1296 XinJiangDaYe-specific PAV genes, mainly involved in defence response, cell growth, and photosynthesis. We have verified the positive roles of MsNCR1 (a XinJiangDaYe-specific PAV gene) in nodulation using an Agrobacterium rhizobia-mediated transgenic method. We also demonstrated that MsSKIP23_1 and MsFBL23_1 (two XinJiangDaYe-specific PAV genes) regulated leaf size by transient overexpression and virus-induced gene silencing analysis. Our study provides a high-quality reference genome of an important diploid alfalfa germplasm and a valuable resource of variation landscape between diploid and autotetraploid, which will facilitate the functional gene discovery and molecular-based breeding for the cultivars in the future.

Genetic insights into the dissolution of dioecy in diploid persimmon Diospyros oleifera Cheng.

BACKGROUND: Dioecy, a sexual system of single-sexual (gynoecious/androecious) individuals, is rare in flowering plants. This rarity may be a result of the frequent transition from dioecy into systems with co-sexual individuals. RESULTS: In this study, co-sexual expression (monoecy and hermaphroditic development), previously thought to be polyploid-specific in Diospyros species, was identified in the diploid D. oleifeara historically. We characterized potential genetic mechanisms that underlie the dissolution of dioecy to monoecy and andro(gyno)monoecy, based on multiscale genome-wide investigations of 150 accessions of Diospyros oleifera. We found all co-sexual plants, including monoecious and andro(gyno)monoecious individuals, possessed the male determinant gene OGI, implying the presence of genetic factors controlling gynoecia development in genetically male D. oleifera. Importantly, discrepancies in the OGI/MeGI module were found in diploid monoecious D. oleifera compared with polyploid monoecious D. kaki, including no Kali insertion on the promoter of OGI, no different abundance of smRNAs targeting MeGI (a counterpart of OGI), and no different expression of MeGI between female and male floral buds. On the contrary, in both single- and co-sexual plants, female function was expressed in the presence of a genome-wide decrease in methylation levels, along with sexually distinct regulatory networks of smRNAs and their targets. Furthermore, a genome-wide association study (GWAS) identified a genomic region and a DUF247 gene cluster strongly associated with the monoecious phenotype and several regions that may contribute to andromonoecy. CONCLUSIONS: Collectively, our findings demonstrate stable breakdown of the dioecious system in D. oleifera, presumably also a result of genomic features of the Y-linked region.

Open reduction and internal fixation of irreducible displaced femoral neck fracture with femoral Neck System: a preliminary study.

BACKGROUND: Most displaced femoral neck fractures can achieve satisfactory anatomical reduction by closed reduction, but there are still some that cannot reset satisfactorily after closed reduction, and open reduction are required. Such fractures that cannot be repositioned successfully by closed reduction are called irreducible displaced femoral neck fractures in this study. The objective of our study was to evaluate the efficacy of direct anterior incision with the Femoral Neck System in the treatment of irreducible displaced femoral fractures. METHODS: A total of 16 young and middle-aged patients with irreducible displaced femoral neck fractures involving Garden type III and IV were treated using Femoral Neck System fixation by open reduction through Direct Anterior Approach between January 2020 to September 2021. Functional outcomes and postoperative complications were assessed during follow-up. Clinical outcomes were evaluated by the Hip Harris score. The postoperative reduction was evaluated by the Garden Index. Observe postoperative complications. RESULTS: All patients were followed up with a mean follow-up time of 21.1(12-30) months, and according to radiological results, all patients achieved fracture healing, with a mean healing time of 4.25 months. All 16 patients received grade Garden I and II reductions, and there was no significant difference in the anteroposterior Garden reduction index between the first day after surgery (166.13 ± 5.61) and the 12th month after surgery(164.94 ± 4.49) (P>0.05) and no significant difference in lateral Garden index between the first day after surgery(171.06 ± 4.46) and the 12th month after surgery(169.38 ± 3.98) (P<0.05). According to the Hip Harris score scale, 13 patients received excellent and 3 patients received good. The postoperative Hip Harris Score(17.19 ± 4.8) was significantly higher than the preoperative score(92.19 ± 3.4), and the difference was statistically significant (P < 0.05). No or mild femoral neck shortness occurred in 12 (75%) patients, moderate shortening occurred in 3 (18.75%) patients, and severe shortening occurred in 1 (6.25%) patient. None of the patients experienced femoral head necrosis, fracture nonunion, or incision infection. One patient developed deep venous thrombosis of the lower extremity. CONCLUSIONS: The Direct Anterior Approach combined with Femoral Neck System is an excellent treatment for irreducible displaced femoral neck fracture and achieved good functional outcomes and anatomical reduction with low complications.

Genomic insight into domestication of rubber tree.

Understanding the genetic basis of rubber tree (Hevea brasiliensis) domestication is crucial for further improving natural rubber production to meet its increasing demand worldwide. Here we provide a high-quality H. brasiliensis genome assembly (1.58 Gb, contig N50 of 11.21 megabases), present a map of genome variations by resequencing 335 accessions and reveal domestication-related molecular signals and a major domestication trait, the higher number of laticifer rings. We further show that HbPSK5, encoding the small-peptide hormone phytosulfokine (PSK), is a key domestication gene and closely correlated with the major domestication trait. The transcriptional activation of HbPSK5 by myelocytomatosis (MYC) members links PSK signaling to jasmonates in regulating the laticifer differentiation in rubber tree. Heterologous overexpression of HbPSK5 in Russian dandelion (Taraxacum kok-saghyz) can increase rubber content by promoting laticifer formation. Our results provide an insight into target genes for improving rubber tree and accelerating the domestication of other rubber-producing plants.

Adaptive evolution of the enigmatic Takakia now facing climate change in Tibet.

The most extreme environments are the most vulnerable to transformation under a rapidly changing climate. These ecosystems harbor some of the most specialized species, which will likely suffer the highest extinction rates. We document the steepest temperature increase (2010-2021) on record at altitudes of above 4,000 m, triggering a decline of the relictual and highly adapted moss Takakia lepidozioides. Its de-novo-sequenced genome with 27,467 protein-coding genes includes distinct adaptations to abiotic stresses and comprises the largest number of fast-evolving genes under positive selection. The uplift of the study site in the last 65 million years has resulted in life-threatening UV-B radiation and drastically reduced temperatures, and we detected several of the molecular adaptations of Takakia to these environmental changes. Surprisingly, specific morphological features likely occurred earlier than 165 mya in much warmer environments. Following nearly 400 million years of evolution and resilience, this species is now facing extinction.

Advances in the application of Mxene nanoparticles in wound healing.

Skin is the largest organ of the human body. It plays a vital role as the body's first barrier: stopping chemical, radiological damage and microbial invasion. The importance of skin to the human body can never be overstated. Delayed wound healing after a skin injury has become a huge challenge in healthcare. In some situations, this can have very serious and even life-threatening effects on people's health. Various wound dressings have been developed to promote quicker wound healing, including hydrogels, gelatin sponges, films, and bandages, all work to prevent the invasion of microbial pathogens. Some of them are also packed with bioactive agents, such as antibiotics, nanoparticles, and growth factors, that help to improve the performance of the dressing it is added to. Recently, bioactive nanoparticles as the bioactive agent have become widely used in wound dressings. Among these, functional inorganic nanoparticles are favored due to their ability to effectively improve the tissue-repairing properties of biomaterials. MXene nanoparticles have attracted the interest of scholars due to their unique properties of electrical conductivity, hydrophilicity, antibacterial properties, and biocompatibility. The potential for its application is very promising as an effective functional component of wound dressings. In this paper, we will review MXene nanoparticles in skin injury repair, particularly its synthesis method, functional properties, biocompatibility, and application.

WeiTsing, a pericycle-expressed ion channel, safeguards the stele to confer clubroot resistance.

Plant roots encounter numerous pathogenic microbes that often cause devastating diseases. One such pathogen, Plasmodiophora brassicae (Pb), causes clubroot disease and severe yield losses on cruciferous crops worldwide. Here, we report the isolation and characterization of WeiTsing (WTS), a broad-spectrum clubroot resistance gene from Arabidopsis. WTS is transcriptionally activated in the pericycle upon Pb infection to prevent pathogen colonization in the stele. Brassica napus carrying the WTS transgene displayed strong resistance to Pb. WTS encodes a small protein localized in the endoplasmic reticulum (ER), and its expression in plants induces immune responses. The cryoelectron microscopy (cryo-EM) structure of WTS revealed a previously unknown pentameric architecture with a central pore. Electrophysiology analyses demonstrated that WTS is a calcium-permeable cation-selective channel. Structure-guided mutagenesis indicated that channel activity is strictly required for triggering defenses. The findings uncover an ion channel analogous to resistosomes that triggers immune signaling in the pericycle.

Genomic decoding of breeding history to guide breeding-by-design in rice.

Deciphering the intrinsic molecular logic of empirical crop breeding from a genomic perspective is a decisive prerequisite for breeding-by-design (BbD), but remains not well established. Here, we decoded the historical features of past rice breeding by phenotyping and haplotyping 546 accessions covering the majority of cultivars bred in the history of Northeast China (NEC). We revealed that three groups founded the genetic diversities in NEC rice with distinct evolution patterns and traced and verified the breeding footprints to known or genome-wide association study (GWAS)-detected quantitative trait loci (QTLs), or introgressions from indica sub-species with chronological changes in allele frequencies. Then we summarized a rice breeding trend/principle in NEC, and combined with the successful example in breeding and application of Zhongkefa5 to demonstrate the guiding value of our conclusion for BbD in practice. Our study provides a paradigm for decoding the breeding history of a specific crop to guide BbD, which may have implications in different crop breeding.

High-quality Fagopyrum esculentum genome provides insights into the flavonoid accumulation among different tissues and self-incompatibility.

Common buckwheat (Fagopyrum esculentum) and Tartary buckwheat (Fagopyrum tataricum), the two most widely cultivated buckwheat species, differ greatly in flavonoid content and reproductive mode. Here, we report the first high-quality and chromosome-level genome assembly of common buckwheat with 1.2 Gb. Comparative genomic analysis revealed that common buckwheat underwent a burst of long terminal repeat retrotransposons insertion accompanied by numerous large chromosome rearrangements after divergence from Tartary buckwheat. Moreover, multiple gene families involved in stress tolerance and flavonoid biosynthesis such as multidrug and toxic compound extrusion (MATE) and chalcone synthase (CHS) underwent significant expansion in buckwheat, especially in common buckwheat. Integrated multi-omics analysis identified high expression of catechin biosynthesis-related genes in flower and seed in common buckwheat and high expression of rutin biosynthesis-related genes in seed in Tartary buckwheat as being important for the differences in flavonoid type and content between these buckwheat species. We also identified a candidate key rutin-degrading enzyme gene (Ft8.2377) that was highly expressed in Tartary buckwheat seed. In addition, we identified a haplotype-resolved candidate locus containing many genes reportedly associated with the development of flower and pollen, which was potentially related to self-incompatibility in common buckwheat. Our study provides important resources facilitating future functional genomics-related research of flavonoid biosynthesis and self-incompatibility in buckwheat.

Biased mutations and gene losses underlying diploidization of the tetraploid broomcorn millet genome.

Broomcorn millet (Panicum miliaceum L.) is one of the earliest domesticated crops, and is a valuable resource to secure food diversity and combat drought stresses under the global warming scenario. However, due to the absence of extant diploid progenitors, the polyploidy genome of broomcorn millet remains poorly understood. Here, we report the chromosome-scale genome assembly of broomcorn millet. We divided the broomcorn millet genome into two subgenomes using the genome sequence of Panicum hallii, a diploid relative of broomcorn millet. Our analyses revealed that the two subgenomes diverged at ~4.8 million years ago (Mya), while the allotetraploidization of broomcorn millet may have occurred about ~0.48 Mya, suggesting that broomcorn millet is a relatively recent allotetraploid. Comparative analyses showed that subgenome B was larger than subgenome A in size, which was caused by the biased accumulation of long terminal repeat retrotransposons in the progenitor of subgenome B before polyploidization. Notably, the accumulation of biased mutations in the transposable element-rich subgenome B led to more gene losses. Although no significant dominance of either subgenome was observed in the expression profiles of broomcorn millet, we found the minimally expressed genes in P. hallii tended to be lost during diploidization of broomcorn millet. These results suggest that broomcorn millet is at the early stage of diploidization and that mutations likely occurred more on genes that were marked with lower expression levels.

Identification of immune-associated prognostic biomarkers in lung adenocarcinoma on the basis of gene co-expression network.

Objective: We aimed to explore immune-related prognosis genes of lung adenocarcinoma (LUAD).Materials and methods: TCGA-LUAD and GSE31210 data sets were accessed from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) respectively. By using "WGCNA" R package, we established a gene co-expression network and clustered genes into various modules. The correlation between immune scores and module eigengenes by using Pearson analysis. Screened hub genes and constructed prognostic model by using LASSO and Cox regression analysis. Evaluated model by survival analysis and receiver operating characteristic (ROC) curves. Hub genes expression in clinical tissues of LUAD patients by qRT-PCR analysis. ssGSEA and TIMER (a website tool for examination of different immune cells in different cancers) analyzed immune correlation of hub genes. Gene set variation analysis (GSVA) uncovered difference of signal pathway between high- and low-risk score group.Results: We found that brown module significantly correlated with the immune scores of immune cells. Therefore, we constructed a 7-gene prognostic model based on brown module genes, and indicated that this model possessed good predictive performance. Patients in training and validation sets were stratified into the high- and low-risk group using this model. Also, hub genes CDCP1, PLSCR1 and CD79A were highly expressed in clinical tissues of LUAD patients, while ID1, CLEC7A, KIAA1324 and CMTM7 were lowly expressed. Both ssGSEA and TIMER revealed a significant negative correlation between risk score and B cell infiltration. Additionally, some signal pathways were suppressed in the high-risk group.Conclusion: We identified 7 immune-associated prognostic markers, which may play vital roles in LUAD and could be used as hopeful targets for immunotherapy of LUAD.

Augmented wound healing potential of photosensitive GelMA hydrogel incorporating antimicrobial peptides and MXene nanoparticles.

Introduction: Wound healing is a delicate and complex process influenced by many factors. The treatment of skin wounds commonly involves the use of wound dressings, which remain a routine approach. An ideal dressing can provide protection and a suitable environment for wound surfaces by maintaining moisture and exhibiting good biocompatibility, mechanical strength, and antibacterial properties to promote healing and prevent infection. Methods: We encapsulated tick-derived antibacterial polypeptides (Os) as a model drug within a methylacrylyl gelatin (GelMA) hydrogel containing MXene nanoparticles. The prepared composite hydrogels were evaluated for their wound dressing potential by analyzing surface morphology, mechanical properties, swelling behavior, degradation properties, antibacterial activity, and cytocompatibility. Results: The results demonstrated excellent mechanical strength, swelling performance, degradation behavior, and antibacterial activity of the prepared composite hydrogels, effectively promoting cell growth, adhesion, and expression of antibacterial peptide activity. A full-thickness rat wound model then observed the wound healing process and surface interactions between the composite hydrogels and wounds. The composite hydrogel significantly accelerated wound closure, reduced inflammation, and sped epithelial formation and maturation. Discussion: Incorporating antibacterial peptides into GelMA provides a feasible strategy for developing excellent antibacterial wound dressings capable of tissue repair. In conclusion, this study presents a GelMA-based approach for designing antibacterial dressings with strong tissue regenerative ability.

Comparative genomic and transcriptomic analyses uncover the molecular basis of high nitrogen-use efficiency in the wheat cultivar Kenong 9204.

Studying the regulatory mechanisms that drive nitrogen-use efficiency (NUE) in crops is important for sustainable agriculture and environmental protection. In this study, we generated a high-quality genome assembly for the high-NUE wheat cultivar Kenong 9204 and systematically analyzed genes related to nitrogen uptake and metabolism. By comparative analyses, we found that the high-affinity nitrate transporter gene family had expanded in Triticeae. Further studies showed that subsequent functional differentiation endowed the expanded family members with saline inducibility, providing a genetic basis for improving the adaptability of wheat to nitrogen deficiency in various habitats. To explore the genetic and molecular mechanisms of high NUE, we compared genomic and transcriptomic data from the high-NUE cultivar Kenong 9204 (KN9204) and the low-NUE cultivar Jing 411 and quantified their nitrogen accumulation under high- and low-nitrogen conditions. Compared with Jing 411, KN9204 absorbed significantly more nitrogen at the reproductive stage after shooting and accumulated it in the shoots and seeds. Transcriptome data analysis revealed that nitrogen deficiency clearly suppressed the expression of genes related to cell division in the young spike of Jing 411, whereas this suppression of gene expression was much lower in KN9204. In addition, KN9204 maintained relatively high expression of NPF genes for a longer time than Jing 411 during seed maturity. Physiological and transcriptome data revealed that KN9204 was more tolerant of nitrogen deficiency than Jing 411, especially at the reproductive stage. The high NUE of KN9204 is an integrated effect controlled at different levels. Taken together, our data provide new insights into the molecular mechanisms of NUE and important gene resources for improving wheat cultivars with a higher NUE trait.

The Pararectus approach in acetabular fractures treatment: functional and radiologcial results.

BACKGROUND: The surgical treatment of complex acetabular fractures is one of the most challenging procedures for orthopedic surgeons. The Pararectus approach, as a reasonable alternative to the existing surgical procedures, was performed for the treatment of acetabular fractures involving the anterior column. This study aimed to evaluate outcome using the Pararectus approach for acetabular fractures involving anterior columns. METHODS: Thirty-seven with displaced acetabular fractures involving anterior columns were treated between July 2016 and October 2019 using the Pararectus approach. The functional outcomes (using the Merle d Aubigné and Postel scoring system, WOMAC and modified Harris scoring), the quality of surgical reduction (using the Matta criteria), and postoperative complications were assessed during approximately 26 months follow-up period. RESULTS: Thirty-seven patients (mean age 53 years, range: 30-71; 28 male) underwent surgery. Mean intraoperative blood loss was 840 ml (rang: 400-2000 ml) and mean operating time was 210 min (rang: 140-500 min). The modified Merle d Aubigné score was excellent and good in 27 cases (73%), fair in 6 cases (16%), and poor in 3 cases (11%). The mean score was 88.5 (range:77-96) for the modified Harris Hip scores, and 22 (range:7-35) for the WOMAC scores after operation. Postoperative functional outcomes were significantly improved compared with preoperative outcomes (P < 0.0001). The quality of reduction was anatomical in 21 cases (57%), satisfactory in 9 cases (24%), and unsatisfactory in 7 cases (19%). At follow-up, four patients developed a DVT, and heterotopic bone formation was observed in one patient. The hip osteoarthritis was not observed. CONCLUSION: The Pararectus approach achieved good functional outcomes and anatomical reduction in the treatment of acetabular fractures involving anterior column with minimal access morbidity.

Omics-based interdisciplinarity is accelerating plant breeding.

Plant breeding is one of the oldest and most important activities accompanying human civilization. During the past thousand years, plant breeding has achieved three significant innovations, each of which derives from introgression of new theories or technologies. These innovations have significantly increased the food supply and allowed for population development. However, with population increases and resource shortages, the world is continuously facing the challenge of food security, which calls for next innovation in plant breeding. Recent technological advances in multiple disciplines have boosted the development of omics, which is accelerating plant breeding. Here, we review the recent advances in omics and discuss our understanding of how interdisciplinary researches will prompt new innovations in plant breeding.

Extensive sequence divergence between the reference genomes of Taraxacum kok-saghyz and Taraxacum mongolicum.

Plants belonging to the genus Taraxacum are widespread all over the world, which contain rubber-producing and non-rubber-producing species. However, the genomic basis underlying natural rubber (NR) biosynthesis still needs more investigation. Here, we presented high-quality genome assemblies of rubber-producing T. kok-saghyz TK1151 and non-rubber-producing T. mongolicum TM5. Comparative analyses uncovered a large number of genetic variations, including inversions, translocations, presence/absence variations, as well as considerable protein divergences between the two species. Two gene duplication events were found in these two Taraxacum species, including one common ancestral whole-genome triplication and one subsequent round of gene amplification. In genomes of both TK1151 and TM5, we identified the genes encoding for each step in the NR biosynthesis pathway and found that the SRPP and CPT gene families have experienced a more obvious expansion in TK1151 compared to TM5. This study will have large-ranging implications for the mechanism of NR biosynthesis and genetic improvement of NR-producing crops.