Optimization of the Preparation Process of Glucuronomannan Oligosaccharides and Their Effects on the Gut Microbiota in MPTP-Induced PD Model Mice.

Parkinson's disease (PD) is a prevalent neurodegenerative disorder, and accumulating evidence suggests a link between dysbiosis of the gut microbiota and the onset and progression of PD. In our previous investigations, we discovered that intraperitoneal administration of glucuronomannan oligosaccharides (GMn) derived from Saccharina japonica exhibited neuroprotective effects in a 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP)-induced PD mouse model. However, the complicated preparation process, difficulties in isolation, and remarkably low yield have constrained further exploration of GMn. In this study, we optimized the degradation conditions in the preparation process of GMn through orthogonal experiments. Subsequently, an MPTP-induced PD model was established, followed by oral administration of GMn. Through a stepwise optimization, we successfully increased the yield of GMn, separated from crude fucoidan, from 1~2/10,000 to 4~8/1000 and indicated the effects on the amelioration of MPTP-induced motor deficits, preservation of dopamine neurons, and elevation in striatal neurotransmitter levels. Importantly, GMn mitigated gut microbiota dysbiosis induced by MPTP in mice. In particular, GM2 significantly reduced the levels of Akkermansia, Verrucomicrobiota, and Lactobacillus, while promoting the abundance of Roseburia and Prevotella compared to the model group. These findings suggest that GM2 can potentially suppress PD by modulating the gut microbiota, providing a foundation for the development of a novel and effective anti-PD marine drug.

Identification and genome-wide analysis provide insights into the genetic diversity and biotechnological potentials of novel cold-adapted Acinetobacter strain.

Extreme cold environments, such as polar regions or high-altitude mountains, are known for their challenging conditions including low temperatures, high salinity, and limited nutrient availability. Microbes that thrive in these environments have evolved specialized strategies to survive and function under such harsh conditions. The study aims to identify, sequence the genome, perform genome assembly, and conduct a comparative genome-wide analysis of Acinetobacter sp. strain P1, which was isolated from the Batura glacier regions of Pakistan. A basic local alignment search tool of NCBI using 16 s RNA gene sequence confirmed the strain Acinetobacter following phylogenetic analysis revealed that strain P1 clustered with Acinetobacter sp. strain AcBz01. The high-throughput Genome sequencing was done by the NovaSeq 6000 sequencing system following de novo genome assembly reported 23 contigs, a genome size of 3,732,502 bp containing approximately 3489 genes and 63 RNAs (60 tRNA, 3 rRNA). The comparative genome analysis revealed that Acinetobacter sp. strain P1 exhibited the highest homology with the Acinetobacter baumannii ATCC 17978 genome and encompassed 1668 indispensable genes, 1280 conserved genes 1821 specific genes suggesting high genomic plasticity and evolutionary diversity. The genes with functional assignments include exopolysaccharide phosphotransferase enzyme, cold-shock proteins, T6SS, membrane modifications, antibiotic resistance, and set of genes related to a wide range of metabolic characteristics such as exopolysaccharides were also present. Moreover, the structural prediction analysis of EPS proteins reveals that structural flexibility allows for conformational modifications during catalysis, which boosts or increases the catalytic effectiveness at lower temperatures. Overall, the identification of Acinetobacter, a cold-adapted bacterium, offers promising applications in bioremediation, enzyme production, food preservation, pharmaceutical development, and astrobiology. Further research and exploration of these microorganisms can unlock their full biotechnological potential and contribute to various industries and scientific endeavors.

Proteomics unite traditional toxicological assessment methods to evaluate the toxicity of iron oxide nanoparticles.

Iron oxide nanoparticles (IONPs) are the first generation of nanomaterials approved by the Food and Drug Administration for use as imaging agents and for the treatment of iron deficiency in chronic kidney disease. However, several IONPs-based imaging agents have been withdrawn because of toxic effects and the poor understanding of the underlying mechanisms. This study aimed to evaluate IONPs toxicity and to elucidate the underlying mechanism after intravenous administration in rats. Seven-week-old rats were intravenously administered IONPs at doses of 0, 10, 30, and 90 mg/kg body weight for 14 consecutive days. Toxicity and molecular perturbations were evaluated using traditional toxicological assessment methods and proteomics approaches, respectively. The administration of 90 mg/kg IONPs induced mild toxic effects, including abnormal clinical signs, lower body weight gain, changes in serum biochemical and hematological parameters, and increased organ coefficients in the spleen, liver, heart, and kidneys. Toxicokinetics, tissue distribution, histopathological, and transmission electron microscopy analyses revealed that the spleen was the primary organ for IONPs elimination from the systemic circulation and that the macrophage lysosomes were the main organelles of IONPs accumulation after intravenous administration. We identified 197 upregulated and 75 downregulated proteins in the spleen following IONPs administration by proteomics. Mechanically, the AKT/mTOR/TFEB signaling pathway facilitated autophagy and lysosomal activation in splenic macrophages. This is the first study to elucidate the mechanism of IONPs toxicity by combining proteomics with traditional methods for toxicity assessment.

Impacts of reference selection on the assembly of suspicious coronavirus genome.

The pandemic caused by SARS-CoV-2 has had a significant impact on the whole world. In a theory of the origin of SARS-CoV-2, pangolins are considered as a potential intermediate host. To assemble the genome of suspicious coronavirus (CoV) found in pangolins, SARS-CoV-2 was used as a reference in most of the previous studies, implicitly assuming the pangolin CoV and SARS-CoV-2 are the closest neighbors in evolution. However, this assumption may not be true. We investigated how the choice of reference genome affected the resulting CoV genome assembly. We explored various representative CoVs as the reference genome, and found significant differences in the resulting assemblies. The assembly obtained using RaTG13 as a reference showed better statistics in total length, N50, and pairwise distance reconstruction (PDR) scores than the assembly guided by SARS-CoV-2, indicating that RaTG13 may be a better reference. Therefore, RaTG13 should also be considered as a reference for assembling suspicious CoV found in pangolins and other potential intermediate hosts.

A genome-wide genetic linkage map and reference quality genome sequence for a new race in the wheat pathogen Pyrenophora tritici-repentis.

Pyrenophora tritici-repentis is an ascomycete fungus that causes tan spot of wheat. The disease has a worldwide distribution and can cause significant yield and quality losses in wheat production. The fungal pathogen is homothallic in nature, which means it can undergo sexual reproduction by selfing to produce pseudothecia on wheat stubble for seasonal survival. Since homothallism precludes the development of bi-parental fungal populations, no genetic linkage map has been developed for P. tritici-repentis for mapping and map-based cloning of fungal virulence genes. In this work, we created two heterothallic strains by deleting one of the mating type genes in each of two parental isolates 86-124 (race 2) and AR CrossB10 (a new race) and developed a bi-parental fungal population between them. The draft genome sequences of the two parental isolates were aligned to the Pt-1C-BFP reference sequence to mine single nucleotide polymorphisms (SNPs). A total of 225 SNP markers were developed for genotyping the entire population. Additionally, 75 simple sequence repeat, and two gene markers were also developed and used in the genotyping. The resulting linkage map consisted of 13 linkage groups spanning 5,075.83 cM in genetic distance. Because the parental isolate AR CrossB10 is a new race and produces Ptr ToxC, it was sequenced using long-read sequencing platforms and de novo assembled into contigs. The majority of the contigs were further anchored into chromosomes with the aid of the linkage maps. The whole genome comparison of AR CrossB10 to the reference genome of M4 revealed a few chromosomal rearrangements. The genetic linkage map and the new AR CrossB10 genome sequence are valuable tools for gene cloning in P. tritici-repentis.

Dissection and physical mapping of wheat chromosome 7B by inducing meiotic recombination with its homoeologues in Aegilops speltoides and Thinopyrum elongatum.

KEY MESSAGE: We constructed a homoeologous recombination-based bin map of wheat chromosome 7B, providing a unique physical framework for further study of chromosome 7B and its homoeologues in wheat and its relatives. Homoeologous recombination leads to the dissection and diversification of the wheat genome. Advances in genome sequencing and genotyping have dramatically improved the efficacy and throughput of homoeologous recombination-based genome studies and alien introgression in wheat and its relatives. In this study, we aimed to physically dissect and map wheat chromosome 7B by inducing meiotic recombination of chromosome 7B with its homoeologues 7E in Thinopyrum elongatum and 7S in Aegilops speltoides. The special genotypes, which were double monosomic for chromosomes 7B' + 7E' or 7B' + 7S' and homozygous for the ph1b mutant, were produced to enhance 7B - 7E and 7B - 7S recombination. Chromosome-specific DNA markers were developed and used to pre-screen the large recombination populations for 7B - 7E and 7B - 7S recombinants. The DNA marker-mediated preselections were verified by fluorescent genomic in situ hybridization (GISH). In total, 29 7B - 7E and 61 7B - 7S recombinants and multiple chromosome aberrations were recovered and delineated by GISH and the wheat 90 K SNP assay. Integrated GISH and SNP analysis of the recombinants physically mapped the recombination breakpoints and partitioned wheat chromosome 7B into 44 bins with 523 SNPs assigned within. A composite bin map was constructed for chromosome 7B, showing the bin size and physical distribution of SNPs. This provides a unique physical framework for further study of chromosome 7B and its homoeologues. In addition, the 7B - 7E and 7B - 7S recombinants extend the genetic variability of wheat chromosome 7B and represent useful germplasm for wheat breeding. Thereby, this genomics-enabled chromosome engineering approach facilitates wheat genome study and enriches the gene pool of wheat improvement.

Partitioning and physical mapping of wheat chromosome 3B and its homoeologue 3E in Thinopyrum elongatum by inducing homoeologous recombination.

KEY MESSAGE: We performed homoeologous recombination-based partitioning and physical mapping of wheat chromosome 3B and Th. elongatum chromosome 3E, providing a unique physical framework of this homoeologous pair for genome studies. The wheat (Triticum aestivum, 2n = 6x = 42, AABBDD) and Thinopyrum elongatum (2n = 2x = 14, EE) genomes can be differentiated from each other by fluorescent genomic in situ hybridization (FGISH) as well as molecular markers. This has facilitated homoeologous recombination-based partitioning and engineering of their genomes for physical mapping and alien introgression. Here, we constructed a special wheat genotype, which was double monosomic for wheat chromosome 3B and Th. elongatum chromosome 3E and homozygous for the ph1b mutant, to induce 3B-3E homoeologous recombination. Totally, 81 3B-3E recombinants were recovered and detected in the primary, secondary, and tertiary homoeologous recombination cycles by FGISH. Comparing to the primary recombination, the secondary and tertiary recombination shifted toward the proximal regions due to the increase in homology between the pairing partners. The 3B-3E recombinants were genotyped by high-throughput wheat 90-K single nucleotide polymorphism (SNP) arrays and their recombination breakpoints physically mapped based on the FGISH patterns and SNP results. The 3B-3E recombination physically partitioned chromosome 3B into 38 bins, and 429 SNPs were assigned to the distinct bins. Integrative analysis of FGISH and SNP results led to the construction of a composite bin map for chromosome 3B. Additionally, we developed 22 SNP-derived semi-thermal asymmetric reverse PCR markers specific for chromosome 3E and constructed a comparative map of homoeologous chromosomes 3E, 3B, 3A, and 3D. In summary, this work provides a unique physical framework for further studies of the 3B-3E homoeologous pair and diversifies the wheat genome for wheat improvement.

Secondary formation damage of low-pressure layer during commingled production in multilayered tight gas reservoirs.

This paper experimentally investigates fluid back-flow behavior and formation damage during commingled production in multilayered tight gas reservoirs. The development of fluid back-flow in commingled tight gas reservoirs was simulated using a newly designed experimental platform. The results indicate that when there is a pressure difference between different layers during commingled production from tight gas reservoir, water produced from the high-pressure layer will invade the low-pressure layer along with gas back-flow and will accumulate in the near-wellbore area. This will lead to an increase in water saturation and a decline in permeability in the low-pressure layer and result in a significant reduction in ultimate recovery. The outcomes of these experiments demonstrate that as well as the formation damage caused by the working fluid during drilling and fracturing, "Secondary Formation Damage" also occurs during commingled production in multilayered tight gas reservoirs. This secondary formation damage mainly occurs in the near-wellbore area of low-pressure layers and is more severe with greater proximity to the wellbore. Through further experimentation to assess the factors influencing secondary formation damage, it is shown that the degree of secondary formation damage increases with decreasing original formation pressure, original water saturation, and permeability in the lower-pressure layer.

An enrichment method for mapping ambiguous reads to the reference genome for NGS analysis.

Mapping short reads to a reference genome is an essential step in many next-generation sequencing (NGS) analyses. In plants with large genomes, a large fraction of the reads can align to multiple locations of the genome with equally good alignment scores. How to map these ambiguous reads to the genome is a challenging problem with big impacts on the downstream analysis. Traditionally, the default method is to assign an ambiguous read randomly to one of the many potential locations. In this study, we explore two alternative methods that are based on the hypothesis that the possibility of an ambiguous read being generated by a location is proportional to the total number of reads produced by that location: (1) the enrichment method that assigns an ambiguous read to the location that has produced the most reads among all the potential locations, (2) the probability method that assigns an ambiguous read to a location based on a probability proportional to the number of reads the location produces. We systematically compared the performance of the proposed methods with that of the default random method. Our results showed that the enrichment method produced better results than the default random method and the probability method in the discovery of single nucleotide polymorphisms (SNPs). Not only did it produce more SNP markers, but it also produced SNP markers with better quality, which was demonstrated using multiple mainstay genomic analyses, including genome-wide association studies (GWAS), minor allele distribution, population structure, and genomic prediction.

Genome sequence and analysis of Mycobacterium tuberculosis strain SWLPK.

OBJECTIVES: Multidrug-resistant Mycobacterium tuberculosis poses a global threat, particularly in developing countries such as Pakistan. Genome sequencing, comparative genomic analysis and drug resistance gene analysis could be beneficial for understanding and monitoring M. tuberculosis disease severity in Pakistan by elucidating the biology of M. tuberculosis. METHODS: Here the draft genome of M. tuberculosis strain SWLPK was sequenced, assembled and annotated using an Illumina MiSeq system. De novo genomic assembly was conducted using Geneious Pro™ v.10. The assembled genome of strain SWLPK was annotated using the Rapid Annotation using Subsystem Technology (RAST) server, tRNAscan-SE 1.21 and RNAmmer v.1.2, which provide high-quality functional annotation. RESULTS: Mycobacterium tuberculosis strain SWLPK yielded an average read depth of 68.5-fold, which covered 97% of the genome of reference strain H37Rv. The genome contains 4305 protein-coding genes, including key drug resistance and virulence-associated genes such as type seven secretion systems. Additionally, it has a 65.6% GC content and contains 48 RNAs and 12 contigs. We determined that all proteins encoded by this strain contain conserved domains, except OxyR, which is associated with first-line antituberculosis drugs such as ethambutol, rifampicin, streptomycin, pyrazinamide and isoniazid. CONCLUSIONS: This genome sequence provides information regarding the drug resistance genes and virulence propensity of M. tuberculosis strain SWLPK. Strain SWLPK appears to be multidrug-resistant, similar to the Beijing genotype, as it clusters in the same group. These findings will pave the way for genomic characterisation, which will provide further insights into adaptation and evolution in human hosts by transcriptome studies and gene manipulation.

Molecular cytogenetic and genomic analyses reveal new insights into the origin of the wheat B genome.

KEY MESSAGE: This work pinpointed the goatgrass chromosomal segment in the wheat B genome using modern cytogenetic and genomic technologies, and provided novel insights into the origin of the wheat B genome. Wheat is a typical allopolyploid with three homoeologous subgenomes (A, B, and D). The donors of the subgenomes A and D had been identified, but not for the subgenome B. The goatgrass Aegilops speltoides (genome SS) has been controversially considered a possible candidate for the donor of the wheat B genome. However, the relationship of the Ae. speltoides S genome with the wheat B genome remains largely obscure. The present study assessed the homology of the B and S genomes using an integrative cytogenetic and genomic approach, and revealed the contribution of Ae. speltoides to the origin of the wheat B genome. We discovered noticeable homology between wheat chromosome 1B and Ae. speltoides chromosome 1S, but not between other chromosomes in the B and S genomes. An Ae. speltoides-originated segment spanning a genomic region of approximately 10.46 Mb was detected on the long arm of wheat chromosome 1B (1BL). The Ae. speltoides-originated segment on 1BL was found to co-evolve with the rest of the B genome. Evidently, Ae. speltoides had been involved in the origin of the wheat B genome, but should not be considered an exclusive donor of this genome. The wheat B genome might have a polyphyletic origin with multiple ancestors involved, including Ae. speltoides. These novel findings will facilitate genome studies in wheat and other polyploids.

Characterization of Methane Excess and Absolute Adsorption in Various Clay Nanopores from Molecular Simulation.

In this work, we use grand canonical Monte Carlo (GCMC) simulation to study methane adsorption in various clay nanopores and analyze different approaches to characterize the absolute adsorption. As an important constituent of shale, clay minerals can have significant amount of nanopores, which greatly contribute to the gas-in-place in shale. In previous works, absolute adsorption is often calculated from the excess adsorption and bulk liquid phase density of absorbate. We find that methane adsorbed phase density keeps increasing with pressure up to 80 MPa. Even with updated adsorbed phase density from GCMC, there is a significant error in absolute adsorption calculation. Thus, we propose to use the excess adsorption and adsorbed phase volume to calculate absolute adsorption and reduce the discrepancy to less than 3% at high pressure conditions. We also find that the supercritical Dubinin-Radushkevich (SDR) fitting method which is commonly used in experiments to convert the excess adsorption to absolute adsorption may not have a solid physical foundation for methane adsorption. The methane excess and absolute adsorptions per specific surface area are similar for different clay minerals in line with previous experimental data. In mesopores, the excess and absolute adsorptions per specific surface area become insensitive to pore size. Our work should provide important fundamental understandings and insights into accurate estimation of gas-in-place in shale reservoirs.

Dbf4 recruitment by forkhead transcription factors defines an upstream rate-limiting step in determining origin firing timing.

Initiation of eukaryotic chromosome replication follows a spatiotemporal program. The current model suggests that replication origins compete for a limited pool of initiation factors. However, it remains to be answered how these limiting factors are preferentially recruited to early origins. Here, we report that Dbf4 is enriched at early origins through its interaction with forkhead transcription factors Fkh1 and Fkh2. This interaction is mediated by the Dbf4 C terminus and was successfully reconstituted in vitro. An interaction-defective mutant, dbf4ΔC, phenocopies fkh alleles in terms of origin firing. Remarkably, genome-wide replication profiles reveal that the direct fusion of the DNA-binding domain (DBD) of Fkh1 to Dbf4 restores the Fkh-dependent origin firing but interferes specifically with the pericentromeric origin activation. Furthermore, Dbf4 interacts directly with Sld3 and promotes the recruitment of downstream limiting factors. These data suggest that Fkh1 targets Dbf4 to a subset of noncentromeric origins to promote early replication in a manner that is reminiscent of the recruitment of Dbf4 to pericentromeric origins by Ctf19.

A Graph Approach to Mining Biological Patterns in the Binding Interfaces.

Protein-RNA interactions play important roles in the biological systems. Searching for regular patterns in the Protein-RNA binding interfaces is important for understanding how protein and RNA recognize each other and bind to form a complex. Herein, we present a graph-mining method for discovering biological patterns in the protein-RNA interfaces. We represented known protein-RNA interfaces using graphs and then discovered graph patterns enriched in the interfaces. Comparison of the discovered graph patterns with UniProt annotations showed that the graph patterns had a significant overlap with residue sites that had been proven crucial for the RNA binding by experimental methods. Using 200 patterns as input features, a support vector machine method was able to classify protein surface patches into RNA-binding sites and non-RNA-binding sites with 84.0% accuracy and 88.9% precision. We built a simple scoring function that calculated the total number of the graph patterns that occurred in a protein-RNA interface. That scoring function was able to discriminate near-native protein-RNA complexes from docking decoys with a performance comparable with that of a state-of-the-art complex scoring function. Our work also revealed possible patterns that might be important for binding affinity.

RNA-Seq Revealed Differences in Transcriptomes between 3ADON and 15ADON Populations of Fusarium graminearum In Vitro and In Planta.

Fusarium graminearum is the major causal agent of Fusarium head blight (FHB) in barley and wheat in North America. The fungus not only causes yield loss of the crops but also produces harmful trichothecene mycotoxins [Deoxynivalenol (DON) and its derivatives-3-acetyldeoxynivalenol (3ADON) and 15-acetyldeoxynivalenol (15ADON), and nivalenol (NIV)] that contaminate grains. Previous studies showed a dramatic increase of 3ADON-producing isolates with higher aggressiveness and DON production than the 15ADON-producing isolates in North America. However, the genetic and molecular basis of differences between the two types of isolates is unclear. In this study, we compared transcriptomes of the 3ADON and 15ADON isolates in vitro (in culture media) and in planta (during infection on the susceptible wheat cultivar 'Briggs') using RNA-sequencing. The in vitro gene expression comparison identified 479 up-regulated and 801 down-regulated genes in the 3ADON isolates; the up-regulated genes were mainly involved in C-compound and carbohydrate metabolism (18.6%), polysaccharide metabolism (7.7%) or were of unknown functions (57.6%). The in planta gene expression analysis revealed that 185, 89, and 62 genes were up-regulated in the 3ADON population at 48, 96, and 144 hours after inoculation (HAI), respectively. The up-regulated genes were significantly enriched in functions for cellular import, C-compound and carbohydrate metabolism, allantoin and allantoate transport at 48 HAI, for detoxification and virulence at 96 HAI, and for metabolism of acetic acid derivatives, detoxification, and cellular import at 144 HAI. Comparative analyses of in planta versus in vitro gene expression further revealed 2,159, 1,981 and 2,095 genes up-regulated in the 3ADON isolates, and 2,415, 2,059 and 1,777 genes up-regulated in the 15ADON isolates at the three time points after inoculation. Collectively, our data provides a foundation for further understanding of molecular mechanisms involved in aggressiveness and DON production of the two chemotype isolates of F. graminearum.

Knockdown of delta-5-desaturase promotes the anti-cancer activity of dihomo-γ-linolenic acid and enhances the efficacy of chemotherapy in colon cancer cells expressing COX-2.

Cyclooxygenase (COX), commonly overexpressed in cancer cells, is a major lipid peroxidizing enzyme that metabolizes polyunsaturated fatty acids (ω-3s and ω-6s). The COX-catalyzed free radical peroxidation of arachidonic acid (ω-6) can produce deleterious metabolites (e.g. 2-series prostaglandins) that are implicated in cancer development. Thus, COX inhibition has been intensively investigated as a complementary therapeutic strategy for cancer. However, our previous study has demonstrated that a free radical-derived byproduct (8-hydroxyoctanoic acid) formed from COX-catalyzed peroxidation of dihomo-γ-linolenic acid (DGLA, the precursor of arachidonic acid) can inhibit colon cancer cell growth. We thus hypothesize that the commonly overexpressed COX in cancer (~90% of colon cancer patients) can be taken advantage to suppress cell growth by knocking down delta-5-desaturase (D5D, a key enzyme that converts DGLA to arachidonic acid). In addition, D5D knockdown along with DGLA supplement may enhance the efficacy of chemotherapeutic drugs. After knocking down D5D in HCA-7 colony 29 cells and HT-29 cells (human colon cancer cell lines with high and low COX levels, respectively), the antitumor activity of DGLA was significantly enhanced along with the formation of a threshold range (~0.5-1.0µM) of 8-hydroxyoctanoic acid. In contrast, DGLA treatment did not inhibit cell growth when D5D was not knocked down and only limited amount of 8-hydroxyoctanoic acid was formed. D5D knockdown along with DGLA treatment also enhanced the cytotoxicities of various chemotherapeutic drugs, including 5-fluorouracil, regorafenib, and irinotecan, potentially through the activation of pro-apoptotic proteins, e.g. p53 and caspase 9. For the first time, we have demonstrated that the overexpressed COX in cancer cells can be utilized in suppressing cancer cell growth. This finding may provide a new option besides COX inhibition to optimize cancer therapy. The outcome of this translational research will guide us to develop a novel ω-6-based diet-care strategy in combination with current chemotherapy for colon cancer prevention and treatment.

RNAseq reveals weed-induced PIF3-like as a candidate target to manipulate weed stress response in soybean.

Weeds reduce yield in soybeans (Glycine max) through incompletely defined mechanisms. The effects of weeds on the soybean transcriptome were evaluated in field conditions during four separate growing seasons. RNASeq data were collected from six biological samples of soybeans growing with or without weeds. Weed species and the methods to maintain weed-free controls varied between years to mitigate treatment effects, and to allow detection of general soybean weed responses. Soybean plants were not visibly nutrient- or water-stressed. We identified 55 consistently downregulated genes in weedy plots. Many of the downregulated genes were heat shock genes. Fourteen genes were consistently upregulated. Several transcription factors including a PHYTOCHROME INTERACTING FACTOR 3-like gene (PIF3) were included among the upregulated genes. Gene set enrichment analysis indicated roles for increased oxidative stress and jasmonic acid signaling responses during weed stress. The relationship of this weed-induced PIF3 gene to genes involved in shade avoidance responses in Arabidopsis provide evidence that this gene may be important in the response of soybean to weeds. These results suggest that the weed-induced PIF3 gene will be a target for manipulating weed tolerance in soybean.

Mechanism of N-Acylthiourea-mediated activation of human histone deacetylase 8 (HDAC8) at molecular and cellular levels.

We reported previously that an N-acylthiourea derivative (TM-2-51) serves as a potent and isozyme-selective activator for human histone deacetylase 8 (HDAC8). To probe the molecular mechanism of the enzyme activation, we performed a detailed account of the steady-state kinetics, thermodynamics, molecular modeling, and cell biology studies. The steady-state kinetic data revealed that TM-2-51 binds to HDAC8 at two sites in a positive cooperative manner. Isothermal titration calorimetric and molecular modeling data conformed to the two-site binding model of the enzyme-activator complex. We evaluated the efficacy of TM-2-51 on SH-SY5Y and BE(2)-C neuroblastoma cells, wherein the HDAC8 expression has been correlated with cellular malignancy. Whereas TM-2-51 selectively induced cell growth inhibition and apoptosis in SH-SY5Y cells, it showed no such effects in BE(2)-C cells, and this discriminatory feature appears to be encoded in the p53 genotype of the above cells. Our mechanistic and cellular studies on HDAC8 activation have the potential to provide insight into the development of novel anticancer drugs.

Genotype-by-sequencing of the plant-pathogenic fungi Pyrenophora teres and Sphaerulina musiva utilizing Ion Torrent sequence technology.

Genetic and genomics tools to characterize host-pathogen interactions are disproportionately directed to the host because of the focus on resistance. However, understanding the genetics of pathogen virulence is equally important and has been limited by the high cost of de novo genotyping of species with limited marker data. Non-resource-prohibitive methods that overcome the limitation of genotyping are now available through genotype-by-sequencing (GBS). The use of a two-enzyme restriction-associated DNA (RAD)-GBS method adapted for Ion Torrent sequencing technology provided robust and reproducible high-density genotyping of several fungal species. A total of 5783 and 2373 unique loci, 'sequence tags', containing 16,441 and 9992 single nucleotide polymorphisms (SNPs) were identified and characterized from natural populations of Pyrenophora teres f. maculata and Sphaerulina musiva, respectively. The data generated from the P. teres f. maculata natural population were used in association mapping analysis to map the mating-type gene to high resolution. To further validate the methodology, a biparental population of P. teres f. teres, previously used to develop a genetic map utilizing simple sequence repeat (SSR) and amplified fragment length polymorphism (AFLP) markers, was re-analysed using the SNP markers generated from this protocol. A robust genetic map containing 1393 SNPs on 997 sequence tags spread across 15 linkage groups with anchored reference markers was generated from the P. teres f. teres biparental population. The robust high-density markers generated using this protocol will allow positional cloning in biparental fungal populations, association mapping of natural fungal populations and population genetics studies.

A graph kernel method for DNA-binding site prediction.

BACKGROUND: Protein-DNA interactions play important roles in many biological processes. Computational methods that can accurately predict DNA-binding sites on proteins will greatly expedite research on problems involving protein-DNA interactions. RESULTS: This paper presents a method for predicting DNA-binding sites on protein structures. The method represents protein surface patches using labeled graphs and uses a graph kernel method to calculate the similarities between graphs. A new surface patch is predicted to be interface or non-interface patch based on its similarities to known DNA-binding patches and non-DNA-binding patches. The proposed method achieved high accuracy when tested on a representative set of 146 protein-DNA complexes using leave-one-out cross-validation. Then, the method was applied to identify DNA-binding sites on 13 unbound structures of DNA-binding proteins. In each of the unbound structure, the top 1 patch predicted by the proposed method precisely indicated the location of the DNA-binding site. Comparisons with other methods showed that the proposed method was competitive in predicting DNA-binding sites on unbound proteins. CONCLUSIONS: The proposed method uses graphs to encode the feature's distribution in the 3-dimensional (3D) space. Thus, compared with other vector-based methods, it has the advantage of taking into account the spatial distribution of features on the proteins. Using an efficient kernel method to compare graphs the proposed method also avoids the demanding computations required for 3D objects comparison. It provides a competitive method for predicting DNA-binding sites without requiring structure alignment.