Recent studies have highlighted the significance of cellular metabolism in the initiation of clonal expansion and effector differentiation of T cells. Upon exposure to antigens, naïve CD4+ T cells undergo metabolic reprogramming to meet their metabolic requirements. However, only few studies have simultaneously evaluated the changes in protein and metabolite levels during T cell differentiation. Our research seeks to fill the gap by conducting a comprehensive analysis of changes in levels of metabolites, including sugars, amino acids, intermediates of the TCA cycle, fatty acids, and lipids. By integrating metabolomics and proteomics data, we discovered that the quantity and composition of cellular lipids underwent significant changes in different effector Th cell subsets. Especially, we found that the sphingolipid biosynthesis pathway was commonly activated in Th1, Th2, Th17, and iTreg cells and that inhibition of this pathway led to the suppression of Th17 and iTreg cells differentiation. Additionally, we discovered that Th17 and iTreg cells enhance glycosphingolipid metabolism, and inhibition of this pathway also results in the suppression of Th17 and iTreg cell generation. These findings demonstrate that the utility of our combined metabolomics and proteomics analysis in furthering the understanding of metabolic transition during Th cell differentiation.

Introduction Comprehensive genomic resources offer insights into biological features, including traits/disease-related genetic loci. The current reference genome assembly for the domestic cat (Felis catus), Felis_Catus_9.0 (felCat9), derived from sequences of the Abyssinian cat, may inadequately represent the general cat population, limiting the extent of deducible genetic variations. Objectives The goal was to develop Anicom American Shorthair 1.0 (AnAms1.0), a reference-grade chromosome-scale cat genome assembly. Methods In contrast to prior assemblies relying on Abyssinian cat sequences, AnAms1.0 was constructed from the sequences of more popular American Shorthair breed, which is related to more breeds than the Abyssinian cat. By combining advanced genomics technologies, including PacBio long-read sequencing and Hi-C- and optical mapping data-based sequence scaffolding, we compared AnAms1.0 to existing Felidae genome assemblies (20 scaffolds, scaffolds N50 > 150 Mbp). Homology-based and ab initio gene annotation through Iso-Seq and RNA-Seq was used to identify new coding genes and splice variants. Results AnAms1.0 demonstrated superior contiguity and accuracy than existing Felidae genome assemblies. Using AnAms1.0, we identified over 1.5 thousand structural variants and 29 million repetitions compared to felCat9. Additionally, we identified > 1,600 novel protein-coding genes. Notably, olfactory receptor structural variants and cardiomyopathy-related variants were identified. Conclusion AnAms1.0 facilitates the discovery of novel genes related to normal and disease phenotypes in domestic cats. The analyzed data are publicly accessible on Cats-I (https://cat.annotation.jp/), which we established as a platform for accumulating and sharing genomic resources to discover novel genetic traits and advance veterinary medicine.

Functionally distinct CD4+ helper T (Th) cell subsets, including Th1, Th2, Th17, and regulatory T cells (Treg), play a pivotal role in the regulation of acquired immunity. Although the key proteins involved in the regulation of Th cell differentiation have already been identified how the proteogenomic landscape changes during the Th cell activation remains unclear. To address this issue, we characterized proteogenomic signatures of differentiation to each Th cell subsets by RNA sequencing and liquid chromatography-assisted mass spectrometry, which enabled us to simultaneously quantify more than 10,000 protein-coding transcripts and 8,000 proteins in a single-shot. The results indicated that T cell receptor activation affected almost half of the transcript and protein levels in a low correlative and gene-specific manner, and specific cytokine treatments modified the transcript and protein profiles in a manner specific to each Th cell subsets: Th17 and Tregs particularly exhibited unique proteogenomic signatures compared to other Th cell subsets. Interestingly, the in-depth proteome data revealed that mRNA profiles alone were not enough to delineate functional changes during Th cell activation, suggesting that the proteogenomic dataset obtained in this study serves as a unique and indispensable data resource for understanding the comprehensive molecular mechanisms underlying effector Th cell differentiation.

Metabolic fluxes involving fatty acid biosynthesis play essential roles in controlling the differentiation of T helper 17 (T H 17) cells. However, the exact enzymes and lipid metabolites involved, as well as their link to promoting the core gene transcriptional signature required for the differentiation of T H 17 cells, remain largely unknown. From a pooled CRISPR-based screen and unbiased lipidomics analyses, we identified that 1-oleoyl-lysophosphatidylethanolamine could act as a lipid modulator of retinoid-related orphan receptor gamma t (RORγt) activity in T H 17 cells. In addition, we specified five enzymes, including Gpam , Gpat3 , Lplat1 , Pla2g12a , and Scd2 , suggestive of the requirement of glycerophospholipids with monounsaturated fatty acids being required for the transcription of Il17a . 1-Oleoyl-lysophosphatidylethanolamine was reduced in Pla2g12a -deficient T H 17 cells, leading to the abolition of interleukin-17 (IL-17) production and disruption to the core transcriptional program required for the differentiation of T H 17 cells. Furthermore, mice with T cell–specific deficiency of Pla2g12a failed to develop disease in an experimental autoimmune encephalomyelitis model of multiple sclerosis. Thus, our data indicate that 1-oleoyl-lysophosphatidylethanolamine is a lipid metabolite that promotes RORγt-induced T H 17 cell differentiation and the pathogenicity of T H 17 cells.

Type I interferons (type I-IFN) are critical for the host defense to viral infection, and at the same time, the dysregulation of type I-IFN responses leads to autoinflammation or autoimmunity. Recently, we reported that the decrease in monounsaturated fatty acid caused by the genetic deletion of Scd2 is essential for the activation of type I-IFN signaling in CD4+ Th1 cells. Although interferon regulatory factor (IRF) is a family of homologous proteins that control the transcription of type I-IFN and interferon stimulated genes (ISGs), the member of the IRF family that is responsible for the type I-IFN responses induced by targeting of SCD2 remains unclear. Here, we report that the deletion of Scd2 triggered IRF3 activation for type I-IFN production, resulting in the nuclear translocation of IRF9 to induce ISG transcriptome in Th1 cells. These data led us to hypothesize that IRF9 plays an essential role in the transcriptional regulation of ISGs in Scd2-deleted (sgScd2) Th1 cells. By employing ChIP-seq analyses, we found a substantial percentage of the IRF9 target genes were shared by sgScd2 and IFNβ-treated Th1 cells. Importantly, our detailed analyses identify a unique feature of IRF9 binding in sgScd2 Th1 cells that were not observed in IFNβ-treated Th1 cells. In addition, our combined analyses of transcriptome and IRF9 ChIP-seq revealed that the autoimmunity related genes, which increase in patient with SLE, were selectively increased in sgScd2 Th1 cells. Thus, our findings provide novel mechanistic insights into the process of fatty acid metabolism that is essential for the type I-IFN response and the activation of the IRF family in CD4+ T cells.

Host lipid metabolism and viral responses are intimately connected. However, the process by which the acquired immune systems adapts lipid metabolism to meet demands, and whether or not the metabolic rewiring confers a selective advantage to host immunity, remains unclear. Here we show that viral infection attenuates the expression of genes related to lipid metabolism in murine CD4+ T cells, which in turn increases the expression of antiviral genes. Inhibition of the fatty acid synthesis pathway substantially increases the basal expression of antiviral genes via the spontaneous production of type I interferon (IFN). Using a combination of CRISPR/Cas9-mediated genome editing technology and a global lipidomics analysis, we found that the decrease in monounsaturated fatty acid caused by genetic deletion of Scd2 in mice was crucial for the induction of an antiviral response through activation of the cGAS-STING pathway. These findings demonstrate the important relationship between fatty acid biosynthesis and type I IFN responses that enhances the antiviral response.

Cultivated strawberry ( Fragaria × ananassa ) is an octoploid species (2n = 8x= 56) that is widely consumed around the world as both fresh and processed fruit. In this study, we report a chromosome-scale strawberry genome assembly of a Japanese variety, Reikou. The Illumina short reads derived from paired-end, mate-pair, and 10X Genomics libraries were assembled using Denovo MAGIC 3.0. The generated phased scaffolds consisted of 32,715 sequences with a total length of 1.4 Gb and an N50 length of 3.9 Mb. A total of 63 pseudomolecules including chr0 were created by aligning the scaffolds onto the Reikou S1 linkage maps with the IStraw90 Axiom SNP array and ddRAD-Seq. Meanwhile, genomes of diploid Fragaria species were resequenced and compared with the most similar chromosome-scale scaffolds to investigate the possible progenitor of each subgenome. Clustering analysis suggested that the most likely progenitors were F. vesca and F. iinumae . The phased pseudomolecules were assigned the scaffolds names with Av, Bi, and X, representing sequence similarity with F. vesca (Av), F. iinumae (Bi), and others (X), respectively. The result of a comparison with the Camerosa genome suggested the possibility of subgenome structure differences between the two varieties.

T cells possess distinguishing effector functions and drive inflammatory disorders. We have previously identified IL-5–producing Th2 cells as the pathogenic population predominantly involved in the pathology of allergic inflammation. However, the cell-intrinsic signaling pathways that control the pathogenic Th2 cell function are still unclear. We herein report the high expression of acetyl-CoA carboxylase 1 (ACC1) in the pathogenic CD4+ T cell population in the lung and skin. The genetic deletion of CD4+ T cell–intrinsic ACC1 dampened eosinophilic and basophilic inflammation in the lung and skin by constraining IL-5 or IL-3 production. Mechanistically, ACC1-dependent fatty acid biosynthesis induces the pathogenic cytokine production of CD4+ T cells via metabolic reprogramming and the availability of acetyl-CoA for epigenetic regulation. We thus identified a distinct phenotype of the pathogenic T cell population in the lung and skin, and ACC1 was shown to be an essential regulator controlling the pathogenic function of these populations to promote type 2 inflammation.

Regulatory T (Treg) cells are critical for immunological tolerance and immune homeostasis. Treg cells strongly rely on mitochondrial metabolism and show a lower level of glycolysis. However, little is known about the role of lipid metabolism in the regulation of Treg cell homeostasis. Some members of the ACSL family of acyl-coenzyme A (CoA) synthases are expressed in T cells, but their function remains unclear. A combination of RNA-sequencing and proteome analyses shows that Acsbg1, a member of ACSL, is selectively expressed in Treg cells. We show that the genetic deletion of Acsbg1 not only causes mitochondrial dysfunction, but it also dampens other metabolic pathways. The extrinsic supplementation of Acsbg1-deficient Treg cells with oleoyl-CoA restores the phenotype of the Treg metabolic signature. Furthermore, this pathway in ST2+ effector Treg cells enhances immunosuppressive capacity in airway inflammation. Thus, Acsbg1 serves as a metabolic checkpoint governing Treg cell homeostasis and the resolution of lung inflammation.

The domestic cat ( Felis catus ) is one of the most popular companion animals in the world. Comprehensive genomic resources will aid the development and application of veterinary medicine including to improve feline health, in particular, to enable precision medicine which is promising in human application. However, currently available cat genome assemblies were mostly built based on the Abyssinian cat breed which is highly inbred and has limited power in representing the vast diversity of the cat population. Moreover, the current reference assembly remains fragmented with sequences contained in thousands of scaffolds. We constructed a reference-grade chromosome-scale genome assembly of a domestic cat, Felis catus genome of American Shorthair breed, Anicom American shorthair 1.0 (AnAms1.0) with high contiguity (scaffold N50 > 120 Mb), by combining multiple advanced genomic technologies, including PacBio long-read sequencing as well as sequence scaffolding by long-range genomic information obtained from Hi-C and optical mapping data. Homology-based and ab initio gene annotation was performed with the Iso-Seq data. Analyzed data is be publicly accessible on Cats genome informatics (Cats-I, https://cat.annotation.jp/ ), a cat genome database established as a platform to facilitate the accumulation and sharing of genomic resources to improve veterinary care.

Ipomoea trifida (H. B. K.) G. Don. is the most likely diploid ancestor of the hexaploid sweet potato, I. batatas (L.) Lam. To assist in analysis of the sweet potato genome, de novo whole-genome sequencing was performed with two lines of I. trifida, namely the selfed line Mx23Hm and the highly heterozygous line 0431-1, using the Illumina HiSeq platform. We classified the sequences thus obtained as either ‘core candidates’ (common to the two lines) or ‘line specific’. The total lengths of the assembled sequences of Mx23Hm (ITR_r1.0) was 513 Mb, while that of 0431-1 (ITRk_r1.0) was 712 Mb. Of the assembled sequences, 240 Mb (Mx23Hm) and 353 Mb (0431-1) were classified into core candidate sequences. A total of 62,407 (62.4 Mb) and 109,449 (87.2 Mb) putative genes were identified, respectively, in the genomes of Mx23Hm and 0431-1, of which 11,823 were derived from core sequences of Mx23Hm, while 28,831 were from the core candidate sequence of 0431-1. There were a total of 1,464,173 single-nucleotide polymorphisms and 16,682 copy number variations (CNVs) in the two assembled genomic sequences (under the condition of log2 ratio of >1 and CNV size >1,000 bases). The results presented here are expected to contribute to the progress of genomic and genetic studies of I. trifida, as well as studies of the sweet potato and the genus Ipomoea in general.

To develop a high density linkage map in faba bean, a total of 1,363 FBES (Faba bean expressed sequence tag [EST]-derived simple sequence repeat [SSR]) markers were designed based on 5,090 non-redundant ESTs developed in this study. A total of 109 plants of a 'Nubaria 2' × 'Misr 3' F2 mapping population were used for map construction. Because the parents were not pure homozygous lines, the 109 F2 plants were divided into three subpopulations according to the original F1 plants. Linkage groups (LGs) generated in each subpopulation were integrated by commonly mapped markers. The integrated 'Nubaria 2' × 'Misr 3' map consisted of six LGs, representing a total length of 684.7 cM, with 552 loci. Of the mapped loci, 47% were generated from multi-loci diagnostic (MLD) markers. Alignment of homologous sequence pairs along each linkage group revealed obvious syntenic relationships between LGs in faba bean and the genomes of two model legumes, Lotus japonicus and Medicago truncatula. In a polymorphic analysis with ten Egyptian faba bean varieties, 78.9% (384/487) of the FBES markers showed polymorphisms. Along with the EST-SSR markers, the dense map developed in this study is expected to accelerate marker assisted breeding in faba bean.

Cultivated strawberry (Fragaria x ananassa) is octoploid and shows allogamous behaviour. The present study aims at dissecting this octoploid genome through comparison with its wild relatives, F. iinumae, F. nipponica, F. nubicola, and F. orientalis by de novo whole-genome sequencing on an Illumina and Roche 454 platforms. The total length of the assembled Illumina genome sequences obtained was 698 Mb for F. x ananassa, and ∼200 Mb each for the four wild species. Subsequently, a virtual reference genome termed FANhybrid_r1.2 was constructed by integrating the sequences of the four homoeologous subgenomes of F. x ananassa, from which heterozygous regions in the Roche 454 and Illumina genome sequences were eliminated. The total length of FANhybrid_r1.2 thus created was 173.2 Mb with the N50 length of 5137 bp. The Illumina-assembled genome sequences of F. x ananassa and the four wild species were then mapped onto the reference genome, along with the previously published F. vesca genome sequence to establish the subgenomic structure of F. x ananassa. The strategy adopted in this study has turned out to be successful in dissecting the genome of octoploid F. x ananassa and appears promising when applied to the analysis of other polyploid plant species.

The cultivated strawberry (Fragaria × ananassa) is an octoploid (2n = 8x = 56) of the Rosaceae family whose genomic architecture is still controversial. Several recent studies support the AAA'A'BBB'B' model, but its complexity has hindered genetic and genomic analysis of this important crop. To overcome this difficulty and to assist genome-wide analysis of F. × ananassa, we constructed an integrated linkage map by organizing a total of 4474 of simple sequence repeat (SSR) markers collected from published Fragaria sequences, including 3746 SSR markers [Fragaria vesca expressed sequence tag (EST)-derived SSR markers] derived from F. vesca ESTs, 603 markers (F. × ananassa EST-derived SSR markers) from F. × ananassa ESTs, and 125 markers (F. × ananassa transcriptome-derived SSR markers) from F. × ananassa transcripts. Along with the previously published SSR markers, these markers were mapped onto five parent-specific linkage maps derived from three mapping populations, which were then assembled into an integrated linkage map. The constructed map consists of 1856 loci in 28 linkage groups (LGs) that total 2364.1 cM in length. Macrosynteny at the chromosome level was observed between the LGs of F. × ananassa and the genome of F. vesca. Variety distinction on 129 F. × ananassa lines was demonstrated using 45 selected SSR markers.

The quality of peanut oil depends on the properties of the oil, one of which is the oleic/linoleic acid ratio (O/L ratio). The two genes related to the O/L ratio, ahFAD2A and ahFAD2B, were previously identified in peanut. This chapter demonstrates marker-assisted backcrossing selection (MABS) for a high O/L ratio in peanut with a breeding line showing a high O/L ratio and an elite line with medium O/L ratio as donor and backcross parents, respectively. The number of DNA markers developed was considerably lower at that time compared to the present. Polymorphic DNA markers and a linkage map were developed in parallel with the MABS program. The ideal scheme of MABS for high O/L would be integration of the schemes used in the two studies: two backcrossings per year with genotyping of genome-wide (GW) and FAD2 loci markers in each generation, then selfing with the genotyping of GW markers.

Background: Peanut (Arachis hypogaea) is an autogamous allotetraploid legume (2n = 4x = 40) that is widely cultivated as a food and oil crop. More than 6,000 DNA markers have been developed in Arachis spp., but high-density linkage maps useful for genetics, genomics, and breeding have not been constructed due to extremely low genetic diversity. Polymorphic marker loci are useful for the construction of such high-density linkage maps. The present study used in silico analysis to develop simple sequence repeat-based and transposon-based markers.Results: The use of in silico analysis increased the efficiency of polymorphic marker development by more than 3-fold. In total, 926 (34.2%) of 2,702 markers showed polymorphisms between parental lines of the mapping population. Linkage analysis of the 926 markers along with 253 polymorphic markers selected from 4,449 published markers generated 21 linkage groups covering 2,166.4 cM with 1,114 loci. Based on the map thus produced, 23 quantitative trait loci (QTLs) for 15 agronomical traits were detected. Another linkage map with 326 loci was also constructed and revealed a relationship between the genotypes of the FAD2 genes and the ratio of oleic/linoleic acid in peanut seed.Conclusions: In silico analysis of polymorphisms increased the efficiency of polymorphic marker development, and contributed to the construction of high-density linkage maps in cultivated peanut. The resultant maps were applicable to QTL analysis. Marker subsets and linkage maps developed in this study should be useful for genetics, genomics, and breeding in Arachis. The data are available at the Kazusa DNA Marker Database (http://marker.kazusa.or.jp).

Large-scale development of expressed sequence tag simple sequence repeat (EST-SSR) markers was performed in peanut (Arachis hypogaea L.) to obtain more informative genetic markers. A total of 10,102 potential non-redundant EST sequences, including 3,445 contigs and 6,657 singletons, were generated from cDNA libraries of the gynophore, roots, leaves and seedlings. A total of 3,187 primer pairs were designed on flanking regions of SSRs, some of which allowed one and two base mismatches. Among the 3,187 markers generated, 2,540 (80%) were trinucleotide repeats, 302 (9%) were dinucleotide repeats, and 345 (11%) were tetranucleotide repeats. Pre-polymorphic analyses of 24 Arachis accessions were performed using 10% polyacrylamide gels. A total of 1,571 EST-SSR markers showing clear polymorphisms were selected for further polymorphic analysis with a Fluoro-fragment Analyzer. The 16 Arachis accessions examined included cultivated peanut varieties as well as diploid species with the A or B genome. Altogether 1,281 (81.5%) of the 1,571 markers were polymorphic among the 16 accessions, and 366 (23.3%) were polymorphic among the 12 cultivated varieties. Diversity analysis was performed and the genotypes of all 16 Arachis accessions showed similarity coefficients ranging from 0.37 to 0.97.

White clover (Trifolium repens L.) is an allotetraploid species (2n - 4X - 32) that is widely distributed in temperate regions and cultivated as a forage legume. In this study, we developed expressed sequence tag (EST)-derived simple sequence repeat (SSR) markers, constructed linkage maps, and performed comparative mapping with other legume species. A total of 7982 ESTs that could be assembled into 5400 contigs and 2582 singletons were generated. Using the EST sequences that were obtained, 1973 primer pairs to amplify EST-derived SSR markers were designed and used for linkage analysis of 188 F-1 progenies, which were generated by a cross between two Japanese plants, '273-7' and 'T17-349,' with previously published SSR markers. An integrated linkage map was constructed by combining parental-specific maps, which consisted of 1743 SSR loci on 16 homeologous linkage groups with a total length of 2511 cM. The primer sequences of the developed EST-SSR markers and their map positions are available on http://clovergarden.jp/. Linkage disequilibrium (LD) was observed on 9 of 16 linkage groups of a parental-specific map. The genome structures were compared among white clover, red clover (T. pratense L.), Medicago truncatula, and Lotus japonicus. Macrosynteny was observed across the four legume species. Surprisingly, the comparative genome structure between white clover and M. truncatula had a higher degree of conservation than that of the two clover species.

The genetic information in the genome of Eucalyptus camaldulensis was investigated by sequencing the genome and the cDNA using a combination of the conventional Sanger method and next-generation sequencing methods, followed by intensive bioinformatics analyses. The total length of the non-redundant genomic sequences thus obtained was 654,922,307 bp consisting of 81,246 scaffolds and 121,194 singlets. These sequences accounted for approximately 92% of the gene-containing regions with an average G + C content of 33.6%. A total of 77,121 complete and partial structures of protein-encoding genes have been deduced. Comparison of the genes mapped on the KEGG pathways or located in the KOG classification with those in other plant species revealed the characteristics of the E. camaldulensis genome, and it was found that 23 pathways contained enzymes present only in the E. camaldulensis genome. Polymorphism analysis using microsatellite markers developed from the genomic sequence data obtained was performed with six Eucalyptus species collected from various parts of the world to estimate their genetic diversity, and the usefulness of these markers was demonstrated. The genomic sequence and accompanying information presented here are expected to serve as valuable resources for the acceleration of fundamental and applied research with Eucalyptus, especially in the fields of paper production and industrial materials. Further information on the genomic and cDNA sequences and microsatellite markers is available at http://www.kazusa.or.jp/eucaly/.

Simple sequence repeat (SSR) markers provide a powerful tool for genetic linkage map construction that can be applied for identification of quantitative trait loci (QTL). In this study, a total of 640 new SSR markers were developed from an enriched genomic DNA library of the cassava variety ‘Huay Bong 60’ and 1,500 novel expressed sequence tag-simple sequence repeat (EST-SSR) loci were developed from the Genbank database. To construct a genetic linkage map of cassava, a 100 F1 line mapping population was developed from the cross Huay Bong 60 by ‘Hanatee’. Polymorphism screening between the parental lines revealed that 199 SSRs and 168 EST-SSRs were identified as novel polymorphic markers. Combining with previously developed SSRs, we report a linkage map consisted of 510 markers encompassing 1,420.3 cM, distributed on 23 linkage groups with a mean distance between markers of 4.54 cM. Comparison analysis of the SSR order on the cassava linkage map and the cassava genome sequences allowed us to locate 284 scaffolds on the genetic map. Although the number of linkage groups reported here revealed that this F1 genetic linkage map is not yet a saturated map, it encompassed around 88% of the cassava genome indicating that the map was almost complete. Therefore, sufficient markers now exist to encompass most of the genomes and efficiently map traits in cassava.