Gene expression profiles associated with early relapse during first remission induction in canine multicentric high-grade B-cell lymphoma

Although chemotherapy using CHOP-based protocol induces remission in most cases of canine multicentric high-grade B-cell lymphoma (mhBCL), some cases develop early relapse during the first induction protocol. In this study, we examined the gene expression profiles of canine mhBCL before chemotherapy and investigated their associations with early relapse during the first whole CHOP-based protocol. Twenty-five cases of mhBCL treated with CHOP-based protocol as first induction chemotherapy were included in this study. Sixteen cases completed the first whole CHOP-based protocol without relapse (S-group), and nine developed relapse during the chemotherapy (R-group). RNA-seq was performed on samples from neoplastic lymph nodes. Differentially expressed genes (DEGs) were extracted by the comparison of gene expression profiles between S- and R-groups, and the differences in the expression levels of these genes were validated by RT-qPCR. Extracted 179 DEGs included the genes related to chemokine CC motif ligand, T-cell receptor signaling pathway, and PD-L1 expression and PD-1 checkpoint pathway. We focused on chemokine CC motif ligand, and CCL4 was confirmed to be significantly downregulated in the R-group (P=0.039). We also focused on the genes related to T-cell signaling pathway, and CD3E (P=0.039), ITK (P=0.023), and LAT (P=0.023) genes were confirmed to be significantly upregulated in the R-group. The current results suggest that both changes in tumor cells and the interactions between tumor cells and immune cells are associated with the efficacy of the chemotherapy for first remission induction.

A practical assembly guideline for genomes with various levels of heterozygosity

Although current long-read sequencing technologies have a long-read length that facilitates assembly for genome reconstruction, they have high sequence errors. While various assemblers with different perspectives have been developed, no systematic evaluation of assemblers with long reads for diploid genomes with varying heterozygosity has been performed. Here, we evaluated a series of processes, including the estimation of genome characteristics such as genome size and heterozygosity, de novo assembly, polishing, and removal of allelic contigs, using six genomes with various heterozygosity levels. We evaluated five long-read-only assemblers (Canu, Flye, miniasm, NextDenovo and Redbean) and five hybrid assemblers that combine short and long reads (HASLR, MaSuRCA, Platanus-allee, SPAdes and WENGAN) and proposed a concrete guideline for the construction of haplotype representation according to the degree of heterozygosity, followed by polishing and purging haplotigs, using stable and high-performance assemblers: Redbean, Flye and MaSuRCA.

Best Practices for Comprehensive Annotation of Neuropeptides of Gryllus bimaculatus

Genome annotation is critically important data that can support research. Draft genome annotations cover representative genes; however, they often do not include genes that are expressed only in limited tissues and stages, or genes with low expression levels. Neuropeptides are responsible for regulation of various physiological and biological processes. A recent study disclosed the genome draft of the two-spotted cricket Gryllus bimaculatus, which was utilized to understand the intriguing physiology and biology of crickets. Thus far, only two of the nine reported neuropeptides in G. bimaculatus were annotated in the draft genome. Even though de novo assembly using transcriptomic analyses can comprehensively identify neuropeptides, this method does not follow those annotations on the genome locus. In this study, we performed the annotations based on the reference mapping, de novo transcriptome assembly, and manual curation. Consequently, we identified 41 neuropeptides out of 43 neuropeptides, which were reported in the insects. Further, 32 of the identified neuropeptides on the genomic loci in G. bimaculatus were annotated. The present annotation methods can be applicable for the neuropeptide annotation of other insects. Furthermore, the methods will help to generate useful infrastructures for studies relevant to neuropeptides.

MarpolBase Expression: A Web-Based, Comprehensive Platform for Visualization and Analysis of Transcriptomes in the Liverwort Marchantia polymorpha

The liverwort Marchantia polymorpha is equipped with a wide range of molecular and genetic tools and resources that have led to its wide use to explore the evo-devo aspects of land plants. Although its diverse transcriptome data are rapidly accumulating, there is no extensive yet user-friendly tool to exploit such a compilation of data and to summarize results with the latest annotations. Here, we have developed a web-based suite of tools, MarpolBase Expression (MBEX, https://marchantia.info/mbex/), where users can visualize gene expression profiles, identify differentially expressed genes, perform co-expression and functional enrichment analyses and summarize their comprehensive output in various portable formats. Using oil body biogenesis as an example, we demonstrated that the results generated by MBEX were consistent with the published experimental evidence and also revealed a novel transcriptional network in this process. MBEX should facilitate the exploration and discovery of the genetic and functional networks behind various biological processes in M. polymorpha and promote our understanding of the evolution of land plants.

A highly contiguous genome assembly of red perilla (Perilla frutescens) domesticated in Japan

Perilla frutescens (Lamiaceae) is an important herbal plant with hundreds of bioactive chemicals, among which perillaldehyde and rosmarinic acid are the two major bioactive compounds in the plant. The leaves of red perilla are used as traditional Kampo medicine or food ingredients. However, the medicinal and nutritional uses of this plant could be improved by enhancing the production of valuable metabolites through the manipulation of key enzymes or regulatory genes using genome editing technology. Here, we generated a high-quality genome assembly of red perilla domesticated in Japan. A near-complete chromosome-level assembly of P. frutescens was generated contigs with N50 of 41.5 Mb from PacBio HiFi reads. 99.2% of the assembly was anchored into 20 pseudochromosomes, among which seven pseudochromosomes consisted of one contig, while the rest consisted of less than six contigs. Gene annotation and prediction of the sequences successfully predicted 86,258 gene models, including 76,825 protein-coding genes. Further analysis showed that potential targets of genome editing for the engineering of anthocyanin pathways in P. frutescens are located on the late-stage pathways. Overall, our genome assembly could serve as a valuable reference for selecting target genes for genome editing of P. frutescens.

Genome evolution of a nonparasitic secondary heterotroph, the diatom Nitzschia putrida

Secondary loss of photosynthesis is observed across almost all plastid-bearing branches of the eukaryotic tree of life. However, genome-based insights into the transition from a phototroph into a secondary heterotroph have so far only been revealed for parasitic species. Free-living organisms can yield unique insights into the evolutionary consequence of the loss of photosynthesis, as the parasitic lifestyle requires specific adaptations to host environments. Here, we report on the diploid genome of the free-living diatom Nitzschia putrida (35 Mbp), a nonphotosynthetic osmotroph whose photosynthetic relatives contribute ca. 40% of net oceanic primary production. Comparative analyses with photosynthetic diatoms and heterotrophic algae with parasitic lifestyle revealed that a combination of gene loss, the accumulation of genes involved in organic carbon degradation, a unique secretome, and the rapid divergence of conserved gene families involved in cell wall and extracellular metabolism appear to have facilitated the lifestyle of a free-living secondary heterotroph.

Identification of the sex-determining factor in the liverwort Marchantia polymorpha reveals unique evolution of sex chromosomes in a haploid system

Sex determination is a central process for sexual reproduction and is often regulated by a sex determinant encoded on a sex chromosome. Rules that govern the evolution of sex chromosomes via specialization and degeneration following the evolution of a sex determinant have been well studied in diploid organisms. However, distinct predictions apply to sex chromosomes in organisms where sex is determined in the haploid phase of the life cycle: both sex chromosomes, female U and male V, are expected to maintain their gene functions, even though both are non-recombining. This is in contrast to the X-Y (or Z-W) asymmetry and Y (W) chromosome degeneration in XY (ZW) systems of diploids. Here, we provide evidence that sex chromosomes diverged early during the evolution of haploid liverworts and identify the sex determinant on the Marchantia polymorpha U chromosome. This gene, Feminizer, encodes a member of the plant-specific BASIC PENTACYSTEINE transcription factor family. It triggers female differentiation via regulation of the autosomal sex-determining locus of FEMALE GAMETOPHYTE MYB and SUPPRESSOR OF FEMINIZATION. Phylogenetic analyses of Feminizer and other sex chromosome genes indicate dimorphic sex chromosomes had already been established 430 mya in the ancestral liverwort. Feminizer also plays a role in reproductive induction that is shared with its gametolog on the V chromosome, suggesting an ancestral function, distinct from sex determination, was retained by the gametologs. This implies ancestral functions can be preserved after the acquisition of a sex determination mechanism during the evolution of a dominant haploid sex chromosome system.

Determination of Lycopene Concentration in Fresh Tomatoes by Spectrophotometry: A Collaborative Study

BACKGROUND

Lycopene has been the object of considerable research attention recently, and the effects of the intake of lycopene, or of tomato products, have been studied in various ways. In Japan, interest in the health-promoting function of food components has increased.

OBJECTIVE

Developing a method to determine lycopene contents in tomato that meets the Japanese Agricultural Standard (JAS).

METHOD

In the proposed JAS method, the test sample consists of fresh tomatoes; a hexane-acetone mixture is utilized as the extraction solvent. A collaborative study was conducted to evaluate the interlaboratory performance of the method.

RESULTS

Ten laboratories participated and analyzed six test materials characterized by a lycopene content between 39 and 170 mg/kg as blind duplicates. After removing statistical outliers, RSDr ranged from 1.2 to 3.0% and RSDR ranged from 2.4 to 4.2%. The HorRat values were calculated and found to be in the 0.26-0.49 range.

CONCLUSIONS

The method for determining the lycopene content in tomato was evaluated by means of a collaborative study, and the reproducibility of this method was found to be acceptable.

HIGHLIGHTS

Intended for standardization in Japan, a method to determine lycopene content in tomato has been developed and shown to have acceptable precision in a collaborative study.

Over-expression of Caj1, a plasma membrane associated J-domain protein in Saccharomyces cerevisiae, stabilizes amino acid permeases

Hsp70: J-domain protein (JDP) machines, along with the cellular protein degradation systems play a central role in regulating cellular proteostasis. An equally robust surveillance system operates at the plasma membrane too that affects proper sorting, stability as well as the turnover of membrane proteins. Although plausible, a definitive role of the Hsp70: JDP machine in regulating the stability of plasma membrane proteins is not well understood in Saccharomyces cerevisiae. Here we show that a moderate over-expression of Caj1, one of the thirteen JDPs residing in the nucleo-cytosolic compartment of S. cerevisiae reduced the cold sensitivity of tryptophan auxotrophic yeast cells by stabilizing tryptophan permeases, Tat1 and Tat2 in a J-domain dependent manner. Concomitantly, higher Caj1 levels also caused slow growth and increased plasma membrane damage at elevated temperatures possibly due to the stabilization of thermolabile plasma membrane proteins. Finally, we show that although majorly cytosolic, Caj1 also co-localizes with the membrane dye FM4-64 at the cellular periphery suggesting that Caj1 might interact with the plasma membrane. Based on the results presented in this study, we implicate the Hsp70: Caj1 chaperone machine in regulating the stability or turnover of plasma membrane proteins in budding yeast.

North Atlantic climate far more predictable than models imply

Quantifying signals and uncertainties in climate models is essential for the detection, attribution, prediction and projection of climate change^1-3. Although inter-model agreement is high for large-scale temperature signals, dynamical changes in atmospheric circulation are very uncertain⁴. This leads to low confidence in regional projections, especially for precipitation, over the coming decades^5,6. The chaotic nature of the climate system^7-9 may also mean that signal uncertainties are largely irreducible. However, climate projections are difficult to verify until further observations become available. Here we assess retrospective climate model predictions of the past six decades and show that decadal variations in North Atlantic winter climate are highly predictable, despite a lack of agreement between individual model simulations and the poor predictive ability of raw model outputs. Crucially, current models underestimate the predictable signal (the predictable fraction of the total variability) of the North Atlantic Oscillation (the leading mode of variability in North Atlantic atmospheric circulation) by an order of magnitude. Consequently, compared to perfect models, 100 times as many ensemble members are needed in current models to extract this signal, and its effects on the climate are underestimated relative to other factors. To address these limitations, we implement a two-stage post-processing technique. We first adjust the variance of the ensemble-mean North Atlantic Oscillation forecast to match the observed variance of the predictable signal. We then select and use only the ensemble members with a North Atlantic Oscillation sufficiently close to the variance-adjusted ensemble-mean forecast North Atlantic Oscillation. This approach greatly improves decadal predictions of winter climate for Europe and eastern North America. Predictions of Atlantic multidecadal variability are also improved, suggesting that the North Atlantic Oscillation is not driven solely by Atlantic multidecadal variability. Our results highlight the need to understand why the signal-to-noise ratio is too small in current climate models¹⁰, and the extent to which correcting this model error would reduce uncertainties in regional climate change projections on timescales beyond a decade.

Chromatin Organization in Early Land Plants Reveals an Ancestral Association between H3K27me3, Transposons, and Constitutive Heterochromatin

Genome packaging by nucleosomes is a hallmark of eukaryotes. Histones and the pathways that deposit, remove, and read histone modifications are deeply conserved. Yet, we lack information regarding chromatin landscapes in extant representatives of ancestors of the main groups of eukaryotes, and our knowledge of the evolution of chromatin-related processes is limited. We used the bryophyte Marchantia polymorpha, which diverged from vascular plants circa 400 mya, to obtain a whole chromosome genome assembly and explore the chromatin landscape and three-dimensional genome organization in an early diverging land plant lineage. Based on genomic profiles of ten chromatin marks, we conclude that the relationship between active marks and gene expression is conserved across land plants. In contrast, we observed distinctive features of transposons and other repetitive sequences in Marchantia compared with flowering plants. Silenced transposons and repeats did not accumulate around centromeres. Although a large fraction of constitutive heterochromatin was marked by H3K9 methylation as in flowering plants, a significant proportion of transposons were marked by H3K27me3, which is otherwise dedicated to the transcriptional repression of protein-coding genes in flowering plants. Chromatin compartmentalization analyses of Hi-C data revealed that repressed B compartments were densely decorated with H3K27me3 but not H3K9 or DNA methylation as reported in flowering plants. We conclude that, in early plants, H3K27me3 played an essential role in heterochromatin function, suggesting an ancestral role of this mark in transposon silencing.

Outcome of maintenance systemic chemotherapy with drug-free interval for metastatic urothelial carcinoma

OBJECTIVE

Aiming to achieve long-term disease control, maintenance systemic chemotherapy (MSC) with a 1-3-month drug-free interval is continued in selected patients. We report our experience of MSC for metastatic urothelial carcinoma (UC).

METHODS

Of 228 metastatic UC patients treated with systemic chemotherapy, 40 (17.5%, 40/228) had continuously undergone MSC. Data on the regimen, cycle number, and reason for the discontinuation of MSC were also collected. We analyzed OS from the initiation of MSC until death or the last follow-up, using the log-rank test to assess the significance of differences.

RESULTS

The median number of cycles of chemotherapy was 6, and the responses were CR in 6, PR in 20, SD in 13, and PD in 1 before MSC. Gemcitabine plus CDDP or carboplatin was mainly performed as MSC (70%, 28/40). MSC was repeated quarterly in 30 (75%, 30/40), every two months in 8 (20%, 8/40), and with other intervals in 2 (5%, 2/40). Overall, a median of 3.5 cycles (range: 1-29) of MSC was performed. The reason for the discontinuation of MSC was PD in 24 (60%, 24/40), favorable disease control in 9 (22.5%, 9/40), and myelosuppression in 3 (7.5%, 3/40), and for other reasons in 2 (5%, 2/40). MSC was ongoing in 2 (5%, 2/40). The median OS was 27 months from the initiation of MSC. PS0 (P = 0.0169), the absence of lung metastasis (P = 0.0387), and resection of the primary site (P = 0.0495) were associated with long-term survival after MSC.

CONCLUSIONS

In selected patients, long-term systemic chemotherapy could be performed with a drug-free interval. Our maintenance strategy with cytotoxic drugs may become one of the treatment options for long-term disease control.

Utility of synthetic MRI in predicting the Ki-67 status of oestrogen receptor-positive breast cancer: a feasibility study

AIM

To evaluate the utility of synthetic magnetic resonance imaging (MRI) of the breast in predicting the Ki-67 status in patients with oestrogen receptor (ER)-positive breast cancer.

MATERIALS AND METHODS

Forty-nine patients with 50 histopathologically proven breast cancers who underwent additional synthetic MRI were enrolled in the present study. Using synthetic MRI images, T1 and T2 relaxation times and their standard deviations (SD) in the breast lesions before (T1-Pre, T2-Pre, PD-Pre, SD of T1-Pre, SD of T2-Pre, SD of PD-Pre) and after (T1-Gd, T2-Gd, PD-Gd, SD of T1-Gd, SD of T2-Gd, SD of PD-Gd) contrast agent injection were obtained. These quantitative values were compared between the low Ki-67 expression (<14%) lesions (low-proliferation group: n=23) and high Ki-67 expression (≥14%) lesions (high-proliferation group: n=27).

RESULTS

The univariate analysis showed that the SD of T1-Gd (p<0.001) and T2-Gd (p=0.042) were significantly higher in the high-proliferation group than in the low-proliferation group. Multivariate analysis further showed that the SD of T1-Gd was a significant and independent predictor of Ki-67 expression, with an area under the receiver operating characteristic (AUROC) curve of 0.885. The sensitivity, specificity, and accuracy of the SD of T1-Gd with an optimal cut-off value of 98.5 were 77.8%, 87%, and 82%, respectively.

CONCLUSION

The SD of T1-Gd obtained from synthetic MRI was useful to predict Ki-67 status.

Long-term Survival in a Kidney Transplantation Patient With Progressive Multifocal Leukoencephalopathy: A Case Report

Post-kidney transplantation progressive multifocal leukoencephalopathy (PML) is a rare disease on which there are very few published reports on record. PML is a demyelinating disease caused by a destructive infection of the oligodendrocytes by the JC polyomavirus. No effective therapeutic protocol has been established other than measures to revive the immune function by reducing or discontinuing the administration of immunosuppressive agents. Most cases are progressive and show a poor prognosis. We herein report a case in which renal function has been maintained for 2 years following the onset of PML, which was initially diagnosed 3 years after kidney transplantation.

Itching in a trichophytin contact dermatitis mouse model and the antipruritic effect of antifungal agents

BACKGROUND

Tinea is an infectious disease by dermatophytes, of which Trichophyton species accounts for the overwhelming majority of case. Tinea often causes itching with inflammation. In terms of pruritus by fungal infection, however, tinea has not been investigated sufficiently to date.

AIM

To evaluate itch caused by Trichophyton infection and the effect of antifungal agents on the infection, by measuring scratch behaviour and profiles of inflammatory cytokines and chemokines.

METHODS

We used a previously established mouse model of contact hypersensitivity induced by trichophytin, a crude extract from Trichophyton mentagrophytes. Scratching behaviour was recorded using a counting device that measured an electric current induced in a coil by movement of magnets that had been inserted into the hind paws of each animal. We investigated expression of various genes in lesional skin of mice and in normal human epidermal keratinocytes. We also investigated the antipruritic effects of the corticosteroid dexamethasone (DEX) and three antifungal agents: ketoconazole (KCZ), terbinafine (TBF) and liranaftate (LNF).

RESULTS

Biphasic peaks of scratching were observed at 1 h and at 6-7 h during an observation period of 14 h after trichophytin induction. For lesional skin, RNA was extracted 24 h after trichophytin challenge, and increased expression was seen in the genes for interleukin (IL)-17A, interferon-γ, tumour necrosis factor (TNF)-α, macrophage inflammatory protein (MIP)-2 and Dectin-1, whereas there was no obvious change in the genes for IL-31 and prostaglandin (PG)E2. Furthermore, KCZ inhibited histidine decarboxylase (HDC) expression in vitro and in vivo, and inhibited scratching in the very early phase. LNF inhibited expression of thymic stromal lymphopoietin (TSLP) and IL-8 in vitro, and TSLP, TNF-α, IL-1α and MIP2 in vivo, and also scratching in the early phase. TBF did not induce any significant alterations in either gene expression or scratching. DEX suppressed expression of all the chemical mediators except HDC in vitro and in vivo, and inhibited scratching.

CONCLUSION

Antifungals can inhibit itching induced by fungal infection through different mechanisms.

Draft Sequencing of the Heterozygous Diploid Genome of Satsuma (Citrus unshiu Marc.) Using a Hybrid Assembly Approach

Satsuma (Citrus unshiu Marc.) is one of the most abundantly produced mandarin varieties of citrus, known for its seedless fruit production and as a breeding parent of citrus. De novo assembly of the heterozygous diploid genome of Satsuma ("Miyagawa Wase") was conducted by a hybrid assembly approach using short-read sequences, three mate-pair libraries, and a long-read sequence of PacBio by the PLATANUS assembler. The assembled sequence, with a total size of 359.7 Mb at the N₅₀ length of 386,404 bp, consisted of 20,876 scaffolds. Pseudomolecules of Satsuma constructed by aligning the scaffolds to three genetic maps showed genome-wide synteny to the genomes of Clementine, pummelo, and sweet orange. Gene prediction by modeling with MAKER-P proposed 29,024 genes and 37,970 mRNA; additionally, gene prediction analysis found candidates for novel genes in several biosynthesis pathways for gibberellin and violaxanthin catabolism. BUSCO scores for the assembled scaffold and predicted transcripts, and another analysis by BAC end sequence mapping indicated the assembled genome consistency was close to those of the haploid Clementine, pummel, and sweet orange genomes. The number of repeat elements and long terminal repeat retrotransposon were comparable to those of the seven citrus genomes; this suggested no significant failure in the assembly at the repeat region. A resequencing application using the assembled sequence confirmed that both kunenbo-A and Satsuma are offsprings of Kishu, and Satsuma is a back-crossed offspring of Kishu. These results illustrated the performance of the hybrid assembly approach and its ability to construct an accurate heterozygous diploid genome.

FUT2 non-secretor status is associated with Type 1 diabetes susceptibility in Japanese children

AIM

To examine the contribution of the FUT2 gene and ABO blood type to the development of Type 1 diabetes in Japanese children.

METHODS

We analysed FUT2 variants and ABO genotypes in a total of 531 Japanese children diagnosed with Type 1 diabetes and 448 control subjects. The possible association of FUT2 variants and ABO genotypes with the onset of Type 1 diabetes was statistically examined.

RESULTS

The se2 genotype (c.385A>T) of the FUT2 gene was found to confer susceptibility to Type 1A diabetes in a recessive effects model [odds ratio for se2/se2, 1.68 (95% CI 1.20-2.35); corrected P value = 0.0075].

CONCLUSIONS

The FUT2 gene contributed to the development of Type 1 diabetes in the present cohort of Japanese children.

DNApod: DNA polymorphism annotation database from next-generation sequence read archives

With the rapid advances in next-generation sequencing (NGS), datasets for DNA polymorphisms among various species and strains have been produced, stored, and distributed. However, reliability varies among these datasets because the experimental and analytical conditions used differ among assays. Furthermore, such datasets have been frequently distributed from the websites of individual sequencing projects. It is desirable to integrate DNA polymorphism data into one database featuring uniform quality control that is distributed from a single platform at a single place. DNA polymorphism annotation database (DNApod; http://tga.nig.ac.jp/dnapod/) is an integrated database that stores genome-wide DNA polymorphism datasets acquired under uniform analytical conditions, and this includes uniformity in the quality of the raw data, the reference genome version, and evaluation algorithms. DNApod genotypic data are re-analyzed whole-genome shotgun datasets extracted from sequence read archives, and DNApod distributes genome-wide DNA polymorphism datasets and known-gene annotations for each DNA polymorphism. This new database was developed for storing genome-wide DNA polymorphism datasets of plants, with crops being the first priority. Here, we describe our analyzed data for 679, 404, and 66 strains of rice, maize, and sorghum, respectively. The analytical methods are available as a DNApod workflow in an NGS annotation system of the DNA Data Bank of Japan and a virtual machine image. Furthermore, DNApod provides tables of links of identifiers between DNApod genotypic data and public phenotypic data. To advance the sharing of organism knowledge, DNApod offers basic and ubiquitous functions for multiple alignment and phylogenetic tree construction by using orthologous gene information.

Subcycle Optical Response Caused by a Terahertz Dressed State with Phase-Locked Wave Functions

The coherent interaction of light with matter imprints the phase information of the light field on the wave function of the photon-dressed electronic state. A driving electric field, together with a stable phase that is associated with the optical probe pulses, enables the role of the dressed state in the optical response to be investigated. We observed optical absorption strengths modulated on a subcycle time scale in a GaAs quantum well in the presence of a multicycle terahertz driving pulse using a near-infrared probe pulse. The measurements were in good agreement with the analytical formula that accounts for the optical susceptibilities caused by the dressed state of the excitons, which indicates that the output probe intensity was coherently reshaped by the excitonic sideband emissions.

Variants associated with autoimmune Type 1 diabetes in Japanese children: implications for age-specific effects of cis-regulatory haplotypes at 17q12-q21

AIMS

The aim of this study was to clarify the significance of previously reported susceptibility variants in the development of autoimmune Type 1 diabetes in non-white children. Tested variants included rs2290400, which has been linked to Type 1 diabetes only in one study on white people. Haplotypes at 17q12-q21 encompassing rs2290400 are known to determine the susceptibility of early-onset asthma by affecting the expression of flanking genes.

METHODS

We genotyped 63 variants in 428 Japanese people with childhood-onset autoimmune Type 1 diabetes and 457 individuals without diabetes. Possible association between variants and age at diabetes onset was examined using age-specific quantitative trait locus analysis and ordered-subset regression analysis.

RESULTS

Ten variants, including rs2290400 in GSDMB, were more frequent among the people with Type 1 diabetes than those without diabetes. Of these, rs689 in INS and rs231775 in CTLA4 yielded particularly high odds ratios of 5.58 (corrected P value 0.001; 95% CI 2.15-14.47) and 1.64 (corrected P value 5.3 × 10^-5 ; 95% CI 1.34-2.01), respectively. Age-specific effects on diabetes susceptibility were suggested for rs2290400; heterozygosity of the risk alleles was associated with relatively early onset of diabetes, and the allele was linked to the phenotype exclusively in the subgroup of age at onset ≤ 5.0 years.

CONCLUSIONS

The results indicate that rs2290400 in GSDMB and polymorphisms in INS and CTLA4 are associated with the risk of Type 1 diabetes in Japanese children. Importantly, cis-regulatory haplotypes at 17q12-q21 encompassing rs2290400 probably determine the risk of autoimmune Type 1 diabetes predominantly in early childhood.

Hybrid Origins of Citrus Varieties Inferred from DNA Marker Analysis of Nuclear and Organelle Genomes

Most indigenous citrus varieties are assumed to be natural hybrids, but their parentage has so far been determined in only a few cases because of their wide genetic diversity and the low transferability of DNA markers. Here we infer the parentage of indigenous citrus varieties using simple sequence repeat and indel markers developed from various citrus genome sequence resources. Parentage tests with 122 known hybrids using the selected DNA markers certify their transferability among those hybrids. Identity tests confirm that most variant strains are selected mutants, but we find four types of kunenbo (Citrus nobilis) and three types of tachibana (Citrus tachibana) for which we suggest different origins. Structure analysis with DNA markers that are in Hardy-Weinberg equilibrium deduce three basic taxa coinciding with the current understanding of citrus ancestors. Genotyping analysis of 101 indigenous citrus varieties with 123 selected DNA markers infers the parentages of 22 indigenous citrus varieties including Satsuma, Temple, and iyo, and single parents of 45 indigenous citrus varieties, including kunenbo, C. ichangensis, and Ichang lemon by allele-sharing and parentage tests. Genotyping analysis of chloroplast and mitochondrial genomes using 11 DNA markers classifies their cytoplasmic genotypes into 18 categories and deduces the combination of seed and pollen parents. Likelihood ratio analysis verifies the inferred parentages with significant scores. The reconstructed genealogy identifies 12 types of varieties consisting of Kishu, kunenbo, yuzu, koji, sour orange, dancy, kobeni mikan, sweet orange, tachibana, Cleopatra, willowleaf mandarin, and pummelo, which have played pivotal roles in the occurrence of these indigenous varieties. The inferred parentage of the indigenous varieties confirms their hybrid origins, as found by recent studies.

Effective Treatment With Everolimus for Recurrent Granulomatous Interstitial Nephritis in a Renal Transplant Recipient: A Case Report

BACKGROUND

Granulomatous interstitial nephritis (GIN) is a rare renal disease, and its etiology remains unknown. We report recurrent GIN in renal allograft successfully treated with everolimus (EVR).

CASE REPORT

A 22-year-old man with GIN received a kidney from his mother. On follow-up 8 months later, his serum creatinine level was increased, from 1.3 mg/dL to 1.7 mg/dL, and he had microhematuria and proteinuria. A protocol graft biopsy at 1 year after transplantation showed epithelioid granuloma with multinucleated giant cells. He received steroid pulse therapy for recurrent GIN twice, but he developed allograft dysfunction, hematuria, and proteinuria. EVR was started in combination with maintenance immunosuppressants at 28 months after transplantation. Thereafter, the serum creatinine level decreased, from 2.1 mg/dL to 1.6 mg/dL, and microhematuria and proteinuria were stable despite reduction of steroid dose.

CONCLUSIONS

Maintenance immunosuppressive therapy combined with EVR may be effective for the recurrence of idiopathic GIN in renal allograft.

OryzaGenome: Genome Diversity Database of Wild Oryza Species

The species in the genus Oryza, encompassing nine genome types and 23 species, are a rich genetic resource and may have applications in deeper genomic analyses aiming to understand the evolution of plant genomes. With the advancement of next-generation sequencing (NGS) technology, a flood of Oryza species reference genomes and genomic variation information has become available in recent years. This genomic information, combined with the comprehensive phenotypic information that we are accumulating in our Oryzabase, can serve as an excellent genotype-phenotype association resource for analyzing rice functional and structural evolution, and the associated diversity of the Oryza genus. Here we integrate our previous and future phenotypic/habitat information and newly determined genotype information into a united repository, named OryzaGenome, providing the variant information with hyperlinks to Oryzabase. The current version of OryzaGenome includes genotype information of 446 O. rufipogon accessions derived by imputation and of 17 accessions derived by imputation-free deep sequencing. Two variant viewers are implemented: SNP Viewer as a conventional genome browser interface and Variant Table as a text-based browser for precise inspection of each variant one by one. Portable VCF (variant call format) file or tab-delimited file download is also available. Following these SNP (single nucleotide polymorphism) data, reference pseudomolecules/scaffolds/contigs and genome-wide variation information for almost all of the closely and distantly related wild Oryza species from the NIG Wild Rice Collection will be available in future releases. All of the resources can be accessed through http://viewer.shigen.info/oryzagenome/.

Global expression differences and tissue specific expression differences in rice evolution result in two contrasting types of differentially expressed genes

BACKGROUND

Since the development of transcriptome analysis systems, many expression evolution studies characterized evolutionary forces acting on gene expression, without explicit discrimination between global expression differences and tissue specific expression differences. However, different types of gene expression alteration should have different effects on an organism, the evolutionary forces that act on them might be different, and different types of genes might show different types of differential expression between species. To confirm this, we studied differentially expressed (DE) genes among closely related groups that have extensive gene expression atlases, and clarified characteristics of different types of DE genes including the identification of regulating loci for differential expression using expression quantitative loci (eQTL) analysis data.

RESULTS

We detected differentially expressed (DE) genes between rice subspecies in five homologous tissues that were verified using japonica and indica transcriptome atlases in public databases. Using the transcriptome atlases, we classified DE genes into two types, global DE genes and changed-tissues DE genes. Global type DE genes were not expressed in any tissues in the atlas of one subspecies, however changed-tissues type DE genes were expressed in both subspecies with different tissue specificity. For the five tissues in the two japonica-indica combinations, 4.6 ± 0.8 and 5.9 ± 1.5 % of highly expressed genes were global and changed-tissues DE genes, respectively. Changed-tissues DE genes varied in number between tissues, increasing linearly with the abundance of tissue specifically expressed genes in the tissue. Molecular evolution of global DE genes was rapid, unlike that of changed-tissues DE genes. Based on gene ontology, global and changed-tissues DE genes were different, having no common GO terms. Expression differences of most global DE genes were regulated by cis-eQTLs. Expression evolution of changed-tissues DE genes was rapid in tissue specifically expressed genes and those rapidly evolved changed-tissues DE genes were regulated not by cis-eQTLs, but by complicated trans-eQTLs.

CONCLUSIONS

Global DE genes and changed-tissues DE genes had contrasting characteristics. The two contrasting types of DE genes provide possible explanations for the previous controversial conclusions about the relationships between molecular evolution and expression evolution of genes in different species, and the relationship between expression breadth and expression conservation in evolution.