Whole-genome Transcriptome Analysis Research Articles

Role of the TLR4-androgen receptor axis and genistein in taxol-resistant ovarian cancer cells

Toll-like receptor 4 (TLR4) is often overexpressed in taxol-resistant cancer cells. Here we used whole-genome transcriptomic analysis to identify 787 upregulated genes in SKOV3 ovarian carcinoma cells that ectopically express TLR4. Using chromatin immunoprecipitation enrichment analysis, we observed that 27.8% of the TLR4-upregulated genes identified were androgen receptor (AR)-regulated genes. Accordingly, AR expression was induced in taxol-resistant SKOV3 cells overexpressing TLR4, whereas depletion of TLR4 by shRNA repressed AR expression. Activation of AR by androgens or silencing of AR using shRNA also regulated expression of AR-related genes. We found that expression of DCDC2, ANKRD18B, ALDH1A1, c14orf105, ITGBL1 and NEB was overexpressed in taxol-resistant cells, suggesting the involvement of these AR-related genes in taxol resistance. Pathway enrichment analysis confirmed that the expression of several upregulated genes enriched in steroid biosynthesis pathways was inducible by androgens, supporting the results of previous studies. We also observed that genistein inhibits AR activation, leading to suppression of AR-driven genes and reduced taxol resistance in ovarian cancer cells. Overall, we identified six TLR4- and AR-regulated genes involved in taxol resistance. Our results reveal that the TLR4/AR axis plays a critical role in taxol resistance and that genistein is a candidate compound to limit chemoresistance and improve cancer treatment in ovarian cancer.

New Insights From Imputed Whole-Genome Sequence-Based Genome-Wide Association Analysis and Transcriptome Analysis: The Genetic Mechanisms Underlying Residual Feed Intake in Chickens.

Poultry feed constitutes the largest cost in poultry production, estimated to be up to 70% of the total cost. Moreover, there is pressure on the poultry industry to increase production to meet the protein demand of humans and simultaneously reduce emissions to protect the environment. Therefore, improving feed efficiency plays an important role to improve profits and the environmental footprint in broiler production. In this study, using imputed whole-genome sequencing data, genome-wide association analysis (GWAS) was performed to identify single-nucleotide polymorphisms (SNPs) and genes associated with residual feed intake (RFI) and its component traits. Furthermore, a transcriptomic analysis between the high-RFI and the low-RFI groups was performed to validate the candidate genes from GWAS. The results showed that the heritability estimates of average daily gain (ADG), average daily feed intake (ADFI), and RFI were 0.29 (0.004), 0.37 (0.005), and 0.38 (0.004), respectively. Using imputed sequence-based GWAS, we identified seven significant SNPs and five candidate genes [MTSS I-BAR domain containing 1, folliculin, COP9 signalosome subunit 3, 5′,3′-nucleotidase (mitochondrial), and gametocyte-specific factor 1] associated with RFI, 20 significant SNPs and one candidate gene (inositol polyphosphate multikinase) associated with ADG, and one significant SNP and one candidate gene (coatomer protein complex subunit alpha) associated with ADFI. After performing a transcriptomic analysis between the high-RFI and the low-RFI groups, both 38 up-regulated and 26 down-regulated genes were identified in the high-RFI group. Furthermore, integrating regional conditional GWAS and transcriptome analysis, ras-related dexamethasone induced 1 was the only overlapped gene associated with RFI, which also suggested that the region (GGA14: 4767015–4882318) is a new quantitative trait locus associated with RFI. In conclusion, using imputed sequence-based GWAS is an efficient method to identify significant SNPs and candidate genes in chicken. Our results provide valuable insights into the genetic mechanisms of RFI and its component traits, which would further improve the genetic gain of feed efficiency rapidly and cost-effectively in the context of marker-assisted breeding selection.

MIXnorm: normalizing RNA-seq data from formalin-fixed paraffin-embedded samples.

Recent studies have shown that RNA-sequencing (RNA-seq) can be used to measure mRNA of sufficient quality extracted from formalin-fixed paraffin-embedded (FFPE) tissues to provide whole-genome transcriptome analysis. However, little attention has been given to the normalization of FFPE RNA-seq data, a key step that adjusts for unwanted biological and technical effects that can bias the signal of interest. Existing methods, developed based on fresh-frozen or similar-type samples, may cause suboptimal performance. We proposed a new normalization method, labeled MIXnorm, for FFPE RNA-seq data. MIXnorm relies on a two-component mixture model, which models non-expressed genes by zero-inflated Poisson distributions and models expressed genes by truncated normal distributions. To obtain maximum likelihood estimates, we developed a nested EM algorithm, in which closed-form updates are available in each iteration. By eliminating the need for numerical optimization in the M-step, the algorithm is easy to implement and computationally efficient. We evaluated MIXnorm through simulations and cancer studies. MIXnorm makes a significant improvement over commonly used methods for RNA-seq expression data. R code available at https://github.com/S-YIN/MIXnorm. swang@smu.edu. Supplementary data are available at Bioinformatics online.

Whole-genome sequencing and comparative transcriptome analysis of Bombyx mori nucleopolyhedrovirus La strain.

The Bombyx mori nucleopolyhedrovirus (BmNPV) La is a variant BmNPV strain isolated in Laos. La has different features from BmNPV type strain T3 in virulence, production of the polyhedrin protein, and the formation of multicapsid occlusion-derived viruses. Here, the whole-genome sequence of La was compared to the sequences of nine BmNPV and two Bombyx mandarina nucleopolyhedrovirus strains. The complete La genome consisted of 127,618 base pairs with a G + C content of 40.3% and contained putative 136 open reading frames encoding more than 60 amino acids. The La genome lacked the bro-b gene and had the highest identity with that of the T3 strain. A comparison of the transcriptomes of La- and T3-infected cells showed that the expression levels of the polyhedrin and cathepsin genes were greater in cells infected with La as compared to those infected with T3. Interestingly, the virus genes with different RNA levels between the two BmNPV strains were assembled into five clusters in the genome of La. Also, the RNA levels of host ribosomal protein genes were significantly decreased in cells infected with La as compared to those infected with T3.

Characterization of immune response against Mycobacterium marinum infection in the main hematopoietic organ of adult zebrafish (Danio rerio)

Tuberculosis remains a major global health challenge. To gain information about genes important for defense against tuberculosis, we used a well-established tuberculosis model; Mycobacterium marinum infection in adult zebrafish. To characterize the immunological response to mycobacterial infection at 14 days post infection, we performed a whole-genome level transcriptome analysis using cells from kidney, the main hematopoietic organ of adult zebrafish. Among the upregulated genes, those associated with immune signaling and regulation formed the largest category, whereas the largest group of downregulated genes had a metabolic role. We also performed a forward genetic screen in adult zebrafish and identified a fish line with severely impaired survival during chronic mycobacterial infection. Based on transcriptome analysis, these fish have decreased expression of several immunological genes. Taken together, these results give new information about the genes involved in the defense against mycobacterial infection in zebrafish.

Genetic Signatures of Drug Response Variability in Drosophila melanogaster.

Knowledge of the genetic basis underlying variation in response to environmental exposures or treatments is important in many research areas. For example, knowing the set of causal genetic variants for drug responses could revolutionize personalized medicine. We used Drosophila melanogaster to investigate the genetic signature underlying behavioral variability in response to methylphenidate (MPH), a drug used in the treatment of attention-deficit/hyperactivity disorder. We exposed a wild-type D. melanogaster population to MPH and a control treatment, and observed an increase in locomotor activity in MPH-exposed individuals. Whole-genome transcriptomic analyses revealed that the behavioral response to MPH was associated with abundant gene expression alterations. To confirm these patterns in a different genetic background and to further advance knowledge on the genetic signature of drug response variability, we used a system of inbred lines, the Drosophila Genetic Reference Panel (DGRP). Based on the DGRP, we showed that the behavioral response to MPH was strongly genotype-dependent. Using an integrative genomic approach, we incorporated known gene interactions into the genomic analyses of the DGRP, and identified putative candidate genes for variability in drug response. We successfully validated 71% of the investigated candidate genes by gene expression knockdown. Furthermore, we showed that MPH has cross-generational behavioral and transcriptomic effects. Our findings establish a foundation for understanding the genetic mechanisms driving genotype-specific responses to medical treatment, and highlight the opportunities that integrative genomic approaches have in optimizing medical treatment of complex diseases.

PECTIN ACETYLESTERASE9 Affects the Transcriptome and Metabolome and Delays Aphid Feeding.

The plant cell wall plays an important role in damage-associated molecular pattern-induced resistance to pathogens and herbivorous insects. Our current understanding of cell wall-mediated resistance is largely based on the degree of pectin methylesterification. However, little is known about the role of pectin acetylesterification in plant immunity. This study describes how one pectin-modifying enzyme, PECTIN ACETYLESTERASE 9 (PAE9), affects the Arabidopsis (Arabidopsis thaliana) transcriptome, secondary metabolome, and aphid performance. Electro-penetration graphs showed that Myzus persicae aphids established phloem feeding earlier on pae9 mutants. Whole-genome transcriptome analysis revealed a set of 56 differentially expressed genes (DEGs) between uninfested pae9-2 mutants and wild-type plants. The majority of the DEGs were enriched for biotic stress responses and down-regulated in the pae9-2 mutant, including PAD3 and IGMT2, involved in camalexin and indole glucosinolate biosynthesis, respectively. Relative quantification of more than 100 secondary metabolites revealed decreased levels of several compounds, including camalexin and oxylipins, in two independent pae9 mutants. In addition, absolute quantification of phytohormones showed that jasmonic acid (JA), jasmonoyl-Ile, salicylic acid, abscisic acid, and indole-3-acetic acid were compromised due to PAE9 loss of function. After aphid infestation, however, pae9 mutants increased their levels of camalexin, glucosinolates, and JA, and no long-term effects were observed on aphid fitness. Overall, these data show that PAE9 is required for constitutive up-regulation of defense-related compounds, but that it is not required for aphid-induced defenses. The signatures of phenolic antioxidants, phytoprostanes, and oxidative stress-related transcripts indicate that the processes underlying PAE9 activity involve oxidation-reduction reactions.

Arabidopsis PP6 phosphatases dephosphorylate PIF proteins to repress photomorphogenesis

The PHYTOCHROME-INTERACTING FACTORs (PIFs) play a central role in repressing photomorphogenesis, and phosphorylation mediates the stability of PIF proteins. Although the kinases responsible for PIF phosphorylation have been extensively studied, the phosphatases that dephosphorylate PIFs remain largely unknown. Here, we report that seedlings with mutations in FyPP1 and FyPP3, 2 genes encoding the catalytic subunits of protein phosphatase 6 (PP6), exhibited short hypocotyls and opened cotyledons in the dark, which resembled the photomorphogenic development of dark-grown pifq mutants. The hypocotyls of dark-grown sextuple mutant fypp1 fypp3 (f1 f3) pifq were shorter than those of parental mutants f1 f3 and pifq, indicating that PP6 phosphatases and PIFs function synergistically to repress photomorphogenesis in the dark. We showed that FyPPs directly interacted with PIF3 and PIF4, and PIF3 and PIF4 proteins exhibited mobility shifts in f1 f3 mutants, consistent with their hyperphosphorylation. Moreover, PIF4 was more rapidly degraded in f1 f3 mutants than in wild type after light exposure. Whole-genome transcriptomic analyses indicated that PP6 and PIFs coregulated many genes, and PP6 proteins may positively regulate PIF transcriptional activity. These data suggest that PP6 phosphatases may repress photomorphogenesis by controlling the stability and transcriptional activity of PIF proteins via regulating PIF phosphorylation.

Sesterterpene MHO7 suppresses breast cancer cells as a novel estrogen receptor degrader

Breast cancer, the most prevalent cancer in women, remains the second in the list of cancer mortality, the majority of these fatalities resulted from estrogen receptor alpha (ERα) positive disease. ERα is well known for its function on breast cancer initiation and development and has become the most successful biomarker in breast cancers. Ophiobolins are sesterterpene compounds with a distinct tricyclic 5-8-5 ring and have presented anti-cancer activities. MHO7(6-epi-ophiobolin G)was isolated from products of a mangrove fungus in our previous research and demonstrated robust activity against breast cancer cells (BCCs). The investigation on the precise mechanism of MHO7 shows that MHO7 acts as a novel ERα down regulator different from the known molecules in ER + breast cancer cells. A whole-genome transcriptomic analysis on MCF-7 cells treated with MHO7 revealed the estrogen signaling pathway was the most affected pathway, and further evidence showed the de novo synthesis of ESR1 mRNA was inhibited. In addition, MHO7 down-regulated ERα at the protein level through multiple approaches. It not only bound to ERα, pushing helix 11 away in the agonist conformation but also increased the ERα degradation through the ubiquitin-proteasome system. These effects consequently caused decreasing of the transcriptional activity of ER modulation which was confirmed by the decreasing of estrogen receptor element (ERE) activity as well as downstream genes expressions like GREB1, BRCA1, MUC1 and CCND1. Combination of tamoxifen and MHO7 yield a synergistic effect on the inhibition of MCF-7 cells when treated around the IC50 values. Our results suggest that MHO7 is a very promising drug candidate and provides a novel drug version on ERα down-regulation to fighting with breast cancer.

Disruption of the NlpD lipoprotein of the plague pathogen Yersinia pestis affects iron acquisition and the activity of the twin-arginine translocation system.

We have previously shown that the cell morphogenesis NlpD lipoprotein is essential for virulence of the plague bacteria, Yersinia pestis. To elucidate the role of NlpD in Y. pestis pathogenicity, we conducted a whole-genome comparative transcriptome analysis of the wild-type Y. pestis strain and an nlpD mutant under conditions mimicking early stages of infection. The analysis suggested that NlpD is involved in three phenomena: (i) Envelope stability/integrity evidenced by compensatory up-regulation of the Cpx and Psp membrane stress-response systems in the mutant; (ii) iron acquisition, supported by modulation of iron metabolism genes and by limited growth in iron-deprived medium; (iii) activity of the twin-arginine (Tat) system, which translocates folded proteins across the cytoplasmic membrane. Virulence studies of Y. pestis strains mutated in individual Tat components clearly indicated that the Tat system is central in Y. pestis pathogenicity and substantiated the assumption that NlpD essentiality in iron utilization involves the activity of the Tat system. This study reveals a new role for NlpD in Tat system activity and iron assimilation suggesting a modality by which this lipoprotein is involved in Y. pestis pathogenesis.

Physiological and transcriptomic analyses of roots from Malus sieversii under drought stress

Water deficit is one of the main limiting factors for apple growth and production. Root architecture plays an important role in drought tolerance of plants. However, little is known about the molecular basis of root system in apple trees under drought. In this study, we compared root morphology of two widely used rootstocks of apple (R3 and Malus sieversii) under drought. Our results suggested that M. sieversii is more tolerant to drought than R3, since M. sieversii had a higher ratio of root to shoot as well as root hydraulic conductivity under long-term drought conditions. We then performed whole-genome transcriptomic analysis to figure out the molecular basis of root responses in M. sieversii under drought. It was found that genes involved in transcription regulation, signaling or biosynthesis of hormones, and oxidative stress were differentially expressed under drought. Consistent with the gene expression profile, roots of M. sieversii had higher activities of peroxidase (POD) and superoxide dismutase (SOD) under drought, as well as higher content of abscisic acid (ABA) and lower content of auxin. Taken together, our results revealed the physiological and transcriptomic analyses of M. sieversii roots in response to drought.

IL-17–high asthma with features of a psoriasis immunophenotype

The role of IL-17 immunity is well established in patients with inflammatory diseases, such as psoriasis and inflammatory bowel disease, but not in asthmatic patients, in whom further study is required. We sought to undertake a deep phenotyping study of asthmatic patients with upregulated IL-17 immunity. Whole-genome transcriptomic analysis was performed by using epithelial brushings, bronchial biopsy specimens (91 asthmatic patients and 46 healthy control subjects), and whole blood samples (n=498) from the Unbiased Biomarkers for the Prediction of Respiratory Disease Outcomes (U-BIOPRED) cohort. Gene signatures induced invitro by IL-17 and IL-13 in bronchial epithelial cells were used to identify patients with IL-17-high and IL-13-high asthma phenotypes. Twenty-two of 91 patients were identified with IL-17, and 9 patients were identified with IL-13 gene signatures. The patients with IL-17-high asthma were characterized by risk of frequent exacerbations, airway (sputum and mucosal) neutrophilia, decreased lung microbiota diversity, and urinary biomarker evidence of activation of the thromboxane B2 pathway. In pathway analysis the differentially expressed genes in patients with IL-17-high asthma were shared with those reported as altered in psoriasis lesions and included genes regulating epithelial barrier function and defense mechanisms, such as IL1B, IL6, IL8, and β-defensin. The IL-17-high asthma phenotype, characterized by bronchial epithelial dysfunction and upregulated antimicrobial and inflammatory response, resembles the immunophenotype of psoriasis, including activation of the thromboxane B2 pathway, which should be considered a biomarker for this phenotype in further studies, including clinical trials targeting IL-17.

Microglial nodules provide the environment for pathogenic T cells in human encephalitis

Microglia nodule formation is a common feature in inflammatory brain diseases mediated by T lymphocytes such as viral and paraneoplastic encephalitis, multiple sclerosis, and Rasmussen encephalitis (RE). However, its role has not been fully understood yet. We hypothesized that, in RE, microglial nodules provide an environment for the initiation of the later dominating T-cell cytotoxicity. In RE stage 0, small primary microglia nodules could be identified in the absence of T cells. These primary nodules showed inflammasome activation and endosomal Toll-like receptor upregulation. In stage 1, T cells migrate into the parenchyma and intermingle with microglial cells, thereby forming secondary nodules in which neurons are destroyed. Whole-genome transcriptome analysis at this point showed upregulation of several inflammatory pathways including interferon signaling and major histocompatibility complex-I signaling. Inflammatory profiles, like the ones observed in RE, could be induced upon TLR3 stimulation in neonatal microglial cell cultures. Taken together, our results point towards activation of endosomal TLRs, resulting in increased interferon signaling, inflammasome activation, and chemokine upregulation as early steps in RE pathogenesis. This activity sets the scene for subsequent infiltration of T cells and destruction of neurons. Similar to RE, this microglial microenvironment might be a crucial step in other T-cell-mediated inflammatory brain diseases.

Tolerance of Transplastomic Tobacco Plants Overexpressing a Theta Class Glutathione Transferase to Abiotic and Oxidative Stresses

Chloroplasts are organelles subjected to extreme oxidative stress conditions. Biomolecules produced in the chloroplasts act as signals guiding plant metabolism toward stress tolerance and play a major role in regulating gene expression in the nucleus. Herein, we used transplastomic plants as an alternative approach to expression of transgenes in the nucleus for conferring stress tolerance to abiotic stresses and herbicides. To investigate the morphophysiological and molecular mechanisms and the role of plastid expressed GSTs in tobacco stress detoxification and stress tolerance, we used transplastomic tobacco lines overexpressing a theta class glutathione transferase (GST) in chloroplasts. The transplastomic plants were tested under drought (0, 100, and 200 mM mannitol) and salinity (0, 150, and 300 mM NaCl) in vitro, and under herbicide stress (Diquat). Our results suggest that ptAtGSTT lines were tolerant to herbicide-induced oxidative and salinity stresses and showed enhanced response tolerance to mannitol-induced osmotic stress compared to WT plants. Overexpression of the Arabidopsis thaliana AtGSTT in the chloroplasts resulted in enhanced photo-tolerance and turgor maintenance under stress. Whole-genome transcriptome analysis revealed that genes related to stress tolerance, were upregulated in ptAtGSTT2a line under both control and high mannitol stress conditions. Transplastomic plants overexpressing the ptAtGSTT2a in the chloroplast showed a state of acclimation to stress, as only limited number of genes were upregulated in the ptAtGSTT2a transplastomic line compared to WT under stress conditions while at the same time genes related to stress tolerance were upregulated in ptAtGSTT2a plants compared to WT in stress-free conditions. In parallel, the metabolic profile indicated limited perturbations of the metabolic homeostasis in the transplastomic lines and greater accumulation of mannitol, and soluble sugars under high mannitol stress. Therefore, transplastomic lines seem to be in a state of acclimation to stress under stress-free conditions, which was maintained even under high mannitol stress. The results help to elucidate the role of GSTs in plant abiotic stress tolerance and the underlying mechanisms of the GSTs expressed in the chloroplast, toward environmental resilience of cultivated crops.

MicroRNAs and their targets in cucumber shoot apices in response to temperature and photoperiod

BackgroundThe cucumber is one of the most important vegetables worldwide and is used as a research model for study of phloem transport, sex determination and temperature-photoperiod physiology. The shoot apex is the most important plant tissue in which the cell fate and organ meristems have been determined. In this study, a series of whole-genome small RNA, degradome and transcriptome analyses were performed on cucumber shoot apical tissues treated with high vs. low temperature and long vs. short photoperiod.ResultsA total of 164 known miRNAs derived from 68 families and 203 novel miRNAs from 182 families were identified. Their 4611 targets were predicted using psRobot and TargetFinder, amongst which 349 were validated by degradome sequencing. Fourteen targets of six miRNAs were differentially expressed between the treatments. A total of eight known and 16 novel miRNAs were affected by temperature and photoperiod. Functional annotations revealed that “Plant hormone signal transduction” pathway was significantly over-represented in the miRNA targets. The miR156/157/SBP-Boxes and novel-mir153/ethylene-responsive transcription factor/senescence-related protein/aminotransferase/acyl-CoA thioesterase are the two most credible miRNA/targets combinations modulating the plant’s responsive processes to the temperature-photoperiod changes. Moreover, the newly evolved, cucumber-specific novel miRNA (novel-mir153) was found to target 2087 mRNAs by prediction and has 232 targets proven by degradome analysis, accounting for 45.26–58.88% of the total miRNA targets in this plant. This is the largest sum of genes targeted by a single miRNA to the best of our knowledge.ConclusionsThese results contribute to a better understanding of the miRNAs mediating plant adaptation to combinations of temperature and photoperiod and sheds light on the recent evolution of new miRNAs in cucumber.

Tandem 3' UTR Patterns and Gene Expression Profiles of Marc-145 Cells During PRRSV Infection.

Porcine reproductive and respiratory syndrome virus (PRRSV) causes substantial economic losses to the global pig industry. Alternative polyadenylation (APA) is a mechanism that diversifies gene expression, which is important for tumorigenesis, development, and cell differentiation. However, it is unclear whether APA plays a role in the course of PRRSV infection. To address this issue, in this study we carried out a whole-genome transcriptome analysis of PRRSV-infected Marc-145 African green monkey kidney cells and identified 185 APA switching genes and 393 differentially expressed genes (DEGs). Most of these genes were involved in cellular process, metabolism, and biological regulation, and there was some overlap between the two gene sets. DEGs were found to be more directly involved in the antiviral response than APA genes. These findings provide insight into the dynamics of host gene regulation during PRRSV infection and a basis for elucidating the pathogenesis of PRRSV.

A whole-genome transcriptome analysis of articular chondrocytes in secondary osteoarthritis of the hip.

ObjectiveTo date, exhaustive gene expression analyses of chondrocytes in hip osteoarthritis (OA) have yielded specific gene expression patterns. No study has reported on the exhaustive transcriptome of secondary hip OA based on acetabular dysplasia in a Japanese population, while previous reports have focused on primary or idiopathic hip OA in Caucasian populations. This study aims to search for specific gene expression patterns of secondary hip OA chondrocytes by transcriptome analysis.DesignHuman articular cartilage was obtained from femoral heads following hemiarthroplasty for femoral neck fracture (N = 8; non-OA) and total hip arthroplasty for secondary hip OA (N = 12). Total RNA was extracted from the articular cartilage and submitted for microarray analysis. The obtained data were used to perform gene expression analysis, GO enrichment analysis and pathway analysis and were compared with data from primary hip OA in Caucasian populations in the literature.ResultsWe identified 888 upregulated (fold change: FC ≥ 2) and 732 downregulated (FC ≤ 0.5) genes in hip OA versus non-OA chondrocytes, respectively. Only 10% of upregulated genes were common between the secondary and primary OA. The newly found genes prominently overexpressed in the secondary hip OA chondrocytes were DPT, IGFBP7, and KLF2. Pathway analysis revealed extracellular matrix (ECM)-receptor interaction as an OA-related pathway, which was similar to previous reports in primary hip OA.ConclusionsThis is the first study to report the genome-wide transcriptome of secondary hip OA chondrocytes and demonstrates new potential OA-related genes. Gene expression patterns were different between secondary and primary hip OA, although the results of pathway and functional analysis were similar.

2,4-D and dicamba resistance mechanisms in wild radish: subtle, complex and population specific?

Resistance to synthetic auxin herbicides such as 2,4-dichlorophenoxyacetic acid (2,4-D) is increasing in weed populations worldwide, which is of concern given the recent introduction of synthetic auxin-resistant transgenic crops. Due to the complex mode of action of the auxinic herbicides, the mechanisms of evolved resistance remain largely uncharacterized. The aims of this study were to assess the level of diversity in resistance mechanisms in 11 populations of the problem weed Raphanus raphanistrum, and to use a high-throughput, whole-genome transcriptomic analysis on one resistant and one susceptible population to identify important changes in gene expression in response to 2,4-D. Levels of 2,4-D and dicamba (3,6-dichloro-2-methoxybenzoic acid) resistance were quantified in a dose-response study and the populations were further screened for auxin selectivity, 2,4-D translocation and metabolism, expression of key 2,4-D-responsive genes and activation of the mitogen-activated proein kinase (MAPK) pathway. Potential links between resistance levels and mechanisms were assessed using correlation analysis. The transcriptomic study revealed early deployment of the plant defence response in the 2,4-D-treated resistant population, and there was a corresponding positive relationship between auxinic herbicide resistance and constitutive MAPK phosphorylation across all populations. Populations with shoot-wide translocation of 2,4-D had similar resistance levels to those with restricted translocation, suggesting that reduced translocation may not be as strong a resistance mechanism as originally thought. Differences in auxin selectivity between populations point to the likelihood of different resistance-conferring alterations in auxin signalling and/or perception in the different populations. 2,4-D resistance in wild radish appears to result from subtly different auxin signalling alterations in different populations, supplemented by an enhanced defence response and, in some cases, reduced 2,4-D translocation. This study highlights the dangers of applying knowledge generated from a few populations of a weed species to the species as a whole.

A Pelota-like gene regulates root development and defence responses in rice.

Pelota (Pelo) are evolutionarily conserved genes reported to be involved in ribosome rescue, cell cycle control and meiotic cell division. However, there is little known about their function in plants. The aim of this study was to elucidate the function of an ethylmethane sulphonate (EMS)-derived mutation of a Pelo-like gene in rice (named Ospelo). A dysfunctional mutant was used to characterize the function of OsPelo. Analyses of its expression and sub-cellular localization were performed. The whole-genome transcriptomic change in leaves of Ospelo was also investigated by RNA sequencing. The Ospelo mutant showed defects in root system development and spotted leaves at early seedling stages. Map-based cloning revealed that the mutation occurred in the putative Pelo gene. OsPelo was found to be expressed in various tissues throughout the plant, and the protein was located in mitochondria. Defence responses were induced in the Ospelo mutant, as shown by enhanced resistance to the bacterial pathogen Xanthomonas oryzae pv. oryzae, coupled with upregulation of three pathogenesis-related marker genes. In addition, whole-genome transcriptome analysis showed that OsPelo was significantly associated with a number of biological processes, including translation, metabolism and biotic stress response. Detailed analysis showed that activation of a number of innate immunity-related genes might be responsible for the enhanced disease resistance in the Ospelo mutant. These results demonstrate that OsPelo positively regulates root development while its loss of function enhances pathogen resistance by pre-activation of defence responses in rice.

Rapid regulatory evolution of a nonrecombining autosome linked to divergent behavioral phenotypes

In the white-throated sparrow (Zonotrichia albicollis), the second chromosome bears a striking resemblance to sex chromosomes. First, within each breeding pair of birds, one bird is homozygous for the standard arrangement of the chromosome (ZAL2/ZAL2) and its mate is heterozygous for a different version (ZAL2/ZAL2m). Second, recombination is profoundly suppressed between the two versions, leading to genetic differentiation between them. Third, the ZAL2m version is linked with phenotypic traits, such as bright plumage, high aggression, and low parental behavior, which are usually associated with males. These similarities to sex chromosomes suggest that the evolutionary mechanisms that shape sex chromosomes, in particular genetic degeneration of the heterogametic version due to the suppression of recombination, are likely important in this system as well. Here, we investigated patterns of protein sequence evolution and gene expression evolution between the ZAL2 and ZAL2m chromosomes by whole-genome sequencing and transcriptome analyses. Patterns of protein evolution exhibited only weak signals of genetic degeneration, and few genes harbored signatures of positive selection. We found substantial evidence of transcriptome evolution, such as significant expression divergence between ZAL2 and ZAL2m alleles and signatures of dosage compensation for highly expressed genes. These results suggest that, early in the evolution of heteromorphic chromosomes, gene expression divergence and dosage compensation can prevail before large-scale genetic degeneration. Our results show further that suppression of recombination between heteromorphic chromosomes can lead to the evolution of alternative (sex-like) behavioral phenotypes before substantial genetic degeneration.