Enhancing grapevine breeding efficiency through genomic prediction and selection index

Grapevine (Vitis vinifera) breeding reaches a critical point. New cultivars are released every year with resistance to powdery and downy mildews. However, the traditional process remains time-consuming, taking 20-25 years, and demands the evaluation of new traits to enhance grapevine adaptation to climate change. Until now, the selection process has relied on phenotypic data and a limited number of molecular markers for simple genetic traits such as resistance to pathogens, without a clearly defined ideotype, and was carried out on a large scale. To accelerate the breeding process and address these challenges, we investigated the use of genomic prediction, a methodology using molecular markers to predict genotypic values. In our study, we focused on 2 existing grapevine breeding programs: Rosé wine and Cognac production. In these programs, several families were created through crosses of emblematic and interspecific resistant varieties to powdery and downy mildews. Thirty traits were evaluated for each program, using 2 genomic prediction methods: Genomic Best Linear Unbiased Predictor and Least Absolute Shrinkage Selection Operator. The results revealed substantial variability in predictive abilities across traits, ranging from 0 to 0.9. These discrepancies could be attributed to factors such as trait heritability and trait characteristics. Moreover, we explored the potential of across-population genomic prediction by leveraging other grapevine populations as training sets. Integrating genomic prediction allowed us to identify superior individuals for each program, using multivariate selection index method. The ideotype for each breeding program was defined collaboratively with representatives from the wine-growing sector.

Target enrichment sequencing coupled with GWAS identifies MdPRX10 as a candidate gene in the control of budbreak in apple

The timing of floral budbreak in apple has a significant effect on fruit production and quality. Budbreak occurs as a result of a complex molecular mechanism that relies on accurate integration of external environmental cues, principally temperature. In the pursuit of understanding this mechanism, especially with respect to aiding adaptation to climate change, a QTL at the top of linkage group (LG) 9 has been identified by many studies on budbreak, but the genes underlying it remain elusive. Here, together with a dessert apple core collection of 239 cultivars, we used a targeted capture sequencing approach to increase SNP resolution in apple orthologues of known or suspected A. thaliana flowering time-related genes, as well as approximately 200 genes within the LG9 QTL interval. This increased the 275 223 SNP Axiom® Apple 480 K array dataset by an additional 40 857 markers. Robust GWAS analyses identified MdPRX10, a peroxidase superfamily gene, as a strong candidate that demonstrated a dormancy-related expression pattern and down-regulation in response to chilling. In-silico analyses also predicted the residue change resulting from the SNP allele associated with late budbreak could alter protein conformation and likely function. Late budbreak cultivars homozygous for this SNP allele also showed significantly up-regulated expression of C-REPEAT BINDING FACTOR (CBF) genes, which are involved in cold tolerance and perception, compared to reference cultivars, such as Gala. Taken together, these results indicate a role for MdPRX10 in budbreak, potentially via redox-mediated signaling and CBF gene regulation. Moving forward, this provides a focus for developing our understanding of the effects of temperature on flowering time and how redox processes may influence integration of external cues in dormancy pathways.

RNAseq based variant dataset in a black poplar association panel

OBJECTIVE

Black poplar (Populus nigra L.) is a species native to Eurasia with a wide distribution area. It is an ecologically important species from riparian ecosystems, that is used as a parent of interspecific (P. deltoides x P. nigra) cultivated poplar hybrids. Variant detection from transcriptomics sequences of 241 P. nigra individuals, sampled in natural populations from 11 river catchments (in four European countries) is described here. These data provide new valuable resources for population structure analysis, population genomics and genome-wide association studies.

DATA DESCRIPTION

We generated transcriptomics data from a mixture of young differentiating xylem and cambium tissues of 480 Populus nigra trees sampled in a common garden experiment located at Orléans (France), corresponding to 241 genotypes (2 clonal replicates per genotype, at maximum) by using RNAseq technology. We launched on the resulting sequences an in-silico pipeline that allowed us to obtain 878,957 biallelic polymorphisms without missing data. More than 99% of these positions are annotated and 98.8% are located on the 19 chromosomes of the P. trichocarpa reference genome. The raw RNAseq sequences are available at the NCBI Sequence Read Archive SPR188754 and the variant dataset at the Recherche Data Gouv repository under https://doi.org/10.15454/8DQXK5 .

Genetic structure and first genome-wide insights into the adaptation of a wild relative of grapevine, Vitis berlandieri

In grafted plants, such as grapevine, increasing the diversity of rootstocks available to growers is an ideal strategy for helping plants to adapt to climate change. The rootstocks used for grapevine are hybrids of various American Vitis, including V. berlandieri. The rootstocks currently use in vineyards are derived from breeding programs involving very small numbers of parental individuals. We investigated the structure of a natural population of V. berlandieri and the association of genetic diversity with environmental variables. In this study, we collected seeds from 78 wild V. berlandieri plants in Texas after open fertilization. We genotyped 286 individuals to describe the structure of the population, and environmental information collected at the sampling site made it possible to perform genome-environment association analysis (GEA). De novo long-read whole-genome sequencing was performed on V. berlandieri and a STRUCTURE analysis was performed. We identified and filtered 104,378 SNPs. We found that there were two subpopulations associated with differences in elevation, temperature, and rainfall between sampling sites. GEA identified three QTL for elevation and 15 QTL for PCA coordinates based on environmental parameter variability. This original study is the first GEA study to be performed on a population of grapevines sampled in natural conditions. Our results shed new light on rootstock genetics and could open up possibilities for introducing greater diversity into genetic improvement programs for grapevine rootstocks.

Use of near-infrared spectroscopy to estimate physical, anatomical and hydraulic properties of Eucalyptus wood

Tree breeding programs and wood industries require simple, time- and cost-effective techniques to process large volumes of samples. In recent decades, near-infrared spectroscopy (NIRS) has been acknowledged as one of the most powerful techniques for wood analysis, making it the most used tool for high-throughput phenotyping. Previous studies have shown that a significant number of anatomical, physical, chemical and mechanical wood properties can be estimated through NIRS, both for angiosperm and gymnosperm species. However, the ability of this technique to predict functional traits related to drought resistance has been poorly explored, especially in angiosperm species. This is particularly relevant since determining xylem hydraulic properties by conventional techniques is complex and time-consuming, clearly limiting its use in studies and applications that demand large amounts of samples. In this study, we measured several wood anatomical and hydraulic traits and collected NIR spectra in branches of two Eucalyptus L'Hér species. We developed NIRS calibration models and discussed their ability to accurately predict the studied traits. The models generated allowed us to adequately calibrate the reference traits, with high R2 (≥0.75) for traits such as P12, P88, the slope of the vulnerability curves to xylem embolism or the fiber wall fraction, and with lower R2 (0.39-0.52) for P50, maximum hydraulic conductivity or frequency of ray parenchyma. We found that certain wavenumbers improve models' calibration, with those in the range of 4000-5500 cm-1 predicting the highest number of both anatomical and functional traits. We concluded that the use of NIRS allows calibrating models with potential predictive value not only for wood structural and chemical variables but also for anatomical and functional traits related to drought resistance in wood types with complex structure as eucalypts. These results are promising in light of the required knowledge about species and genotypes adaptability to global climatic change.

Regulation of CO2 by the sea in areas around Latin America in a context of climate change

Anthropogenic activities are increasing the atmospheric carbon dioxide (CO₂); around a third of the CO₂ emitted by these activities has been taken up by the ocean. Nevertheless, this marine ecosystem service of regulation remains largely invisible to society, and not enough is known about regional differences and trends in sea-air CO₂ fluxes (FCO2), especially in the Southern Hemisphere. The objectives of this work were as follows: first to put values of FCO2 integrated over the exclusive economic zones (EEZ) of five Latin-American countries (Argentina, Brazil, Mexico, Peru, and Venezuela) into perspective regarding total country-level greenhouse gases (GHG) emissions. Second, to assess the variability of two main biological factors affecting FCO2 at marine ecological time series (METS) in these areas. FCO2 over the EEZs were estimated using the NEMO model, and GHG emissions were taken from reports to the UN Framework Convention on Climate Change. For each METS, the variability in phytoplankton biomass (indexed by chlorophyll-a concentration, Chla) and abundance of different cell sizes (phy-size) were analyzed at two time periods (2000-2015 and 2007-2015). Estimates of FCO2 at the analyzed EEZs showed high variability among each other and non-negligible values in the context of greenhouse gas emissions. The trends observed at the METS indicated, in some cases, an increase in Chla (e.g., EPEA-Argentina) and a decrease in others (e.g., IMARPE-Peru). Evidence of increasing populations of small size-phytoplankton was observed (e.g., EPEA-Argentina, Ensenada-Mexico), which would affect the carbon export to the deep ocean. These results highlight the relevance of ocean health and its ecosystem service of regulation when discussing carbon net emissions and budgets.

Interest of phenomic prediction as an alternative to genomic prediction in grapevine

BACKGROUND

Phenomic prediction has been defined as an alternative to genomic prediction by using spectra instead of molecular markers. A reflectance spectrum provides information on the biochemical composition within a tissue, itself being under genetic determinism. Thus, a relationship matrix built from spectra could potentially capture genetic signal. This new methodology has been mainly applied in several annual crop species but little is known so far about its interest in perennial species. Besides, phenomic prediction has only been tested for a restricted set of traits, mainly related to yield or phenology. This study aims at applying phenomic prediction for the first time in grapevine, using spectra collected on two tissues and over two consecutive years, on two populations and for 15 traits, related to berry composition, phenology, morphological and vigour. A major novelty of this study was to collect spectra and phenotypes several years apart from each other. First, we characterized the genetic signal in spectra and under which condition it could be maximized, then phenomic predictive ability was compared to genomic predictive ability.

RESULTS

For the first time, we showed that the similarity between spectra and genomic relationship matrices was stable across tissues or years, but variable across populations, with co-inertia around 0.3 and 0.6 for diversity panel and half-diallel populations, respectively. Applying a mixed model on spectra data increased phenomic predictive ability, while using spectra collected on wood or leaves from one year or another had less impact. Differences between populations were also observed for predictive ability of phenomic prediction, with an average of 0.27 for the diversity panel and 0.35 for the half-diallel. For both populations, a significant positive correlation was found across traits between predictive ability of genomic and phenomic predictions.

CONCLUSION

NIRS is a new low-cost alternative to genotyping for predicting complex traits in perennial species such as grapevine. Having spectra and phenotypes from different years allowed us to exclude genotype-by-environment interactions and confirms that phenomic prediction can rely only on genetics.

eQTLs are key players in the integration of genomic and transcriptomic data for phenotype prediction

Background

Multi-omics represent a promising link between phenotypes and genome variation. Few studies yet address their integration to understand genetic architecture and improve predictability.

Results

Our study used 241 poplar genotypes, phenotyped in two common gardens, with xylem and cambium RNA sequenced at one site, yielding large phenotypic, genomic (SNP), and transcriptomic datasets. Prediction models for each trait were built separately for SNPs and transcripts, and compared to a third model integrated by concatenation of both omics. The advantage of integration varied across traits and, to understand such differences, an eQTL analysis was performed to characterize the interplay between the genome and transcriptome and classify the predicting features into cis or trans relationships. A strong, significant negative correlation was found between the change in predictability and the change in predictor ranking for trans eQTLs for traits evaluated in the site of transcriptomic sampling.

Conclusions

Consequently, beneficial integration happens when the redundancy of predictors is decreased, likely leaving the stage to other less prominent but complementary predictors. An additional gene ontology (GO) enrichment analysis appeared to corroborate such statistical output. To our knowledge, this is a novel finding delineating a promising method to explore data integration.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08690-7.

Tree architecture, light interception and water-use related traits are controlled by different genomic regions in an apple tree core collection

Tree architecture shows large genotypic variability, but how this affects water-deficit responses is poorly understood. To assess the possibility of reaching ideotypes with adequate combinations of architectural and functional traits in the face of climate change, we combined high-throughput field phenotyping and genome-wide association studies (GWAS) on an apple tree (Malus domestica) core-collection. We used terrestrial light detection and ranging (T-LiDAR) scanning and airborne multispectral and thermal imagery to monitor tree architecture, canopy shape, light interception, vegetation indices and transpiration on 241 apple cultivars submitted to progressive field soil drying. GWAS was performed with single nucleotide polymorphism (SNP)-by-SNP and multi-SNP methods. Large phenotypic and genetic variability was observed for all traits examined within the collection, especially canopy surface temperature in both well-watered and water deficit conditions, suggesting control of water loss was largely genotype-dependent. Robust genomic associations revealed independent genetic control for the architectural and functional traits. Screening associated genomic regions revealed candidate genes involved in relevant pathways for each trait. We show that multiple allelic combinations exist for all studied traits within this collection. This opens promising avenues to jointly optimize tree architecture, light interception and water use in breeding strategies. Genotypes carrying favourable alleles depending on environmental scenarios and production objectives could thus be targeted.

Across-population genomic prediction in grapevine opens up promising prospects for breeding

Crop breeding involves two selection steps: choosing progenitors and selecting individuals within progenies. Genomic prediction, based on genome-wide marker estimation of genetic values, could facilitate these steps. However, its potential usefulness in grapevine (Vitis vinifera L.) has only been evaluated in non-breeding contexts mainly through cross-validation within a single population. We tested across-population genomic prediction in a more realistic breeding configuration, from a diversity panel to ten bi-parental crosses connected within a half-diallel mating design. Prediction quality was evaluated over 15 traits of interest (related to yield, berry composition, phenology and vigour), for both the average genetic value of each cross (cross mean) and the genetic values of individuals within each cross (individual values). Genomic prediction in these conditions was found useful: for cross mean, average per-trait predictive ability was 0.6, while per-cross predictive ability was halved on average, but reached a maximum of 0.7. Mean predictive ability for individual values within crosses was 0.26, about half the within-half-diallel value taken as a reference. For some traits and/or crosses, these across-population predictive ability values are promising for implementing genomic selection in grapevine breeding. This study also provided key insights on variables affecting predictive ability. Per-cross predictive ability was well predicted by genetic distance between parents and when this predictive ability was below 0.6, it was improved by training set optimization. For individual values, predictive ability mostly depended on trait-related variables (magnitude of the cross effect and heritability). These results will greatly help designing grapevine breeding programs assisted by genomic prediction.

Genomic signatures of a major adaptive event in the pathogenic fungus Melampsora larici-populina

The recent availability of genome-wide sequencing techniques has allowed systematic screening for molecular signatures of adaptation, including in nonmodel organisms. Host–pathogen interactions constitute good models due to the strong selective pressures that they entail. We focused on an adaptive event which affected the poplar rust fungus Melampsora larici-populina when it overcame a resistance gene borne by its host, cultivated poplar. Based on 76 virulent and avirulent isolates framing narrowly the estimated date of the adaptive event, we examined the molecular signatures of selection. Using an array of genome scan methods based on different features of nucleotide diversity, we detected a single locus exhibiting a consistent pattern suggestive of a selective sweep in virulent individuals (excess of differentiation between virulent and avirulent samples, linkage disequilibrium, genotype–phenotype statistical association, and long-range haplotypes). Our study pinpoints a single gene and further a single amino acid replacement which may have allowed the adaptive event. Although our samples are nearly contemporary to the selective sweep, it does not seem to have affected genome diversity further than the immediate vicinity of the causal locus, which can be explained by a soft selective sweep (where selection acts on standing variation) and by the impact of recombination in mitigating the impact of selection. Therefore, it seems that properties of the life cycle of M. larici-populina, which entails both high genetic diversity and outbreeding, has facilitated its adaptation.

RNAi suppression of DNA methylation affects the drought stress response and genome integrity in transgenic poplar

Trees are long-lived organisms that continuously adapt to their environments, a process in which epigenetic mechanisms are likely to play a key role. Via downregulation of the chromatin remodeler DECREASED IN DNA METHYLATION 1 (DDM1) in poplar (Populus tremula × Populus alba) RNAi lines, we examined how DNA methylation coordinates genomic and physiological responses to moderate water deficit. We compared the growth and drought response of two RNAi-ddm1 lines to wild-type (WT) trees under well-watered and water deficit/rewatering conditions, and analyzed their methylomes, transcriptomes, mobilomes and phytohormone contents in the shoot apical meristem. The RNAi-ddm1 lines were more tolerant to drought-induced cavitation but did not differ in height or stem diameter growth. About 5000 differentially methylated regions were consistently detected in both RNAi-ddm1 lines, colocalizing with 910 genes and 89 active transposable elements. Under water deficit conditions, 136 differentially expressed genes were found, including many involved in phytohormone pathways; changes in phytohormone concentrations were also detected. Finally, the combination of hypomethylation and drought led to the mobility of two transposable elements. Our findings suggest major roles for DNA methylation in regulation of genes involved in hormone-related stress responses, and the maintenance of genome integrity through repression of transposable elements.

Singular patterns of skull shape and brain size change in the domestication of South American camelids

Patterns of selection in South American camelids (Lamini) and their unique demographic history establish the llama and alpaca as remarkable cases of domestication among large herd animals. Skull shape is implicated in many changes reported between wild and domestic taxa. We apply 3D geometric morphometric methods to describe skull shape, form, and size, differences among the four species of Lamini. In so doing, we test if domesticated Lamini exhibit changes similar to those in other domesticated groups: not only in the skull, but also in brain and body size. In contrast to other domesticated artiodactyls, very little change has occurred in domestic alpacas and llamas compared to their wild counterparts. Nevertheless, their differences are statistically significant and include a flatter cranium, inclined palate and increased airorhynchy in the domestics. Selection pressures that contrast with those on other herd animals, as well as recent population bottlenecks, likely have influenced the morphological patterns we note in Lamini. High-resolution 3D morphospace allows skull size, shape, and form (shape + size), to discriminate all four species, with form providing the greatest separation. These results help differentiate morphologically the Lamini, which in nature are distinguished mainly by body size, and provide an additional tool to archaeologists for distinction of wild and domestic remains. Most of our shape analyses suggest a marginally closer relationship between the alpaca and vicuña, to the exclusion of the guanaco, supporting the genetic relationships for this group. The expected brain size change between wild and domestic populations is lower than previously thought, with a 15.4% reduction in llama, and 6.8% reduction in alpaca. This is the lowest reduction in brain size thus far reported among domesticated Artiodactyla.

Favorable Conditions for Genomic Evaluation to Outperform Classical Pedigree Evaluation Highlighted by a Proof-of-Concept Study in Poplar

Forest trees like poplar are particular in many ways compared to other domesticated species. They have long juvenile phases, ongoing crop-wild gene flow, extensive outcrossing, and slow growth. All these particularities tend to make the conduction of breeding programs and evaluation stages costly both in time and resources. Perennials like trees are therefore good candidates for the implementation of genomic selection (GS) which is a good way to accelerate the breeding process, by unchaining selection from phenotypic evaluation without affecting precision. In this study, we tried to compare GS to pedigree-based traditional evaluation, and evaluated under which conditions genomic evaluation outperforms classical pedigree evaluation. Several conditions were evaluated as the constitution of the training population by cross-validation, the implementation of multi-trait, single trait, additive and non-additive models with different estimation methods (G-BLUP or weighted G-BLUP). Finally, the impact of the marker densification was tested through four marker density sets. The population under study corresponds to a pedigree of 24 parents and 1,011 offspring, structured into 35 full-sib families. Four evaluation batches were planted in the same location and seven traits were evaluated on 1 and 2 years old trees. The quality of prediction was reported by the accuracy, the Spearman rank correlation and prediction bias and tested with a cross-validation and an independent individual test set. Our results show that genomic evaluation performance could be comparable to the already well-optimized pedigree-based evaluation under certain conditions. Genomic evaluation appeared to be advantageous when using an independent test set and a set of less precise phenotypes. Genome-based methods showed advantages over pedigree counterparts when ranking candidates at the within-family levels, for most of the families. Our study also showed that looking at ranking criteria as Spearman rank correlation can reveal benefits to genomic selection hidden by biased predictions.

Distribution of segmental chromosomal alterations in neuroblastoma

BACKGROUND

Neuroblastoma (NB) is a heterogeneous tumor with extremely diverse prognosis according to clinical and genetic factors such as specific combinations of chromosomal imbalances.

METHODS

Molecular karyotyping data from a national neuroblastic tumor database of 155 NB samples were analyzed and related to clinical data.

RESULTS

Segmental chromosomal alterations (SCA) were detected in 102 NB, whereas 45 only displayed numerical alterations. Incidence of SCA was higher in stage M (92%) and MYCN amplified (MNA) NB (96%). Presence of SCA was associated with older age, especially 1q gain and 3p deletion. 96% of the deaths were observed in the SCA group and 85% of the relapsed NB contained SCA. The alteration most commonly associated with a higher number of other segmental rearrangements was 11q deletion, followed by 4p deletion. Whole-chromosome 19 gain was associated with lower stages, absence of SCA and better outcome.

CONCLUSIONS

SCA are not randomly distributed and are concentrated on recurrent chromosomes. The most frequently affected chromosomes identify prognostic factors in specific risk groups. SCA are associated with older age and MNA. We have identified a small subset of patients with better outcome that share whole-chromosome 19 numeric gain, suggesting its use as a prognostic biomarker in NB.

Gene expression predictions and networks in natural populations supports the omnigenic theory

Background

Recent literature on the differential role of genes within networks distinguishes core from peripheral genes. If previous works have shown contrasting features between them, whether such categorization matters for phenotype prediction remains to be studied.

Results

We measured 17 phenotypic traits for 241 cloned genotypes from a Populus nigra collection, covering growth, phenology, chemical and physical properties. We also sequenced RNA for each genotype and built co-expression networks to define core and peripheral genes. We found that cores were more differentiated between populations than peripherals while being less variable, suggesting that they have been constrained through potentially divergent selection. We also showed that while cores were overrepresented in a subset of genes statistically selected for their capacity to predict the phenotypes (by Boruta algorithm), they did not systematically predict better than peripherals or even random genes.

Conclusion

Our work is the first attempt to assess the importance of co-expression network connectivity in phenotype prediction. While highly connected core genes appear to be important, they do not bear enough information to systematically predict better quantitative traits than other gene sets.

Li-Fraumeni syndrome heterogeneity

Clinical variability is commonly seen in Li-Fraumeni syndrome. Phenotypic heterogeneity is present among different families affected by the same pathogenic variant in TP53 gene and among members of the same family. However, causes of this huge clinical spectrum have not been studied in depth. TP53 type mutation, polymorphic variants in TP53 gene or in TP53-related genes, copy number variations in particular regions, and/or epigenetic deregulation of TP53 expression might be responsible for clinical heterogeneity. In this review, recent advances in the understanding of genetic and epigenetic aspects influencing Li-Fraumeni phenotype are discussed.

Accuracy of RNAseq based SNP discovery and genotyping in Populusnigra

Backgroud

Populus nigra is a major tree species of ecological and economic importance for which several initiatives have been set up to create genomic resources. In order to access the large number of Single Nucleotide Polymorphisms (SNPs) typically needed to carry out a genome scan, the present study aimed at evaluating RNA sequencing as a tool to discover and type SNPs in genes within natural populations of P. nigra.

Results

We have devised a bioinformatics pipeline to call and type SNPs from RNAseq reads and applied it to P. nigra transcriptomic data. The accuracy of the resulting RNAseq-based SNP calling and typing has been evaluated by (i) comparing their position and alleles to those previously reported in candidate genes, (ii) assessing their genotyping accuracy with respect to a previously available SNP chip and (iii) evaluating their inter-annual repeatability. We found that a combination of several callers yields a good compromise between the number of variants type and the accuracy of genotyping. We further used the resulting genotypic data to carry out basic genetic analyses whose results confirm the quality of the RNAseq-based SNP dataset.

Conclusions

We demonstrated the potential and accuracy of RNAseq as an efficient way to genotype SNPs in P. nigra. These results open prospects towards the use of this technology for quantitative and population genomics studies.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5239-z) contains supplementary material, which is available to authorized users.

Phenomic Selection Is a Low-Cost and High-Throughput Method Based on Indirect Predictions: Proof of Concept on Wheat and Poplar

Genomic selection - the prediction of breeding values using DNA polymorphisms - is a disruptive method that has widely been adopted by animal and plant breeders to increase productivity. It was recently shown that other sources of molecular variations such as those resulting from transcripts or metabolites could be used to accurately predict complex traits. These endophenotypes have the advantage of capturing the expressed genotypes and consequently the complex regulatory networks that occur in the different layers between the genome and the phenotype. However, obtaining such omics data at very large scales, such as those typically experienced in breeding, remains challenging. As an alternative, we proposed using near-infrared spectroscopy (NIRS) as a high-throughput, low cost and non-destructive tool to indirectly capture endophenotypic variants and compute relationship matrices for predicting complex traits, and coined this new approach "phenomic selection" (PS). We tested PS on two species of economic interest (Triticum aestivum L. and Populus nigra L.) using NIRS on various tissues (grains, leaves, wood). We showed that one could reach predictions as accurate as with molecular markers, for developmental, tolerance and productivity traits, even in environments radically different from the one in which NIRS were collected. Our work constitutes a proof of concept and provides new perspectives for the breeding community, as PS is theoretically applicable to any organism at low cost and does not require any molecular information.

Winter-dormant shoot apical meristem in poplar trees shows environmental epigenetic memory

Trees have a long lifespan and must continually adapt to environmental pressures, notably in the context of climate change. Epigenetic mechanisms are doubtless involved in phenotypic plasticity and in stress memory; however, little evidence of the role of epigenetic processes is available for trees growing in fields. Here, we analyzed the possible involvement of epigenetic mechanisms in the winter-dormant shoot apical meristem of Populus × euramericana clones in memory of the growing conditions faced during the vegetative period. We aimed to estimate the range of genetic and environmentally induced variations in global DNA methylation and to evaluate their correlation with changes in biomass production, identify differentially methylated regions (DMRs), and characterize common DMRs between experiments. We showed that the variations in global DNA methylation between conditions were genotype dependent and correlated with biomass production capacity. Microarray chip analysis allowed detection of DMRs 6 months after the stressful summer period. The 161 DMRs identified as common to three independent experiments most notably targeted abiotic stress and developmental response genes. Results are consistent with a winter-dormant shoot apical meristem epigenetic memory of stressful environmental conditions that occurred during the preceding summer period. This memory may facilitate tree acclimation.

Retinoblastoma: towards an earlier diagnosis

OBJECTIVES

The number of enucleations and visual sequels due to retinoblastoma is high. The aim of this study was to evaluate the different diagnostic aspects and propose strategies that might improve the clinical management of this condition.

METHOD

A retrospective study was conducted on 38 patients with retinoblastoma studied genetically (29 unilateral, 9 bilateral). The evaluation included: age of onset, clinical signs, and time since onset, number of enucleations, time to diagnosis, and survival at 5 years.

RESULTS

Leukocoria was the main clinical sign (present in 90% of cases). The mean diagnostic delay was 3.2 months. Among the unilateral cases, the eyes were enucleated in 76%, and 55% in the bilateral forms. Only one death was found among the 25 patients followed-up for at least 5 years.

CONCLUSIONS

Retinoblastoma diagnostic and treatment strategies need to be updated. Good coordination between paediatricians and ophthalmologists is essential for this. Its management in reference centres, which have the necessary technology and experience, should contribute to increase the rate of organ preservation.

Changes in the epigenome and transcriptome of the poplar shoot apical meristem in response to water availability affect preferentially hormone pathways

The adaptive capacity of long-lived organisms such as trees to the predicted climate changes, including severe and successive drought episodes, will depend on the presence of genetic diversity and phenotypic plasticity. Here, the involvement of epigenetic mechanisms in phenotypic plasticity toward soil water availability was examined in Populus×euramericana. This work aimed at characterizing (i) the transcriptome plasticity, (ii) the genome-wide plasticity of DNA methylation, and (iii) the function of genes affected by a drought-rewatering cycle in the shoot apical meristem. Using microarray chips, differentially expressed genes (DEGs) and differentially methylated regions (DMRs) were identified for each water regime. The rewatering condition was associated with the highest variations of both gene expression and DNA methylation. Changes in methylation were observed particularly in the body of expressed genes and to a lesser extent in transposable elements. Together, DEGs and DMRs were significantly enriched in genes related to phytohormone metabolism or signaling pathways. Altogether, shoot apical meristem responses to changes in water availability involved coordinated variations in DNA methylation, as well as in gene expression, with a specific targeting of genes involved in hormone pathways, a factor that may enable phenotypic plasticity.

Identification of major loci and genomic regions controlling acid and volatile content in tomato fruit: implications for flavor improvement

Plant metabolites are important to world food security due to their roles in crop yield and nutritional quality. Here we report the metabolic profile of 300 tomato accessions (Solanum lycopersicum and related wild species) by quantifying 60 primary and secondary metabolites, including volatile organic compounds, over a period of 2 yr. Metabolite content and genetic inheritance of metabolites varied broadly, both within and between different genetic groups. Using genotype information gained from 10 000 single nucleotide polymorphism markers, we performed a metabolite genome-wide association mapping (GWAS) study. We identified 79 associations influencing 13 primary and 19 secondary metabolites with large effects at high resolution. Four genome regions were detected, highlighting clusters of associations controlling the variation of several metabolites. Local linkage disequilibrium analysis and allele mining identified possible candidate genes which may modulate the content of metabolites that are of significant importance for human diet and fruit consumption. We precisely characterized two associations involved in fruit acidity and phenylpropanoid volatile production. Taken together, this study reveals complex and distinct metabolite regulation in tomato subspecies and demonstrates that GWAS is a powerful tool for gene-metabolite annotation and identification, pathways elucidation, and further crop improvement.

ODZ1 allows glioblastoma to sustain invasiveness through a Myc-dependent transcriptional upregulation of RhoA

Long-term survival remains low for most patients with glioblastoma (GBM), which reveals the need for markers of disease outcome and novel therapeutic targets. We describe that ODZ1 (also known as TENM1), a type II transmembrane protein involved in fetal brain development, plays a crucial role in the invasion of GBM cells. Differentiation of glioblastoma stem-like cells drives the nuclear translocation of an intracellular fragment of ODZ1 through proteolytic cleavage by signal peptide peptidase-like 2a. The intracellular fragment of ODZ1 promotes cytoskeletal remodelling of GBM cells and invasion of the surrounding environment both in vitro and in vivo. Absence of ODZ1 by gene deletion or downregulation of ODZ1 by small interfering RNAs drastically reduces the invasive capacity of GBM cells. This activity is mediated by an ODZ1-triggered transcriptional pathway, through the E-box binding Myc protein, that promotes the expression and activation of Ras homolog family member A (RhoA) and subsequent activation of Rho-associated, coiled-coil containing protein kinase (ROCK). Overexpression of ODZ1 in GBM cells reduced survival of xenografted mice. Consistently, analysis of 122 GBM tumour samples revealed that the number of ODZ1-positive cells inversely correlated with overall and progression-free survival. Our findings establish a novel marker of invading GBM cells and consequently a potential marker of disease progression and a therapeutic target in GBM.

Association mapping reveals the genetic architecture of tomato response to water deficit: focus on major fruit quality traits

Water scarcity constitutes a crucial constraint for agriculture productivity. High-throughput approaches in model plant species identified hundreds of genes potentially involved in survival under drought, but few having beneficial effects on quality and yield. Nonetheless, controlled water deficit may improve fruit quality through higher concentration of flavor compounds. The underlying genetic determinants are still poorly known. In this study, we phenotyped 141 highly diverse small fruit tomato accessions for 27 traits under two contrasting watering conditions. A subset of 55 accessions exhibited increased metabolite contents and maintained yield under water deficit. Using 6100 single nucleotide polymorphisms (SNPs), association mapping revealed 31, 41, and 44 quantitative trait loci (QTLs) under drought, control, and both conditions, respectively. Twenty-five additional QTLs were interactive between conditions, emphasizing the interest in accounting for QTLs by watering regime interactions in fruit quality improvement. Combining our results with the loci previously identified in a biparental progeny resulted in 11 common QTLs and contributed to a first detailed characterization of the genetic determinants of response to water deficit in tomato. Major QTLs for fruit quality traits were dissected and candidate genes were proposed using expression and polymorphism data. The outcomes provide a basis for fruit quality improvement under deficit irrigation while limiting yield losses.

Safe Reduction in CD4 Cell Count Monitoring in Stable, Virally Suppressed Patients With HIV Infection or HIV/Hepatitis C Virus Coinfection

BACKGROUND

It has been suggested that routine CD4 cell count monitoring in human immunodeficiency virus (HIV)-monoinfected patients with suppressed viral loads and CD4 cell counts >300 cell/μL could be reduced to annual. HIV/hepatitis C virus (HCV) coinfection is frequent, but evidence supporting similar reductions in CD4 cell count monitoring is lacking for this population. We determined whether CD4 cell count monitoring could be reduced in monoinfected and coinfected patients by estimating the probability of maintaining CD4 cell counts ≥200 cells/µL during continuous HIV suppression.

METHODS

The PISCIS Cohort study included data from 14 539 patients aged ≥16 years from 10 hospitals in Catalonia and 2 in the Balearic Islands (Spain) since January 1998. All patients who had at least one period of 6 months of continuous HIV suppression were included in this analysis. Cumulative probabilities with 95% confidence intervals were calculated using the Kaplan-Meier estimator stratified by the initial CD4 cell count at the period of continuous suppression initiation.

RESULTS

A total of 8695 patients were included. CD4 cell counts fell to <200 cells/µL in 7.4% patients, and the proportion was lower in patients with an initial count >350 cells/µL (1.8%) and higher in those with an initial count of 200-249 cells/µL (23.1%). CD4 cell counts fell to <200 cells/µL in 5.7% of monoinfected and 11.1% of coinfected patients. Of monoinfected patients with an initial CD4 cell count of 300-349 cells/µL, 95.6% maintained counts ≥200 cells/µL. In the coinfected group with the same initial count, this rate was lower, but 97.6% of coinfected patients with initial counts >350 cells/µL maintained counts ≥200 cells/µL.

CONCLUSIONS

From our data, it can be inferred that CD4 cell count monitoring can be safely performed annually in HIV-monoinfected patients with CD4 cell counts >300 cells/µL and HIV/HCV-coinfected patients with counts >350 cells/µL.

New resources for genetic studies in Populus nigra: genome-wide SNP discovery and development of a 12k Infinium array

Whole genome resequencing of 51 Populus nigra (L.) individuals from across Western Europe was performed using Illumina platforms. A total number of 1 878 727 SNPs distributed along the P. nigra reference sequence were identified. The SNP calling accuracy was validated with Sanger sequencing. SNPs were selected within 14 previously identified QTL regions, 2916 expressional candidate genes related to rust resistance, wood properties, water-use efficiency and bud phenology and 1732 genes randomly spread across the genome. Over 10 000 SNPs were selected for the construction of a 12k Infinium Bead-Chip array dedicated to association mapping. The SNP genotyping assay was performed with 888 P. nigra individuals. The genotyping success rate was 91%. Our high success rate was due to the discovery panel design and the stringent parameters applied for SNP calling and selection. In the same set of P. nigra genotypes, linkage disequilibrium throughout the genome decayed on average within 5-7 kb to half of its maximum value. As an application test, ADMIXTURE analysis was performed with a selection of 600 SNPs spread throughout the genome and 706 individuals collected along 12 river basins. The admixture pattern was consistent with genetic diversity revealed by neutral markers and the geographical distribution of the populations. These newly developed SNP resources and genotyping array provide a valuable tool for population genetic studies and identification of QTLs through natural-population based genetic association studies in P. nigra.

Paediatric tumour boards in Spain: a national survey

PURPOSE

Multidisciplinary tumour boards (MDTs) are conducted worldwide for the management of patients with cancer, and they deliver a higher standard of care by simultaneously involving different specialists in diagnosis and treatment planning. However, information of paediatric MDTs functioning is scarce. A pilot study was conducted in Spain in the frame of the European Expert Paediatric Oncology Reference Network for Diagnostics and Treatment (ExPO-r-Net).

METHODS

A specific questionnaire was designed regarding various features of MDT practice. Data collected included information on the centres and the team, infrastructure for meetings, MDT organization/logistics and clinical decision-making. The survey was distributed to all Paediatric Oncology Units that register patients in the Spanish Registry of Childhood Tumours (RETI-SEHOP).

RESULTS

32 out of 43 contacted centres responded the questionnaire (74 % response rate; 88 % response rate for centres with >25 new patients/year). All units with >25 new patients/year have a dedicated Paediatric MDT compared to 76 % of units with ≤25 new patients/year. MDTs should be improved at institutional level by clear protected time in service planning for all specialists involved, incentives for attendance and attendance registration. Clinical decision-making process and follow-up of recommendation adherence should be assessed and potential legal responsibilities for physicians participating in Tumour Board defined. Network collaboration through virtual MDTs, using available videoconferencing tools, is an opportunity to share expertise among centres.

Concomitant histone deacetylase and phosphodiesterase 5 inhibition synergistically prevents the disruption in synaptic plasticity and it reverses cognitive impairment in a mouse model of Alzheimer's disease

BACKGROUND

Given the implication of histone acetylation in memory processes, histone deacetylase inhibitors (HDACIs) have been postulated as potential modulators of cognitive impairment in Alzheimer's disease (AD). However, dose-dependent side effects have been described in patients with the currently available broad-spectrum HDACIs, explaining why their therapeutic potential has not been realized for chronic diseases. Here, by simultaneously targeting two independent enzyme activities, histone deacetylase (HDAC) and phosphodiesterase-5 (PDE5), we propose a novel mode of inhibitory action that might increase the therapeutic specificity of HDACIs.

RESULTS

The combination of vorinostat, a pan-HDACI, and tadalafil, a PDE5 inhibitor, rescued the long-term potentiation impaired in slices from APP/PS1 mice. When administered in vivo, the combination of these drugs alleviated the cognitive deficits in AD mice, as well as the amyloid and tau pathology, and it reversed the reduced dendritic spine density on hippocampal neurons. Significantly, the combination of vorinostat and tadalafil was more effective than each drug alone, both against the symptoms and in terms of disease modification, and importantly, these effects persisted after a 4-week washout period.

CONCLUSIONS

The results highlight the pharmacological potential of a combination of molecules that inhibit HDAC and PDE5 as a therapeutic approach for AD treatment.

Genome-wide prediction methods in highly diverse and heterozygous species: proof-of-concept through simulation in grapevine

Nowadays, genome-wide association studies (GWAS) and genomic selection (GS) methods which use genome-wide marker data for phenotype prediction are of much potential interest in plant breeding. However, to our knowledge, no studies have been performed yet on the predictive ability of these methods for structured traits when using training populations with high levels of genetic diversity. Such an example of a highly heterozygous, perennial species is grapevine. The present study compares the accuracy of models based on GWAS or GS alone, or in combination, for predicting simple or complex traits, linked or not with population structure. In order to explore the relevance of these methods in this context, we performed simulations using approx 90,000 SNPs on a population of 3,000 individuals structured into three groups and corresponding to published diversity grapevine data. To estimate the parameters of the prediction models, we defined four training populations of 1,000 individuals, corresponding to these three groups and a core collection. Finally, to estimate the accuracy of the models, we also simulated four breeding populations of 200 individuals. Although prediction accuracy was low when breeding populations were too distant from the training populations, high accuracy levels were obtained using the sole core-collection as training population. The highest prediction accuracy was obtained (up to 0.9) using the combined GWAS-GS model. We thus recommend using the combined prediction model and a core-collection as training population for grapevine breeding or for other important economic crops with the same characteristics.

Different β-adrenoceptor subtypes coupling to cAMP or NO/cGMP pathways: implications in the relaxant response of rat conductance and resistance vessels

BACKGROUND AND PURPOSE

To analyse the relative contribution of β1 -, β2 - and β3 -adrenoceptors (Adrb) to vasodilatation in conductance and resistance vessels, assessing the role of cAMP and/or NO/cGMP signalling pathways.

EXPERIMENTAL APPROACH

Rat mesenteric resistance artery (MRA) and aorta were used to analyse the Adrb expression by real-time-PCR and immunohistochemistry, and for the pharmacological characterization of Adrb-mediated activity by wire myography and tissue nucleotide accumulation.

KEY RESULTS

The mRNAs and protein for all Adrb were identified in endothelium and/or smooth muscle cells (SMCs) in both vessels. In MRA, Adrb1 signalled through cAMP, Adrb3 through both cAMP and cGMP, but Adrb2, did not activate nucleotide formation; isoprenaline relaxation was inhibited by propranolol (β1 , β2 ), CGP20712A (β1 ), and SQ22536 (adenylyl cyclase inhibitor), but not by ICI118,551 (β2 ), SR59230A (β3 ), ODQ (soluble guanylyl cyclase inhibitor), L-NAME or endothelium removal. In aorta, Adrb1 signalled through cAMP, while β2 - and β3 -subtypes through cGMP; isoprenaline relaxation was inhibited by propranolol, ICI118,551, ODQ, L-NAME, and to a lesser extent, by endothelium removal. CL316243 (β3 -agonist) relaxed aorta, but not MRA.

CONCLUSION AND IMPLICATION

Despite all three Adrb subtypes being found in both vessels, Adrb1, located in SMCs and acting through the adenylyl cyclase/cAMP pathway, are primarily responsible for vasodilatation in MRA. However, Adrb-mediated vasodilatation in aorta is driven by endothelial Adrb2 and Adrb3, but also by the Adrb2 present in SMCs, and is coupled to the NO/cGMP pathway. These results could help to understand the different physiological roles played by Adrb signalling in regulating conductance and resistance vessels.

Genome-Wide Association in Tomato Reveals 44 Candidate Loci for Fruit Metabolic Traits

Genome-wide association studies have been successful in identifying genes involved in polygenic traits and are valuable for crop improvement. Tomato (Solanum lycopersicum) is a major crop and is highly appreciated worldwide for its health value. We used a core collection of 163 tomato accessions composed of S. lycopersicum, S. lycopersicum var cerasiforme, and Solanum pimpinellifolium to map loci controlling variation in fruit metabolites. Fruits were phenotyped for a broad range of metabolites, including amino acids, sugars, and ascorbate. In parallel, the accessions were genotyped with 5,995 single-nucleotide polymorphism markers spread over the whole genome. Genome-wide association analysis was conducted on a large set of metabolic traits that were stable over 2 years using a multilocus mixed model as a general method for mapping complex traits in structured populations and applied to tomato. We detected a total of 44 loci that were significantly associated with a total of 19 traits, including sucrose, ascorbate, malate, and citrate levels. These results not only provide a list of candidate loci to be functionally validated but also a powerful analytical approach for finding genetic variants that can be directly used for crop improvement and deciphering the genetic architecture of complex traits.

The three α1-adrenoceptor subtypes show different spatio-temporal mechanisms of internalization and ERK1/2 phosphorylation

We analyzed the kinetic and spatial patterns characterizing activation of the MAP kinases ERK 1 and 2 (ERK1/2) by the three α1-adrenoceptor (α1-AR) subtypes in HEK293 cells and the contribution of two different pathways to ERK1/2 phosphorylation: protein kinase C (PKC)-dependent ERK1/2 activation and internalization-dependent ERK1/2 activation. The different pathways of phenylephrine induced ERK phosphorylation were determined by western blot, using the PKC inhibitor Ro 31-8425, the receptor internalization inhibitor concanavalin A and the siRNA targeting β-arrestin 2. Receptor internalization properties were studied using CypHer5 technology and VSV-G epitope-tagged receptors. Activation of α1A- and α1B-ARs by phenylephrine elicited rapid ERK1/2 phosphorylation that was directed to the nucleus and inhibited by Ro 31-8425. Concomitant with phenylephrine induced receptor internalization α1A-AR, but not α1B-AR, produced a maintained and PKC-independent ERK phosphorylation, which was restricted to the cytosol and inhibited by β-arrestin 2 knockdown or concanavalin A treatment. α1D-AR displayed constitutive ERK phosphorylation, which was reduced by incubation with prazosin or the selective α1D antagonist BMY7378. Following activation by phenylephrine, α1D-AR elicited rapid, transient ERK1/2 phosphorylation that was restricted to the cytosol and not inhibited by Ro 31-8425. Internalization of the α1D-AR subtype was not observed via CypHer5 technology. The three α1-AR subtypes present different spatio-temporal patterns of receptor internalization, and only α1A-AR stimulation translates to a late, sustained ERK1/2 phosphorylation that is restricted to the cytosol and dependent on β-arrestin 2 mediated internalization.

The use of combining ability analysis to identify elite parents for Artemisia annua F1 hybrid production

Artemisia annua is an important medicinal crop used for the production of the anti-malarial compound artemisinin. In order to assist in the production of affordable high quality artemisinin we have carried out an A. annua breeding programme aimed at improving artemisinin concentration and biomass. Here we report on a combining ability analysis of a diallel cross to identify robust parental lines for hybrid breeding. The parental lines were selected based on a range of phenotypic traits to encourage heterosis. The general combining ability (GCA) values for the diallel parental lines correlated to the positive alleles of quantitative trait loci (QTL) in the same parents indicating the presence of beneficial alleles that contribute to parental performance. Hybrids generated from crossing specific parental lines with good GCA were identified as having an increase in both artemisinin concentration and biomass when grown either in glasshouse or experimental field trials and compared to controls. This study demonstrates that combining ability as determined by a diallel cross can be used to identify elite parents for the production of improved A. annua hybrids. Furthermore, the selection of material for breeding using this approach was found to be consistent with our QTL-based molecular breeding approach.

Spanish human proteome project: dissection of chromosome 16

The Chromosome 16 Consortium forms part of the Human Proteome Project that aims to develop an entire map of the proteins encoded by the human genome following a chromosome-centric strategy (C-HPP) to make progress in the understanding of human biology in health and disease (B/D-HPP). A Spanish consortium of 16 laboratories was organized into five working groups: Protein/Antibody microarrays, protein expression and Peptide Standard, S/MRM, Protein Sequencing, Bioinformatics and Clinical healthcare, and Biobanking. The project is conceived on a multicenter configuration, assuming the standards and integration procedures already available in ProteoRed-ISCIII, which is encompassed within HUPO initiatives. The products of the 870 protein coding genes in chromosome 16 were analyzed in Jurkat T lymphocyte cells, MCF-7 epithelial cells, and the CCD18 fibroblast cell line as it is theoretically expected that most chromosome 16 protein coding genes are expressed in at least one of these. The transcriptome and proteome of these cell lines was studied using gene expression microarray and shotgun proteomics approaches, indicating an ample coverage of chromosome 16. With regard to the B/D section, the main research areas have been adopted and a biobanking initiative has been designed to optimize methods for sample collection, management, and storage under normalized conditions and to define QC standards. The general strategy of the Chr-16 HPP and the current state of the different initiatives are discussed.

GWAPP: a web application for genome-wide association mapping in Arabidopsis

Arabidopsis thaliana is an important model organism for understanding the genetics and molecular biology of plants. Its highly selfing nature, small size, short generation time, small genome size, and wide geographic distribution make it an ideal model organism for understanding natural variation. Genome-wide association studies (GWAS) have proven a useful technique for identifying genetic loci responsible for natural variation in A. thaliana. Previously genotyped accessions (natural inbred lines) can be grown in replicate under different conditions and phenotyped for different traits. These important features greatly simplify association mapping of traits and allow for systematic dissection of the genetics of natural variation by the entire A. thaliana community. To facilitate this, we present GWAPP, an interactive Web-based application for conducting GWAS in A. thaliana. Using an efficient implementation of a linear mixed model, traits measured for a subset of 1386 publicly available ecotypes can be uploaded and mapped with a mixed model and other methods in just a couple of minutes. GWAPP features an extensive, interactive, and user-friendly interface that includes interactive Manhattan plots and linkage disequilibrium plots. It also facilitates exploratory data analysis by implementing features such as the inclusion of candidate polymorphisms in the model as cofactors.

A mixed-model approach for genome-wide association studies of correlated traits in structured populations

Genome-wide association studies (GWAS) are a standard approach for studying the genetics of natural variation. A major concern in GWAS is the need to account for the complicated dependence-structure of the data both between loci as well as between individuals. Mixed models have emerged as a general and flexible approach for correcting for population structure in GWAS. Here we extend this linear mixed model approach to carry out GWAS of correlated phenotypes, deriving a fully parameterized multi-trait mixed model (MTMM) that considers both the within-trait and between-trait variance components simultaneously for multiple traits. We apply this to human cohort data for correlated blood lipid traits from the Northern Finland Birth Cohort 1966, and demonstrate greatly increased power to detect pleiotropic loci that affect more than one blood lipid trait. We also apply this to an Arabidopsis dataset for flowering measurements in two different locations, identifying loci whose effect depends on the environment.

An efficient multi-locus mixed-model approach for genome-wide association studies in structured populations

Population structure causes genome-wide linkage disequilibrium between unlinked loci, leading to statistical confounding in genome-wide association studies. Mixed models have been shown to handle the confounding effects of a diffuse background of large numbers of loci of small effect well, but they do not always account for loci of larger effect. Here we propose a multi-locus mixed model as a general method for mapping complex traits in structured populations. Simulations suggest that our method outperforms existing methods in terms of power as well as false discovery rate. We apply our method to human and Arabidopsis thaliana data, identifying new associations and evidence for allelic heterogeneity. We also show how a priori knowledge from an A. thaliana linkage mapping study can be integrated into our method using a Bayesian approach. Our implementation is computationally efficient, making the analysis of large data sets (n > 10,000) practicable.

Preclinical activity of LBH589 alone or in combination with chemotherapy in a xenogeneic mouse model of human acute lymphoblastic leukemia

Histone deacetylases (HDACs) have been identified as therapeutic targets due to their regulatory function in chromatin structure and organization. Here, we analyzed the therapeutic effect of LBH589, a class I-II HDAC inhibitor, in acute lymphoblastic leukemia (ALL). In vitro, LBH589 induced dose-dependent antiproliferative and apoptotic effects, which were associated with increased H3 and H4 histone acetylation. Intravenous administration of LBH589 in immunodeficient BALB/c-RAG2(-/-)γc(-/-) mice in which human-derived T and B-ALL cell lines were injected induced a significant reduction in tumor growth. Using primary ALL cells, a xenograft model of human leukemia in BALB/c-RAG2(-/-)γc(-/-) mice was established, allowing continuous passages of transplanted cells to several mouse generations. Treatment of mice engrafted with T or B-ALL cells with LBH589 induced an in vivo increase in the acetylation of H3 and H4, which was accompanied with prolonged survival of LBH589-treated mice in comparison with those receiving vincristine and dexamethasone. Notably, the therapeutic efficacy of LBH589 was significantly enhanced in combination with vincristine and dexamethasone. Our results show the therapeutic activity of LBH589 in combination with standard chemotherapy in pre-clinical models of ALL and suggest that this combination may be of clinical value in the treatment of patients with ALL.

Blockade of the NFκB pathway drives differentiating glioblastoma-initiating cells into senescence both in vitro and in vivo

Glioblastoma multiforme is one of the most devastating cancers and presents unique challenges to therapy because of its aggressive behavior. Cancer-initiating or progenitor cells have been described to be the only cell population with tumorigenic capacity in glioblastoma. Therefore, effective therapeutic strategies targeting these cells or the early precursors may be beneficial. We have established different cultures of glioblastoma-initiating cells (GICs) derived from surgical specimens and found that, after induction of differentiation, the NFκB transcriptional pathway was activated, as determined by analyzing key proteins such as p65 and IκB and the upregulation of a number of target genes. We also showed that blockade of nuclear factor (NF)κB signaling in differentiating GICs by different genetic strategies or treatment with small-molecule inhibitors, promoted replication arrest and senescence. This effect was partly mediated by reduced levels of the NFκB target gene cyclin D1, because its downregulation by RNA interference reproduced a similar phenotype. Furthermore, these results were confirmed in a xenograft model. Intravenous treatment of immunodeficient mice bearing human GIC-derived tumors with a novel small-molecule inhibitor of the NFκB pathway induced senescence of tumor cells but no ultrastructural alterations of the brain parenchyma were detected. These findings reveal that activation of NFκB may keep differentiating GICs from acquiring a mature postmitotic phenotype, thus allowing cell proliferation, and support the rationale for therapeutic strategies aimed to promote premature senescence of differentiating GICs by blocking key factors within the NFκB pathway.

Chronic stress and impaired glutamate function elicit a depressive-like phenotype and common changes in gene expression in the mouse frontal cortex

Major depression might originate from both environmental and genetic risk factors. The environmental chronic mild stress (CMS) model mimics some environmental factors contributing to human depression and induces anhedonia and helplessness. Mice heterozygous for the synaptic vesicle protein (SVP) vesicular glutamate transporter 1 (VGLUT1) have been proposed as a genetic model of deficient glutamate function linked to depressive-like behaviour. Here, we aimed to identify, in these two experimental models, gene expression changes in the frontal cortex, common to stress and impaired glutamate function. Both VGLUT1(+/-) and CMS mice showed helpless and anhedonic-like behavior. Microarray studies in VGLUT1(+/-) mice revealed regulation of genes involved in apoptosis, neurogenesis, synaptic transmission, protein metabolic process or learning and memory. In addition, RT-PCR studies confirmed gene expression changes in several glutamate, GABA, dopamine and serotonin neurotransmitter receptors. On the other hand, CMS affected the regulation of 147 transcripts, some of them involved in response to stress and oxidoreductase activity. Interestingly, 52 genes were similarly regulated in both models. Specifically, a dowregulation in genes that promote cell proliferation (Anapc7), cell growth (CsnK1g1), cell survival (Hdac3), and inhibition of apoptosis (Dido1) was observed. Genes linked to cytoskeleton (Hspg2, Invs), psychiatric disorders (Grin1, MapK12) or an antioxidant enzyme (Gpx2) were also downregulated. Moreover, genes that inhibit the MAPK pathways (Dusp14), stimulate oxidative metabolism (Eif4a2) and enhance glutamate transmission (Rab8b) were upregulated. We suggest that these genes could form part of the altered "molecular context" underlying depressive-like behaviour in animal models. The clinical relevance of these findings is discussed.