Unveiling interatomic distances influencing the reaction coordinates in alanine dipeptide isomerization: An explainable deep learning approach.

The present work shows that the free energy landscape associated with alanine dipeptide isomerization can be effectively represented by specific interatomic distances without explicit reference to dihedral angles. Conventionally, two stable states of alanine dipeptide in vacuum, i.e., C7eq (ß-sheet structure) and C7ax (left handed α-helix structure), have been primarily characterized using the main chain dihedral angles, φ (C-N-Cα-C) and ψ (N-Cα-C-N). However, our recent deep learning combined with the "Explainable AI" (XAI) framework has shown that the transition state can be adequately captured by a free energy landscape using φ and θ (O-C-N-Cα) [Kikutsuji et al., J. Chem. Phys. 156, 154108 (2022)]. In the perspective of extending these insights to other collective variables, a more detailed characterization of the transition state is required. In this work, we employ interatomic distances and bond angles as input variables for deep learning rather than the conventional and more elaborate dihedral angles. Our approach utilizes deep learning to investigate whether changes in the main chain dihedral angle can be expressed in terms of interatomic distances and bond angles. Furthermore, by incorporating XAI into our predictive analysis, we quantified the importance of each input variable and succeeded in clarifying the specific interatomic distance that affects the transition state. The results indicate that constructing a free energy landscape based on the identified interatomic distance can clearly distinguish between the two stable states and provide a comprehensive explanation for the energy barrier crossing.

Widespread and transgenerational retrotransposon activation in inter- and intraspecies recombinant inbred populations of Lotus japonicus.

Transposable elements (TEs) constitute a large proportion of genomes of multicellular eukaryotes, including flowering plants. TEs are normally maintained in a silenced state and their transpositions rarely occur. Hybridization between distant species has been regarded as a 'shock' that stimulates genome reorganization, including TE mobilization. However, whether crosses between genetically close parents that result in viable and fertile offspring can induce TE transpositions has remained unclear. Here, we investigated the activation of long terminal repeat (LTR) retrotransposons in three Lotus japonicus recombinant inbred line (RIL) populations. We found that at least six LTR retrotransposon families were activated and transposed in 78% of the RILs investigated. LORE1a, one of the transposed LTR retrotransposons, showed transgenerational epigenetic activation, indicating the long-term effects of epigenetic instability induced by hybridization. Our study highlights TE activation as an unexpectedly common event in plant reproduction.

Characterization of FLOWERING LOCUS C 5 in Brassica rapa L.

Brassica rapa L., which includes Chinese cabbage, turnip, and pak choi, has more complex flowering time regulation than does Arabidopsis thaliana due to the presence of multiple paralogous flowering time genes. FLOWERING LOCUS C (FLC) is one of the key genes regulating the flowering time, and B. rapa has four FLC paralogs. BrFLC5 on the reference genome is deemed a pseudogene because of a mutation (from G to A) in the splice site of the third intron, but there are some accessions with a G nucleotide in the splice site. In this study, we genotyped 310 B. rapa accessions and found that 19 had homozygous and 81 had heterozygous putative functional BrFLC5 alleles. Accessions of turnip showed the highest proportion with a functional BrFLC5 allele. BrFLC5 acts as a floral repressor when overexpressed in A. thaliana. The BrFLC5 expression level varied in pre-vernalized plants, and this transcriptional variation was not associated with the G/A polymorphism in the third intron. Three accessions having a higher BrFLC5 expression in pre-vernalized plants had a 584-bp insertion in the promoter region. Many regions homologous to this 584-bp sequence are present in the B. rapa genome, and this 584-bp inserted region has tandem duplications of an AT-rich sequence in its central region. The possibility that a high expression of a functional BrFLC5 could contribute to producing premature bolting-resistant lines in B. rapa vegetables is discussed. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-023-01405-0.

Quantitative trait locus (QTL) analysis and fine-mapping for Fusarium oxysporum disease resistance in Raphanus sativus using GRAS-Di technology.

Fusarium wilt is a significant disease in radish, but the genetic mechanisms controlling yellows resistance (YR) are not well understood. This study aimed to identify YR-QTLs and to fine-map one of them using F2:3 populations developed from resistant and susceptible radish parents. In this study, two high-density genetic maps each containing shared co-dominant markers and either female or male dominant markers that spanned 988.6 and 1127.5 cM with average marker densities of 1.40 and 1.53 cM, respectively, were generated using Genotyping by Random Amplicon Sequencing-Direct (GRAS-Di) technology. We identified two YR-QTLs on chromosome R2 and R7, and designated the latter as ForRs1 as the major QTL. Fine mapping narrowed down the ForRs1 locus to a 195 kb region. Among the 16 predicted genes in the delimited region, 4 genes including two receptor-like protein and -kinase genes (RLP/RLK) were identified as prime candidates for ForRs1 based on the nucleotide sequence comparisons between the parents and their predicted functions. This study is the first to use a GRAS-Di for genetic map construction of cruciferous crops and fine map the YR-QTL on the R7 chromosome of radish. These findings will provide groundbreaking insights into radish YR breeding and understanding the genetics of YR mechanism.

Molecular Design of FRET Probes Based on Domain Rearrangement of Protein Disulfide Isomerase for Monitoring Intracellular Redox Status.

Multidomain proteins can exhibit sophisticated functions based on cooperative interactions and allosteric regulation through spatial rearrangements of the multiple domains. This study explored the potential of using multidomain proteins as a basis for Förster resonance energy transfer (FRET) biosensors, focusing on protein disulfide isomerase (PDI) as a representative example. PDI, a well-studied multidomain protein, undergoes redox-dependent conformational changes, enabling the exposure of a hydrophobic surface extending across the b' and a' domains that serves as the primary binding site for substrates. Taking advantage of the dynamic domain rearrangements of PDI, we developed FRET-based biosensors by fusing the b' and a' domains of thermophilic fungal PDI with fluorescent proteins as the FRET acceptor and donor, respectively. Both experimental and computational approaches were used to characterize FRET efficiency in different redox states. In vitro and in vivo evaluations demonstrated higher FRET efficiency of this biosensor in the oxidized form, reflecting the domain rearrangement and its responsiveness to intracellular redox environments. This novel approach of exploiting redox-dependent domain dynamics in multidomain proteins offers promising opportunities for designing innovative FRET-based biosensors with potential applications in studying cellular redox regulation and beyond.

Comparative transcriptome analysis during tuberous stem formation in Kohlrabi (B. oleracea var. gongylodes) at early growth periods (seedling stages).

Tuberous stem of kohlrabi is an important agronomic trait, however, the molecular basis of tuberization is poorly understood. To elucidate the tuberization mechanism, we conducted a comparative transcriptomic analysis between kohlrabi and broccoli at 10 and 20 days after germination (DAG) as tuberous stem initiated between these time points. A total of 5580 and 2866 differentially expressed transcripts (DETs) were identified between genotypes (kohlrabi vs. broccoli) and growth stages (10 DAG vs. 20 DAG), respectively, and most of the DETs were down-regulated in kohlrabi. Gene ontology (GO) and KEGG pathway enrichment analyses showed that the DETs between genotypes are involved in cell wall loosening and expansion, cell cycle and division, carbohydrate metabolism, hormone transport, hormone signal transduction and in several transcription factors. The DETs identified in those categories may directly/indirectly relate to the initiation and development of tuberous stem in kohlrabi. In addition, the expression pattern of the hormone synthesis related DETs coincided with the endogenous IAA, IAAsp, GA, ABA, and tZ profiles in kohlrabi and broccoli seedlings, that were revealed in our phytohormone analysis. This is the first report on comparative transcriptome analysis for tuberous stem formation in kohlrabi at early growth periods. The resulting data could provide significant insights into the molecular mechanism underlying tuberous stem development in kohlrabi as well as in other tuberous organ forming crops.

Explaining reaction coordinates of alanine dipeptide isomerization obtained from deep neural networks using Explainable Artificial Intelligence (XAI).

A method for obtaining appropriate reaction coordinates is required to identify transition states distinguishing the product and reactant in complex molecular systems. Recently, abundant research has been devoted to obtaining reaction coordinates using artificial neural networks from deep learning literature, where many collective variables are typically utilized in the input layer. However, it is difficult to explain the details of which collective variables contribute to the predicted reaction coordinates owing to the complexity of the nonlinear functions in deep neural networks. To overcome this limitation, we used Explainable Artificial Intelligence (XAI) methods of the Local Interpretable Model-agnostic Explanation (LIME) and the game theory-based framework known as Shapley Additive exPlanations (SHAP). We demonstrated that XAI enables us to obtain the degree of contribution of each collective variable to reaction coordinates that is determined by nonlinear regressions with deep learning for the committor of the alanine dipeptide isomerization in vacuum. In particular, both LIME and SHAP provide important features to the predicted reaction coordinates, which are characterized by appropriate dihedral angles consistent with those previously reported from the committor test analysis. The present study offers an AI-aided framework to explain the appropriate reaction coordinates, which acquires considerable significance when the number of degrees of freedom increases.

Learning reaction coordinates via cross-entropy minimization: Application to alanine dipeptide.

We propose a cross-entropy minimization method for finding the reaction coordinate from a large number of collective variables in complex molecular systems. This method is an extension of the likelihood maximization approach describing the committor function with a sigmoid. By design, the reaction coordinate as a function of various collective variables is optimized such that the distribution of the committor pB * values generated from molecular dynamics simulations can be described in a sigmoidal manner. We also introduce the L2-norm regularization used in the machine learning field to prevent overfitting when the number of considered collective variables is large. The current method is applied to study the isomerization of alanine dipeptide in vacuum, where 45 dihedral angles are used as candidate variables. The regularization parameter is determined by cross-validation using training and test datasets. It is demonstrated that the optimal reaction coordinate involves important dihedral angles, which are consistent with the previously reported results. Furthermore, the points with pB *∼0.5 clearly indicate a separatrix distinguishing reactant and product states on the potential of mean force using the extracted dihedral angles.

F1-ATPase conformational cycle from simultaneous single-molecule FRET and rotation measurements.

Despite extensive studies, the structural basis for the mechanochemical coupling in the rotary molecular motor F1-ATPase (F1) is still incomplete. We performed single-molecule FRET measurements to monitor conformational changes in the stator ring-α3ß3, while simultaneously monitoring rotations of the central shaft-γ. In the ATP waiting dwell, two of three ß-subunits simultaneously adopt low FRET nonclosed forms. By contrast, in the catalytic intermediate dwell, two ß-subunits are simultaneously in a high FRET closed form. These differences allow us to assign crystal structures directly to both major dwell states, thus resolving a long-standing issue and establishing a firm connection between F1 structure and the rotation angle of the motor. Remarkably, a structure of F1 in an ε-inhibited state is consistent with the unique FRET signature of the ATP waiting dwell, while most crystal structures capture the structure in the catalytic dwell. Principal component analysis of the available crystal structures further clarifies the five-step conformational transitions of the αß-dimer in the ATPase cycle, highlighting the two dominant modes: the opening/closing motions of ß and the loosening/tightening motions at the αß-interface. These results provide a new view of tripartite coupling among chemical reactions, stator conformations, and rotary angles in F1-ATPase.

Allele specific DNA marker for fusarium resistance gene FocBo1 in Brassica oleracea.

The fusarium yellows resistance (YR) gene FocBo1 was previously identified and the DNA markers were developed to assist the breeding of YR cultivars in Brassica oleracea. However, the further analysis revealed discrepancies between the phenotypes and the genotypes predicted by those DNA markers in cabbage commercial cultivars. Since this discrepancy seemed to be due to unknown susceptible alleles of focbo1, we sequenced the gene in 19 accessions to determine the sequence variations between alleles and found that there were two resistant FocBo1 alleles and six susceptible alleles in the investigated population. The newly designed PCR markers detected three mutations in the susceptible alleles that generate premature termination codons. These were shown to accurately distinguish resistant and susceptible alleles in more than 200 accessions of B. oleracea inbred lines and cultivars. The study revealed that the locus is represented by 37.2% resistant and 62.8% susceptible alleles within seventy-eight commercial cultivars. Structural analysis of the gene revealed that a part of the allelic variation comes from intragenic recombination between alleles. Our results enable a more precise prediction of the phenotype by marker assisted selection, promoting the production of YR cultivars in B. oleracea.

Elasticity, friction, and pathway of γ-subunit rotation in FoF1-ATP synthase.

We combine molecular simulations and mechanical modeling to explore the mechanism of energy conversion in the coupled rotary motors of FoF1-ATP synthase. A torsional viscoelastic model with frictional dissipation quantitatively reproduces the dynamics and energetics seen in atomistic molecular dynamics simulations of torque-driven γ-subunit rotation in the F1-ATPase rotary motor. The torsional elastic coefficients determined from the simulations agree with results from independent single-molecule experiments probing different segments of the γ-subunit, which resolves a long-lasting controversy. At steady rotational speeds of ∼ 1 kHz corresponding to experimental turnover, the calculated frictional dissipation of less than k(B)T per rotation is consistent with the high thermodynamic efficiency of the fully reversible motor. Without load, the maximum rotational speed during transitions between dwells is reached at ∼ 1 MHz. Energetic constraints dictate a unique pathway for the coupled rotations of the Fo and F1 rotary motors in ATP synthase, and explain the need for the finer stepping of the F1 motor in the mammalian system, as seen in recent experiments. Compensating for incommensurate eightfold and threefold rotational symmetries in Fo and F1, respectively, a significant fraction of the external mechanical work is transiently stored as elastic energy in the γ-subunit. The general framework developed here should be applicable to other molecular machines.

Breeding of lilies and tulips-Interspecific hybridization and genetic background.

Lilies and tulips (Liliaceae family) are economically very important ornamental bulbous plants. Here, we summarize major breeding goals, the role of an integrated method of cut-style pollination and fertilization followed by embryo rescue and mitotic and meiotic polyploidization involved in new assortment development. Both crops have been subjected to extensive interspecific hybridization followed by selection. Additionally, spontaneous polyploidization has played a role in their evolution. In lilies, there is a tendency to replace diploids with polyploid cultivars, whereas in tulip a majority of the cultivars that exist today are still diploid except for triploid Darwin hybrid tulips. The introduction of molecular cytogenetic techniques such as genomic in situ hybridization (GISH) permitted the detailed studies of genome composition in lily and tulip interspecific hybrids and to follow the chromosome inheritance in interspecific crosses. In addition, this review presents the latest information on phylogenetic relationship in lily and tulip and recent developments in molecular mapping using different DNA molecular techniques.

Genome re-sequencing, SNP analysis, and genetic mapping of the parental lines of a commercial F1 hybrid cultivar of Chinese cabbage.

The genome-wide characterization of single nucleotide polymorphism (SNP) between cultivars or between inbred lines contributes to the creation of genetic markers that are important for plant breeding. Functional markers derived from polymorphisms within genes that affect phenotypic variation are especially valuable in plant breeding. Here, we report on the genome re-sequencing and analysis of the two parental inbred lines of the commercial F1 hybrid Chinese cabbage cultivar "W77". Through the genome-wide identification and classification of the SNPs and indels present in each parental line, we identified about 1,500 putative non-functional genes in each parent. We designed cleaved amplified polymorphic sequence (CAPS) markers using specific mutations found at Eco RI restriction sites in the parental lines and confirmed their Mendelian segregation by constructing a linkage map using 96 F2 plants derived from the F1 hybrid cultivar, "W77". Our results and data will be a useful genomic resource for future studies of gene function and metagenomic studies in Chinese cabbage.

The production and characterization of a BoFLC2 introgressed Brassica rapa by repeated backcrossing to an F1.

Flowering time is an important agronomic trait for Brassica rapa crops, and previous breeding work in Brassica has successfully transmitted other important agronomic traits from donor species. However, there has been no previous attempts to produce hybrids replacing the original Brassica FLC alleles with alien FLC alleles. In this paper, we introduce the creation of a chromosome substitution line (CSSL) containing a homozygous introgression of Flowering Locus C from Brassica oleracea (BoFLC2) into a B. rapa genomic background, and characterize the CSSL line with respect to the parental cultivars. The preferential transmission of alien chromosome inheritance and the pattern of transmission observed during the production of the CSSLs are also discussed.

IntroMap: a signal analysis based method for the detection of genomic introgressions.

BACKGROUND: Breeding programs often rely on marker-assisted tests or variant calling of next generation sequence (NGS) data to identify regions of genomic introgression arising from the hybridization of two plant species. In this paper we present IntroMap, a bioinformatics pipeline for the screening of candidate plants through the application of signal processing techniques to NGS data, using alignment to a reference genome sequence (annotation is not required) that shares homology with the recurrent parental cultivar, and without the need for de novo assembly of the read data or variant calling. RESULTS: We show the accurate identification of introgressed genomic regions using both in silico simulated genomes, and a hybridized cultivar data set using our pipeline. Additionally we show, through targeted marker-based assays, validation of the IntroMap predicted regions for the hybrid cultivar. CONCLUSIONS: This approach can be used to automate the screening of large populations, reducing the time and labor required, and can improve the accuracy of the detection of introgressed regions in comparison to a marker-based approach. In contrast to other approaches that generally rely upon a variant calling step, our method achieves accurate identification of introgressed regions without variant calling, relying solely upon alignment.

Transition path sampling of rare events by shooting from the top.

Transition path sampling is a powerful tool in the study of rare events. Shooting trial trajectories from configurations along existing transition paths proved particularly efficient in the sampling of reactive trajectories. However, most shooting attempts tend not to result in transition paths, in particular in cases where the transition dynamics has diffusive character. To overcome the resulting efficiency problem, we developed an algorithm for "shooting from the top." We first define a shooting range through which all paths have to pass and then shoot off trial trajectories only from within this range. For a well chosen shooting range, nearly every shot is successful, resulting in an accepted transition path. To deal with multiple mechanisms, weighted shooting ranges can be used. To cope with the problem of unsuitably placed shooting ranges, we developed an algorithm that iteratively improves the location of the shooting range. The transition path sampling procedure is illustrated for models of diffusive and Langevin dynamics. The method should be particularly useful in cases where the transition paths are long so that only relatively few shots are possible, yet reasonable order parameters are known.

Comparison of transcriptome profiles by Fusarium oxysporum inoculation between Fusarium yellows resistant and susceptible lines in Brassica rapa L.

KEY MESSAGE: Resistant and susceptible lines in Brassica rapa have different immune responses against Fusarium oxysporum inoculation. Fusarium yellows caused by Fusarium oxysporum f. sp. conglutinans (Foc) is an important disease of Brassicaceae; however, the mechanism of how host plants respond to Foc is still unknown. By comparing with and without Foc inoculation in both resistant and susceptible lines of Chinese cabbage (Brassica rapa var. pekinensis), we identified differentially expressed genes (DEGs) between the bulked inoculated (6, 12, 24, and 72 h after inoculation (HAI)) and non-inoculated samples. Most of the DEGs were up-regulated by Foc inoculation. Quantitative real-time RT-PCR showed that most up-regulated genes increased their expression levels from 24 HAI. An independent transcriptome analysis at 24 and 72 HAI was performed in resistant and susceptible lines. GO analysis using up-regulated genes at 24 HAI indicated that Foc inoculation activated systemic acquired resistance (SAR) in resistant lines and tryptophan biosynthetic process and responses to chitin and ethylene in susceptible lines. By contrast, GO analysis using up-regulated genes at 72 HAI showed the overrepresentation of some categories for the defense response in susceptible lines but not in the resistant lines. We also compared DEGs between B. rapa and Arabidopsis thaliana after F. oxysporum inoculation at the same time point, and identified genes related to defense response that were up-regulated in the resistant lines of Chinese cabbage and A. thaliana. Particular genes that changed expression levels overlapped between the two species, suggesting that they are candidates for genes involved in the resistance mechanisms against F. oxysporum.

Tube feeding decreases pneumonia rate in patients with severe dementia: comparison between pre- and post-intervention.

BACKGROUND: It is widely supposed that there is no benefit, including extended survival and decreased rate of pneumonia, in patients with severe dementia receiving enteral tube feeding (TF). However, there have been few studies comparing the frequency of pneumonia before and after TF in severe dementia. METHODS: Nine psychiatric hospitals in Okayama Prefecture participated in this retrospective survey. All inpatients fulfilling the entry criteria were evaluated. All subjects suffered from difficulty in oral intake. Attending physicians thought that the patients could not live without long-term artificial nutrition, and they decided whether or not to make use of long-term artificial nutrition from January 1, 2014 to December 31, 2014. RESULTS: We evaluated 58 patients including 46 with TF and 12 without. The mean age of all patients was 79.6 ± 9.0 years old. Patients with probable Alzheimer's disease (n = 38) formed the biggest group, and those with vascular dementia the second (n = 14). Median survival times were 23 months among patients with TF and two months among patients without TF. The start of TF decreased the frequency of pneumonia and the use of intravenous antibiotics. CONCLUSIONS: TF decreased pneumonia and antibiotic use, even in patients with severe dementia. The results of this study do not necessarily indicate that we should administer TF to patients with severe dementia. We should consider the quality of life of patients carefully before deciding the use or disuse of TF for patients with severe dementia.

Molecular and cellular characteristics of hybrid vigour in a commercial hybrid of Chinese cabbage.

BACKGROUND: Heterosis or hybrid vigour is a phenomenon in which hybrid progeny exhibit superior performance compared to their parental inbred lines. Most commercial Chinese cabbage cultivars are F1 hybrids and their level of hybrid vigour is of critical importance and is a key selection criterion in the breeding system. RESULTS: We have characterized the heterotic phenotype of one F1 hybrid cultivar of Chinese cabbage and its parental lines from early- to late-developmental stages of the plants. Hybrid cotyledons are larger than those of the parents at 4 days after sowing and biomass in the hybrid, determined by the fresh weight of leaves, is greater than that of the larger parent line by approximately 20% at 14 days after sowing. The final yield of the hybrid harvested at 63 days after sowing is 25% greater than the yield of the better parent. The larger leaves of the hybrid are a consequence of increased cell size and number of the photosynthetic palisade mesophyll cells and other leaf cells. The accumulation of plant hormones in the F1 was within the range of the parental levels at both 2 and 10 days after sowing. Two days after sowing, the expression levels of chloroplast-targeted genes in the cotyledon cells were upregulated in the F1 hybrid relative to their mid parent values. Shutdown of chlorophyll biosynthesis in the cotyledon by norflurazon prevented the increased leaf area in the F1 hybrid. CONCLUSIONS: In the cotyledons of F1 hybrids, chloroplast-targeted genes were upregulated at 2 days after sowing. The increased activity levels of this group of genes suggested that their differential transcription levels could be important for establishing early heterosis but the increased transcription levels were transient. Inhibition of the photosynthetic process in the cotyledon reduced heterosis in later seedling stages. These observations suggest early developmental events in the germinating seedling of the hybrid may be important for later developmental vigour and yield advantage.