BackgroundIn Japan and East Asia, endemic frogs appear to be tolerant or not susceptible to chytridiomycosis, a deadly amphibian disease caused by the chytrid fungus Batrachochytridium dendrobatidis (Bd). Japanese frogs may have evolved mechanisms of immune resistance to pathogens such as Bd. This study characterizes immune genes expressed in various tissues of healthy Japanese Rana frogs.ResultsWe generated transcriptome data sets of skin, spleen and blood from three adult Japanese Ranidae frogs (Japanese brown frog Rana japonica, the montane brown frog Rana ornativentris, and Tago’s brown frog Rana tagoi tagoi) as well as whole body of R. japonica and R. ornativentris tadpoles. From this, we identified tissue- and stage-specific differentially expressed genes; in particular, the spleen was most enriched for immune-related genes. A specific immune gene, major histocompatibility complex class IIB (MHC-IIB), was further characterized due to its role in pathogen recognition. We identified a total of 33 MHC-IIB variants from the three focal species (n = 7 individuals each), which displayed evolutionary signatures related to increased MHC variation, including balancing selection. Our supertyping analyses of MHC-IIB variants from Japanese frogs and previously studied frog species identified potential physiochemical properties of MHC-II that may be important for recognizing and binding chytrid-related antigens.ConclusionsThis is one of the first studies to generate transcriptomic resources for Japanese frogs, and contributes to further understanding the immunogenetic factors associated with resistance to infectious diseases in amphibians such as chytridiomycosis. Notably, MHC-IIB supertyping analyses identified unique functional properties of specific MHC-IIB alleles that may partially contribute to Bd resistance, and such properties provide a springboard for future experimental validation.

BackgroundDrought is a lifestyle disease. Plant metabolomics has been exercised for understanding the fine-tuning of the potential pathways to surmount the adverse effects of drought stress. A broad spectrum of morphological and metabolic responses from seven Triticeae species including wild types with different drought tolerance/susceptibility level was investigated under control and water scarcity conditions.ResultsSignificant morphological parameters measured were root length, surface area, average root diameter and overall root development. Principal Component Analysis, Partial Least-Squares-Discriminant Analysis and Hierarchical Cluster Analysis were applied to the metabolomic data obtained by Gas Chromatography-Mass Spectrometry technique in order to determine the important metabolites of the drought tolerance across seven different Triticeae species. The metabolites showing significant accumulation under the drought stress were considered as the key metabolites and correlated with potential biochemical pathways, enzymes or gene locations for a better understanding of the tolerance mechanisms. In all tested species, 45 significantly active metabolites with possible roles in drought stress were identified. Twenty-one metabolites out of forty-five including sugars, amino acids, organic acids and low molecular weight compounds increased in both leaf and root samples of TR39477, IG132864 and Bolal under the drought stress, contrasting to TTD-22, Tosunbey, Ligustica and Meyeri samples. Three metabolites including succinate, aspartate and trehalose were selected for further genome analysis due to their increased levels in TR39477, IG132864, and Bolal upon drought stress treatment as well as their significant role in energy producing biochemical pathways.ConclusionThese results demonstrated that the genotypes with high drought tolerance skills, especially wild emmer wheat, have a great potential to be a genetic model system for experiments aiming to validate metabolomics–genomics networks.

BackgroundIn 2016, it is estimated that there will be 62,700 new cases of kidney cancer in the United States, and 14,240 patients will die from the disease. Because the incidence of kidney renal clear cell carcinoma (KIRC), the most common type of kidney cancer, is expected to continue to increase in the US, there is an urgent need to find effective diagnostic biomarkers for KIRC that could help earlier detection of and customized treatment strategies for the disease. Accordingly, in this study we systematically investigated KIRC’s prognostic biomarkers for survival using the reverse phase protein array (RPPA) data and the high throughput sequencing data from The Cancer Genome Atlas (TCGA).ResultsWith comprehensive data available in TCGA, we systematically screened protein expression based survival biomarkers in 10 major cancer types, among which KIRC presented many protein prognostic biomarkers of survival time. This is in agreement with a previous report that expression level changes (mRNAs, microRNA and protein) may have a better performance for prognosis of KIRC. In this study, we also identified 52 prognostic genes for KIRC, many of which are involved in cell-cycle and cancer signaling, as well as 15 tumor-stage-specific prognostic biomarkers. Notably, we found fewer prognostic biomarkers for early-stage than for late-stage KIRC. Four biomarkers (the RPPA protein IDs: FASN, ACC1, Cyclin_B1 and Rad51) were found to be prognostic for survival based on both protein and mRNA expression data.ConclusionsThrough pan-cancer screening, we found that many protein biomarkers were prognostic for patients’ survival in KIRC. Stage-specific survival biomarkers in KIRC were also identified. Our study indicated that these protein biomarkers might have potential clinical value in terms of predicting survival in KIRC patients and developing individualized treatment strategies. Importantly, we found many biomarkers in KIRC at both the mRNA expression level and the protein expression level. These biomarkers shared a significant overlap, indicating that they were technically replicable.

BackgroundFatty acid composition contributes importantly to meat quality and is essential to the nutritional value of the meat. Identification of genetic factors underlying levels of fatty acids can be used to breed for pigs with healthier meat. The aim of this study was to conduct genome-wide association studies (GWAS) to identify QTL regions affecting fatty acid composition in backfat from the pig breeds Duroc and Landrace.ResultsUsing data from the Axiom porcine 660 K array, we performed GWAS on 454 Duroc and 659 Landrace boars for fatty acid phenotypes measured by near-infrared spectroscopy (NIRS) technology (C16:0, C16:1n-7, C18:0, C18:1n-9, C18:2n-6, C18:3n-3, total saturated fatty acids, monounsaturated fatty acids and polyunsaturated fatty acids). Two QTL regions on SSC4 and SSC14 were identified in Duroc for the de novo synthesized fatty acids traits, whereas one QTL on SSC8 was detected in Landrace for C16:1n-7. The QTL region on SSC14 has been reported in previous studies and a putative causative mutation has been suggested in the promoter region of the SCD gene. Whole genome re-sequencing data was used for genotype imputation and to fine map the SSC14 QTL region in Norwegian Duroc. This effort confirms the location of the QTL on this chromosome as well as suggesting other putative candidate genes in the region. The most significant single nucleotide polymorphisms (SNPs) located on SSC14 explain between 55 and 76% of the genetic variance and between 27 and 54% of the phenotypic variance for the de novo synthesized fatty acid traits in Norwegian Duroc. For the QTL region on SSC8 in Landrace, the most significant SNP explained 19% of the genetic variance and 5% of the phenotypic variance for C16:1n-7.ConclusionsThis study confirms a major QTL affecting fatty acid composition on SSC14 in Duroc, which can be used in genetic selection to increase the level of fatty acid desaturation. The SSC14 QTL was not segregating in the Landrace population, but another QTL on SSC8 affecting C16:1n-7 was identified and might be used to increase the level of desaturation in meat products from this breed.

BackgroundThe family Brassicaceae encompasses diverse species, many of which have high scientific and economic importance. Early diversifications and phylogenetic relationships between major lineages or clades remain unclear. Here we re-investigate Brassicaceae phylogeny with complete plastomes from 51 species representing all four lineages or 5 of 6 major clades (A, B, C, E and F) as identified in earlier studies.ResultsBayesian and maximum likelihood phylogenetic analyses using a partitioned supermatrix of 77 protein coding genes resulted in nearly identical tree topologies exemplified by highly supported relationships between clades. All four lineages were well identified and interrelationships between them were resolved. The previously defined Clade C was found to be paraphyletic (the genus Megadenia formed a separate lineage), while the remaining clades were monophyletic. Clade E (lineage III) was sister to clades B + C rather than to all core Brassicaceae (clades A + B + C or lineages I + II), as suggested by a previous transcriptome study. Molecular dating based on plastome phylogeny supported the origin of major lineages or clades between late Oligocene and early Miocene, and the following radiative diversification across the family took place within a short timescale. In addition, gene losses in the plastomes occurred multiple times during the evolutionary diversification of the family.ConclusionsPlastome phylogeny illustrates the early diversification of cruciferous species. This phylogeny will facilitate our further understanding of evolution and adaptation of numerous species in the model family Brassicaceae.

BackgroundOur main purpose was to evaluate the expression of plastic and evolved genes involved in ecological speciation in the noctuid moth Spodoptera frugiperda, the fall armyworm (FAW); and to demonstrate how host plants might influence lineage differentiation in this polyphagous insect. FAW is an important pest of several crops worldwide, and it is differentiated into host plant-related strains, corn (CS) and rice strains (RS). RNA-Seq and transcriptome characterization were applied to evaluate unbiased genetic expression differences in larvae from the two strains, fed on primary (corn) and alternative (rice) host plants. We consider that genes that are differently regulated by the same FAW strain, as a response to different hosts, are “plastic”. Otherwise, differences in gene expression between the two strains fed on the same host are considered constitutive differences.ResultsIndividual performance parameters (larval and pupal weight) varied among conditions (strains vs. hosts). A total of 3657 contigs was related to plastic response, and 2395 contigs were differentially regulated in the two strains feeding on preferential and alternative hosts (constitutive contigs). Three molecular functions were present in all comparisons, both down- and up-regulated: oxidoreductase activity, metal-ion binding, and hydrolase activity.ConclusionsMetabolization of foreign chemicals is among the key functions involved in the phenotypic variation of FAW strains. From an agricultural perspective, high plasticity in families of detoxifying genes indicates the capacity for a rapid response to control compounds such as insecticides.