Plateau-linear reaction norm model analysis of number born alive in purebred Landrace pigs using meteorological data in Japan.

We performed a plateau-linear reaction norm model (RNM) analysis of number born alive (NBA) in purebred Landrace pigs, where breeding value changes according to maximum temperature at mating day, using public meteorological observation data in Japan. We analysed 52,668 NBA records obtained from 10,320 Landrace sows. Pedigree data contained 99,201 animals. Off-farm daily temperature data at the nearest weather station from each of the farms were downloaded from the Japan Meteorological Agency website. A plateau-linear RNM analysis based on daily maximum temperature on mating day (threshold temperature of 16.6°C) was performed. The percentage of the records with daily maximum temperatures at mating days of ≤16.6, ≥25.0, ≥30.0 and ≥35.0°C were 34.3%, 33.6%, 14.0% and 0.8%, respectively. The value of Akaike's information criterion for the plateau-linear RNM was lower than that for a simple repeatability model (RM). With the plateau-linear RNM, estimated value of heritability ranged from 0.14 to 0.15, while that from the RM analysis was 0.15. Additive genetic correlation between intercept and slope terms was estimated to be -0.52 from the plateau-linear RNM analysis. Estimated additive genetic correlations were >0.9 between NBA at different temperatures ranging from 16.6 to 37.6°C. For the 10,320 sows, average values of prediction reliability of the intercept and slope terms for breeding values in the plateau-linear RNM were 0.47 and 0.16, respectively. Increasing weight for slope term in linear selection index could bring positive genetic gain in the slope part, but prediction accuracy would decrease. Our results imply that genetically improving heat tolerance in sows reared in Japan focusing on NBA using RNM is possible, while RNM is more complex to implement and interpret. Therefore, further study should be encouraged to make genetic improvement for heat tolerance in sows more efficient.

Reaction-norm animal model analysis of average daily gain heat tolerance in purebred Duroc pigs.

The present study aimed to genetically improve growth performance under high-heat environments by specifically designing a reaction-norm animal model (RNAM) for purebred Duroc pigs in Japan. A total of 54,750 records of average daily gain (ADG) measured for pigs reared at four farms in different prefectures were analyzed. Estimated maximum daily temperatures at the respective farm locations were used to calculate the average cumulative thermal load (TL). The TL values served as an indicator of high-heat environments for pigs. The plausible cumulative period length and threshold temperature for calculating TL were determined to be 8 weeks until just before shipping and 25°C, respectively. Variance components were estimated via RNAM analysis using TL as a linear covariate. The estimated additive genetic variances under both responsive and non-responsive to TL were found to be significant. Moreover, the estimated heritability of ADG ranged from 0.38 to 0.73 for TL values of 0-8. These results suggest that the RNAM developed holds the potential for improving the genetic ability of growth under high-heat environments in pigs.

Application of linear and machine learning models to genomic prediction of fatty acid composition in Japanese Black cattle.

We collected 3180 records of oleic acid (C18:1) and monounsaturated fatty acid (MUFA) measured using gas chromatography (GC) and 6960 records of C18:1 and MUFA measured using near-infrared spectroscopy (NIRS) in intermuscular fat samples of Japanese Black cattle. We compared genomic prediction performance for four linear models (genomic best linear unbiased prediction [GBLUP], kinship-adjusted multiple loci [KAML], BayesC, and BayesLASSO) and five machine learning models (Gaussian kernel [GK], deep kernel [DK], random forest [RF], extreme gradient boost [XGB], and convolutional neural network [CNN]). For GC-based C18:1 and MUFA, KAML showed the highest accuracies, followed by BayesC, XGB, DK, GK, and BayesLASSO, with more than 6% gain of accuracy by KAML over GBLUP. Meanwhile, DK had the highest prediction accuracy for NIRS-based C18:1 and MUFA, but the difference in accuracies between DK and KAML was slight. For all traits, accuracies of RF and CNN were lower than those of GBLUP. The KAML extends GBLUP methods, of which marker effects are weighted, and involves only additive genetic effects; whereas machine learning methods capture non-additive genetic effects. Thus, KAML is the most suitable method for breeding of fatty acid composition in Japanese Black cattle.

Analysis of the trehalose synthesis pathway of Physarum polycehalum.

Desiccation is a severe survival problem for organisms. We have been studying the desiccation tolerance mechanisms in the true slime mold Physarum polycephalum. We measured the trehalose content of P. polycephalum vegetative cells (plasmodia) and drought cells (sclerotia). Surprisingly, we found that the content in sclerotia was about 473-fold greater than in the plasmodia. We then examined trehalose metabolism-related genes via RNAseq, and consequently found that trehalose 6-phosphate phosphorylase (T6pp) expression levels increased following desiccation. Next, we cloned and expressed the genes for T6pp, trehalose 6-phosphate synthase/phosphatase (Tps/Tpp), maltooligosyltrehalose trehalohydrolase (TreZ), and maltooligosyltrehalose synthase (TreY) in E. coli. Incidentally, TreY and TreZ clones have been reported in several prokaryotes, but not in eukaryotes. This report in P. polycephalum is the first evidence of their presence in a eukaryote species. Recombinant T6pp, TreY, and TreZ were purified and confirmed to be active. Our results showed that these enzymes catalyze reactions related to trehalose production, and their reaction kinetics follow the Michaelis-Menten equation. The t6pp mRNA levels of the sclerotia were about 15-fold higher than in the plasmodia. In contrast, the expression levels of TreZ and TreY showed no significant change between the sclerotia and plasmodia. Thus, T6pp is probably related to desiccation tolerance, whereas the contribution of TreY and TreZ is insufficient to account for the considerable accumulation of trehalose in sclerotia.

Annual Increase of Acute Inpatients with Both Cancer and Cardiovascular Diseases in Japan 2011-2015: Analysis From National Database of Health Insurance Claims and Specific Health Checkups of Japan.

BACKGROUND: Patients with cancer were able to live longer due to improvements in cancer treatment. Additionally, cardiovascular disease (CVD) is the second leading cause of mortality in cancer survivors. However, epidemiological data on onco-CVD have not been sufficiently provided. We aimed to investigate the clinical characteristics of cancer in CVD patients using the National Database of Health Insurance Claims and Specific Health Checkups of Japan (NDB). METHOD AND RESULTS: The NDB sampling dataset used in this study was randomly sampled 10% from the whole Diagnosis Procedure Combination (DPC) records from every January, April, July, and October from 2011 to 2015. The significance of the increase trend in the percentage of records in each disease group to the total number of all DPC records from 2011 to 2015 was checked with Chi-square test with a Bonferroni correction. The percentage of records in cancer with the CVD group to the total number of all DPC records significantly increased over time, and their average age also increased since 2011. Their proportion over 75 years was approximately 56 % in 2015. There was no difference in the cancer sites. However, the prevalence of heart failure dramatically elevated. CONCLUSION: We were able to assess the increase in cancer among CVD patients using DPC inpatient records obtained from the NDB. Both cardiologists and oncologists should be more aware of this phenomenon.

Comparing pedigree and genomic inbreeding coefficients, and inbreeding depression of reproductive traits in Japanese Black cattle.

BACKGROUND: Pedigree-based inbreeding coefficients have been generally included in statistical models for genetic evaluation of Japanese Black cattle. The use of genomic data is expected to provide precise assessment of inbreeding level and depression. Recently, many measures have been used for genome-based inbreeding coefficients; however, with no consensus on which is the most appropriate. Therefore, we compared the pedigree- ([Formula: see text]) and multiple genome-based inbreeding coefficients, which were calculated from the genomic relationship matrix with observed allele frequencies ([Formula: see text]), correlation between uniting gametes ([Formula: see text]), the observed vs expected number of homozygous genotypes ([Formula: see text]), runs of homozygosity (ROH) segments ([Formula: see text]) and heterozygosity by descent segments ([Formula: see text]). We quantified inbreeding depression from estimating regression coefficients of inbreeding coefficients on three reproductive traits: age at first calving (AFC), calving difficulty (CD) and gestation length (GL) in Japanese Black cattle. RESULTS: The highest correlations with [Formula: see text] were for [Formula: see text] (0.86) and [Formula: see text] (0.85) whereas [Formula: see text] and [Formula: see text] provided weak correlations with [Formula: see text], with range 0.33-0.55. Except for [Formula: see text] and [Formula: see text], there were strong correlations among genome-based inbreeding coefficients ([Formula: see text] 0.94). The estimates of regression coefficients of inbreeding depression for [Formula: see text] was 2.1 for AFC, 0.63 for CD and -1.21 for GL, respectively, but [Formula: see text] had no significant effects on all traits. Genome-based inbreeding coefficients provided larger effects on all reproductive traits than [Formula: see text]. In particular, for CD, all estimated regression coefficients for genome-based inbreeding coefficients were significant, and for GL, that for [Formula: see text] had a significant.. Although there were no significant effects when using overall genome-level inbreeding coefficients for AFC and GL, [Formula: see text] provided significant effects at chromosomal level in four chromosomes for AFC, three chromosomes for CD, and two chromosomes for GL. In addition, similar results were obtained for [Formula: see text]. CONCLUSIONS: Genome-based inbreeding coefficients can capture more phenotypic variation than [Formula: see text]. In particular, [Formula: see text] and [Formula: see text] can be considered good estimators for quantifying inbreeding level and identifying inbreeding depression at the chromosome level. These findings might improve the quantification of inbreeding and breeding programs using genome-based inbreeding coefficients.

Evaluating the performance of genomic prediction accounting for effects of single nucleotide polymorphism markers in reproductive traits of Japanese Black cattle.

We examined the prediction accuracies of genomic best linear unbiased prediction (GBLUP), various weighted GBLUP according to the degrees of marker effects (WGBLUP) and machine learning (ML) methods, and compared them with traditional BLUP for age at first calving (AFC), calving difficulty (CD), and gestation length in Japanese Black cattle. For WGBLUP, firstly, BayesC and FarmCPU were used to estimate marker effects. Then, we constructed three weighted genomic relationship matrices from information of estimated marker effects in the first step: absolute value of the estimated marker-effect WGBLUP, estimated marker-variance WGBLUP, and genomic-feature WGBLUP. For ML, we applied Gaussian kernel, random forest, extreme gradient boost, and support vector regression. We collected a total of 2583 animals having both phenotypic records and genotypes with 30,105 markers and 16,406 pedigree records. For AFC, prediction accuracies of WGBLUP methods using FarmCPU exceeded BLUP by 25.7%-39.5%. For CD, estimated marker-variance WGBLUP using BayesC achieved the highest prediction accuracy. Among ML methods, extreme gradient boost, support vector regression, and Gaussian kernel increased prediction accuracies by 28.4%, 19.0%, and 36.4% for AFC, CD, and gestation length compared with BLUP, respectively. Thus, prediction performance could be improved using suitable WGBLUP and ML methods according to target reproductive traits for the population used.

Combined actions of bacteriophage-encoded genes in Wolbachia-induced male lethality.

Some Wolbachia endosymbionts induce male killing, whereby male offspring of infected females are killed during development; however, the origin and diversity of the underlying mechanisms remain unclear. In this study, we identified a 76 kbp prophage region specific to male-killing Wolbachia hosted by the moth Homona magnanima. The prophage encoded a homolog of the male-killing gene oscar in Ostrinia moths and the wmk gene that induces various toxicities in Drosophila melanogaster. Upon overexpressing these genes in D. melanogaster, wmk-1 and wmk-3 killed all males and most females, whereas Hm-oscar, wmk-2, and wmk-4 had no impact on insect survival. Strikingly, co-expression of tandemly arrayed wmk-3 and wmk-4 killed 90% of males and restored 70% of females, suggesting their conjugated functions for male-specific lethality. While the male-killing gene in the native host remains unknown, our findings highlight the role of bacteriophages in male-killing evolution and differences in male-killing mechanisms among insects.

Gut microbiota of Parkinson's disease in an appendectomy cohort: a preliminary study.

In patients with Parkinson's disease (PD), α-synuclein pathology is thought to spread to the brain via the dorsal motor nucleus of the vagus nerve. The link between the gut microbiome and PD has been explored in various studies. The appendix might play an important role in immunity by maintaining the microbiota as a reservoir. In recent times, appendectomy has been linked to a lower risk of PD, possibly owing to the role of the appendix in altering the gut microbiome. We aimed to elucidate whether the gut microbiota affects PD development in the appendectomy cohort. We analyzed the fecal microbial composition in patients with PD and healthy controls with and without a history of appendectomy. The abundance of microbes from the family Enterobacteriaceae was higher in feces samples from patients with Parkinson's disease compared to that in samples collected from healthy controls. Furthermore, there was a significant phylogenetic difference between patients with PD and healthy controls who had undergone appendectomy. There was a significant phylogenetic difference between patients with PD and HCs who had undergone APP. These results suggest the correlation between gut microbiota and PD in patients who have undergone APP.

Risk factor structure of heart failure in patients with cancer after treatment with anticancer agents' assessment by big data from a Japanese electronic health record.

As the prognosis of cancer patients has been improved, comorbidity of heart failure (HF) in cancer survivors is a serious concern, especially in the aged population. This study aimed to examine the risk factors of HF development after treatment by anticancer agents, using a machine learning-based analysis of a massive dataset obtained from the electronic health record (EHR) in Japan. This retrospective, cohort study, using a dataset from 2008 to 2017 in the Diagnosis Procedure Combination (DPC) database in Japan, enrolled 140,327 patients. The structure of risk factors was determined using multivariable analysis and classification and regression tree (CART) algorithm for time-to-event data. The mean follow-up period was 1.55 years. The prevalence of HF after anticancer agent administration were 4.0%. HF was more prevalent in the older than the younger. As the presence of cardiovascular diseases and various risk factors predicted HF, CART analysis of the risk factors revealed that the risk factor structures complicatedly differed among different age groups. The highest risk combination was hypertension, diabetes mellitus, and atrial fibrillation in the group aged ≤ 64 years, and the presence of ischemic heart disease was a key in both groups aged 65-74 years and 75 ≤ years. The machine learning-based approach was able to develop complicated HF risk structures in cancer patients after anticancer agents in different age population, of which knowledge would be essential for realizing precision medicine to improve the prognosis of cancer patients.

Whole-Genome Sequencing Revealed a Late-Maturing Isogenic Rice Koshihikari Integrated with Hd16 Gene Derived from an Ise Shrine Mutant.

We identified the key genes controlling the late maturation of the Japonica cultivar Isehikari, which was found at Ise Jingu Shrine and matures 6 days later than Koshihikari. We conducted a genetics-based approach through this study. First, the latest mature plants, which flowered later than Isehikari, were segregated in the F2 and F3 generations of Koshihikari×Isehikari. Next, the linkage relationship of a single late-maturing gene with the SSR markers on the long arm of chromosome 3 was inferred by using late-maturing homozygous F2 segregants. Moreover, genetic analyses of late maturity were conducted through the process of six times of continuous backcross with Koshihikari as a recurrent parent by using the late-maturing homozygous F3 line as a nonrecurrent parent, thus developing a late-maturing isogenic Koshihikari (BC6F2). As a result, we elucidated a single late-maturing gene with incomplete dominance that caused the 14-day maturation delay of Koshihikari. The whole-genome sequencing was conducted on both of Koshihikari and the late-maturing isogenic Koshihikari. Then, the SNP call was conducted as the reference genome of Koshihikari. Finally, a single SNP was identified in the key gene Hd16 of the late-maturing isogenic Koshihikari.

Genetic analysis for sow stayability at different parities in purebred Landrace and Large White pigs.

Genetic parameters for sow stayability were estimated from farrowing records of 10,295 Landrace sows and 8192 Large White sows. The record for sow stayability from parity k to parity k + 1 (k = 1, …, 6) was 0 when a sow had a farrowing record at parity k but not at parity k + 1, and 1 when a sow had both records. Heritability was estimated by using single-trait linear and threshold animal models. Genetic correlations among parities were estimated by using two-trait linear-linear and single-trait random regression linear animal models. Genetic correlations with litter traits at birth were estimated by using a two-trait linear-linear animal model. Heritability estimates by linear model analysis were low (0.065-0.119 in Landrace & 0.061-0.157 in Large White); those by threshold model analysis were higher (0.136-0.200 & 0.110-0.283). Genetic correlations among parities differed between breeds and models. Genetic correlation between sow stayability and number born alive was positive in many cases, implying that selection for number born alive does not reduce sow stayability. The results seem to be affected by decisions on culling made by farmers.

Utilizing Cattle Manure Compost Increases Ammonia Monooxygenase A Gene Expression and Ammonia-oxidizing Activity of Both Bacteria and Archaea in Biofiltration Media for Ammonia Deodorization.

Malodorous emissions are a crucial and inevitable issue during the decomposition of biological waste and contain a high concentration of ammonia. Biofiltration technology is a feasible, low-cost, energy-saving method that reduces and eliminates malodors without environmental impact. In the present study, we evaluated the effectiveness of compost from cattle manure and food waste as deodorizing media based on their removal of ammonia and the expression of ammonia-oxidizing genes, and identified the bacterial and archaeal communities in these media. Ammonia was removed by cattle manure compost, but not by food waste compost. The next-generation sequencing of 16S ribosomal RNA obtained from cattle manure compost revealed the presence of ammonia-oxidizing bacteria (AOB), including Cytophagia, Alphaproteobacteria, and Gammaproteobacteria, and ammonia-oxidizing archaea (AOA), such as Thaumarchaeota. In cattle manure compost, the bacterial and archaeal ammonia monooxygenase A (amoA) genes were both up-regulated after exposure to ammonia (fold ratio of 14.2±11.8 after/before), and the bacterial and archaeal communities were more homologous after than before exposure to ammonia, which indicates the adaptation of these communities to ammonia. These results suggest the potential of cattle manure compost as an efficient biological deodorization medium due to the activation of ammonia-oxidizing microbes, such as AOB and AOA, and the up-regulation of their amoA genes.

Maternal effect on body measurement and meat production traits in purebred Duroc pigs.

We investigated maternal effect on nine body measurement traits (body height, body length, front width (FW), chest width (CW), hind width (HW), chest depth, chest girth (CHG), front cannon circumference (FCC) and rear cannon circumference (RCC)) measured at the end of performance testing and five meat production traits (ages at the start and end of performance testing (D30 and D105), average daily gain (ADG), backfat thickness and loin muscle area) in purebred Duroc pigs. Genetic parameters for each trait were estimated by using six single-trait models with and without common litter environmental effect, maternal genetic effect and direct-maternal genetic correlation. The value of Akaike's information criterion was lowest with the model including direct additive genetic and common litter environmental effects for 10 traits. The estimated proportion of common litter environmental variance to phenotypic variance was approximately ≥0.1 for D30, D105, ADG, FW, CW, HW, CHG, FCC and RCC. Using a model without common litter environmental effect would overestimate the direct heritability of most traits. Standard errors of estimated genetic parameters tended to be larger in models including maternal genetic effect. The results indicate that a compromise could be made for accurate genetic parameter estimation for body measurement traits, as well as meat production traits, in pigs by considering common litter environmental effect.

Genetic relationship between litter size traits at birth and body measurement and production traits in purebred Duroc pigs.

Heritabilities of litter size traits at birth (total number born (TNB), number born alive (NBA), and number still born (NSB)) and their genetic correlations with body measurement (body height, body length, front width (FW), chest width (CW), hind width, chest depth, chest girth, front cannon circumference, and rear cannon circumference) and production traits (ages at the start and end of performance testing (D30 and D105), average daily gain (ADG), backfat thickness, and loin muscle area) in purebred Duroc pigs were estimated. Records of performance testing for 2,835 animals and farrowing records of 1,168 litters from 437 dams were used. Genetic parameters were estimated using single-trait and two-trait animal models. Permanent environment effect was considered for litter size traits and common litter environmental effect was considered for body measurement and production traits. The estimated heritability was 0.10 ± 0.06 for TNB, 0.16 ± 0.06 for NBA, and 0.08 ± 0.05 for NSB. Positive genetic correlation of NBA was estimated with D30, D105, and ADG (0.51, 0.11, and 0.39). The estimated genetic correlation of NBA was 0.47 ± 0.17 with FW and 0.55 ± 0.18 with CW, implying that FW and CW could be promising indicator traits for efficiently improving NBA.

Detection of candidate proteins in the indican biosynthetic pathway of Persicaria tinctoria (Polygonum tinctorium) using protein-protein interactions and transcriptome analyses.

Persicaria tinctoria (Polygonum tinctorium) synthesizes indican (indoxyl-ß-D-glucoside) as a specialized metabolite. Indican is synthesized in the cytosol of leaf cells from indoxyl and UDP-glucose by the catalysis of indoxyl-ß-D-glucoside synthase (PtIGS), then transported into vacuoles. As a portion of PtIGS is found on the microsomal membrane, we assume that it is present on the ER membrane as a large complex involving other indican metabolism-related proteins. Based on this hypothesis, the existence of such a complex was investigated using two separate approaches: a protein-protein interaction assay and transcriptome analysis. We first performed a co-immunoprecipitation using the anti-PtIGS antibody and a pull-down assay using recombinant PtIGS, then identified the candidate proteins through MS/MS analysis. Secondly, we performed a transcriptome analysis to examine the differential gene expression between the first and the second leaves. The expressions of candidate genes detected by protein-protein interaction analyses were collated with transcriptome data and validated by quantitative reverse transcription polymerase chain reaction, showing that the expression of sucrose synthase and cytochrome P450 genes decreased in the second leaves compared with the first leaves. Furthermore, we detected several additional proteins, such as heat shock and cytoskeletal proteins, suggesting that PtIGS may form a large complex, a metabolon.

Differential Metabolomics Profiles Identified by CE-TOFMS between High and Low Intramuscular Fat Amount in Fattening Pigs.

The amount of intramuscular fat (IMF) present in the loin eye area is one of the most important characteristics of high-quality pork. IMF measurements are currently impractical without a labor-intensive process. Metabolomic profiling could be used as an IMF indicator to avoid this process; however, no studies have investigated their use during the fattening period of pigs. This study examined the metabolite profiles in the plasma of two groups of pigs derived from the same Duroc genetic line and fed the same diet. Five plasma samples were collected from each individual the day before slaughter. Capillary electrophoresis-time of flight mass spectrometry (CE-TOFMS) was used to analyze the purified plasma from each sample. Principle component analysis (PCA) and partial least squares (PLS) were used to find the semi-quantitative values of the compounds. The results indicate that branched-chain amino acids are significantly associated with high IMF content, while amino acids are associated with low IMF content. These differences were validated using the quantification analyses by high-performance liquid chromatograph, which supported our results. These results suggest that the concentration of branched-chain amino acids in plasma could be an indicative biomarker for the IMF content in the loin eye area.

Comparison of two models to estimate genetic parameters for number of born alive in pigs.

The performance of the two-trait animal model that regards the first parity and later parities as two different traits in estimating genetic parameters for number of born alive (NBA) was examined using real and simulated data. Genetic parameters for NBA were estimated in purebred Landrace and Large White pigs using a single-trait repeatability model (Model 1) that regards all parities as the same trait and a two-trait animal model (Model 2) that regards the first and the later parities as different traits. For Model 2, the permanent environmental effect was fitted to only the records of the later parities. Heritability for NBA estimated using Model 1 was 0.12 for Landrace and 0.11 for Large White. Estimated heritability for NBA of the first parity and the later parities was 0.21 and 0.16, respectively, for Landrace; 0.18 and 0.16, respectively, for Large White obtained using Model 2, and higher than those in both breeds obtained using Model 1. Further results based on data simulated using the Monte Carlo method suggest that estimated additive genetic variance could be more biased using Model 2 than Model 1.

Inferring phenotypic causal structure among farrowing and weaning traits in pigs.

Direct selection for litter size or weight at weaning in pigs is often hindered by external interventions such as cross-fostering. The objective of this study was to infer the causal structure among phenotypes of reproductive traits in pigs to enable subsequent direct selection for these traits. Examined traits included: number born alive (NBA), litter size on day 21 (LS21), and litter weight on day 21 (LW21). The study included 6,240 litters from 1,673 Landrace dams and 5,393 litters from 1,484 Large White dams. The inductive causation (IC) algorithm was used to infer the causal structure, which was then fitted to a structural equation model (SEM) to estimate causal coefficients and genetic parameters. Based on the IC algorithm and temporal and biological information, the causal structure among traits was identified as: NBA â†’ LS21 â†’ LW21 and NBA â†’ LW21. Owing to the causal effect of NBA on LS21 and LW21, the genetic, permanent environmental, and residual variances of LS21 and LW21were much lower in the SEM than in the multiple-trait model for both breeds. Given the strong effect of NBA on LS21 and LW21, the SEM and causal information might assist with selective breeding for LS21 and LW21 when cross-fostering occurs.

Comprehensive analysis of DNA methylation and gene expression in orally tolerized T cells.

T cell anergy is known to be a crucial mechanism for various types of immune tolerance, including oral tolerance. The expression of several anergy-specific genes was reportedly up-regulated in anergic T cells, and played important roles in the cells. However, how the genes were up-regulated has not been understood. In this study, we comprehensively analyzed the altered gene expression and DNA methylation status in T cells tolerized by oral antigen in vivo. Our results showed that many genes were significantly up-regulated in the orally tolerized T cells, and most of the genes found in this study have not been reported previously as anergy related genes; for example, ribosomal protein L41 (FC = 3.54E06, p = 3.70E-09: Fisher's exact test; the same applies hereinafter) and CD52 (FC = 2.18E05, p = 3.44E-06). Furthermore, we showed that the DNA methylation statuses of many genes; for example, enoyl-coenzyme A delta isomerase 3 (FC = 3.62E-01, p = 3.01E-02) and leucine zipper protein 1 (FC = 4.80E-01, p = 3.25E-02), including the ones distinctly expressed in tolerized T cells; for example, latexin (FC = 3.85E03, p = 4.06E-02 for expression; FC = 7.75E-01, p = 4.13E-01 for DNA methylation) and small nuclear ribonucleoprotein polypeptide F (FC = 3.12E04, p = 4.46E-04 for expression; FC = 8.56E-01, p = 5.15E-01 for DNA methylation), changed during tolerization, suggesting that the distinct expression of some genes was epigenetically regulated in the tolerized T cells. This study would contribute to providing a novel clue to the fine understanding of the mechanism for T cell anergy and oral tolerance.