Development of EST-SSR markers and construction of a linkage map in faba bean (Vicia faba).

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.

The integration of metabolic and proteomic data uncovers an augmentation of the sphingolipid biosynthesis pathway during T-cell differentiation.

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.

1-Oleoyl-lysophosphatidylethanolamine stimulates RORγt activity in TH17 cells.

Metabolic fluxes involving fatty acid biosynthesis play essential roles in controlling the differentiation of T helper 17 (TH17) 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 TH17 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 TH17 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 TH17 cells, leading to the abolition of interleukin-17 (IL-17) production and disruption to the core transcriptional program required for the differentiation of TH17 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 TH17 cell differentiation and the pathogenicity of TH17 cells.

Characterization of proteogenomic signatures of differentiation of CD4+ T cell subsets.

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.

In silico polymorphism analysis for the development of simple sequence repeat and transposon markers and construction of linkage map in cultivated peanut.

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).

SCD2-mediated cooperative activation of IRF3-IRF9 regulatory circuit controls type I interferon transcriptome in CD4+ T 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.

SSR and EST-SSR-based genetic linkage map of cassava (Manihot esculenta Crantz).

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 F(1) 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 F(1) 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.

ACC1-expressing pathogenic T helper 2 cell populations facilitate lung and skin inflammation in mice.

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.

Acsbg1-dependent mitochondrial fitness is a metabolic checkpoint for tissue Treg cell homeostasis.

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.

SCD2-mediated monounsaturated fatty acid metabolism regulates cGAS-STING-dependent type I IFN responses 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.

An interspecific linkage map of SSR and intronic polymorphism markers in tomato.

Despite the collection and availability of abundant tomato genome sequences, PCR-based markers adapted to large scale analysis have not been developed in tomato species. Therefore, using public genome sequence data in tomato, we developed three types of DNA markers: expressed sequence tag (EST)-derived simple sequence repeat (SSR) markers (TES markers), genome-derived SSR markers (TGS markers) and EST-derived intronic polymorphism markers (TEI markers). A total of 2,047 TES, 3,510 TGS and 674 TEI markers were established and used in the polymorphic analysis of a cultivated tomato (Solanum lycopersicum) 'LA925' and its wild relative Solanum pennellii 'LA716', parents of the Tomato-EXPEN 2000 mapping population. The polymorphic ratios between parents revealed by the TES, TGS and TEI markers were 37.3, 22.6 and 80.0%, respectively. Those showing polymorphisms were used to genotype the Tomato-EXPEN 2000 mapping population, and a high-density genetic linkage map composed of 1,433 new and 683 existing marker loci was constructed on 12 chromosomes, covering 1,503.1 cM. In the present map, 48% of the mapped TGS loci were located within heterochromatic regions, while 18 and 21% of TES and TEI loci, respectively, were located in heterochromatin. The large number of SSR and SNP markers developed in this study provide easily handling genomic tools for molecular breeding in tomato. Information on the DNA markers developed in this study is available at http://www.kazusa.or.jp/tomato/.

Construction of a consensus linkage map for red clover (Trifolium pratense L.).

BACKGROUND: Red clover (Trifolium pratense L.) is a major forage legume that has a strong self-incompatibility system and exhibits high genetic diversity within populations. For several crop species, integrated consensus linkage maps that combine information from multiple mapping populations have been developed. For red clover, three genetic linkage maps have been published, but the information in these existing maps has not been integrated. RESULTS: A consensus linkage map was constructed using six mapping populations originating from eight parental accessions. Three of the six mapping populations were established for this study. The integrated red clover map was composed of 1804 loci, including 1414 microsatellite loci, 181 amplified fragment length polymorphism (AFLP) loci and 204 restriction fragment length polymorphism (RFLP) loci, in seven linkage groups. The average distance between loci and the total length of the consensus map were 0.46 cM and 836.6 cM, respectively. The locus order on the consensus map correlated highly with that of accession-specific maps. Segregation distortion was observed across linkage groups. We investigated genome-wide allele frequency in 1144 red clover individuals using 462 microsatellite loci randomly chosen from the consensus map. The average number of alleles and polymorphism information content (PIC) were 9.17 and 0.69, respectively. CONCLUSION: A consensus genetic linkage map for red clover was constructed for the first time based on six mapping populations. The locus order on the consensus map was highly conserved among linkage maps and was sufficiently reliable for use as a reference for genetic analysis of random red clover germplasms.

Survey of genome sequences in a wild sweet potato, Ipomoea trifida (H. B. K.) G. Don.

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.

Structural analysis of four large plasmids harboring in a unicellular cyanobacterium, Synechocystis sp. PCC 6803.

The genome of the unicellular cyanobacterium Synechocystis sp. PCC 6803 consists of a single chromosome and several plasmids of different sizes, and the nucleotide sequences of the chromosome and three small plasmids (5.2 kb, 2.4 kb, and 2.3 kb) have already been sequenced. We newly determined the nucleotide sequences of four large plasmids, which have been identified in our laboratory (pSYSM:120 kb, pSYSX:106 kb, pSYSA:103 kb, and pSYSG:44 kb). Computer-aided analysis was performed to explore the genetic information carried by these plasmids. A total of 397 potential protein-encoding genes were predicted, but little information was obtained about the functional relationship of plasmids to host cell, as a large portion of the predicted genes (77%) were of unknown function. The occurrence of the potential genes on plasmids was divergent, and parA was the only gene common to all four large plasmids. The distribution data of a Cyanobacterium-specific sequence (HIP1: 5'-GCGATCGC-3') suggested that respective plasmids could have originated from different cyanobacterial strains.

Complete genome structure of the thermophilic cyanobacterium Thermosynechococcus elongatus BP-1.

The entire genome of a thermophilic unicellular cyanobacterium, Thermosynechococcus elongatus BP-1, was sequenced. The genome consisted of a circular chromosome 2,593,857 bp long, and no plasmid was detected. A total of 2475 potential protein-encoding genes, one set of rRNA genes, 42 tRNA genes representing 42 tRNA species and 4 genes for small structural RNAs were assigned to the chromosome by similarity search and computer prediction. The translated products of 56% of the potential protein-encoding genes showed sequence similarity to experimentally identified and predicted proteins of known function, and the products of 34% of these genes showed sequence similarity to the translated products of hypothetical genes. The remaining 10% lacked significant similarity to genes for predicted proteins in the public DNA databases. Sixty-three percent of the T. elongatus genes showed significant sequence similarity to those of both Synechocystis sp. PCC 6803 and Anabaena sp. PCC 7120, while 22% of the genes were unique to this species, indicating a high degree of divergence of the gene information among cyanobacterial strains. The lack of genes for typical fatty acid desaturases and the presence of more genes for heat-shock proteins in comparison with other mesophilic cyanobacteria may be genomic features of thermophilic strains. A remarkable feature of the genome is the presence of 28 copies of group II introns, 8 of which contained a presumptive gene for maturase/reverse transcriptase. A trace of genome rearrangement mediated by the group II introns was also observed.

Complete genomic sequence of nitrogen-fixing symbiotic bacterium Bradyrhizobium japonicum USDA110.

The complete nucleotide sequence of the genome of a symbiotic bacterium Bradyrhizobium japonicum USDA110 was determined. The genome of B. japonicum was a single circular chromosome 9,105,828 bp in length with an average GC content of 64.1%. No plasmid was detected. The chromosome comprises 8317 potential protein-coding genes, one set of rRNA genes and 50 tRNA genes. Fifty-two percent of the potential protein genes showed sequence similarity to genes of known function and 30% to hypothetical genes. The remaining 18% had no apparent similarity to reported genes. Thirty-four percent of the B. japonicum genes showed significant sequence similarity to those of both Mesorhizobium loti and Sinorhizobium meliloti, while 23% were unique to this species. A presumptive symbiosis island 681 kb in length, which includes a 410-kb symbiotic region previously reported by Göttfert et al., was identified. Six hundred fifty-five putative protein-coding genes were assigned in this region, and the functions of 301 genes, including those related to symbiotic nitrogen fixation and DNA transmission, were deduced. A total of 167 genes for transposases/104 copies of insertion sequences were identified in the genome. It was remarkable that 100 out of 167 transposase genes are located in the presumptive symbiotic island. DNA segments of 4 to 97 kb inserted into tRNA genes were found at 14 locations in the genome, which generates partial duplication of the target tRNA genes. These observations suggest plasticity of the B. japonicum genome, which is probably due to complex genome rearrangements such as horizontal transfer and insertion of various DNA elements, and to homologous recombination.

Complete genome structure of Gloeobacter violaceus PCC 7421, a cyanobacterium that lacks thylakoids.

The nucleotide sequence of the entire genome of a cyanobacterium Gloeobacter violaceus PCC 7421 was determined. The genome of G. violaceus was a single circular chromosome 4,659,019 bp long with an average GC content of 62%. No plasmid was detected. The chromosome comprises 4430 potential protein-encoding genes, one set of rRNA genes, 45 tRNA genes representing 44 tRNA species and genes for tmRNA, B subunit of RNase P, SRP RNA and 6Sa RNA. Forty-one percent of the potential protein-encoding genes showed sequence similarity to genes of known function, 37% to hypothetical genes, and the remaining 22% had no apparent similarity to reported genes. Comparison of the assigned gene components with those of other cyanobacteria has unveiled distinctive features of the G. violaceus genome. Genes for PsaI, PsaJ, PsaK, and PsaX for Photosystem I and PsbY, PsbZ and Psb27 for Photosystem II were missing, and those for PsaF, PsbO, PsbU, and PsbV were poorly conserved. cpcG for a rod core linker peptide for phycobilisomes and nblA related to the degradation of phycobilisomes were also missing. Potential signal peptides of the presumptive products of petJ and petE for soluble electron transfer catalysts were less conserved than the remaining portions. These observations may be related to the fact that photosynthesis in G. violaceus takes place not in thylakoid membranes but in the cytoplasmic membrane. A large number of genes for sigma factors and transcription factors in the LuxR, LysR, PadR, TetR, and MarR families could be identified, while those for major elements for circadian clock, kaiABC were not found. These differences may reflect the phylogenetic distance between G. violaceus and other cyanobacteria.

Dissection of the octoploid strawberry genome by deep sequencing of the genomes of Fragaria species.

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.

Construction of an integrated high density simple sequence repeat linkage map in cultivated strawberry (Fragaria × ananassa) and its applicability.

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.

Comparative Genetic Mapping and Discovery of Linkage Disequilibrium Across Linkage Groups in White Clover (Trifolium repens L.).

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.