Coupling clearing and hybridization chain reaction approaches to investigate gene expression in organs inside intact insect heads

Detecting RNA molecules within their natural environment inside intact arthropods has long been challenging, particularly in small organisms covered by a tanned and pigmented cuticle. Here, we have developed a methodology that enables high-resolution analysis of the spatial distribution of transcripts of interest without having to dissect tiny organs or tissues, thereby preserving their integrity. We have combined an in situ amplification approach based on hybridization chain reaction, which enhances the signal-to-noise ratio, and a clearing approach that allows the visualization of inner organs beneath the cuticle. We have implemented this methodology for the first time in Hemiptera, mapping two salivary aphid (Acyrthosiphon pisum) transcripts, the effector c002 and the salivary sheath protein SHP. With a multiplex approach, we could simultaneously detect different mRNAs in mounted pea aphid head-thorax samples and show that they were distributed in distinct secretory cells of salivary glands. RESEARCH HIGHLIGHTS: Combining hybridisation chain reaction and clearing allows the detection of transcripts in intact aphids heads. The transcripts of the two salivary proteins c002 and SHP are compartmentalized in distinct secretory cells of the principal glands.

Localized osmotic stress activates systemic responses to N limitation in Medicago truncatula-Sinorhizobium symbiotic plants

In mature symbiotic root nodules, differentiated rhizobia fix atmospheric dinitrogen and provide ammonium to fulfill the plant nitrogen (N) demand. The plant enables this process by providing photosynthates to the nodules. The symbiosis is adjusted to the whole plant N demand thanks to systemic N signaling controlling nodule development. Symbiotic plants under N deficit stimulate nodule expansion and activate nodule senescence under N satiety. Besides, nodules are highly sensitive to drought. Here, we used split-root systems to characterize the systemic responses of symbiotic plants to a localized osmotic stress. We showed that polyéthylène glycol (PEG) application rapidly inhibited the symbiotic dinitrogen fixation activity of nodules locally exposed to the treatment, resulting to the N limitation of the plant supplied exclusively by symbiotic dinitrogen fixation. The localized PEG treatment triggered systemic signaling stimulating nodule development in the distant untreated roots. This response was associated with an enhancement of the sucrose allocation. Our analyses showed that transcriptomic reprogramming associated with PEG and N deficit systemic signaling(s) shared many targets transcripts. Altogether, our study suggests that systemic N signaling is a component of the adaptation of the symbiotic plant to the local variations of its edaphic environment.

Evidence of Covid-19 lockdown effects on riverine dissolved organic matter dynamics provides a proof-of-concept for needed regulations of anthropogenic emissions

The fast spread of SARS-CoV-2 virus in Italy resulted in a 3-months lockdown of the entire country. During this period, the effect of the relieved anthropogenic activities on the environment was plainly clear all over the country. Herein, we provide the first evidence of the lockdown effects on riverine dissolved organic matter (DOM) dynamics. The strong reduction in anthropogenic activities resulted in a marked decrease in dissolved organic carbon (DOC) concentration in the Arno River (-44%) and the coastal area affected by its input (-15%), compared to previous conditions. The DOM optical properties (absorption and fluorescence) showed a change in its quality, with a shift toward smaller and less aromatic molecules during the lockdown. The reduced human activity and the consequent change in DOM dynamics affected the abundance and annual dynamics of heterotrophic prokaryotes. The results of this study highlight the extent to which DOM dynamics in small rivers is affected by secondary and tertiary human activities as well as the quite short time scales to return to the impacted conditions. Our work also supports the importance of long-term research to disentangle the effects of casual events from the natural variability.

Systemic control of nodule formation by plant nitrogen demand requires autoregulation-dependent and independent mechanisms

In legumes interacting with rhizobia, the formation of symbiotic organs involved in the acquisition of atmospheric nitrogen gas (N2) is dependent on the plant nitrogen (N) demand. We used Medicago truncatula plants cultivated in split-root systems to discriminate between responses to local and systemic N signaling. We evidenced a strong control of nodule formation by systemic N signaling but obtained no clear evidence of a local control by mineral nitrogen. Systemic signaling of the plant N demand controls numerous transcripts involved in root transcriptome reprogramming associated with early rhizobia interaction and nodule formation. SUPER NUMERIC NODULES (SUNN) has an important role in this control, but we found that major systemic N signaling responses remained active in the sunn mutant. Genes involved in the activation of nitrogen fixation are regulated by systemic N signaling in the mutant, explaining why its hypernodulation phenotype is not associated with higher nitrogen fixation of the whole plant. We show that the control of transcriptome reprogramming of nodule formation by systemic N signaling requires other pathway(s) that parallel the SUNN/CLE (CLAVATA3/EMBRYO SURROUNDING REGION-LIKE PEPTIDES) pathway.

The transposable element-rich genome of the cereal pest Sitophilus oryzae

BACKGROUND

The rice weevil Sitophilus oryzae is one of the most important agricultural pests, causing extensive damage to cereal in fields and to stored grains. S. oryzae has an intracellular symbiotic relationship (endosymbiosis) with the Gram-negative bacterium Sodalis pierantonius and is a valuable model to decipher host-symbiont molecular interactions.

RESULTS

We sequenced the Sitophilus oryzae genome using a combination of short and long reads to produce the best assembly for a Curculionidae species to date. We show that S. oryzae has undergone successive bursts of transposable element (TE) amplification, representing 72% of the genome. In addition, we show that many TE families are transcriptionally active, and changes in their expression are associated with insect endosymbiotic state. S. oryzae has undergone a high gene expansion rate, when compared to other beetles. Reconstruction of host-symbiont metabolic networks revealed that, despite its recent association with cereal weevils (30 kyear), S. pierantonius relies on the host for several amino acids and nucleotides to survive and to produce vitamins and essential amino acids required for insect development and cuticle biosynthesis.

CONCLUSIONS

Here we present the genome of an agricultural pest beetle, which may act as a foundation for pest control. In addition, S. oryzae may be a useful model for endosymbiosis, and studying TE evolution and regulation, along with the impact of TEs on eukaryotic genomes.

Identification and characterization of Vietnamese coffee bacterial endophytes displaying in vitro antifungal and nematicidal activities

The endophytic bacteria were isolated from coffee roots and seeds in Vietnam and identified with 16S rDNA sequencing as belonging to the Actinobacteria, Firmicutes and Proteobacteria phyla with the Nocardia, Bacillus and Burkholderia as dominant genera, respectively. Out of the thirty genera recovered from Coffea canephora and Coffea liberica, twelve were reported for the first time in endophytic association with coffee including members of the genera Brachybacterium, Caballeronia, Kitasatospora, Lechevalieria, Leifsonia, Luteibacter, Lysinibacillus, Mycolicibacterium, Nakamurella, Paracoccus, Sinomonas and Sphingobium. A total of eighty bacterial endophytes were characterized in vitro for several plant growth promoting and biocontrol traits including: the phosphate solubilization, the indolic compounds, siderophores, HCN, esterase, lipase, gelatinase and chitinase production. A subset of fifty selected bacteria were tested for their potential as biocontrol agents with in vitro confrontations with the fungal pathogen Fusarium oxysporum as well as the coffee parasitic nematodes Radopholus duriophilus and Pratylenchus coffeae. The three most efficient isolates on F. oxysporum belonging to the Bacillus, Burkholderia, and Streptomyces genera displayed a growth inhibition rate higher than 40%. Finally, five isolates from the Bacillus genus were able to lead to 100% of mortality in 24 h on both R. duriophilus and P. coffeae.

Responses of mature symbiotic nodules to the whole-plant systemic nitrogen signaling

In symbiotic root nodules of legumes, terminally differentiated rhizobia fix atmospheric N2 producing an NH4+ influx that is assimilated by the plant. The plant, in return, provides photosynthates that fuel the symbiotic nitrogen acquisition. Mechanisms responsible for the adjustment of the symbiotic capacity to the plant N demand remain poorly understood. We have investigated the role of systemic signaling of whole-plant N demand on the mature N2-fixing nodules of the model symbiotic association Medicago truncatula/Sinorhizobium using split-root systems. The whole-plant N-satiety signaling rapidly triggers reductions of both N2 fixation and allocation of sugars to the nodule. These responses are associated with the induction of nodule senescence and the activation of plant defenses against microbes, as well as variations in sugars transport and nodule metabolism. The whole-plant N-deficit responses mirror these changes: a rapid increase of sucrose allocation in response to N-deficit is associated with a stimulation of nodule functioning and development resulting in nodule expansion in the long term. Physiological, transcriptomic, and metabolomic data together provide evidence for strong integration of symbiotic nodules into whole-plant nitrogen demand by systemic signaling and suggest roles for sugar allocation and hormones in the signaling mechanisms.

THU0257 ESTIMATED 10-YEARS CARDIOVASCULAR RISK IN SYSTEMIC LUPUS ERYTHEMATOSUS PATIENTS: PRELIMINARY RESULTS FROM THE “CARDIOVASCULAR OBESITY AND RHEUMATIC DISEASE (CORDIS)” STUDY GROUP OF THE ITALIAN SOCIETY OF RHEUMATOLOGY

Background:

Systemic lupus erythematosus (SLE) patients are at high risk for CV events, and EULAR recommends assessing the 10-year CV-risk using the Systematic Coronary Risk Evaluation (SCORE) [1]. The QRISK3, another score to assess CV-risk in UK population, considers different factors among which also SLE. The Progetto Cuore score (PCS) is validated to estimate CV risk in Italian people and largely replicates the SCORE project [2].

Objectives:

This cross-sectional study aimed to estimate CV-risk using SCORE, QRISK3 and, for the first time, PCS in a multicentric cohort of Italian SLE patients.

Methods:

During 2019 we evaluated 173 SLE patients (87.7% female; age 40±16 years; disease duration 138±105 months), fulfilling the 1997 ACR classification criteria. Clinical and laboratory data were registered, and individual CV-risk was calculated using suitable algorithms for the SCORE, QRISK3 and PCS. Statistical analysis was performed using Graphpad Instat 8.0 (San Diego, CA-USA).

Results:

In 13 (7%) SLE patients a previous CV event was recorded. Hypertension was present in 60 (37.5%) and diabetes in 27 (16.9%) patients. Mean total cholesterol was 184±39 mg/dL, HDLc 58±18 mg/dL, LDLc 124±37 mg/dL, triglycerides 105±63 mg/dL; dyslipidaemia was reported in 58 (36.2%) patients and 29 (18.1%) were on statin. Mean BMI was 24.9±5.3 Kg/sm, 60 (37.5%) and 23 (14.3%) patients were overweight and obese, while 25 (15.6%) patients were smokers. 87 (54.3%) SLE patients had a SLEDAI<4, 91% of patients were taken HCQ and 65% were on prednisone (mean dose 5.4±5.9 mg/day), but only 7.5% took >7.5 mg/day. The CV-risk of SLE patients according to SCORE, QRISK3 and PCS was 1.1±2.1%, 10.5±12.3% and 3.7±5.4%, respectively. Stratifying patients at low, moderate or high CV risk according to the PCS and SCORE a double proportion of patients was at moderate (8% vs 3.9%) or high (1.9% vs 0.9%) CV risk (p=0.03). Finally, CV-risk according to QRISK3 was higher than 20% (high risk) in 32/160 (20%) patients.

Conclusion:

This multicentre study demonstrated that the mean estimated CV-risk in SLE patients is globally low using the SCORE, QRISK3 and PCS. The PCS seems to better intercept those patients at moderate/high risk, at least in Italian SLE patients, while QRISK3 predicts the highest CV risk. The lack of disease-specific CV-risk factors (such as autoantibodies profile or organ involvement) probably account for the underestimation of CV risk using the SCORE and PCS.

References:

[1]ARD 2019;78(6):736-745.

[2]ARD 2019;0:1–2.doi:10.1136/annrheumdis-2019-215715

Disclosure of Interests:

Fabio Cacciapaglia Speakers bureau: BMS; Roche; Pfizer; Abbvie, Andreina Manfredi: None declared, Gianluca Erre: None declared, Elena Bartoloni Bocci: None declared, Garifallia Sakellariou Speakers bureau: Abbvie, Novartis, MSD, Ombretta Viapiana: None declared, Sergio Colella: None declared, Anna Abbruzzese: None declared, Marco Fornaro: None declared, Giacomo Cafaro: None declared, Maria Antonietta Fenu: None declared, Bianca Lucia Palermo: None declared, Martina Dessì: None declared, Adalgisa Palermo: None declared, Alessandro Giollo: None declared, Elisa Gremese Consultant of: AbbVie, Bristol-Myers Squibb, Celgene, Eli Lilly, Janssen, Merck Sharp & Dohme, Novartis, Sanofi, UCB, Roche, Pfizer, Speakers bureau: AbbVie, Bristol-Myers Squibb, Celgene, Eli Lilly, Janssen, Merck Sharp & Dohme, Novartis, Sanofi, UCB, Roche, Pfizer, Francesca Romana Spinelli Grant/research support from: Pfizer, Speakers bureau: Lilly, BMS, Celgene, Fabiola Atzeni: None declared, Matteo Piga: None declared

THU0127 Estimated cardiovascular risk in a large cohort of rheumatoid arthritis patients from the “Cardiovascular Obesity and Rheumatic DISease (CORDIS)” Study Group of the Italian Society of Rheumatology

Background:

Rheumatoid arthritis (RA) patients present high cardiovascular (CV) morbidity and mortality and EULAR recommends estimating their CV-risk [1]. The Systematic Coronary Risk Evaluation (SCORE) algorithm is suggested if National Guidelines are lack, but few data are available about different strategies.

Objectives:

To estimate the 10-years CV-risk using different algorithms in RA compared to osteoarthritis (OA) patients, as control group.

Methods:

A total of 1467 RA patients (78.3% female; mean age 59.8±11.5 years; mean disease duration 131±109 months), fulfilling the 2010 EULAR/ACR classification criteria, and 342 age and sex matched patients with OA (79.8% female; mean age 58.7±11.5 years) were enrolled in this multicentre cross-sectional study during 2019. Clinical and laboratory data were registered, and individual CV-risk was calculated using: SCORE chart, “Progetto Cuore” model (PCM), QRisk3, Reynolds Risk Scores (RRS) and Expanded Risk Score in RA (ERS-RA), as stated by suitable algorithms. Statistical analysis was performed using the Statistical System Graphpad Instat 8.0 (San Diego, CA-USA).

Results:

In 46 (3%) RA patients a previous CV event was observed. Among traditional CV-risk factors, RA patients presented higher frequency of diabetes (9.9% vs 6.4%; p=0.04) and lower prevalence of dyslipidaemia (21.7% vs 32.5%; p<0.0001) compared to OA patients. Prevalence of hypertension was similar in both groups (40% vs 39.2%). Mean BMI (25.6±4.8 vs 26.6±4.4; p<0.0001) and prevalence of obesity (15% vs 21%; p=0.003) were significantly lower in RA patients. Finally, RA patients were more frequently smokers (20.4% vs 12.5% - p=0.002). 441 (30%) RA patients were in CDAI remission, 998 (68%) patients were on csDMARDs while a biologic agent was used in 617 (42%) patients. About 43% of RA patients were on a mean prednisone-dose of 4.5±3.5 mg/day. The 10-years CV-risk resulted 2 to 3-fold higher in RA compared to OA patients using the different algorithms. The QRisk3 estimated the highest CV risk in our cohort of patients, while the ERS-RA and RRS were significantly higher than PCM and SCORE.

Conclusion:

Our study demonstrates a higher estimated CV-risk in RA compared to OA patients. The commonly used algorithms to estimate CV-risk in clinical practice perform differently, evaluating different traditional CV-risk factors and disease specific characteristic, as for QRisk3 or ERS-RA. Rheumatologist should impact on both traditional and RA related modifiable CV-risk factors.

References:

[1]Agca R, et al. Ann Rheum Dis 2017;76:17–28.

Disclosure of Interests:

Fabio Cacciapaglia Speakers bureau: BMS; Roche; Pfizer; Abbvie, Matteo Piga: None declared, Gianluca Erre: None declared, Andreina Manfredi: None declared, Elena Bartoloni Bocci: None declared, Garifallia Sakellariou Speakers bureau: Abbvie, Novartis, MSD, Ombretta Viapiana: None declared, Sergio Colella: None declared, Anna Abbruzzese: None declared, Martina Dessì: None declared, Caterina Vacchi: None declared, Floriana Castagna: None declared, Giacomo Cafaro: None declared, Bianca Lucia Palermo: None declared, Alessandro GIollo: None declared, Marco Fornaro: None declared, Elisa Gremese Consultant of: AbbVie, Bristol-Myers Squibb, Celgene, Eli Lilly, Janssen, Merck Sharp & Dohme, Novartis, Sanofi, UCB, Roche, Pfizer, Speakers bureau: AbbVie, Bristol-Myers Squibb, Celgene, Eli Lilly, Janssen, Merck Sharp & Dohme, Novartis, Sanofi, UCB, Roche, Pfizer, Francesca Romana Spinelli Grant/research support from: Pfizer, Speakers bureau: Lilly, BMS, Celgene, Fabiola Atzeni: None declared

Host-specific competitiveness to form nodules in Rhizobium leguminosarum symbiovar viciae

Fabeae legumes such as pea and faba bean form symbiotic nodules with a large diversity of soil Rhizobium leguminosarum symbiovar viciae (Rlv) bacteria. However, bacteria competitive to form root nodules (CFN) are generally not the most efficient to fix dinitrogen, resulting in a decrease in legume crop yields. Here, we investigate differential selection by host plants on the diversity of Rlv. A large collection of Rlv was collected by nodule trapping with pea and faba bean from soils at five European sites. Representative genomes were sequenced. In parallel, diversity and abundance of Rlv were estimated directly in these soils using metabarcoding. The CFN of isolates was measured with both legume hosts. Pea/faba bean CFN were associated to Rlv genomic regions. Variations of bacterial pea and/or faba bean CFN explained the differential abundance of Rlv genotypes in pea and faba bean nodules. No evidence was found for genetic association between CFN and variations in the core genome, but variations in specific regions of the nod locus, as well as in other plasmid loci, were associated with differences in CFN. These findings shed light on the genetic control of CFN in Rlv and emphasise the importance of host plants in controlling Rhizobium diversity.

DiCoExpress: a tool to process multifactorial RNAseq experiments from quality controls to co-expression analysis through differential analysis based on contrasts inside GLM models

BACKGROUND

RNAseq is nowadays the method of choice for transcriptome analysis. In the last decades, a high number of statistical methods, and associated bioinformatics tools, for RNAseq analysis were developed. More recently, statistical studies realised neutral comparison studies using benchmark datasets, shedding light on the most appropriate approaches for RNAseq data analysis.

RESULTS

DiCoExpress is a script-based tool implemented in R that includes methods chosen based on their performance in neutral comparisons studies. DiCoExpress uses pre-existing R packages including FactoMineR, edgeR and coseq, to perform quality control, differential, and co-expression analysis of RNAseq data. Users can perform the full analysis, providing a mapped read expression data file and a file containing the information on the experimental design. Following the quality control step, the user can move on to the differential expression analysis performed using generalized linear models thanks to the automated contrast writing function. A co-expression analysis is implemented using the coseq package. Lists of differentially expressed genes and identified co-expression clusters are automatically analyzed for enrichment of annotations provided by the user. We used DiCoExpress to analyze a publicly available RNAseq dataset on the transcriptional response of Brassica napus L. to silicon treatment in plant roots and mature leaves. This dataset, including two biological factors and three replicates for each condition, allowed us to demonstrate in a tutorial all the features of DiCoExpress.

CONCLUSIONS

DiCoExpress is an R script-based tool allowing users to perform a full RNAseq analysis from quality controls to co-expression analysis through differential analysis based on contrasts inside generalized linear models. DiCoExpress focuses on the statistical modelling of gene expression according to the experimental design and facilitates the data analysis leading the biological interpretation of the results.

Molecular evolutionary trends and feeding ecology diversification in the Hemiptera, anchored by the milkweed bug genome

BACKGROUND

The Hemiptera (aphids, cicadas, and true bugs) are a key insect order, with high diversity for feeding ecology and excellent experimental tractability for molecular genetics. Building upon recent sequencing of hemipteran pests such as phloem-feeding aphids and blood-feeding bed bugs, we present the genome sequence and comparative analyses centered on the milkweed bug Oncopeltus fasciatus, a seed feeder of the family Lygaeidae.

RESULTS

The 926-Mb Oncopeltus genome is well represented by the current assembly and official gene set. We use our genomic and RNA-seq data not only to characterize the protein-coding gene repertoire and perform isoform-specific RNAi, but also to elucidate patterns of molecular evolution and physiology. We find ongoing, lineage-specific expansion and diversification of repressive C2H2 zinc finger proteins. The discovery of intron gain and turnover specific to the Hemiptera also prompted the evaluation of lineage and genome size as predictors of gene structure evolution. Furthermore, we identify enzymatic gains and losses that correlate with feeding biology, particularly for reductions associated with derived, fluid nutrition feeding.

CONCLUSIONS

With the milkweed bug, we now have a critical mass of sequenced species for a hemimetabolous insect order and close outgroup to the Holometabola, substantially improving the diversity of insect genomics. We thereby define commonalities among the Hemiptera and delve into how hemipteran genomes reflect distinct feeding ecologies. Given Oncopeltus's strength as an experimental model, these new sequence resources bolster the foundation for molecular research and highlight technical considerations for the analysis of medium-sized invertebrate genomes.

Molecular evolutionary trends and feeding ecology diversification in the Hemiptera, anchored by the milkweed bug genome

BACKGROUND

The Hemiptera (aphids, cicadas, and true bugs) are a key insect order, with high diversity for feeding ecology and excellent experimental tractability for molecular genetics. Building upon recent sequencing of hemipteran pests such as phloem-feeding aphids and blood-feeding bed bugs, we present the genome sequence and comparative analyses centered on the milkweed bug Oncopeltus fasciatus, a seed feeder of the family Lygaeidae.

RESULTS

The 926-Mb Oncopeltus genome is well represented by the current assembly and official gene set. We use our genomic and RNA-seq data not only to characterize the protein-coding gene repertoire and perform isoform-specific RNAi, but also to elucidate patterns of molecular evolution and physiology. We find ongoing, lineage-specific expansion and diversification of repressive C2H2 zinc finger proteins. The discovery of intron gain and turnover specific to the Hemiptera also prompted the evaluation of lineage and genome size as predictors of gene structure evolution. Furthermore, we identify enzymatic gains and losses that correlate with feeding biology, particularly for reductions associated with derived, fluid nutrition feeding.

CONCLUSIONS

With the milkweed bug, we now have a critical mass of sequenced species for a hemimetabolous insect order and close outgroup to the Holometabola, substantially improving the diversity of insect genomics. We thereby define commonalities among the Hemiptera and delve into how hemipteran genomes reflect distinct feeding ecologies. Given Oncopeltus's strength as an experimental model, these new sequence resources bolster the foundation for molecular research and highlight technical considerations for the analysis of medium-sized invertebrate genomes.

Bacteriocyte Reprogramming to Cope With Nutritional Stress in a Phloem Sap Feeding Hemipteran, the Pea Aphid Acyrthosiphon pisum

Nutritional symbioses play a central role in the ability of insects to thrive on unbalanced diets and in ensuring their evolutionary success. A genomic model for nutritional symbiosis comprises the hemipteran Acyrthosiphon pisum, and the gamma-3-proteobacterium, Buchnera aphidicola, with genomes encoding highly integrated metabolic pathways. A. pisum feeds exclusively on plant phloem sap, a nutritionally unbalanced diet highly variable in composition, thus raising the question of how this symbiotic system responds to nutritional stress. We addressed this by combining transcriptomic, phenotypic and life history trait analyses to determine the organismal impact of deprivation of tyrosine and phenylalanine. These two aromatic amino acids are essential for aphid development, are synthesized in a metabolic pathway for which the aphid host and the endosymbiont are interdependent, and their concentration can be highly variable in plant phloem sap. We found that this nutritional challenge does not have major phenotypic effects on the pea aphid, except for a limited weight reduction and a 2-day delay in onset of nymph laying. Transcriptomic analyses through aphid development showed a prominent response in bacteriocytes (the core symbiotic tissue which houses the symbionts), but not in gut, thus highlighting the role of bacteriocytes as major modulators of this homeostasis. This response does not involve a direct regulation of tyrosine and phenylalanine biosynthetic pathway and transporter genes. Instead, we observed an extensive transcriptional reprogramming of the bacteriocyte with a rapid down-regulation of genes encoding sugar transporters and genes required for sugar metabolism. Consistently, we observed continued overexpression of the A. pisum homolog of RRAD, a small GTPase implicated in repressing aerobic glycolysis. In addition, we found increased transcription of genes involved in proliferation, cell size control and signaling. We experimentally confirmed the significance of these gene expression changes detecting an increase in bacteriocyte number and cell size in vivo under tyrosine and phenylalanine depletion. Our results support a central role of bacteriocytes in the aphid response to amino acid deprivation: their transcriptional and cellular responses fine-tune host physiology providing the host insect with an effective way to cope with the challenges posed by the variability in composition of phloem sap.

Elucidating the effect of the lead iodide complexation degree behind the morphology and performance of perovskite solar cells

The inclusion of iodide additives in hybrid perovskite precursor solutions has been successfully exploited to improve the solar cell efficiency but their impact on perovskite formation, morphology and photovoltaic performance is still not clear. Here an extensive analysis of the effect of iodide additives in the solution-phase and during the perovskite film formation, as well as their effect on device performance is provided. The results demonstrate that in the solution-phase the additives promote the formation of lead poly-iodide species resulting in the disaggregation of the inorganic lead iodide framework and in the formation of smaller nuclei inducing the growth of uniform and smooth perovskite films. Most importantly, the complexation capability of different iodide additives does not only directly affect film morphology but also influences the density of defect states by varying the stoichiometry of precursors. These findings demonstrate that the fine control of the interactions of the chemical species in the solution-phase is essential for the precise control of the morphology at the nanoscale and the growth of the perovskite films with a reduced density of defect states. Therefore, the in-depth understanding of all the processes involved in the solution-phase is the first step for the development of a facile and reproducible approach for the fabrication of hybrid perovskite solar cells with enhanced photovoltaic performance.

Rapid transcriptional plasticity of duplicated gene clusters enables a clonally reproducing aphid to colonise diverse plant species

Background

The prevailing paradigm of host-parasite evolution is that arms races lead to increasing specialisation via genetic adaptation. Insect herbivores are no exception and the majority have evolved to colonise a small number of closely related host species. Remarkably, the green peach aphid, Myzus persicae, colonises plant species across 40 families and single M. persicae clonal lineages can colonise distantly related plants. This remarkable ability makes M. persicae a highly destructive pest of many important crop species.

Results

To investigate the exceptional phenotypic plasticity of M. persicae, we sequenced the M. persicae genome and assessed how one clonal lineage responds to host plant species of different families. We show that genetically identical individuals are able to colonise distantly related host species through the differential regulation of genes belonging to aphid-expanded gene families. Multigene clusters collectively upregulate in single aphids within two days upon host switch. Furthermore, we demonstrate the functional significance of this rapid transcriptional change using RNA interference (RNAi)-mediated knock-down of genes belonging to the cathepsin B gene family. Knock-down of cathepsin B genes reduced aphid fitness, but only on the host that induced upregulation of these genes.

Conclusions

Previous research has focused on the role of genetic adaptation of parasites to their hosts. Here we show that the generalist aphid pest M. persicae is able to colonise diverse host plant species in the absence of genetic specialisation. This is achieved through rapid transcriptional plasticity of genes that have duplicated during aphid evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-016-1145-3) contains supplementary material, which is available to authorized users.

An unusual thoracic localizations of Erdheim- Chester disease: A case report

Erdheim- Chester disease is a rare non- Langerhans cell histiocytosis that usually involves the bones, heart, central nervous system, retroperitoneum, eyes, kidneys, skin and adrenals. Lungs are affected in up to one-half cases; at CT scan various patterns are described: interstitial disease, consolidations, micronodules and microcysts, with or without pleural involvement. We presented a case of a 59 year-old man with unusual intrathoracic manifestation of Erdheim- Chester disease. Singularities of our report are the lonely thoracic involvement at the onset of the disease and a histiocytic lesion in the posterior mediastinum.

Cytokine modulation in patients with idiopathic pulmonary fibrosis undergoing treatment with steroids, immunosuppressants, and IFN-γ 1b

Idiopathic pulmonary fibrosis (IPF) is a chronic lung disease of unknown etiology and pathogenic mechanisms. From an etiopathogenic point of view, alveolar macrophages play a key role in accumulation of fibroblasts and deposition of collagen and extracellular matrix by releasing specific cytokines and inflammatory mediators. IPF seems to be also associated with circulating fibrocytes, which might be involved with an abnormal pulmonary vascular repair and remodeling. Based on its hypothesized pathologic mechanisms, anti-inflammatory, anti-fibrotic and immunosuppressive therapies are often used. For these reasons, Interferon-g (IFN-g) has been used to exploit its activity on macrophages and fibroblasts. The aim of this study was to investigate the response to corticosteroids and/or IFN-g 1b treatments based on pulmonary function tests and on inflammatory cytokine patterns of expression on bronchoalveolar lavage (BAL), at baseline and during and after the therapies. Unlike previous studies, we analyzed a period of therapy longer than 1 year. Our results demonstrated the effectiveness of IFN-γ in a group of IPF patients in whom the treatment was prolonged for over a year. These data suggest a positive role of IFN-γ; treatment in patients in the initial stage of the disease.

Absorption correction and phase function shape effects on the closure of apparent optical properties

We present a closure experiment between new inherent optical properties (IOPs: absorption a, scattering b, backscattering b_b) and apparent optical properties (AOPs: remote-sensing reflectance R_rs, irradiance reflectance R, and anisotropic factor at nadir Q_n) data of Ionian and Adriatic seawaters, from very clear to turbid waters, ranging across one order of magnitude in R_rs. The internal consistency of the IOP-AOP matchups was investigated though radiative transfer closure. Using the in situ IOPs, we predicted the AOPs with the commercial radiative transfer solver Hydrolight. Closure was limited by two unresolved issues, one regarding processing of in situ data and the other related to radiative transfer modeling. First, different correction methods of the absorption data measured by the Wetlabs ac-s produced high variations in simulated reflectances, reaching 40% for the highest reflectances in our dataset. Second, the lack of detailed volume scattering function measurements forces us to adopt analytical functions that are consistent with a given particle backscattering ratio. The analytical phase functions named Fournier-Forand and two-term Kopelevich presented reasonable angular shapes with respect to measurements at a few backward angles. Between these phase functions, induced changes were within 4% for R_rs, within 11% for R, and within 10% for Q_n. Additionally, closure of Q_n was generally not successful considering radiometric uncertainties. Simulated Q_n overestimated low values and underestimated high values, especially at 665 nm, where Hydrolight was unable to predict measured Q_n values greater than 6 sr. The physical nature of Q_n makes this mismatch almost independent of the measured IOPs, thus precluding Q_n tuning by varying the former. The non-closure of Q_n might be caused by an inaccurate phase function and, to a lesser extent, by the modeling of the incoming radiance. For the future, this remains the task of accurate absorption and phase function measurements, especially at red wavelengths.

Disruption of phenylalanine hydroxylase reduces adult lifespan and fecundity, and impairs embryonic development in parthenogenetic pea aphids

Phenylalanine hydroxylase (PAH) is a key tyrosine-biosynthetic enzyme involved in neurological and melanin-associated physiological processes. Despite extensive investigations in holometabolous insects, a PAH contribution to insect embryonic development has never been demonstrated. Here, we have characterized, for the first time, the PAH gene in a hemimetabolous insect, the aphid Acyrthosiphon pisum. Phylogenetic and sequence analyses confirmed that ApPAH is closely related to metazoan PAH, exhibiting the typical ACT regulatory and catalytic domains. Temporal expression patterns suggest that ApPAH has an important role in aphid developmental physiology, its mRNA levels peaking at the end of embryonic development. We used parental dsApPAH treatment to generate successful knockdown in aphid embryos and to study its developmental role. ApPAH inactivation shortens the adult aphid lifespan and considerably affects fecundity by diminishing the number of nymphs laid and impairing embryonic development, with newborn nymphs exhibiting severe morphological defects. Using single nymph HPLC analyses, we demonstrated a significant tyrosine deficiency and a consistent accumulation of the upstream tyrosine precursor, phenylalanine, in defective nymphs, thus confirming the RNAi-mediated disruption of PAH activity. This study provides first insights into the role of PAH in hemimetabolous insects and demonstrates that this metabolic gene is essential for insect embryonic development.

Endobronchial-ultrasound needle aspiration and endoscopic ultrasound-fine-needle aspiration in thoracic diseases

EBUS-TBNA and EUS-FNA are minimally invasive techniques rapidly gaining ground in the non-surgical invasive diagnostic approach to thoracic diseases due to their high accuracy and low morbidity and mortality compared to surgical techniques. Moreover, in the diagnosis and staging of lung cancer the combination of the two techniques is superior to either test alone. In this review we focus on the role of EBUS-TBNA and EUS-FNA in both malignant and non-malignant thoracic diseases.

ArthropodaCyc: a CycADS powered collection of BioCyc databases to analyse and compare metabolism of arthropods

Arthropods interact with humans at different levels with highly beneficial roles (e.g. as pollinators), as well as with a negative impact for example as vectors of human or animal diseases, or as agricultural pests. Several arthropod genomes are available at present and many others will be sequenced in the near future in the context of the i5K initiative, offering opportunities for reconstructing, modelling and comparing their metabolic networks. In-depth analysis of these genomic data through metabolism reconstruction is expected to contribute to a better understanding of the biology of arthropods, thereby allowing the development of new strategies to control harmful species. In this context, we present here ArthropodaCyc, a dedicated BioCyc collection of databases using the Cyc annotation database system (CycADS), allowing researchers to perform reliable metabolism comparisons of fully sequenced arthropods genomes. Since the annotation quality is a key factor when performing such global genome comparisons, all proteins from the genomes included in the ArthropodaCyc database were re-annotated using several annotation tools and orthology information. All functional/domain annotation results and their sources were integrated in the databases for user access. Currently, ArthropodaCyc offers a centralized repository of metabolic pathways, protein sequence domains, Gene Ontology annotations as well as evolutionary information for 28 arthropod species. Such database collection allows metabolism analysis both with integrated tools and through extraction of data in formats suitable for systems biology studies.

Database URL:http://arthropodacyc.cycadsys.org/

New insight into the RNA interference response against cathepsin-L gene in the pea aphid, Acyrthosiphon pisum: molting or gut phenotypes specifically induced by injection or feeding treatments

RNA interference (RNAi) has been widely and successfully used for gene inactivation in insects, including aphids, where dsRNA administration can be performed either by feeding or microinjection. However, several aspects related to the aphid response to RNAi, as well as the influence of the administration method on tissue response, or the mixed success to observe phenotypes specific to the gene targeted, are still unclear in this insect group. In the present study, we made the first direct comparison of two administration methods (injection or feeding) for delivery of dsRNA targeting the cathepsin-L gene in the pea aphid, Acyrthosiphon pisum. In order to maximize the possibility of discovering specific phenotypes, the effect of the treatment was analyzed in single individual aphids at the level of five body compartments: the bacteriocytes, the gut, the embryonic chains, the head and the remaining body carcass. Our analysis revealed that gene expression knockdown effect in each single body compartment was dependent on the administration method used, and allowed us to discover new functions for the cathepsin-L gene in aphids. Injection of cathepsin-L dsRNA was much more effective on carcass and head, inducing body morphology alterations, and suggesting a novel role of this gene in the molting of these insects. Administration by feeding provoked cathepsin-L knockdown in the gut and specific gut epithelial cell alteration, therefore allowing a better characterization of tissue specific role of this gene in aphids.

Hierarchical multiple peeling simulations

We model hierarchical multiple peeling through finite element simulations, quantifying how hierarchy optimizes the strength of the bioinspired attachment.

Tyrosine pathway regulation is host-mediated in the pea aphid symbiosis during late embryonic and early larval development

BACKGROUND

Nutritional symbioses play a central role in insects' adaptation to specialized diets and in their evolutionary success. The obligatory symbiosis between the pea aphid, Acyrthosiphon pisum, and the bacterium, Buchnera aphidicola, is no exception as it enables this important agricultural pest insect to develop on a diet exclusively based on plant phloem sap. The symbiotic bacteria provide the host with essential amino acids lacking in its diet but necessary for the rapid embryonic growth seen in the parthenogenetic viviparous reproduction of aphids. The aphid furnishes, in exchange, non-essential amino acids and other important metabolites. Understanding the regulations acting on this integrated metabolic system during the development of this insect is essential in elucidating aphid biology.

RESULTS

We used a microarray-based approach to analyse gene expression in the late embryonic and the early larval stages of the pea aphid, characterizing, for the first time, the transcriptional profiles in these developmental phases. Our analyses allowed us to identify key genes in the phenylalanine, tyrosine and dopamine pathways and we identified ACYPI004243, one of the four genes encoding for the aspartate transaminase (E.C. 2.6.1.1), as specifically regulated during development. Indeed, the tyrosine biosynthetic pathway is crucial for the symbiotic metabolism as it is shared between the two partners, all the precursors being produced by B. aphidicola. Our microarray data are supported by HPLC amino acid analyses demonstrating an accumulation of tyrosine at the same developmental stages, with an up-regulation of the tyrosine biosynthetic genes. Tyrosine is also essential for the synthesis of cuticular proteins and it is an important precursor for cuticle maturation: together with the up-regulation of tyrosine biosynthesis, we observed an up-regulation of cuticular genes expression. We were also able to identify some amino acid transporter genes which are essential for the switch over to the late embryonic stages in pea aphid development.

CONCLUSIONS

Our data show that, in the development of A. pisum, a specific host gene set regulates the biosynthetic pathways of amino acids, demonstrating how the regulation of gene expression enables an insect to control the production of metabolites crucial for its own development and symbiotic metabolism.

Multimodal dynamic response of the Buchnera aphidicola pLeu plasmid to variations in leucine demand of its host, the pea aphid Acyrthosiphon pisum

Aphids, important agricultural pests, can grow and reproduce thanks to their intimate symbiosis with the γ-proteobacterium Buchnera aphidicola that furnishes them with essential amino acids lacking in their phloem sap diet. To study how B. aphidicola, with its reduced genome containing very few transcriptional regulators, responds to variations in the metabolic requirements of its host, we concentrated on the leucine metabolic pathway. We show that leucine is a limiting factor for aphid growth and it displays a stimulatory feeding effect. Our metabolic analyses demonstrate that symbiotic aphids are able to respond to leucine starvation or excess by modulating the neosynthesis of this amino acid. At a molecular level, this response involves an early important transcriptional regulation (after 12 h of treatment) followed by a moderate change in the pLeu plasmid copy number. Both responses are no longer apparent after 7 days of treatment. These experimental data are discussed in the light of a re-annotation of the pLeu plasmid regulatory elements. Taken together, our data show that the response of B. aphidicola to the leucine demand of its host is multimodal and dynamically regulated, providing new insights concerning the genetic regulation capabilities of this bacterium in relation to its symbiotic functions.

CycADS: an annotation database system to ease the development and update of BioCyc databases

In recent years, genomes from an increasing number of organisms have been sequenced, but their annotation remains a time-consuming process. The BioCyc databases offer a framework for the integrated analysis of metabolic networks. The Pathway tool software suite allows the automated construction of a database starting from an annotated genome, but it requires prior integration of all annotations into a specific summary file or into a GenBank file. To allow the easy creation and update of a BioCyc database starting from the multiple genome annotation resources available over time, we have developed an ad hoc data management system that we called Cyc Annotation Database System (CycADS). CycADS is centred on a specific database model and on a set of Java programs to import, filter and export relevant information. Data from GenBank and other annotation sources (including for example: KAAS, PRIAM, Blast2GO and PhylomeDB) are collected into a database to be subsequently filtered and extracted to generate a complete annotation file. This file is then used to build an enriched BioCyc database using the PathoLogic program of Pathway Tools. The CycADS pipeline for annotation management was used to build the AcypiCyc database for the pea aphid (Acyrthosiphon pisum) whose genome was recently sequenced. The AcypiCyc database webpage includes also, for comparative analyses, two other metabolic reconstruction BioCyc databases generated using CycADS: TricaCyc for Tribolium castaneum and DromeCyc for Drosophila melanogaster. Linked to its flexible design, CycADS offers a powerful software tool for the generation and regular updating of enriched BioCyc databases. The CycADS system is particularly suited for metabolic gene annotation and network reconstruction in newly sequenced genomes. Because of the uniform annotation used for metabolic network reconstruction, CycADS is particularly useful for comparative analysis of the metabolism of different organisms.

Database URL:http://www.cycadsys.org

A statistical approach for detecting genomic aberrations in heterogeneous tumor samples from single nucleotide polymorphism genotyping data

We describe a statistical method for the characterization of genomic aberrations in single nucleotide polymorphism microarray data acquired from cancer genomes. Our approach allows us to model the joint effect of polyploidy, normal DNA contamination and intra-tumour heterogeneity within a single unified Bayesian framework. We demonstrate the efficacy of our method on numerous datasets including laboratory generated mixtures of normal-cancer cell lines and real primary tumours.

Altered intra-nuclear organisation of heterochromatin and genes in ICF syndrome

The ICF syndrome is a rare autosomal recessive disorder, the most common symptoms of which are immunodeficiency, facial anomalies and cytogenetic defects involving decondensation and instability of chromosome 1, 9 and 16 centromeric regions. ICF is also characterised by significant hypomethylation of the classical satellite DNA, the major constituent of the juxtacentromeric heterochromatin. Here we report the first attempt at analysing some of the defining genetic and epigenetic changes of this syndrome from a nuclear architecture perspective. In particular, we have compared in ICF (Type 1 and Type 2) and controls the large-scale organisation of chromosome 1 and 16 juxtacentromeric heterochromatic regions, their intra-nuclear positioning, and co-localisation with five specific genes (BTG2, CNN3, ID3, RGS1, F13A1), on which we have concurrently conducted expression and methylation analysis. Our investigations, carried out by a combination of molecular and cytological techniques, demonstrate the existence of specific and quantifiable differences in the genomic and nuclear organisation of the juxtacentromeric heterochromatin in ICF. DNA hypomethylation, previously reported to correlate with the decondensation of centromeric regions in metaphase described in these patients, appears also to correlate with the heterochromatin spatial configuration in interphase. Finally, our findings on the relative positioning of hypomethylated satellite sequences and abnormally expressed genes suggest a connection between disruption of long-range gene-heterochromatin associations and some of the changes in gene expression in ICF. Beyond its relevance to the ICF syndrome, by addressing fundamental principles of chromosome functional organisation within the cell nucleus, this work aims to contribute to the current debate on the epigenetic impact of nuclear architecture in development and disease.

Systemic analysis of the symbiotic function of Buchnera aphidicola, the primary endosymbiont of the pea aphid Acyrthosiphon pisum

Buchnera aphidicola is the primary obligate intracellular symbiont of most aphid species. B. aphidicola and aphids have been evolving in parallel since their association started, about 150 Myr ago. Both partners have lost their autonomy, and aphid diversification has been confined to smaller ecological niches by this co-evolution. B. aphidicola has undergone major genomic and biochemical changes as a result of adapting to intracellular life. Several genomes of B. aphidicola from different aphid species have been sequenced in the last decade, making it possible to carry out analyses and comparative studies using system-level in silico methods. This review attempts to provide a systemic description of the symbiotic function of aphid endosymbionts, particularly of B. aphidicola from the pea aphid Acyrthosiphon pisum, by analyzing their structural genomic properties, as well as their genetic and metabolic networks.

Genome-wide hypomethylation in head and neck cancer is more pronounced in HPV-negative tumors and is associated with genomic instability

Loss of genome-wide methylation is a common feature of cancer, and the degree of hypomethylation has been correlated with genomic instability. Global methylation of repetitive elements possibly arose as a defense mechanism against parasitic DNA elements, including retrotransposons and viral pathogens. Given the alterations of global methylation in both viral infection and cancer, we examined genome-wide methylation levels in head and neck squamous cell carcinoma (HNSCC), a cancer causally associated with human papilloma virus (HPV). We assayed global hypomethylation levels in 26 HNSCC samples, compared with their matched normal adjacent tissue, using Pyrosequencing-based methylation assays for LINE repeats. In addition, we examined cell lines derived from a variety of solid tumors for LINE and SINE (Alu) repeats. The degree of LINE and Alu hypomethylation varied among different cancer cell lines. There was only moderate correlation between LINE and Alu methylation levels, with the range of variation in methylation levels being greater for the LINE elements. LINE hypomethylation was more pronounced in HPV-negative than in HPV-positive tumors. Moreover, genomic instability, as measured by genome-wide loss-of-heterozygosity (LOH) single nucleotide polymorphism (SNP) analysis, was greater in HNSCC samples with more pronounced LINE hypomethylation. Global hypomethylation was variable in HNSCC. Its correlation with both HPV status and degree of LOH as a surrogate for genomic instability may reflect alternative oncogenic pathways in HPV-positive versus HPV-negative tumors.

Molecular signatures of metastasis in head and neck cancer

BACKGROUND

Metastases are the primary cause of cancer treatment failure and death, yet metastatic mechanisms remain incompletely understood.

METHODS

We studied the molecular basis of head and neck cancer metastasis by transcriptionally profiling 70 samples from 27 patients-matching normal adjacent tissue, primary tumor, and cervical lymph node metastases.

RESULTS

We identified tumor-associated expression signatures common to both primary tumors and metastases. Use of matching metastases revealed an additional 46 dysregulated genes associated solely with head and neck cancer metastasis. However, despite being metastasis-specific in our sample set, these 46 genes are concordant with genes previously discovered in primary tumors that metastasized.

CONCLUSIONS

Although our data and related studies show that most of the metastatic potential appears to be inherent to the primary tumor, they are also consistent with the notion that a limited number of additional clonal changes are necessary to yield the final metastatic cell(s), albeit in a variable temporal order.

GenoSNP: a variational Bayes within-sample SNP genotyping algorithm that does not require a reference population

Current genotyping algorithms typically call genotypes by clustering allele-specific intensity data on a single nucleotide polymorphism (SNP) by SNP basis. This approach assumes the availability of a large number of control samples that have been sampled on the same array and platform. We have developed a SNP genotyping algorithm for the Illumina Infinium SNP genotyping assay that is entirely within-sample and does not require the need for a population of control samples nor parameters derived from such a population. Our algorithm exhibits high concordance with current methods and >99% call accuracy on HapMap samples. The ability to call genotypes using only within-sample information makes the method computationally light and practical for studies involving small sample sizes and provides a valuable independent quality control metric for other population-based approaches.

AVAILABILITY

http://www.stats.ox.ac.uk/~giannoul/GenoSNP/.

hsa-miR-210 Is induced by hypoxia and is an independent prognostic factor in breast cancer

PURPOSE

MicroRNA (miRNA) expression alterations have been described in cancer. Many cancers are characterized by areas of hypoxia, enhanced hypoxia-inducible factor (HIF) levels, and increased expression of hypoxically regulated genes, all of which correlate with patient outcome. We examined hypoxia-induced miRNA expression changes to identify markers of survival in breast cancer.

EXPERIMENTAL DESIGN

We used microarrays to analyze miRNA expression changes induced by hypoxia in MCF7 breast cancer cell lines and validated results by quantitative-PCR (Q-PCR). Small interfering RNA against HIF-1alpha and HIF-2alpha, and RCC4 cells transfected with the von Hippel-Lindau (VHL) protein were used to investigate HIF dependency of miRNA expression. miRNA Q-PCR assays were done on 219 early breast cancer samples with long-term follow-up. Correlation of expression with clinical variables was done using Pearson and Spearman's rank tests, univariate, and Cox multivariate analysis.

RESULTS

hsa-miR-210 induction was the most significant change under hypoxia by microarray analysis (3.4-fold, P < 0.001). hsa-miR-210 expression changes were validated by Q-PCR and detected in other cancer cell lines. Using small interfering RNAs and RCC4 cells transfected with VHL, we showed that the regulation by hypoxia of hsa-miR-210 was mediated by the HIF-1alpha/VHL transcriptional system but not HIF-2alpha. hsa-miR-210 expression levels in breast cancer samples correlated directly with a hypoxia score based on the expression of 99 genes. hsa-miR-210 expression levels showed an inverse correlation with disease-free and overall survival, significant in both univariate and multivariate analyses.

CONCLUSIONS

We show that hsa-miR-210 overexpression is induced by hypoxia in a HIF-1alpha- and VHL-dependent fashion and its expression levels in breast cancer samples are an independent prognostic factor.

LRRTM1 on chromosome 2p12 is a maternally suppressed gene that is associated paternally with handedness and schizophrenia

Left-right asymmetrical brain function underlies much of human cognition, behavior and emotion. Abnormalities of cerebral asymmetry are associated with schizophrenia and other neuropsychiatric disorders. The molecular, developmental and evolutionary origins of human brain asymmetry are unknown. We found significant association of a haplotype upstream of the gene LRRTM1 (Leucine-rich repeat transmembrane neuronal 1) with a quantitative measure of human handedness in a set of dyslexic siblings, when the haplotype was inherited paternally (P=0.00002). While we were unable to find this effect in an epidemiological set of twin-based sibships, we did find that the same haplotype is overtransmitted paternally to individuals with schizophrenia/schizoaffective disorder in a study of 1002 affected families (P=0.0014). We then found direct confirmatory evidence that LRRTM1 is an imprinted gene in humans that shows a variable pattern of maternal downregulation. We also showed that LRRTM1 is expressed during the development of specific forebrain structures, and thus could influence neuronal differentiation and connectivity. This is the first potential genetic influence on human handedness to be identified, and the first putative genetic effect on variability in human brain asymmetry. LRRTM1 is a candidate gene for involvement in several common neurodevelopmental disorders, and may have played a role in human cognitive and behavioral evolution.

QuantiSNP: an Objective Bayes Hidden-Markov Model to detect and accurately map copy number variation using SNP genotyping data

Array-based technologies have been used to detect chromosomal copy number changes (aneuploidies) in the human genome. Recent studies identified numerous copy number variants (CNV) and some are common polymorphisms that may contribute to disease susceptibility. We developed, and experimentally validated, a novel computational framework (QuantiSNP) for detecting regions of copy number variation from BeadArray SNP genotyping data using an Objective Bayes Hidden-Markov Model (OB-HMM). Objective Bayes measures are used to set certain hyperparameters in the priors using a novel re-sampling framework to calibrate the model to a fixed Type I (false positive) error rate. Other parameters are set via maximum marginal likelihood to prior training data of known structure. QuantiSNP provides probabilistic quantification of state classifications and significantly improves the accuracy of segmental aneuploidy identification and mapping, relative to existing analytical tools (Beadstudio, Illumina), as demonstrated by validation of breakpoint boundaries. QuantiSNP identified both novel and validated CNVs. QuantiSNP was developed using BeadArray SNP data but it can be adapted to other platforms and we believe that the OB-HMM framework has widespread applicability in genomic research. In conclusion, QuantiSNP is a novel algorithm for high-resolution CNV/aneuploidy detection with application to clinical genetics, cancer and disease association studies.

Matrix-assisted laser desorption/ionisation, time-of-flight mass spectrometry in genomics research

The beginning of this millennium has seen dramatic advances in genomic research. Milestones such as the complete sequencing of the human genome and of many other species were achieved and complemented by the systematic discovery of variation at the single nucleotide (SNP) and whole segment (copy number polymorphism) level. Currently most genomics research efforts are concentrated on the production of whole genome functional annotations, as well as on mapping the epigenome by identifying the methylation status of CpGs, mainly in CpG islands, in different tissues. These recent advances have a major impact on the way genetic research is conducted and have accelerated the discovery of genetic factors contributing to disease. Technology was the critical driving force behind genomics projects: both the combination of Sanger sequencing with high-throughput capillary electrophoresis and the rapid advances in microarray technologies were keys to success. MALDI-TOF MS–based genome analysis represents a relative newcomer in this field. Can it establish itself as a long-term contributor to genetics research, or is it only suitable for niche areas and for laboratories with a passion for mass spectrometry? In this review, we will highlight the potential of MALDI-TOF MS–based tools for resequencing and for epigenetics research applications, as well as for classical complex genetic studies, allele quantification, and quantitative gene expression analysis. We will also identify the current limitations of this approach and attempt to place it in the context of other genome analysis technologies.

alpha-Haemoglobin stabilising protein is a quantitative trait gene that modifies the phenotype of beta-thalassaemia

It has been suggested that altered levels or function of alpha-haemoglobin stabilising protein (AHSP), an erythroid-specific protein that binds specifically to free alpha-(haemo)globin, might account for some of the clinical variability in beta-thalassaemia. To assess the variation of AHSP expression, mRNA levels in circulating reticulocytes of 103 healthy individuals were measured by quantitative reverse transcription-polymerase chain reaction. AHSP expression varied up to threefold, and did not correlate with age or sex. A systematic survey of the AHSP locus identified eight sequence variants, of which six were common. Four common variants, including the longer homopolymer (T18) in the putative promoter, are strongly associated with AHSP expression. Reporter assays in K562 cells showed that the activity of the shorter (T15) reporter was relatively lower than that of the T18 reporter. In a study of nine anaemic patients who were heterozygous for beta-thalassaemia and also heterozygous for the triplicated alpha-globin gene (alpha alpha alpha/alpha alpha), frequency of the shorter homopolymer was higher than expected. AHSP expression is variable, with cis control accounting for some of its variance. In some families, the subtle altered levels in AHSP related to the AHSP genotype appears to be a relevant contributory factor in the haematological phenotype.

N-(4-Hydroxyphenyl)retinamide and nitric oxide pro-drugs exhibit apoptotic and anti-invasive effects against bone metastatic breast cancer cells

Breast cancer most frequently metastasizes to bone causing decreased quality of life and morbidity. Since current treatments are palliative, strategies to prevent bone metastases in breast cancer patients are required. There is substantial evidence indicating that high levels of nitric oxide (NO) suppress tumor growth and metastasis in vivo. We hypothesize that agents that produce high concentrations of NO could prevent the spread of breast cancer to bone. We previously demonstrated that the synthetic retinoid N-(4-hydroxyphenyl)retinamide (4-HPR) produces high levels of NO via the induction of NO synthases. NO pro-drugs are designed to produce large amounts of NO without inducing NO synthases but upon metabolism by their intracellular targets. The objective of this study was to determine the effectiveness of 4-HPR and an NO pro-drug, diethylamineNONOate/AM (NONO-AM), in inhibiting the growth and invasiveness of bone metastatic breast cancer cells. Parental MDA-MB-231 breast cancer cells were resistant to 4-HPR-induced apoptosis at clinically relevant doses, whereas 4-HPR-induced apoptosis in a dose-dependent manner in MDA-MB-231/F10 bone metastatic breast cancer cells. Unlike 4-HPR, NONO-AM induced apoptosis in a dose-dependent manner in both parental MDA-MB-231 cells and F10 cells. The bone metastatic F10 cells were more sensitive to the anti-invasive effects of 4-HPR and NONO-AM than were MDA-MB-231 cells. Although suppression of matrix metalloprotease-9 activity may be one mechanism by which 4-HPR decreases the invasion of F10 cells, it does not appear to be the anti-invasion mechanism of NONO-AM. These in vitro results suggest that 4-HPR and NO pro-drugs may be effective chemopreventive agents against bone metastatic breast cancer.

Genome-wide loss of heterozygosity analysis of WT1-wild-type and WT1-mutant Wilms tumors

Wilms tumor (WT) is genetically heterogeneous, and the one known WT gene, WT1 at 11p13, is altered in only a subset of WTs. Previous loss of heterozygosity (LOH) analyses have revealed the existence of additional putative WT genes at 11p15, 16q, and 1p, but these analyses examined only one or a handful of chromosomes or looked at LOH at only a few markers per chromosome. We conducted a genome-wide scan for LOH in WT by using 420 markers spaced at an average of 10 cM throughout the genome and analyzed the data for two genetically defined subsets of WTs: those with mutations in WT1 and those with no detectable WT1 alteration. Our findings indicated that the incidence of LOH throughout the genome was significantly lower in our group of WTs with WT1 mutations. In WT1-wild-type tumors, we observed the expected LOH at 11p, 16q, and 1p, and, in addition, we localized a previously unobserved region of LOH at 9q. Using additional 9q markers within this region of interest, we sublocalized the region of 9q LOH to the 12.2 Mb between D9S283 and a simple tandem repeat in BAC RP11-177I8, a region containing several potential tumor-suppressor genes. As a result, we have established for the first time that WT1-mutant and WT1-wild-type WTs differ significantly in their patterns of LOH throughout the genome, suggesting that the genomic regions showing LOH in WT1-wild-type tumors harbor genes whose expression is regulated by the pleiotropic effects of WT1. Our results implicate 9q22.2-q31.1 as a region containing such a gene.

Reduced expression of the Aalpha subunit of protein phosphatase 2A in human gliomas in the absence of mutations in the Aalpha and Abeta subunit genes

Protein phosphatase 2A (PP2A) consists of 3 subunits: the catalytic subunit, C, and the regulatory subunits, A and B. The A and C subunits both exist as 2 isoforms (alpha and beta) and the B subunit as multiple forms subdivided into 3 families, B, B' and B". It has been reported that the genes encoding the Aalpha and Abeta subunits are mutated in various human cancers, suggesting that they may function as tumor suppressors. We investigated whether Aalpha and Abeta mutations occur in human gliomas. Using single strand conformational polymorphism analysis and DNA sequencing, 58 brain tumors were investigated, including 23 glioblastomas, 19 oligodendrogliomas and 16 anaplastic oligodendrogliomas. Only silent mutations were detected in the Aalpha gene and no mutations in the Abeta gene. However, in 43% of the tumors, the level of Aalpha was reduced at least 10-fold. By comparison, the levels of the Balpha and Calpha subunits were mostly normal. Our data indicate that these tumors contain very low levels of core and holoenzyme and high amounts of unregulated catalytic C subunit.

Impaired insulin-like growth factor I vasorelaxant effects in hypertension

Insulin-like growth factor I (IGF-I) can be considered a factor potentially involved in arterial hypertension not only for its growth-promoting features but also for its effects on vascular tone. Nevertheless, the actions of the hormone on vascular reactivity are still unexplored in hypertension. Therefore, the vasodilation induced by increasing doses of IGF-I and the modulation of norepinephrine vasoconstriction induced by low levels of the hormone were tested on aortic rings of spontaneously hypertensive and normotensive rats. The results indicate that the vasodilation evoked by IGF-I is impaired in hypertensive rats (Delta% of maximal vasorelaxation, 30+/-1 versus 41+/-1; P<0.01), and after the removal of endothelium or the inhibition of endothelial NO synthase, the vasodilation evoked by the hormone was blunted in both rat strains and became similar between hypertensive and normotensive rats (Delta% of maximal vasorelaxation, 21+/-1 versus 20+/-1; P=NS). Moreover, IGF-I does not show any effect on norepinephrine vasoconstriction in hypertensive rats, and this alteration may depend on the lack of sensitizing effect exerted by IGF-I on alpha(2)-adrenergic-evoked NO vasorelaxation. The defect in IGF-I vascular action is also present in young spontaneously hypertensive rats (age 5 weeks). In conclusion, our data demonstrate that IGF-I vasorelaxant properties are impaired in spontaneously hypertensive rats, suggesting that such defect may play a causative or permissive role in the development of hypertensive conditions.

Mutation analysis of hBUB1, hBUBR1 and hBUB3 genes in glioblastomas

Glioblastomas, the most malignant human brain tumors, are characterized by marked aneuploidy, suggesting chromosomal instability which may be caused by a defective mitotic spindle checkpoint. We screened 22 glioblastomas for mutations in the mitotic spindle check-point genes hBUB1, hBUBR1 and hBUB3. DNA sequencing revealed a silent mutation at codon 144 of hBUB1 (CAG-->CAA, Gln-->Gln) in one glioblastoma, a silent mutation at codon 952 of hBUBR1 (GAC-->GAT, Asp-->Asp) in another glioblastoma, and a silent mutation at codon 388 of the hBUBR1 gene (GCG-->GCA, Ala-->Ala) in 8 glioblastomas. We also observed a known polymorphism at hBUBR1 codon 349 (CAA/CGA, Gln/Arg), with an allelic frequency of 0.75 for Gln and 0.25 for Arg, which is similar to that among healthy Caucasian individuals (0.73 vs 0.27). The coding sequence of the hBUB3 gene did not contain any mutation, but in 4 glioblastomas (18%), a C-->T point mutation was detected at position -6 (6 nucleotides upstream of the ATG initiator codon). Analysis of blood DNA of these patients showed identical sequence alterations, indicating that this is a polymorphism. Again, the frequency in glioblastomas was similar to that in healthy Caucasians (15%). We further screened hBUB1 in 18 cases of giant cell glioblastoma, a variant characterized by a predominance of bizarre, multinucleated giant cells. There were no changes, except for a silent mutation at codon 144 in two cases. These results suggest that mutations in these mitotic spindle checkpoint genes do not play a significant role in the causation of chromosomal instability in glioblastomas.