Mechanisms of osteoclastogenesis inhibition by a novel class of biphenyl-type cannabinoid CB(2) receptor inverse agonists

The cannabinoid CB(2) receptor is known to modulate osteoclast function by poorly understood mechanisms. Here, we report that the natural biphenyl neolignan 4'-O-methylhonokiol (MH) is a CB(2) receptor-selective antiosteoclastogenic lead structure (K(i) < 50 nM). Intriguingly, MH triggers a simultaneous G(i) inverse agonist response and a strong CB(2) receptor-dependent increase in intracellular calcium. The most active inverse agonists from a library of MH derivatives inhibited osteoclastogenesis in RANK ligand-stimulated RAW264.7 cells and primary human macrophages. Moreover, these ligands potently inhibited the osteoclastogenic action of endocannabinoids. Our data show that CB(2) receptor-mediated cAMP formation, but not intracellular calcium, is crucially involved in the regulation of osteoclastogenesis, primarily by inhibiting macrophage chemotaxis and TNF-α expression. MH is an easily accessible CB(2) receptor-selective scaffold that exhibits a novel type of functional heterogeneity.

Drug Hypersensitivity: How Drugs Stimulate T Cells via Pharmacological Interaction with Immune Receptors.

Small chemicals like drugs tend to bind to proteins via noncovalent bonds, e.g. hydrogen bonds, salt bridges or electrostatic interactions. Some chemicals interact with other molecules than the actual target ligand, representing so-called 'off-target' activities of drugs. Such interactions are a main cause of adverse side effects to drugs and are normally classified as predictable type A reactions. Detailed analysis of drug-induced immune reactions revealed that off-target activities also affect immune receptors, such as highly polymorphic human leukocyte antigens (HLA) or T cell receptors (TCR). Such drug interactions with immune receptors may lead to T cell stimulation, resulting in clinical symptoms of delayed-type hypersensitivity. They are assigned the 'pharmacological interaction with immune receptors' (p-i) concept. Analysis of p-i has revealed that drugs bind preferentially or exclusively to distinct HLA molecules (p-i HLA) or to distinct TCR (p-i TCR). P-i reactions differ from 'conventional' off-target drug reactions as the outcome is not due to the effect on the drug-modified cells themselves, but is the consequence of reactive T cells. Hence, the complex and diverse clinical manifestations of delayed-type hypersensitivity are caused by the functional heterogeneity of T cells. In the abacavir model of p-i HLA, the drug binding to HLA may result in alteration of the presenting peptides. More importantly, the drug binding to HLA generates a drug-modified HLA, which stimulates T cells directly, like an allo-HLA. In the sulfamethoxazole model of p-i TCR, responsive T cells likely require costimulation for full T cell activation. These findings may explain the similarity of delayed-type hypersensitivity reactions to graft-versus-host disease, and how systemic viral infections increase the risk of delayed-type hypersensitivity reactions.

IL-4 and -5 prime human mast cells for different profiles of IgE-dependent cytokine production

Mast cells (MC) are stem cell factor-dependent tissue-based hematopoietic cells with substantial functional heterogeneity. Cord blood-derived human MC (hMC) express functional receptors for IL-5, and IL-5 mediates stem cell factor-dependent comitogenesis of hMC in vitro. Although IL-5 is not required for normal hMC development, we considered that it might prime hMC for their high-affinity Fc receptor for IgE (FcɛRI)-dependent generation of cytokines, as previously demonstrated for IL-4. Compared with hMC maintained in stem cell factor alone, hMC primed with IL-5 expressed 2- to 4-fold higher steady-state levels of TNF-α, IL-5, IL-13, macrophage inflammatory protein 1α, and granulocyte-macrophage colony-stimulating factor transcripts 2 h after FcɛRI crosslinking and secreted 2- to 5-fold greater quantities of the corresponding cytokines, except IL-13, at 6 h. Unlike IL-4, IL-5 priming did not enhance FcɛRI-dependent histamine release. Thus, IL-5 augments cytokine production by hMC by a mechanism distinct from that of IL-4 and with a different resultant profile of cytokine production. These observations suggest a potentially autocrine effect of IL-5 on hMC for amplification of allergic immune responses, in addition to its recognized paracrine effects on eosinophils, and implicate both IL-4 and IL-5 in the modulation of the hMC phenotype.

Primary and secondary somatosensory cortex responses to anticipation and pain: a magnetoencephalography study

Several brain regions, including the primary and secondary somatosensory cortices (SI and SII, respectively), are functionally active during the pain experience. Both of these regions are thought to be involved in the sensory-discriminative processing of pain and recent evidence suggests that SI in particular may also be involved in more affective processing. In this study we used MEG to investigate the hypothesis that frequency-specific oscillatory activity may be differentially associated with the sensory and affective components of pain. In eight healthy participants (four male), MEG was recorded during a visceral pain experiment comprising baseline, anticipation, pain and post-pain phases. Pain was delivered via intraluminal oesophageal balloon distension (four stimuli at 1 Hz). Significant bilateral but asymmetrical changes in neural activity occurred in the beta-band within SI and SII. In SI, a continuous increase in neural activity occurred during the anticipation phase (20-30 Hz), which continued during the pain phase but at a lower frequency (10-15 Hz). In SII, oscillatory changes only occurred during the pain phase, predominantly in the 20-30 Hz beta band, and were coincident with the stimulus. These data provide novel evidence of functional diversity within SI, indicating a role in attentional and sensory aspects of pain processing. In SII, oscillatory changes were predominantly stimulus-related, indicating a role in encoding the characteristics of the stimulus. We therefore provide objective evidence of functional heterogeneity within SI and functional segregation between SI and SII, and suggest that the temporal and frequency dynamics within cortical regions may offer valuable insights into pain processing.

Radiotherapy Treatment Assessment using DCE-MRI

Abstract

The goal of modern radiotherapy is to precisely deliver a prescribed radiation dose to delineated target volumes that contain a significant amount of tumor cells while sparing the surrounding healthy tissues/organs. Precise delineation of treatment and avoidance volumes is the key for the precision radiation therapy. In recent years, considerable clinical and research efforts have been devoted to integrate MRI into radiotherapy workflow motivated by the superior soft tissue contrast and functional imaging possibility. Dynamic contrast-enhanced MRI (DCE-MRI) is a noninvasive technique that measures properties of tissue microvasculature. Its sensitivity to radiation-induced vascular pharmacokinetic (PK) changes has been preliminary demonstrated. In spite of its great potential, two major challenges have limited DCE-MRI’s clinical application in radiotherapy assessment: the technical limitations of accurate DCE-MRI imaging implementation and the need of novel DCE-MRI data analysis methods for richer functional heterogeneity information.

This study aims at improving current DCE-MRI techniques and developing new DCE-MRI analysis methods for particular radiotherapy assessment. Thus, the study is naturally divided into two parts. The first part focuses on DCE-MRI temporal resolution as one of the key DCE-MRI technical factors, and some improvements regarding DCE-MRI temporal resolution are proposed; the second part explores the potential value of image heterogeneity analysis and multiple PK model combination for therapeutic response assessment, and several novel DCE-MRI data analysis methods are developed.

I. Improvement of DCE-MRI temporal resolution. First, the feasibility of improving DCE-MRI temporal resolution via image undersampling was studied. Specifically, a novel MR image iterative reconstruction algorithm was studied for DCE-MRI reconstruction. This algorithm was built on the recently developed compress sensing (CS) theory. By utilizing a limited k-space acquisition with shorter imaging time, images can be reconstructed in an iterative fashion under the regularization of a newly proposed total generalized variation (TGV) penalty term. In the retrospective study of brain radiosurgery patient DCE-MRI scans under IRB-approval, the clinically obtained image data was selected as reference data, and the simulated accelerated k-space acquisition was generated via undersampling the reference image full k-space with designed sampling grids. Two undersampling strategies were proposed: 1) a radial multi-ray grid with a special angular distribution was adopted to sample each slice of the full k-space; 2) a Cartesian random sampling grid series with spatiotemporal constraints from adjacent frames was adopted to sample the dynamic k-space series at a slice location. Two sets of PK parameters’ maps were generated from the undersampled data and from the fully-sampled data, respectively. Multiple quantitative measurements and statistical studies were performed to evaluate the accuracy of PK maps generated from the undersampled data in reference to the PK maps generated from the fully-sampled data. Results showed that at a simulated acceleration factor of four, PK maps could be faithfully calculated from the DCE images that were reconstructed using undersampled data, and no statistically significant differences were found between the regional PK mean values from undersampled and fully-sampled data sets. DCE-MRI acceleration using the investigated image reconstruction method has been suggested as feasible and promising.

Second, for high temporal resolution DCE-MRI, a new PK model fitting method was developed to solve PK parameters for better calculation accuracy and efficiency. This method is based on a derivative-based deformation of the commonly used Tofts PK model, which is presented as an integrative expression. This method also includes an advanced Kolmogorov-Zurbenko (KZ) filter to remove the potential noise effect in data and solve the PK parameter as a linear problem in matrix format. In the computer simulation study, PK parameters representing typical intracranial values were selected as references to simulated DCE-MRI data for different temporal resolution and different data noise level. Results showed that at both high temporal resolutions (<1s) and clinically feasible temporal resolution (~5s), this new method was able to calculate PK parameters more accurate than the current calculation methods at clinically relevant noise levels; at high temporal resolutions, the calculation efficiency of this new method was superior to current methods in an order of 102. In a retrospective of clinical brain DCE-MRI scans, the PK maps derived from the proposed method were comparable with the results from current methods. Based on these results, it can be concluded that this new method can be used for accurate and efficient PK model fitting for high temporal resolution DCE-MRI.

II. Development of DCE-MRI analysis methods for therapeutic response assessment. This part aims at methodology developments in two approaches. The first one is to develop model-free analysis method for DCE-MRI functional heterogeneity evaluation. This approach is inspired by the rationale that radiotherapy-induced functional change could be heterogeneous across the treatment area. The first effort was spent on a translational investigation of classic fractal dimension theory for DCE-MRI therapeutic response assessment. In a small-animal anti-angiogenesis drug therapy experiment, the randomly assigned treatment/control groups received multiple fraction treatments with one pre-treatment and multiple post-treatment high spatiotemporal DCE-MRI scans. In the post-treatment scan two weeks after the start, the investigated Rényi dimensions of the classic PK rate constant map demonstrated significant differences between the treatment and the control groups; when Rényi dimensions were adopted for treatment/control group classification, the achieved accuracy was higher than the accuracy from using conventional PK parameter statistics. Following this pilot work, two novel texture analysis methods were proposed. First, a new technique called Gray Level Local Power Matrix (GLLPM) was developed. It intends to solve the lack of temporal information and poor calculation efficiency of the commonly used Gray Level Co-Occurrence Matrix (GLCOM) techniques. In the same small animal experiment, the dynamic curves of Haralick texture features derived from the GLLPM had an overall better performance than the corresponding curves derived from current GLCOM techniques in treatment/control separation and classification. The second developed method is dynamic Fractal Signature Dissimilarity (FSD) analysis. Inspired by the classic fractal dimension theory, this method measures the dynamics of tumor heterogeneity during the contrast agent uptake in a quantitative fashion on DCE images. In the small animal experiment mentioned before, the selected parameters from dynamic FSD analysis showed significant differences between treatment/control groups as early as after 1 treatment fraction; in contrast, metrics from conventional PK analysis showed significant differences only after 3 treatment fractions. When using dynamic FSD parameters, the treatment/control group classification after 1st treatment fraction was improved than using conventional PK statistics. These results suggest the promising application of this novel method for capturing early therapeutic response.

The second approach of developing novel DCE-MRI methods is to combine PK information from multiple PK models. Currently, the classic Tofts model or its alternative version has been widely adopted for DCE-MRI analysis as a gold-standard approach for therapeutic response assessment. Previously, a shutter-speed (SS) model was proposed to incorporate transcytolemmal water exchange effect into contrast agent concentration quantification. In spite of richer biological assumption, its application in therapeutic response assessment is limited. It might be intriguing to combine the information from the SS model and from the classic Tofts model to explore potential new biological information for treatment assessment. The feasibility of this idea was investigated in the same small animal experiment. The SS model was compared against the Tofts model for therapeutic response assessment using PK parameter regional mean value comparison. Based on the modeled transcytolemmal water exchange rate, a biological subvolume was proposed and was automatically identified using histogram analysis. Within the biological subvolume, the PK rate constant derived from the SS model were proved to be superior to the one from Tofts model in treatment/control separation and classification. Furthermore, novel biomarkers were designed to integrate PK rate constants from these two models. When being evaluated in the biological subvolume, this biomarker was able to reflect significant treatment/control difference in both post-treatment evaluation. These results confirm the potential value of SS model as well as its combination with Tofts model for therapeutic response assessment.

In summary, this study addressed two problems of DCE-MRI application in radiotherapy assessment. In the first part, a method of accelerating DCE-MRI acquisition for better temporal resolution was investigated, and a novel PK model fitting algorithm was proposed for high temporal resolution DCE-MRI. In the second part, two model-free texture analysis methods and a multiple-model analysis method were developed for DCE-MRI therapeutic response assessment. The presented works could benefit the future DCE-MRI routine clinical application in radiotherapy assessment.

The role of habitat heterogeneity for the functional response of the spider Nesticodes rufipes (Araneae : Theridiidae) to houseflies

We investigated whether or not different degrees of refuge for prey influence the characteristic of functional response exhibited by the spider Nesticodes rufipes on Musca domestica, comparing the inherent ability of N. rufipes to kill individual houseflies in such environments at two distinct time intervals. To investigate these questions, two artificial habitats were elaborated in the laboratory. For 168 h of predator-prey interaction, logistic regression analyses revealed a type 11 functional response, and a significant decrease in prey capture in the highest prey density was observed when habitat complexity was increased. Data from habitat 1 (less complex) presented a greater coefficient of determination than those from habitat 2 (more complex), indicating a higher variation of predation of the latter. For a 24 h period of predator-prey interaction, spiders killed significantly fewer prey in habitat 2 than in habitat 1. Although prey capture did not enable data to fit properly in the random predator equation in this case, predation data from habitat 2 presented a higher variation than data from habitat 1, corroborating results from 168 h of interaction. The high variability observed on data from habitat 2 (more complex habitat) is an interesting result because it reinforces the importance of refuge in promoting spatial heterogeneity, which can affect the extent of predator-prey interactions.

Application of density functional theory to analysis of energetic heterogeneity and pore size distribution of activated carbons

A new approach based on the nonlocal density functional theory to determine pore size distribution (PSD) of activated carbons and energetic heterogeneity of the pore wall is proposed. The energetic heterogeneity is modeled with an energy distribution function (EDF), describing the distribution of solid-fluid potential well depth (this distribution is a Dirac delta function for an energetic homogeneous surface). The approach allows simultaneous determining of the PSD (assuming slit shape) and EDF from nitrogen or argon isotherms at their respective boiling points by using a set of local isotherms calculated for a range of pore widths and solid-fluid potential well depths. It is found that the structure of the pore wall surface significantly differs from that of graphitized carbon black. This could be attributed to defects in the crystalline structure of the surface, active oxide centers, finite size of the pore walls (in either wall thickness or pore length), and so forth. Those factors depend on the precursor and the process of carbonization and activation and hence provide a fingerprint for each adsorbent. The approach allows very accurate correlation of the experimental adsorption isotherm and leads to PSDs that are simpler and more realistic than those obtained with the original nonlocal density functional theory.

Kidney side population reveals multilineage potential and renal functional capacity but also cellular heterogeneity

Side population (SP) cells in the adult kidney are proposed to represent a progenitor population. However, the size, origin, phenotype, and potential of the kidney SP has been controversial. In this study, the SP fraction of embryonic and adult kidneys represented 0.1 to 0.2% of the total viable cell population. The immunophenotype and the expression profile of kidney SP cells was distinct from that of bone marrow SP cells, suggesting that they are a resident nonhematopoietic cell population. Affymetrix expression profiling implicated a role for Notch signaling in kidney SP cells and was used to identify markers of kidney SP. Localization by in situ hybridization confirmed a primarily proximal tubule location, supporting the existence of a tubular niche, but also revealed considerable heterogeneity, including the presence of renal macrophages. Adult kidney SP cells demonstrated multilineage differentiation in vitro, whereas microinjection into mouse metanephroi showed that SP cells had a 3.5- to 13-fold greater potential to contribute to developing kidney than non-SP main population cells. However, although reintroduction of SP cells into an Adriamycin-nephropathy model reduced albuminuria:creatinine ratios, this was without significant tubular integration, suggesting a humoral role for SP cells in renal repair. The heterogeneity of the renal SP highlights the need for further fractionation to distinguish the cellular subpopulations that are responsible for the observed multilineage capacity and transdifferentiative and humoral activities.

Greater CD8+ TCR heterogeneity and functional flexibility in HIV-2 compared to HIV-1 infection

Virus-specific CD8+ T cells are known to play an important role in the control of HIV infection. In this study we investigated whether there may be qualitative differences in the CD8+ T cell response in HIV-1- and HIV-2-infected individuals that contribute to the relatively efficient control of the latter infection. A molecular comparison of global TCR heterogeneity showed a more oligoclonal pattern of CD8 cells in HIV-1- than HIV-2-infected patients. This was reflected in restricted and conserved TCR usage by CD8+ T cells recognizing individual HLA-A2- and HLA-B57-restricted viral epitopes in HIV-1, with limited plasticity in their response to amino acid substitutions within these epitopes. The more diverse TCR usage observed for HIV-2-specific CD8 T cells was associated with an enhanced potential for CD8+ expansion and IFN- production on cross-recognition of variant epitopes. Our data suggest a mechanism that could account for any possible cross-protection that may be mediated by HIV-2-specific CD8+ T cells against HIV-1 infection. Furthermore, they have implications for HIV vaccine development, demonstrating an association between a polyclonal, virus-specific CD8+ T cell response and an enhanced capacity to tolerate substitutions within T cell epitopes.

Should athletes return to sport after applying Ice? A systematic review of the effect of local 3 cooling on functional performance

Applying ice or other forms of topical cooling is a popular method of treating sports injuries. It is commonplace for athletes to return to competitive activity, shortly or immediately after the application of a cold treatment. In this article, we examine the effect of local tissue cooling on outcomes relating to functional performance and to discuss their relevance to the sporting environment. A computerized literature search, citation tracking and hand search was performed up to April, 2011. Eligible studies were trials involving healthy human participants, describing the effects of cooling on outcomes relating to functional performance. Two reviewers independently assessed the validity of included trials and calculated effect sizes. Thirty five trials met the inclusion criteria; all had a high risk of bias. The mean sample size was 19. Meta-analyses were not undertaken due to clinical heterogeneity. The majority of studies used cooling durations >20 minutes. Strength (peak torque/force) was reported by 25 studies with approximately 75% recording a decrease in strength immediately following cooling. There was evidence from six studies that cooling adversely affected speed, power and agility-based running tasks; two studies found this was negated with a short rewarming period. There was conflicting evidence on the effect of cooling on isolated muscular endurance. A small number of studies found that cooling decreased upper limb dexterity and accuracy. The current evidence base suggests that athletes will probably be at a performance disadvantage if they return to activity immediately after cooling. This is based on cooling for longer than 20 minutes, which may exceed the durations employed in some sporting environments. In addition, some of the reported changes were clinically small and may only be relevant in elite sport. Until better evidence is available, practitioners should use short cooling applications and/or undertake a progressive warm up prior to returning to play.

Functional analysis of missense variants in the TRESK (KCNK18) K+ channel

A loss of function mutation in the TRESK K2P potassium channel (KCNK18), has recently been linked with typical familial migraine with aura. We now report the functional characterisation of additional TRESK channel missense variants identified in unrelated patients. Several variants either had no apparent functional effect, or they caused a reduction in channel activity. However, the C110R variant was found to cause a complete loss of TRESK function, yet is present in both sporadic migraine and control cohorts, and no variation in KCNK18 copy number was found. Thus despite the previously identified association between loss of TRESK channel activity and migraine in a large multigenerational pedigree, this finding indicates that a single non-functional TRESK variant is not alone sufficient to cause typical migraine and highlights the genetic complexity of this disorder. Migraine is a common, disabling neurological disorder with a genetic, environmental and in some cases hormonal component. It is characterized by attacks of severe, usually unilateral and throbbing headache, can be accompanied by nausea, vomiting and photophobia and is clinically divided into two main subtypes, migraine with aura (MA) when a migraine is accompanied by transient and reversible focal neurological symptoms and migraine without aura (MO)1. The multifactorial and clinical heterogeneity of the disorder have considerably hindered the identification of common migraine susceptibility genes and most of our current understanding comes from the studies of familial hemiplegic migraine (FHM), a rare monogenic autosomal dominant form of MA2. So far, the three susceptibility genes that have been convincingly identified in FHM families all encode ion channels or transporters: CACNA1A encoding the α1 subunit of the Cav2.1 calcium channel3, SCN1A encoding the Nav1.1 sodium channel4 and ATP1A2 encoding the α2 subunit of the Na+/K+ pump5. It is believed that mutations in these genes may lead to increased efflux of glutamate and potassium in the synapse and thereby cause migraine by rendering the brain more susceptible to cortical spreading depression (CSD)6 which is thought to play a role in initiating a migraine attack7,8. However, these genes have not to date been implicated in common forms of migraine9. Nevertheless, current opinion suggests that typical migraine, like FHM, is also disorder of neuronal excitability, ion homeostasis and neurotransmitter release10,11,12. Mutations in the SLC4A4 gene encoding the sodium-bicarbonate cotransporter NBCe1, have recently been implicated in several different forms of migraine13, and a variety of genes involved in glutamate homeostasis (PGCP, MTDH14 and LRP115) and a cation channel (TRPM8)15 have also recently been implicated in migraine via genome-wide association studies. Ion channels are therefore highly likely to play an important role in the pathogenesis of typical migraine. TRESK (KCNK18), is a member of the two-pore domain (K2P) family of potassium channels involved in the control of cellular electrical excitability16. Regulation of TRESK activity by the calcium-dependent phosphatase calcineurin17, as well as its expression in dorsal root ganglia (DRG)18 and trigeminal ganglia (TG)19,20 has led to a proposed role for this channel in a variety of pain pathways. In a recent study, a frameshift mutation (F139Wfsx24) in TRESK was identified in a large multigenerational pedigree where it co-segregated perfectly with typical MA and a significant genome-wide linkage LOD score of 3.0. Furthermore, functional analysis revealed that this mutation caused a complete loss of TRESK function and that the truncated subunit was also capable of down regulating wild-type channel function. This therefore highlighted KCNK18 as potentially important candidate gene and suggested that TRESK dysfunction might play a possible role in the pathogenesis of familial migraine with visual aura20. Additional screening for KCNK18 mutations in unrelated sporadic migraine and control cohorts also identified a number of other missense variants; R10G, A34V, C110R, S231P and A233V20. The A233V variant was found only in the control cohort, whilst A34V was identified in a single Australian migraine proband for which family samples were not available, but it was not detected in controls. By contrast, the R10G, C110R, and S231P variants were found in both migraineurs and controls in both cohorts. In this study, we have investigated the functional effect of these variants to further probe the potential association of TRESK dysfunction with typical migraine.

Predicting clinical response to anticancer drugs using an ex vivo platform that captures tumour heterogeneity

Predicting clinical response to anticancer drugs remains a major challenge in cancer treatment. Emerging reports indicate that the tumour microenvironment and heterogeneity can limit the predictive power of current biomarker-guided strategies for chemotherapy. Here we report the engineering of personalized tumour ecosystems that contextually conserve the tumour heterogeneity, and phenocopy the tumour microenvironment using tumour explants maintained in defined tumour grade-matched matrix support and autologous patient serum. The functional response of tumour ecosystems, engineered from 109 patients, to anticancer drugs, together with the corresponding clinical outcomes, is used to train a machine learning algorithm; the learned model is then applied to predict the clinical response in an independent validation group of 55 patients, where we achieve 100% sensitivity in predictions while keeping specificity in a desired high range. The tumour ecosystem and algorithm, together termed the CANScript technology, can emerge as a powerful platform for enabling personalized medicine.

Functional genomic studies of the structure and regulation of eukaryotic transcriptomes

The main focus of this thesis is the use of high-throughput sequencing technologies in functional genomics (in particular in the form of ChIP-seq, chromatin immunoprecipitation coupled with sequencing, and RNA-seq) and the study of the structure and regulation of transcriptomes. Some parts of it are of a more methodological nature while others describe the application of these functional genomic tools to address various biological problems. A significant part of the research presented here was conducted as part of the ENCODE (ENCyclopedia Of DNA Elements) Project.

The first part of the thesis focuses on the structure and diversity of the human transcriptome. Chapter 1 contains an analysis of the diversity of the human polyadenylated transcriptome based on RNA-seq data generated for the ENCODE Project. Chapter 2 presents a simulation-based examination of the performance of some of the most popular computational tools used to assemble and quantify transcriptomes. Chapter 3 includes a study of variation in gene expression, alternative splicing and allelic expression bias on the single-cell level and on a genome-wide scale in human lymphoblastoid cells; it also brings forward a number of critical to the practice of single-cell RNA-seq measurements methodological considerations.

The second part presents several studies applying functional genomic tools to the study of the regulatory biology of organellar genomes, primarily in mammals but also in plants. Chapter 5 contains an analysis of the occupancy of the human mitochondrial genome by TFAM, an important structural and regulatory protein in mitochondria, using ChIP-seq. In Chapter 6, the mitochondrial DNA occupancy of the TFB2M transcriptional regulator, the MTERF termination factor, and the mitochondrial RNA and DNA polymerases is characterized. Chapter 7 consists of an investigation into the curious phenomenon of the physical association of nuclear transcription factors with mitochondrial DNA, based on the diverse collections of transcription factor ChIP-seq datasets generated by the ENCODE, mouseENCODE and modENCODE consortia. In Chapter 8 this line of research is further extended to existing publicly available ChIP-seq datasets in plants and their mitochondrial and plastid genomes.

The third part is dedicated to the analytical and experimental practice of ChIP-seq. As part of the ENCODE Project, a set of metrics for assessing the quality of ChIP-seq experiments was developed, and the results of this activity are presented in Chapter 9. These metrics were later used to carry out a global analysis of ChIP-seq quality in the published literature (Chapter 10). In Chapter 11, the development and initial application of an automated robotic ChIP-seq (in which these metrics also played a major role) is presented.

The fourth part presents the results of some additional projects the author has been involved in, including the study of the role of the Piwi protein in the transcriptional regulation of transposon expression in Drosophila (Chapter 12), and the use of single-cell RNA-seq to characterize the heterogeneity of gene expression during cellular reprogramming (Chapter 13).

The last part of the thesis provides a review of the results of the ENCODE Project and the interpretation of the complexity of the biochemical activity exhibited by mammalian genomes that they have revealed (Chapters 15 and 16), an overview of the expected in the near future technical developments and their impact on the field of functional genomics (Chapter 14), and a discussion of some so far insufficiently explored research areas, the future study of which will, in the opinion of the author, provide deep insights into many fundamental but not yet completely answered questions about the transcriptional biology of eukaryotes and its regulation.

Temporal and spatial variation in the Nazaré Canyon (Western Iberian margin): Inter-annual and canyon heterogeneity effects on meiofauna biomass and diversity

The Nazaré Canyon on the Portuguese Margin (NE Atlantic) was sampled during spring-summer for three consecutive years (2005–2007), permitting the first inter-annual study of the meiofaunal communities at the Iberian Margin at two abyssal depths (~3500 m and ~4400 m). Using new and already published data, the meiofauna standing stocks (abundance and biomass) and nematode structural and functional diversity were investigated in relation to the sediment biogeochemistry (e.g. organic carbon, nitrogen, chlorophyll a, phaeopigments) and grain size. A conspicuous increase in sand content from 2005 to 2006 and decrease of phytodetritus at both sites, suggested the occurrence of one or more physical disturbance events. Nematode standing stocks and trophic diversity decreased after these events, seemingly followed by a recovery/recolonisation period in 2007, which was strongly correlated with an increase in the quantity and bioavailability of phytodetrital organic matter supplied. Changes in meiofauna assemblages, however, also differed between stations, likely because of the contrasting hydrodynamic and food supply conditions. Higher meiofauna and nematode abundances, biomass and trophic complexity were found at the shallowest canyon station, where the quantity, quality and bioavailability of food material were higher than at the deeper site. The present results suggest that even though inter-annual variations in the sedimentary environment can regulate the meiofauna in the abyssal Nazaré Canyon, heterogeneity between sampling locations in the canyon were more pronounced.

Flatworm nerve-muscle: structural and functional analysis

Platyhelminthes occupy a unique position in nerve-muscle evolution, being the most primitive of metazoan phyla. Essentially, their nervous system consists of an archaic brain and associated pairs of longitudinal nerve cords cross-linked as an orthogon by transverse commissures. Confocal imaging reveals that these central nervous system elements are in continuity with an array of peripheral nerve plexuses which innervate a well-differentiated grid work of somatic muscle as well as a complexity of myofibres associated with organs of attachment, feeding, and reproduction. Electrophysiological studies of flatworm muscles have exposed a diversity of voltage-activated ion channels that influence muscle contractile events. Neuronal cell types are mainly multi- and bi-polar and highly secretory in nature, producing a heterogeneity of vesicular inclusions whose contents have been identified cytochemically to include all three major types of cholinergic, aminergic, and peptidergic messenger molecules. A landmark discovery in flatworm neurobiology was the biochemical isolation and amino acid sequencing of two groups of native neuropeptides: neuropeptide F and FMRFamide-related peptides (FaRPs). Both families of neuropeptide are abundant and broadly distributed in platyhelminths, occurring in neuronal vesicles in representatives of all major flatworm taxa. Dual localization studies have revealed that peptidergic and cholinergic substances occupy neuronal sets separate from those of serotoninergic components. The physiological actions of neuronal messengers in flatworms are beginning to be established, and where examined, FaRPs and 5-HT are myoexcitatory, while cholinomimetic substances are generally inhibitory. There is immunocytochemical evidence that FaRPs and 5-HT have a regulatory role in the mechanism of egg assembly. Use of muscle strips and (or) muscle fibres from free-living and parasitic flatworms has provided baseline information to indicate that muscle responses to FaRPs are mediated by a G-protein-coupled receptor, and that the signal transduction pathway for contraction involves the second messengers cAMP and protein kinase C.