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.

Large-scale epidemiological study on feline autosomal dominant polycystic kidney disease and identification of novel PKD1 gene variants.

OBJECTIVES: Autosomal dominant polycystic kidney disease (ADPKD) is a common inherited disease in cats. In most cases, the responsible abnormality is a nonsense single nucleotide polymorphism in exon 29 of the PKD1 gene (chrE3:g.42858112C>A, the conventional PKD1 variant). The aim of this study was to conduct a large-scale epidemiological study of ADPKD caused by the conventional PKD1 variant in Japan and to search for novel polymorphisms by targeted resequencing of the PKD1 using a next-generation sequencer. METHODS: A total of 1281 cats visiting the Veterinary Medical Center of the University of Tokyo were included in this study. DNA was extracted from the blood of each cat. We established a novel TaqMan real-time PCR genotyping assay for the conventional PKD1 variant, and all cases were examined for the presence of this variant. Targeted resequencing of all exons of the PKD1 was performed on the DNA of 23 cats with the conventional PKD1 variant, six cats diagnosed with cystic kidneys but without this variant, and 61 wild-type normal cats. RESULTS: Among the 1281 cats examined in this study, 23 (1.8%) harboured the conventional PKD1 variant. The odds of having the conventional PKD1 variant were significantly higher in Persian cats, Scottish Folds and Exotic Shorthairs than in the other breeds, although the number of cases in each breed was small. Furthermore, we identified four variants unique to cats with cystic kidneys that were not found in wild-type normal cats, all of which were in exon 15. In particular, two (chrE:g.42848725delC, pGly1641fs and chrE:g.42850283C>T, pArg2162Trp) were candidate variants. CONCLUSIONS AND RELEVANCE: This study revealed that the conventional PKD1 variant was prevalent in Scottish Fold, Persian and Exotic Shorthair breeds in Japan, and variants in exon 15 of PKD1, in addition to the conventional variant in exon 29, would be key factors in the pathogenesis of ADPKD in cats.

Prevalence of pathogenic variants in cancer-predisposing genes in second cancer after childhood solid cancers.

BACKGROUND: Second malignant neoplasms (SMNs) are one of the most severe late complications after pediatric cancer treatment. However, the effect of genetic variation on SMNs remains unclear. In this study, we revealed germline genetic factors that contribute to the development of SMNs after treatment of pediatric solid tumors. METHODS: We performed whole-exome sequencing in 14 pediatric patients with SMNs, including three brain tumors. RESULTS: Our analysis revealed that five of 14 (35.7%) patients had pathogenic germline variants in cancer-predisposing genes (CPGs), which was significantly higher than in the control cohort (p < 0.01). The identified genes with variants were TP53 (n = 2), DICER1 (n = 1), PMS2 (n = 1), and PTCH1 (n = 1). In terms of the type of subsequent cancer, leukemia and multiple episodes of SMN had an exceptionally high rate of CPG pathogenic variants. None of the patients with germline variants had a family history of SMN development. Mutational signature analysis showed that platinum drugs contributed to the development of SMN in three cases, which suggests the role of platinum agents in SMN development. CONCLUSIONS: We highlight that overlapping effects of genetic background and primary cancer treatment contribute to the development of second cancers after treatment of pediatric solid tumors. A comprehensive analysis of germline and tumor samples may be useful to predict the risk of secondary cancers.

Fucosyltransferase 8 (FUT8) and core fucose expression in oxidative stress response.

GlycoMaple is a new tool to predict glycan structures based on the expression levels of 950 genes encoding glycan biosynthesis-related enzymes and proteins using RNA-seq data. The antioxidant response, protecting cells from oxidative stress, has been focused on because its activation may relieve pathological conditions, such as neurodegenerative diseases. Genes involved in the antioxidant response are defined within the GO:0006979 category, including 441 human genes. Fifteen genes overlap between the glycan biosynthesis-related genes defined by GlycoMaple and the antioxidant response genes defined by GO:0006979, one of which is FUT8. 5-Hydroxy-4-phenyl-butenolide (5H4PB) extracted from Chinese aromatic vinegar induces the expression of a series of antioxidant response genes that protect cells from oxidative stress via activation of the nuclear factor erythroid 2-related factor 2-antioxidant response element pathway. Here, we show that FUT8 is upregulated in both our RNA-seq data set of 5H4PB-treated cells and publicly available RNA-seq data set of cells treated with another antioxidant, sulforaphane. Applying our RNA-seq data set to GlycoMaple led to a prediction of an increase in the core fucose of N-glycan that was confirmed by flow cytometry using a fucose-binding lectin. These results suggest that FUT8 and core fucose expression may increase upon the antioxidant response.

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.

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.

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.

Determination of language areas in patients with epilepsy using the super-selective Wada test.

The Wada test is the gold standard for determining language-dominant hemisphere. However, the precise determination of language areas in each patient requires more invasive methods, such as electrocortical stimulation. Some studies have reported the use of anesthetic injection into selective cerebral arteries to predict postoperative function. To assess the function of the anterior and posterior language areas separately, we developed an advanced test named the "super-selective Wada test" (ssWada). The ssWada procedure is as follows: an endovascular neurosurgeon identifies the arterial branches of the middle cerebral artery (MCA) perfusing the anterior language area of the inferior frontal gyrus and the posterior language area of the posterior part of the superior temporal gyrus using angiography. Behavioral neurologists assess language symptoms before and after propofol administration using a microcatheter tip in the selected arterial branch. From 30 serial patients with epilepsy who underwent ssWada test at Tohoku University Hospital, we retrospectively reviewed patients in whom multiple areas in the bilateral MCA region was examined. Eight cases were identified in this study. All eight cases had been considered for resection of the area overlapping the classical language area. Three of the eight cases were left-dominant, and the within-hemisphere distribution was also considered typical. One case was determined to be left-dominant but atypical in the intra-hemispheric functional distribution. Two cases were right-dominant, and the intra-hemispheric functional distribution was considered a mirror image of the typical pattern. The remaining two cases were considered atypical, not only in terms of bilateral language function, but also in terms of anterior-posterior functional distribution. This case series demonstrates the potential utility of ssWada in revealing separate function of the anterior and posterior language areas. The ssWada allows simulation of local surgical brain resection and detailed investigation of language function, which potentially contributes to planning the resection area. Although indications for ssWada are quite limited, it could play a complementary role to noninvasive testing because it provides information related to resection using a different approach.

Abnormal heart rate variability during non-REM sleep and postictal generalized EEG suppression in focal epilepsy.

OBJECTIVE: To identify any relationship between abnormal interictal heart rate variability (HRV) during sleep and the occurrence of postictal generalized EEG suppression (PGES), a potential biomarker of sudden unexpected death in epilepsy (SUDEP), in patients with focal epilepsy. METHODS: We retrospectively reviewed 34 consecutive patients with focal epilepsy, who presented with focal to bilateral tonic-clonic seizures (FBTCSs), 19 with PGES (PGES+) and 15 without PGES (PGES-), and 14 patients without epilepsy as controls. HRV spectrum analysis was performed for periods of 10 minutes during wakefulness, non-REM sleep, and REM sleep. HRV spectra consisted of high frequency (HF) and low frequency (LF) regions. Normalized HF representing parasympathetic activities and normalized LF were compared between the PGES+, PGES-, and control groups for each sleep stage. RESULTS: The PGES+ group showed significantly decreased normalized HF and increased normalized LF during non-REM sleep compared to the PGES- (P < 0.05) and control (P < 0.01) groups, but not during wakefulness or REM sleep. CONCLUSIONS: Abnormal interictal HRV during non-REM sleep, indicating abnormally decreased parasympathetic activities, was associated with PGES in patients with FBTCS. SIGNIFICANCE: This study proposes to further investigate the relation between decreased parasympathetic activities and PGES.

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.

Collection of 2429 constrained headshots of 277 volunteers for deep learning.

Deep learning has rapidly been filtrating many aspects of human lives. In particular, image recognition by convolutional neural networks has inspired numerous studies in this area. Hardware and software technologies as well as large quantities of data have contributed to the drastic development of the field. However, the application of deep learning is often hindered by the need for big data and the laborious manual annotation thereof. To experience deep learning using the data compiled by us, we collected 2429 constrained headshot images of 277 volunteers. The collection of face photographs is challenging in terms of protecting personal information; we therefore established an online procedure in which both the informed consent and image data could be obtained. We did not collect personal information, but issued agreement numbers to deal with withdrawal requests. Gender and smile labels were manually and subjectively annotated only from the appearances, and final labels were determined by majority among our team members. Rotated, trimmed, resolution-reduced, decolorized, and matrix-formed data were allowed to be publicly released. Moreover, simplified feature vectors for data sciences were released. We performed gender and smile recognition by building convolutional neural networks based on the Inception V3 model with pre-trained ImageNet data to demonstrate the usefulness of our dataset.

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.

NUDT15 variants confer high incidence of second malignancies in children with acute lymphoblastic leukemia.

The effect of genetic variation on second malignant neoplasms (SMNs) remains unclear. First, we identified the pathogenic germline variants in cancer-predisposing genes among 15 children with SMNs after childhood leukemia/lymphoma using whole-exome sequencing. Because the prevalence was low, we focused on the association between SMNs and NUDT15 in primary acute lymphoblastic leukemia (ALL) cases. NUDT15 is one of the 6-mercaptopurine (6-MP) metabolic genes, and its variants are common in East Asian individuals. The prevalence of NUDT15 hypomorphic variants was higher in patients with SMNs (n = 14; 42.9%) than in the general population in the gnomAD database (19.7%; P = .042). In the validation study with a cohort of 438 unselected patients with ALL, the cumulative incidence of SMNs was significantly higher among those with (3.0%; 95% confidence interval [CI], 0.6% to 9.4%) than among those without NUDT15 variants (0.3%; 95% CI, 0.0% to 1.5%; P = .045). The 6-MP dose administered to patients with ALL with a NUDT15 variant was higher than that given to those without SMNs (P = .045). The 6-MP-related mutational signature was observed in SMN specimens after 6-MP exposure. In cells exposed to 6-MP, a higher level of 6-MP induced DNA damage in NUDT15-knockdown induced pluripotent stem cells. Our study indicates that NUDT15 variants may confer a risk of SMNs after treatment with 6-MP in patients with ALL.

Robust classification of cell cycle phase and biological feature extraction by image-based deep learning.

Across the cell cycle, the subcellular organization undergoes major spatiotemporal changes that could in principle contain biological features that could potentially represent cell cycle phase. We applied convolutional neural network-based classifiers to extract such putative features from the fluorescence microscope images of cells stained for the nucleus, the Golgi apparatus, and the microtubule cytoskeleton. We demonstrate that cell images can be robustly classified according to G1/S and G2 cell cycle phases without the need for specific cell cycle markers. Grad-CAM analysis of the classification models enabled us to extract several pairs of quantitative parameters of specific subcellular features as good classifiers for the cell cycle phase. These results collectively demonstrate that machine learning-based image processing is useful to extract biological features underlying cellular phenomena of interest in an unbiased and data-driven manner.

A prospective compound screening contest identified broader inhibitors for Sirtuin 1.

Potential inhibitors of a target biomolecule, NAD-dependent deacetylase Sirtuin 1, were identified by a contest-based approach, in which participants were asked to propose a prioritized list of 400 compounds from a designated compound library containing 2.5 million compounds using in silico methods and scoring. Our aim was to identify target enzyme inhibitors and to benchmark computer-aided drug discovery methods under the same experimental conditions. Collecting compound lists derived from various methods is advantageous for aggregating compounds with structurally diversified properties compared with the use of a single method. The inhibitory action on Sirtuin 1 of approximately half of the proposed compounds was experimentally accessed. Ultimately, seven structurally diverse compounds were identified.

Relationship between conformation shift and disease related variation sites in ATP-binding cassette transporter proteins.

Transport of small molecules across the cell membrane is a crucial biological mechanism for the maintenance of the cell activity. ABC transporter family is a huge group in the transporter membrane proteins and actively transports the substrates using the energy derived from ATP hydrolysis. In humans, there are 48 distinct genes for ABC transporters. A variation of a single amino acid in the amino acid sequence of ABC transporter has been known to be linked with certain disease. The mechanism of the onset of the disease by the variation is, however, still unclear. Recent progress in the method to measure the structures of huge membrane proteins has enabled determination of the 3D structures of ABC transporters and the accumulation of coordinate data of ABC transporter has enabled us to obtain clues for the onset of the disease caused by a single variation of amino acid residue. We compared the structures of ABC transporter in apo and ATP-binding forms and found a possible conformation shift around pivot-like residues in the transmembrane domains. When this conformation change in ABC transporter and the location of pathogenic variation were compared, we found a reasonable match between the two, explaining the onset of the disease by the variation. They likely cause impairment of the pivot-like movement, weakening of ATP binding and weakening of membrane surface interactions. These findings will give a new interpretation of the variations on ABC transporter genes and pave a way to analyse the effect of variation on protein structure and function.

Passive cigarette smoking changes the circadian rhythm of clock genes in rat intervertebral discs.

We aimed to elucidate the molecular changes in intervertebral discs (IVDs) caused by passive smoking. Rats were subjected to 8 weeks of passive smoking; thereafter, their lumbar vertebrae were harvested. The annulus fibrosus and cartilage endplate (AF/CEP) were harvested together, and the nucleus pulposus (NP) was isolated separately. The expression of 27,342 rat genes was analyzed. In 3 "nonsmoking" rats, 96 of 112 genes whose expression varied ≥10-fold between the AF/CEP and NP were more highly expressed in the AF/CEP. With these differentially expressed genes, we uncovered novel AF/CEP and NP marker genes and indicated their possible novel functions. Although passive smoking induced less marked alteration in the gene expression profiles of both the AF/CEP and NP, multiple clock-related genes showed altered expression. These genes were expressed with a circadian rhythm in IVD cells, and most genes showed a phase shift of -6 to -9 h induced by passive smoking. Some clock-related genes showed abolished oscillation in the NP. Passive smoking also changed the expression levels of proteases and protease inhibitors and reduced the expression of NP marker genes. Thus, passive smoking induces changes in the circadian rhythm of a peripheral clock (IVD clock) that might be involved in molecular events related to IVD degeneration.

Multiple gamma-glutamylation: a novel type of post-translational modification in a diapausing Artemia cyst protein.

A highly hydrophilic, glutamate-rich protein was identified in the aqueous phenol extract from the cytosolic fraction of brine shrimp (Artemia franciscana) diapausing cysts and termed Artemia phenol soluble protein (PSP). Mass spectrometric analysis revealed the presence of many protein peaks around m/z 11,000, separated by 129 atomic mass units; this value corresponds to that of glutamate, which is strongly suggestive of heterogeneous polyglutamylation. Polyglutamylation has long been known as the functionally important post-translational modification of tubulins, which carry poly(L-glutamic acid) chains of heterogeneous length branching off from the main chain at the gamma-carboxy groups of a few specific glutamate residues. In Artemia PSP, however, Edman degradation of enzymatic peptides revealed that at least 13, and presumably 16, glutamate residues were modified by the attachment of a single L-glutamate, representing a hitherto undescribed type of post-translational modification: namely, multiple gamma-glutamylation or the addition of a large number of glutamate residues along the polypeptide chain. Although biological significance of PSP and its modification is yet to be established, suppression of in vitro thermal aggregation of lactate dehydrogenase by glutamylated PSP was observed.

[The levels of recognition and hand functions of patients with Alzheimer's disease: a pilot study based on the simple test for evaluating hand function].

AIM: The purpose of this study was to analyze the levels of recognition and hand functions of the patients with Alzheimer's disease (AD) based on the Simple Test for Evaluating Hand Function (STEF). METHODS: Forty subjects with AD (AD group) and 20 without AD (control group) in a facility for elderly health care services participated in this study. The AD group consisted of 20 subjects with Clinical Dementia Rating (CDR) 1 and 20 other subjects with CDR2. RESULTS: The results were as follows: Time required for full motion of right hand dexterity was significantly longer in both the CDR1 and CDR2 groups compared to the control group (p < 0.05, p < 0.01). Time required for the control group and CDR1 group were significantly shorter for motion of right hand dexterity compared to the left hand. No significant difference was observed between the times of right and left hand dexterity among subjects of CDR2. The disparity of right-left hand dexterity decreased in the order of control, CDR1 and CDR2. The correlations between the sub-items for CDR scale and the right-left hand disparity were significant in "memory" and "judgment". The right-left hand disparity decreased with the shortened time for left hand motion, as the "memory" or "judgment" scores increased. CONCLUSION: These findings suggest that right hand dexterity speed in patients in the mild AD stage was low when their dominant hand superiority declined. Moreover the decline of dominant hand superiority correlated with the recognition of memory and judgment.