The role of lysosomal phospholipase A2 in the catabolism of bis(monoacylglycerol)phosphate and association with phospholipidosis.

Bis(monoacylglycerol)phosphate (BMP) is an acidic glycerophospholipid localized to late endosomes and lysosomes. However, the metabolism of BMP is poorly understood. Because many drugs that cause phospholipidosis inhibit lysosomal phospholipase A2 (LPLA2, PLA2G15, LYPLA3) activity, we investigated whether this enzyme has a role in BMP catabolism. The incubation of recombinant human LPLA2 (hLPLA2) and liposomes containing the naturally occurring BMP (sn-(2-oleoyl-3-hydroxy)-glycerol-1-phospho-sn-1'-(2'-oleoyl-3'-hydroxy)-glycerol (S,S-(2,2',C18:1)-BMP) resulted in the deacylation of this BMP isomer. The deacylation rate was 70 times lower than that of dioleoyl phosphatidylglycerol (DOPG), an isomer and precursor of BMP. The release rates of oleic acid from DOPG and four BMP stereoisomers by LPLA2 differed. The rank order of the rates of hydrolysis were DOPG>S,S-(3,3',C18:1)-BMP>R,S-(3,1',C18:1)-BMP>R,R-(1,1',C18:1)>S,S-(2,2')-BMP. The cationic amphiphilic drug amiodarone (AMD) inhibited the deacylation of DOPG and BMP isomers by hLPLA2 in a concentration dependent manner. Under these experimental conditions, the IC50s of amiodarone-induced inhibition of the four BMP isomers and DOPG were less than 20 µM and approximately 30 µM, respectively. BMP accumulation was observed in AMD-treated RAW 264.7 cells. The accumulated BMP was significantly reduced by exogenous treatment of cells with active recombinant hLPLA2 but not with diisopropylfluorophosphate-inactivated recombinant hLPLA2. Finally, a series of cationic amphiphilic drugs known to cause phospholipidosis were screened for inhibition of LPLA2 activity as measured by either the transacylation or fatty acid hydrolysis of BMP or phosphatidylcholine as substrates. Fifteen compounds demonstrated significant inhibition with IC50s ranging from 6.8 to 63.3 µM. These results indicate that LPLA2 degrades BMP isomers with different substrate specificities under acidic conditions and may be the key enzyme associated with BMP accumulation in drug-induced phospholipidosis.

Latitudinal adaptation and dispersal pathway suggested by geographical distribution of transposable elements inserted in the SiPRR37 gene in foxtail millet.

We investigated the variation and geographical distribution of the Pseudo-regulator response 37 (Setaria italica PRR37= SiPRR37) gene involved in heading time (photoperiodism) in foxtail millet. An allele of the SiPRR37 gene, in which an approximately 4.9-kb transposable element (designated as TE1 in this article) is inserted (a loss-of-function or reduction-of-function type), is distributed sporadically in East Asia and broadly in Southeast Asia and South Asia, which implies that this gene is important in latitudinal adaptation. In addition, we found a new allele of SiPRR37 with an insertion of a 360-bp TE (designated as TE2) at this locus and investigated the geographical distribution of this new type. This SiPRR37 allele with TE2 is distributed in Japan, Korea, Nepal, Iran, and Turkey. Both TE1 and TE2 are useful markers for tracing foxtail millet dispersal pathways in Asia.

V-primer: software for the efficient design of genome-wide InDel and SNP markers from multi-sample variant call format (VCF) genotyping data.

DNA markers are indispensable tools in genetics and genomics research as well as in crop breeding, particularly for marker-assisted selection. Recent advances in next-generation sequencing technology have made it easier to obtain genome sequences for various crop species, enabling the large-scale identification of DNA polymorphisms among varieties, which in turn has made DNA marker design more accessible. However, existing primer design software is not suitable for designing many types of genome-wide DNA markers from next-generation sequencing data. Here, we describe the development of V-primer, high-throughput software for designing insertion/deletion, cleaved amplified polymorphic sequence, and single-nucleotide polymorphism (SNP) markers. We validated the applicability of these markers in different crops. In addition, we performed multiplex PCR targeted amplicon sequencing using SNP markers designed with V-primer. Our results demonstrate that V-primer facilitates the efficient and accurate design of primers and is thus a useful tool for genetics, genomics, and crop breeding. V-primer is freely available at https://github.com/ncod3/vprimer.

Modulation of hepatic transcription factor EB activity during cold exposure uncovers direct regulation of bis(monoacylglycero)phosphate lipids by Pla2g15.

Cold exposure is an environmental stress that elicits a rapid metabolic shift in endotherms and is required for survival. The liver provides metabolic flexibility through its ability to rewire lipid metabolism to respond to an increased demand in energy for thermogenesis. We leveraged cold exposure to identify novel lipids contributing to energy homeostasis and found that lysosomal bis(monoacylglycero)phosphate (BMP) lipids were significantly increased in the liver during acute cold exposure. BMP lipid changes occurred independently of lysosomal abundance but were dependent on the lysosomal transcriptional regulator transcription factor EB (TFEB). Knockdown of TFEB in hepatocytes decreased BMP lipid levels. Through molecular biology and biochemical assays, we found that TFEB regulates lipid catabolism during cold exposure and that TFEB knockdown mice were cold intolerant. To identify how TFEB regulates BMP lipid levels, we used a combinatorial approach to identify TFEB target Pla2g15 , a lysosomal phospholipase, as capable of degrading BMP lipids in in vitro liposome assays. Knockdown of Pla2g15 in hepatocytes led to a decrease in BMP lipid species. Together, our studies uncover a required role of TFEB in mediating lipid liver remodeling during cold exposure and identified Pla2g15 as an enzyme that regulates BMP lipid catabolism.

Upcycling rice yield trial data using a weather-driven crop growth model.

Efficient plant breeding plays a significant role in increasing crop yields and attaining food security under climate change. Screening new cultivars through yield trials in multi-environments has improved crop yields, but the accumulated data from these trials has not been effectively upcycled. We propose a simple method that quantifies cultivar-specific productivity characteristics using two regression coefficients: yield-ability (ß) and yield-plasticity (α). The recorded yields of each cultivar are expressed as a unique linear regression in response to the theoretical potential yield (Yp) calculated by a weather-driven crop growth model, called as the "YpCGM method". We apply this to 72510 independent datasets from yield trials of rice that used 237 cultivars measured at 110 locations in Japan over 38 years. The YpCGM method can upcycle accumulated yield data for use in genetic-gain analysis and genome-wide-association studies to guide future breeding programs for developing new cultivars suitable for the world's changing climate.

Whole-genome analysis of recombinant inbred rice lines reveals a quantitative trait locus on chromosome 3 with genotype-by-environment interaction effects.

Elucidating genotype-by-environment interactions is fundamental for understanding the interplay between genetic and environmental factors that shape complex traits in crops. Genotype-by-environment interactions are of practical importance, as they determine the performance of cultivars grown in different environments, prompting the need for an efficient approach for evaluating genotype-by-environment interactions. Here, we describe a method for genotype-by-environment detection that involves comparing linear mixed models. This method successfully detected genotype-by-environment interactions in rice (Oryza sativa) recombinant inbred lines grown at 3 locations. We identified a quantitative trait locus (QTL) on chromosome 3 that was associated with heading date, grain number, and leaf length. The effect of this QTL on plant growth-related traits varied with environmental conditions, indicating the presence of genotype-by-environment interactions. Therefore, our method enables a powerful genotype-by-environment detection pipeline that should facilitate the production of high-yielding crops in a given environment.

Isolation of Pikps, an allele of Pik, from the aus rice cultivar Shoni.

Blast disease caused by the filamentous fungus Pyricularia oryzae (syn. Magnaporthe oryzae) is one of the most destructive diseases of rice (Oryza sativa L.) around the globe. An aus cultivar, Shoni, showed resistance against at least four Japanese P. oryzae isolates. To understand Shoni's resistance against the P. oryzae isolate Naga69-150, genetic analysis was carried out using recombinant inbred lines developed by a cross between Shoni and the japonica cultivar Hitomebore, which is susceptible to Naga69-150. The result indicated that the resistance was controlled by a single locus, which was named Pi-Shoni. A QTL analysis identified Pi-Shoni as being located in the telomeric region of chromosome 11. A candidate gene approach in the region indicated that Pi-Shoni corresponds to the previously cloned Pik locus, and we named this allele Pikps. Loss of gene function mediated by RNA interference demonstrated that a head-to-head-orientated pair of NBS-LRR receptor genes (Pikps-1 and Pikps-2) are required for the Pikps-mediated resistance. Amino acid sequence comparison showed that Pikps-1 is 99% identical to Pikp-1, while Pikps-2 is identical to Pikp-2. Pikps-1 had one amino acid substitution (Pro351Ser) in the NBS domain as compared to Pikp-1. The recognition specificity of Pikps against known AVR-Pik alleles is identical to that of Pikp.

Disentangling the complex gene interaction networks between rice and the blast fungus identifies a new pathogen effector.

Studies focused solely on single organisms can fail to identify the networks underlying host-pathogen gene-for-gene interactions. Here, we integrate genetic analyses of rice (Oryza sativa, host) and rice blast fungus (Magnaporthe oryzae, pathogen) and uncover a new pathogen recognition specificity of the rice nucleotide-binding domain and leucine-rich repeat protein (NLR) immune receptor Pik, which mediates resistance to M. oryzae expressing the avirulence effector gene AVR-Pik. Rice Piks-1, encoded by an allele of Pik-1, recognizes a previously unidentified effector encoded by the M. oryzae avirulence gene AVR-Mgk1, which is found on a mini-chromosome. AVR-Mgk1 has no sequence similarity to known AVR-Pik effectors and is prone to deletion from the mini-chromosome mediated by repeated Inago2 retrotransposon sequences. AVR-Mgk1 is detected by Piks-1 and by other Pik-1 alleles known to recognize AVR-Pik effectors; recognition is mediated by AVR-Mgk1 binding to the integrated heavy metal-associated (HMA) domain of Piks-1 and other Pik-1 alleles. Our findings highlight how complex gene-for-gene interaction networks can be disentangled by applying forward genetics approaches simultaneously to the host and pathogen. We demonstrate dynamic coevolution between an NLR integrated domain and multiple families of effector proteins.

De Novo Transcriptome Analysis Reveals Flowering-Related Genes That Potentially Contribute to Flowering-Time Control in the Japanese Cultivated Gentian Gentiana triflora.

Japanese cultivated gentians are perennial plants that flower in early summer to late autumn in Japan, depending on the cultivar. Several flowering-related genes, including GtFT1 and GtTFL1, are known to be involved in regulating flowering time, but many such genes remain unidentified. In this study, we obtained transcriptome profiling data using the Gentiana triflora cultivar 'Maciry', which typically flowers in late July. We conducted deep RNA sequencing analysis using gentian plants grown under natural field conditions for three months before flowering. To investigate diurnal changes, the plants were sampled at 4 h intervals over 24 h. Using these transcriptome data, we determined the expression profiles of leaves based on homology searches against the Flowering-Interactive Database of Arabidopsis. In particular, we focused on transcription factor genes, belonging to the BBX and MADS-box families, and analyzed their developmental and diurnal variation. The expression levels of representative BBX genes were also analyzed under long- and short-day conditions using in-vitro-grown seedlings, and the expression patterns of some BBX genes differed. Clustering analysis revealed that the transcription factor genes were coexpressed with GtFT1. Overall, these expression profiles will facilitate further analysis of the molecular mechanisms underlying the control of flowering time in gentians.

Rice apoplastic CBM1-interacting protein counters blast pathogen invasion by binding conserved carbohydrate binding module 1 motif of fungal proteins.

When infecting plants, fungal pathogens secrete cell wall-degrading enzymes (CWDEs) that break down cellulose and hemicellulose, the primary components of plant cell walls. Some fungal CWDEs contain a unique domain, named the carbohydrate binding module (CBM), that facilitates their access to polysaccharides. However, little is known about how plants counteract pathogen degradation of their cell walls. Here, we show that the rice cysteine-rich repeat secretion protein OsRMC binds to and inhibits xylanase MoCel10A of the blast fungus pathogen Magnaporthe oryzae, interfering with its access to the rice cell wall and degradation of rice xylan. We found binding of OsRMC to various CBM1-containing enzymes, suggesting that it has a general role in inhibiting the action of CBM1. OsRMC is localized to the apoplast, and its expression is strongly induced in leaves infected with M. oryzae. Remarkably, knockdown and overexpression of OsRMC reduced and enhanced rice defense against M. oryzae, respectively, demonstrating that inhibition of CBM1-containing fungal enzymes by OsRMC is crucial for rice defense. We also identified additional CBM-interacting proteins (CBMIPs) from Arabidopsis thaliana and Setaria italica, indicating that a wide range of plants counteract pathogens through this mechanism.

Genome Analysis Revives a Forgotten Hybrid Crop Edo-dokoro in the Genus Dioscorea.

A rhizomatous Dioscorea crop 'Edo-dokoro' was described in old records of Japan, but its botanical identity has not been characterized. We found that Edo-dokoro is still produced by four farmers in Tohoku-machi of the Aomori prefecture, Japan. The rhizomes of Edo-dokoro are a delicacy to the local people and are sold in the markets. Morphological characters of Edo-dokoro suggest its hybrid origin between the two species, Dioscorea tokoro and Dioscorea tenuipes. Genome analysis revealed that Edo-dokoro likely originated by hybridization of a male D. tokoro to a female D. tenuipes, followed by a backcross with a male plant of D. tokoro. Edo-dokoro is a typical minor crop possibly maintained for more than 300 years but now almost forgotten by the public. We hypothesize that there are many such uncharacterized genetic heritages passed over generations by small-scale farmers that await serious scientific investigation for future use and improvement by using modern genomics information.

A genetically linked pair of NLR immune receptors shows contrasting patterns of evolution.

Throughout their evolution, plant nucleotide-binding leucine-rich-repeat receptors (NLRs) have acquired widely divergent unconventional integrated domains that enhance their ability to detect pathogen effectors. However, the functional dynamics that drive the evolution of NLRs with integrated domains (NLR-IDs) remain poorly understood. Here, we reconstructed the evolutionary history of an NLR locus prone to unconventional domain integration and experimentally tested hypotheses about the evolution of NLR-IDs. We show that the rice (Oryza sativa) NLR Pias recognizes the effector AVR-Pias of the blast fungal pathogen Magnaporthe oryzae. Pias consists of a functionally specialized NLR pair, the helper Pias-1 and the sensor Pias-2, that is allelic to the previously characterized Pia pair of NLRs: the helper RGA4 and the sensor RGA5. Remarkably, Pias-2 carries a C-terminal DUF761 domain at a similar position to the heavy metal-associated (HMA) domain of RGA5. Phylogenomic analysis showed that Pias-2/RGA5 sensor NLRs have undergone recurrent genomic recombination within the genus Oryza, resulting in up to six sequence-divergent domain integrations. Allelic NLRs with divergent functions have been maintained transspecies in different Oryza lineages to detect sequence-divergent pathogen effectors. By contrast, Pias-1 has retained its NLR helper activity throughout evolution and is capable of functioning together with the divergent sensor-NLR RGA5 to respond to AVR-Pia. These results suggest that opposite selective forces have driven the evolution of paired NLRs: highly dynamic domain integration events maintained by balancing selection for sensor NLRs, in sharp contrast to purifying selection and functional conservation of immune signaling for helper NLRs.

High-performance pipeline for MutMap and QTL-seq.

Summary: Bulked segregant analysis implemented in MutMap and QTL-seq is a powerful and efficient method to identify loci contributing to important phenotypic traits. However, the previous pipelines were not user-friendly to install and run. Here, we describe new pipelines for MutMap and QTL-seq. These updated pipelines are approximately 5-8 times faster than the previous pipeline, are easier for novice users to use, and can be easily installed through bioconda with all dependencies. Availability: The new pipelines of MutMap and QTL-seq are written in Python and can be installed via bioconda. The source code and manuals are available online (MutMap: https://github.com/YuSugihara/MutMap, QTL-seq: https://github.com/YuSugihara/QTL-seq).

Recombinant inbred lines and next-generation sequencing enable rapid identification of candidate genes involved in morphological and agronomic traits in foxtail millet.

We constructed recombinant inbred lines (RILs) between a Japanese and a Taiwanese landrace of foxtail millet and employed next-generation sequencing, such as flexible ddRAD-seq and Nanopore sequencing to identify the candidate genes involved in the crop evolution of foxtail millet. We successfully constructed a linkage map using flexible ddRAD-seq with parents and RILs and detected major QTLs for each of three traits: leaf sheath colors, spikelet-tipped bristles (stb), and days to heading (DTH). (1) For leaf sheath colors, we identified the C gene on chromosome IV. (2) We identified a homeobox (HOX14) gene for stb on chromosome II, which shows homology with HvVrs1 in barley. (3) Finally, we identified a QTL with a large effect on DTH on chromosome II. A parent of the RILs from Taiwan and Yugu1 had a Harbinger-like TE in intron 3 of this gene. We also investigated the geographical distribution of the TE insertion type of this gene and found that the insertion type is distributed in the northern part of East Asia and intensively in South and Southeast Asia, suggesting that loss/reduction of function of this gene plays an important role in spreading into the northern part of East Asia and subtropical and tropical zones.

Sat-BSA: an NGS-based method using local de novo assembly of long reads for rapid identification of genomic structural variations associated with agronomic traits.

Advances in next generation sequencing (NGS)-based methodologies have accelerated the identifications of simple genetic variants such as point mutations and small insertions/deletions (InDels). Structural variants (SVs) including large InDels and rearrangements provide vital sources of genetic diversity for plant breeding. However, their analysis remains a challenge due to their complex nature. Consequently, novel NGS-based approaches are needed to rapidly and accurately identify SVs. Here, we present an NGS-based bulked-segregant analysis (BSA) technique called Sat-BSA (SVs associated with traits) for identifying SVs controlling traits of interest in crops. Sat-BSA targets allele frequencies at all SNP positions to first identify candidate genomic regions associated with a trait, which is then reconstructed by long reads-based local de novo assembly. Finally, the association between SVs, RNA-seq-based gene expression patterns and trait is evaluated for multiple cultivars to narrow down the candidate genes. We applied Sat-BSA to segregating F2 progeny obtained from crosses between turnip cultivars with different tuber colors and successfully isolated two genes harboring SVs that are responsible for tuber phenotypes. The current study demonstrates the utility of Sat-BSA for the identification of SVs associated with traits of interest in species with large and heterozygous genomes.

Inhibition of lysosomal phospholipase A2 predicts drug-induced phospholipidosis.

Phospholipidosis, the excessive accumulation of phospholipids within lysosomes, is a pathological response observed following exposure to many drugs across multiple therapeutic groups. A clear mechanistic understanding of the causes and implications of this form of drug toxicity has remained elusive. We previously reported the discovery and characterization of a lysosome-specific phospholipase A2 (PLA2G15) and later reported that amiodarone, a known cause of drug-induced phospholipidosis, inhibits this enzyme. Here, we assayed a library of 163 drugs for inhibition of PLA2G15 to determine whether this phospholipase was the cellular target for therapeutics other than amiodarone that cause phospholipidosis. We observed that 144 compounds inhibited PLA2G15 activity. Thirty-six compounds not previously reported to cause phospholipidosis inhibited PLA2G15 with IC50 values less than 1 mM and were confirmed to cause phospholipidosis in an in vitro assay. Within this group, fosinopril was the most potent inhibitor (IC50 0.18 µM). Additional characterization of the inhibition of PLA2G15 by fosinopril was consistent with interference of PLA2G15 binding to liposomes. PLA2G15 inhibition was more accurate in predicting phospholipidosis compared with in silico models based on pKa and ClogP, measures of protonation, and transport-independent distribution in the lysosome, respectively. In summary, PLA2G15 is a primary target for cationic amphiphilic drugs that cause phospholipidosis, and PLA2G15 inhibition by cationic amphiphilic compounds provides a potentially robust screening platform for potential toxicity during drug development.

RIL-StEp: epistasis analysis of rice recombinant inbred lines reveals candidate interacting genes that control seed hull color and leaf chlorophyll content.

Characterizing epistatic gene interactions is fundamental for understanding the genetic architecture of complex traits. However, due to the large number of potential gene combinations, detecting epistatic gene interactions is computationally demanding. A simple, easy-to-perform method for sensitive detection of epistasis is required. Due to their homozygous nature, use of recombinant inbred lines excludes the dominance effect of alleles and interactions involving heterozygous genotypes, thereby allowing detection of epistasis in a simple and interpretable model. Here, we present an approach called RIL-StEp (recombinant inbred lines stepwise epistasis detection) to detect epistasis using single-nucleotide polymorphisms in the genome. We applied the method to reveal epistasis affecting rice (Oryza sativa) seed hull color and leaf chlorophyll content and successfully identified pairs of genomic regions that presumably control these phenotypes. This method has the potential to improve our understanding of the genetic architecture of various traits of crops and other organisms.

Genome analyses reveal the hybrid origin of the staple crop white Guinea yam (Dioscorea rotundata).

White Guinea yam (Dioscorea rotundata) is an important staple tuber crop in West Africa. However, its origin remains unclear. In this study, we resequenced 336 accessions of white Guinea yam and compared them with the sequences of wild Dioscorea species using an improved reference genome sequence of D. rotundata In contrast to a previous study suggesting that D. rotundata originated from a subgroup of Dioscorea praehensilis, our results suggest a hybrid origin of white Guinea yam from crosses between the wild rainforest species D. praehensilis and the savannah-adapted species Dioscorea abyssinica We identified a greater genomic contribution from D. abyssinica in the sex chromosome of Guinea yam and extensive introgression around the SWEETIE gene. Our findings point to a complex domestication scenario for Guinea yam and highlight the importance of wild species as gene donors for improving this crop through molecular breeding.

Optimization of Eliglustat-Based Glucosylceramide Synthase Inhibitors as Substrate Reduction Therapy for Gaucher Disease Type 3.

There remain no approved therapies for rare but devastating neuronopathic glyocosphingolipid storage diseases, such as Sandhoff, Tay-Sachs, and Gaucher disease type 3. We previously reported initial optimization of the scaffold of eliglustat, an approved therapy for the peripheral symptoms of Gaucher disease type 1, to afford 2, which effected modest reductions in brain glucosylceramide (GlcCer) in normal mice at 60 mg/kg. The relatively poor pharmacokinetic properties and high Pgp-mediated efflux of 2 prompted further optimization of the scaffold. With a general objective of reducing topological polar surface area, and guided by multiple metabolite identification studies, we were successful at identifying 17 (CCG-222628), which achieves remarkably greater brain exposure in mice than 2. After demonstrating an over 60-fold improvement in potency over 2 at reducing brain GlcCer in normal mice, we compared 17 with Sanofi clinical candidate venglustat (Genz-682452) in the CBE mouse model of Gaucher disease type 3. At doses of 10 mg/kg, 17 and venglustat effected comparable reductions in both brain GlcCer and glucosylsphingosine. Importantly, 17 achieved these equivalent pharmacodynamic effects at significantly lower brain exposure than venglustat.

A novel WD40-repeat protein involved in formation of epidermal bladder cells in the halophyte quinoa.

Halophytes are plants that grow in high-salt environments and form characteristic epidermal bladder cells (EBCs) that are important for saline tolerance. To date, however, little has been revealed about the formation of these structures. To determine the genetic basis for their formation, we applied ethylmethanesulfonate mutagenesis and obtained two mutants with reduced levels of EBCs (rebc) and abnormal chloroplasts. In silico subtraction experiments revealed that the rebc phenotype was caused by mutation of REBC, which encodes a WD40 protein that localizes to the nucleus and chloroplasts. Phylogenetic and transformant analyses revealed that the REBC protein differs from TTG1, a WD40 protein involved in trichome formation. Furthermore, rebc mutants displayed damage to their shoot apices under abiotic stress, suggesting that EBCs may protect the shoot apex from such stress. These findings will help clarify the mechanisms underlying EBC formation and function.