Genome-wide identification of histone methylation (H3K9me2) and acetylation (H4K12ac) marks in two ecotypes of switchgrass (Panicum virgatum L.).

BACKGROUND: Histone modifications play a significant role in the regulation of transcription and various biological processes, such as development and regeneration. Though a few genomic (including DNA methylation patterns) and transcriptomic studies are currently available in switchgrass, the genome-wide distribution of histone modifications has not yet been studied to help elucidate gene regulation and its application to switchgrass improvement. RESULTS: This study provides a comprehensive epigenomic analyses of two contrasting switchgrass ecotypes, lowland (AP13) and upland (VS16), by employing chromatin immunoprecipitation sequencing (ChIP-Seq) with two histone marks (suppressive- H3K9me2 and active- H4K12ac). In this study, most of the histone binding was in non-genic regions, and the highest enrichment was seen between 0 and 2 kb regions from the transcriptional start site (TSS). Considering the economic importance and potential of switchgrass as a bioenergy crop, we focused on genes, transcription factors (TFs), and pathways that were associated with C4-photosynthesis, biomass, biofuel production, biotic stresses, and abiotic stresses. Using quantitative real-time PCR (qPCR) the relative expression of five genes selected from the phenylpropanoid-monolignol pathway showed preferential binding of acetylation marks in AP13 rather than in VS16. CONCLUSIONS: The genome-wide histone modifications reported here can be utilized in understanding the regulation of genes important in the phenylpropanoid-monolignol biosynthesis pathway, which in turn, may help understand the recalcitrance associated with conversion of biomass to biofuel, a major roadblock in utilizing lignocellulosic feedstocks.

Comparative transcriptome profiling of upland (VS16) and lowland (AP13) ecotypes of switchgrass.

KEY MESSAGE: Transcriptomes of two switchgrass genotypes representing the upland and lowland ecotypes will be key tools in switchgrass genome annotation and biotic and abiotic stress functional genomics. Switchgrass (Panicum virgatum L.) is an important bioenergy feedstock for cellulosic ethanol production. We report genome-wide transcriptome profiling of two contrasting tetraploid switchgrass genotypes, VS16 and AP13, representing the upland and lowland ecotypes, respectively. A total of 268 million Illumina short reads (50 nt) were generated, of which, 133 million were obtained in AP13 and the rest 135 million in VS16. More than 90% of these reads were mapped to the switchgrass reference genome (V1.1). We identified 6619 and 5369 differentially expressed genes in VS16 and AP13, respectively. Gene ontology and KEGG pathway analysis identified key genes that regulate important pathways including C4 photosynthesis, photorespiration and phenylpropanoid metabolism. A series of genes (33) involved in photosynthetic pathway were up-regulated in AP13 but only two genes showed higher expression in VS16. We identified three dicarboxylate transporter homologs that were highly expressed in AP13. Additionally, genes that mediate drought, heat, and salinity tolerance were also identified. Vesicular transport proteins, syntaxin and signal recognition particles were seen to be up-regulated in VS16. Analyses of selected genes involved in biosynthesis of secondary metabolites, plant-pathogen interaction, membrane transporters, heat, drought and salinity stress responses confirmed significant variation in the relative expression reflected in RNA-Seq data between VS16 and AP13 genotypes. The phenylpropanoid pathway genes identified here are potential targets for biofuel conversion.

Metatranscriptomic profiles of Eastern subterranean termites, Reticulitermes flavipes (Kollar) fed on second generation feedstocks.

BACKGROUND: Second generation lignocellulosic feedstocks are being considered as an alternative to first generation biofuels that are derived from grain starches and sugars. However, the current pre-treatment methods for second generation biofuel production are inefficient and expensive due to the recalcitrant nature of lignocellulose. In this study, we used the lower termite Reticulitermes flavipes (Kollar), as a model to identify potential pretreatment genes/enzymes specifically adapted for use against agricultural feedstocks. RESULTS: Metatranscriptomic profiling was performed on worker termite guts after feeding on corn stover (CS), soybean residue (SR), or 98% pure cellulose (paper) to identify (i) microbial community, (ii) pathway level and (iii) gene-level responses. Microbial community profiles after CS and SR feeding were different from the paper feeding profile, and protist symbiont abundance decreased significantly in termites feeding on SR and CS relative to paper. Functional profiles after CS feeding were similar to paper and SR; whereas paper and SR showed different profiles. Amino acid and carbohydrate metabolism pathways were downregulated in termites feeding on SR relative to paper and CS. Gene expression analyses showed more significant down regulation of genes after SR feeding relative to paper and CS. Stereotypical lignocellulase genes/enzymes were not differentially expressed, but rather were among the most abundant/constitutively-expressed genes. CONCLUSIONS: These results suggest that the effect of CS and SR feeding on termite gut lignocellulase composition is minimal and thus, the most abundantly expressed enzymes appear to encode the best candidate catalysts for use in saccharification of these and related second-generation feedstocks. Further, based on these findings we hypothesize that the most abundantly expressed lignocellulases, rather than those that are differentially expressed have the best potential as pretreatment enzymes for CS and SR feedstocks.

Genomics reveals heterogeneous Plasmodium falciparum transmission and population differentiation in Zambia and bordering countries.

Genomic surveillance plays a critical role in monitoring malaria transmission and understanding how the parasite adapts in response to interventions. We conducted genomic surveillance of malaria by sequencing 241 Plasmodium falciparum genomes from regions with varying levels of malaria transmission across Zambia. We found genomic evidence of high levels of within-host polygenomic infections, regardless of epidemiological characteristics, underscoring the extensive and ongoing endemic malaria transmission in the country. We identified country-level clustering of parasites from Zambia and neighboring countries, and distinct clustering of parasites from West Africa. Within Zambia, our identity by descent (IBD) relatedness analysis uncovered spatial clustering of closely related parasite pairs at the local level and rare cases of long-distance sharing. Genomic regions with large shared IBD segments and strong positive selection signatures identified genes involved in sulfadoxine-pyrimethamine and artemisinin combination therapies drug resistance, but no signature related to chloroquine resistance. Together, our findings enhance our understanding of P. falciparum transmission nationwide in Zambia and highlight the urgency of strengthening malaria control programs and surveillance of antimalarial drug resistance.

Genomics reveals heterogeneous Plasmodium falciparum transmission and selection signals in Zambia.

BACKGROUND: Genomic surveillance is crucial for monitoring malaria transmission and understanding parasite adaptation to interventions. Zambia lacks prior nationwide efforts in malaria genomic surveillance among African countries. METHODS: We conducted genomic surveillance of Plasmodium falciparum parasites from the 2018 Malaria Indicator Survey in Zambia, a nationally representative household survey of children under five years of age. We whole-genome sequenced and analyzed 241 P. falciparum genomes from regions with varying levels of malaria transmission across Zambia and estimated genetic metrics that are informative about transmission intensity, genetic relatedness between parasites, and selection. RESULTS: We provide genomic evidence of widespread within-host polygenomic infections, regardless of epidemiological characteristics, underscoring the extensive and ongoing endemic malaria transmission in Zambia. Our analysis reveals country-level clustering of parasites from Zambia and neighboring regions, with distinct separation in West Africa. Within Zambia, identity by descent (IBD) relatedness analysis uncovers local spatial clustering and rare cases of long-distance sharing of closely related parasite pairs. Genomic regions with large shared IBD segments and strong positive selection signatures implicate genes involved in sulfadoxine-pyrimethamine and artemisinin combination therapies drug resistance, but no signature related to chloroquine resistance. Furthermore, differences in selection signatures, including drug resistance loci, are observed between eastern and western Zambian parasite populations, suggesting variable transmission intensity and ongoing drug pressure. CONCLUSIONS: Our findings enhance our understanding of nationwide P. falciparum transmission in Zambia, establishing a baseline for analyzing parasite genetic metrics as they vary over time and space. These insights highlight the urgency of strengthening malaria control programs and surveillance of antimalarial drug resistance.

Malaria is caused by a parasite that is spread to humans via mosquito bites. It is a leading cause of death in children under five years old in sub-Saharan Africa. Analysis of the malaria parasite's complete set of DNA (its genome) can help us to understand transmission of the disease and how this changes in response to different strategies to control the disease. We analyzed the genomes of malaria parasites from children across Zambia. Our study revealed that 77% of children harbored multiple parasite strains, which suggests that local transmission (transmission between people within the same local area) is high. Genetic evidence for long-distance transmission was rarer. Furthermore, our findings suggest parasites are evolving in response to antimalarial drugs. Our study enhances our understanding of malaria dynamics in Zambia and may help to inform strategies for improved surveillance and control.

Phase I/II Trial of Vemurafenib in Dogs with Naturally Occurring, BRAF-mutated Urothelial Carcinoma.

BRAF-targeted therapies including vemurafenib (Zelboraf) induce dramatic cancer remission; however, drug resistance commonly emerges. The purpose was to characterize a naturally occurring canine cancer model harboring complex features of human cancer, to complement experimental models to improve BRAF-targeted therapy. A phase I/II clinical trial of vemurafenib was performed in pet dogs with naturally occurring invasive urothelial carcinoma (InvUC) harboring the canine homologue of human BRAF V600E The safety, MTD, pharmacokinetics, and antitumor activity were determined. Changes in signaling and immune gene expression were assessed by RNA sequencing and phosphoproteomic analyses of cystoscopic biopsies obtained before and during treatment, and at progression. The vemurafenib MTD was 37.5 mg/kg twice daily. Anorexia was the most common adverse event. At the MTD, partial remission occurred in 9 of 24 dogs (38%), with a median progression-free interval of 181 days (range, 53-608 days). In 18% of the dogs, new cutaneous squamous cell carcinoma and papillomas occurred, a known pharmacodynamic effect of vemurafenib in humans. Upregulation of genes in the classical and alternative MAPK-related pathways occurred in subsets of dogs at cancer progression. The most consistent transcriptomic changes were the increase in patterns of T lymphocyte infiltration during the first month of vemurafenib, and of immune failure accompanying cancer progression. In conclusion, the safety, antitumor activity, and cutaneous pharmacodynamic effects of vemurafenib, and the development of drug resistance in dogs closely mimic those reported in humans. This suggests BRAF-mutated canine InvUC offers an important complementary animal model to improve BRAF-targeted therapies in humans.

Proteomic Analysis of 3T3-L1 Adipocytes Treated with Insulin and TNF-α.

Insulin resistance is an indication of early stage Type 2 diabetes (T2D). Insulin resistant adipose tissues contain higher levels of insulin than the physiological level, as well as higher amounts of intracellular tumor necrosis factor-α (TNF-α) and other cytokines. However, the mechanism of insulin resistance remains poorly understood. To better understand the roles played by insulin and TNF-α in insulin resistance, we performed proteomic analysis of differentiated 3T3-L1 adipocytes treated with insulin (Ins), TNF-α (TNF), and both (Ins + TNF). Out of the 693 proteins identified, the abundances of 78 proteins were significantly different (p < 0.05). Carnitine parmitoyltransferase-2 (CPT2), acetyl CoA carboxylase 1 (ACCAC-1), ethylmalonyl CoA decarboxylase (ECHD1), and methylmalonyl CoA isomerase (MCEE), enzymes required for fatty acid ß-oxidation and respiratory electron transport, and ß-glucuronidase, an enzyme responsible for the breakdown of complex carbohydrates, were down-regulated in all the treatment groups, compared to the control group. In contrast, superoxide dismutase 2 (SOD2), protein disulfide isomerase (PDI), and glutathione reductase, which are the proteins responsible for cytoskeletal structure, protein folding, degradation, and oxidative stress responses, were up-regulated. This suggests higher oxidative stress in cells treated with Ins, TNF, or both. We proposed a conceptual metabolic pathway impacted by the treatments and their possible link to insulin resistance or T2D.

Genome-Wide Profiling of Histone Modifications (H3K9me2 and H4K12ac) and Gene Expression in Rust (Uromyces appendiculatus) Inoculated Common Bean (Phaseolus vulgaris L.).

Histone modifications such as methylation and acetylation play a significant role in controlling gene expression in unstressed and stressed plants. Genome-wide analysis of such stress-responsive modifications and genes in non-model crops is limited. We report the genome-wide profiling of histone methylation (H3K9me2) and acetylation (H4K12ac) in common bean (Phaseolus vulgaris L.) under rust (Uromyces appendiculatus) stress using two high-throughput approaches, chromatin immunoprecipitation sequencing (ChIP-Seq) and RNA sequencing (RNA-Seq). ChIP-Seq analysis revealed 1,235 and 556 histone methylation and acetylation responsive genes from common bean leaves treated with the rust pathogen at 0, 12 and 84 hour-after-inoculation (hai), while RNA-Seq analysis identified 145 and 1,763 genes differentially expressed between mock-inoculated and inoculated plants. The combined ChIP-Seq and RNA-Seq analyses identified some key defense responsive genes (calmodulin, cytochrome p450, chitinase, DNA Pol II, and LRR) and transcription factors (WRKY, bZIP, MYB, HSFB3, GRAS, NAC, and NMRA) in bean-rust interaction. Differential methylation and acetylation affected a large proportion of stress-responsive genes including resistant (R) proteins, detoxifying enzymes, and genes involved in ion flux and cell death. The genes identified were functionally classified using Gene Ontology (GO) and EuKaryotic Orthologous Groups (KOGs). The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis identified a putative pathway with ten key genes involved in plant-pathogen interactions. This first report of an integrated analysis of histone modifications and gene expression involved in the bean-rust interaction as reported here provides a comprehensive resource for other epigenomic regulation studies in non-model species under stress.