Application of a rank-based genetic association test to age-at-onset data from the Collaborative Study on the Genetics of Alcoholism study.

Association studies of quantitative traits have often relied on methods in which a normal distribution of the trait is assumed. However, quantitative phenotypes from complex human diseases are often censored, highly skewed, or contaminated with outlying values. We recently developed a rank-based association method that takes into account censoring and makes no distributional assumptions about the trait. In this study, we applied our new method to age-at-onset data on ALDX1 and ALDX2. Both traits are highly skewed (skewness > 1.9) and often censored. We performed a whole genome association study of age at onset of the ALDX1 trait using Illumina single-nucleotide polymorphisms. Only slightly more than 5% of markers were significant. However, we identified two regions on chromosomes 14 and 15, which each have at least four significant markers clustering together. These two regions may harbor genes that regulate age at onset of ALDX1 and ALDX2. Future fine mapping of these two regions with densely spaced markers is warranted.

Computational Methods for RNA Structure Validation and Improvement.

With increasing recognition of the roles RNA molecules and RNA/protein complexes play in an unexpected variety of biological processes, understanding of RNA structure-function relationships is of high current importance. To make clean biological interpretations from three-dimensional structures, it is imperative to have high-quality, accurate RNA crystal structures available, and the community has thoroughly embraced that goal. However, due to the many degrees of freedom inherent in RNA structure (especially for the backbone), it is a significant challenge to succeed in building accurate experimental models for RNA structures. This chapter describes the tools and techniques our research group and our collaborators have developed over the years to help RNA structural biologists both evaluate and achieve better accuracy. Expert analysis of large, high-resolution, quality-conscious RNA datasets provides the fundamental information that enables automated methods for robust and efficient error diagnosis in validating RNA structures at all resolutions. The even more crucial goal of correcting the diagnosed outliers has steadily developed toward highly effective, computationally based techniques. Automation enables solving complex issues in large RNA structures, but cannot circumvent the need for thoughtful examination of local details, and so we also provide some guidance for interpreting and acting on the results of current structure validation for RNA.

Sea Urchin Body Plan Development and Evolution: An Integrative Transcriptomic Approach

My dissertation work integrates comparative transcriptomics and functional analyses to investigate gene expression changes underlying two significant aspects of sea urchin evolution and development: the dramatic developmental changes associated with an ecologically significant shift in life history strategy and the development of the unusual radial body plan of adult sea urchins.

In Chapter 2, I investigate evolutionary changes in gene expression underlying the switch from feeding (planktotrophic) to nonfeeding (lecithotrophic) development in sea urchins. In order to identify these changes, I used Illumina RNA-seq to measure expression dynamics across 7 developmental stages in three sea urchin species: the lecithotroph Heliocidaris erythrogramma, the closely related planktotroph Heliocidaris tuberculata, and an outgroup planktotroph Lytechinus variegatus. My analyses draw on a well-characterized developmental gene regulatory network (GRN) in sea urchins to understand how the ancestral planktotrophic developmental program was altered during the evolution of lecithotrophic development. My results suggest that changes in gene expression profiles occurred more frequently across the transcriptome during the evolution of lecithotrophy than during the persistence of planktotrophy. These changes were even more pronounced within the GRN than across the transcriptome as a whole, and occurred in each network territory (skeletogenic, endomesoderm and ectoderm). I found evidence for both conservation and divergence of regulatory interactions in the network, as well as significant changes in the expression of genes with known roles in larval skeletogenesis, which is dramatically altered in lecithotrophs. I further explored network dynamics between species using coexpression analyses, which allowed me to identify novel players likely involved in sea urchin neurogenesis and endoderm patterning.

In Chapter 3, I investigate developmental changes in gene expression underlying radial body plan development and metamorphosis in H. erythrogramma. Using Illumina RNA-seq, I measured gene expression profiles across larval, metamorphic, and post-metamorphic life cycle phases. My results present a high-resolution view of gene expression dynamics during the complex transition from pre- to post-metamorphic development and suggest that distinct sets of regulatory and effector proteins are used during different life history phases.

Collectively, my investigations provide an important foundation for future, empirical studies to investigate the functional role of gene expression change in the evolution of developmental differences between species and also for the generation of the unusual radial body plan of sea urchins.

Expression of the Ets transcription factor Erm is regulated through a conventional PKC signaling pathway in the Molt4 lymphoblastic cell line.

Erm, a member of the PEA3 group within the Ets family of transcription factors, is expressed in murine and human lymphocytes. Here, we show that in the human Molt4 lymphoblastic cell line, the erm gene expression is regulated by the conventional PKC (cPKC) pathway. To better characterize the molecular mechanism by which cPKC regulates Erm transcription in Molt4 cells, we tested proximal promoter deletions of the human gene, and identified a specific cPKC-regulated region between positions -420 and -115 upstream of the first exon.

Optimal kinetics for quantification of antigen-induced cytokines in human peripheral blood mononuclear cells by real-time PCR and by ELISA.

Real-time polymerase chain reaction (PCR) has recently been described as a new tool to measure and accurately quantify mRNA levels. In this study, we have applied this technique to evaluate cytokine mRNA synthesis induced by antigenic stimulation with purified protein derivative (PPD) or heparin-binding haemagglutinin (HBHA) in human peripheral blood mononuclear cells (PBMC) from Mycobacterium tuberculosis-infected individuals. Whereas PPD and HBHA optimally induced IL-2 mRNA after respectively 8 and 16 to 24 h of in vitro stimulation, longer in vitro stimulation times were necessary for optimal induction of interferon-gamma (IFN-gamma) mRNA, respectively 16 to 24 h for PPD and 24 to 96 h for HBHA. IL-13 mRNA was optimally induced by in vitro stimulation after 16-48 h for PPD and after 48 to 96 h for HBHA. Comparison of antigen-induced Th1 and Th2 cytokines appears, therefore, valuable only if both cytokine types are analysed at their optimal time point of production, which, for a given cytokine, may differ for each antigen tested. Results obtained by real-time PCR for IFN-gamma and IL-13 mRNA correlated well with those obtained by measuring the cytokine concentrations in cell culture supernatants, provided they were high enough to be detected. We conclude that real-time PCR can be successfully applied to the quantification of antigen-induced cytokine mRNA and to the evaluation of the Th1/Th2 balance, only if the kinetics of cytokine mRNA appearance are taken into account and evaluated for each cytokine measured and each antigen analysed.

Characterization of autonomous thyroid adenoma: metabolism, gene expression, and pathology.

Fifty-one in vivo characterized autonomous single adenomas have been studied for functional parameters in vitro, for gene and protein expression and for pathology, and have been systematically compared to the corresponding extratumoral quiescent tissue. The adenomas were characterized by a high level of iodide trapping that corresponds to a high level of Na+ /iodide symporter gene expression, a high thyroperoxidase mRNA and protein content, and a low H2O2 generation. This explains the iodide metabolism characteristics demonstrated before, ie, the main cause of the "hot" character of the adenomas is their increased iodide transport. The adenomas spontaneously secreted higher amounts of thyroid hormone than the quiescent tissue and in agreement with previous in vivo data, this secretion could be further enhanced by thyrotropin (TSH). Inositol uptake was also increased but there was no spontaneous increase of the generation of inositol phosphates and this metabolism could be further activated by TSH. These positive responses to TSH are in agreement with the properties of TSH-stimulated thyroid cells in vitro and in vivo. They are compatible with the characteristics of mutated TSH receptors whose constitutive activation accounts for the majority of autonomous thyroid adenomas in Europe. The number of cycling cells, as evaluated by MIB-1 immunolabeling was low but increased in comparison with the corresponding quiescent tissue or normal tissue. The cycling cells are observed mainly at the periphery; there was very little apoptosis. Both findings account for the slow growth of these established adenomas. On the other hand, by thyroperoxidase immunohistochemistry, the whole lesion appeared hyperfunctional, which demonstrates a dissociation of mitogenic and functional stimulations. Thyroglobulin, TSH receptor, and E-cadherin mRNA accumulations were not modified in a consistent way, which confirms the near-constitutive expression of the corresponding genes in normal differentiated tissue. On the contrary, early immediate genes expressions (c-myc, NGF1B, egr 1, genes of the fos and jun families) were decreased. This may be explained by the proliferative heterogeneity of the lesion and the previously described short, biphasic expression of these genes when induced by mitogenic agents. All the characteristics of the autonomous adenomas can therefore be explained by the effect of the known activating mutations of genes coding for proteins of the TSH cyclic adenosine monophosphate (cAMP) cascade, all cells being functionally activated while only those at the periphery multiply. The reason of this heterogeneity is unknown.

Study of gene expression in thyrotropin-stimulated thyroid cells by cDNA expression array: ID3 transcription modulating factor as an early response protein and tumor marker in thyroid carcinomas.

Induction of cell proliferation by mitogen or growth factor stimulation leads to the specific induction or repression of a large number of genes. To identify genes differentially regulated by the cAMP-dependent transduction pathway, which is poorly characterized so far, we used the cDNA expression array technology. Hybridizations of Atlas human cDNA expression arrays with (32)P-labeled cDNA probes derived from control or thyrotropin (TSH)-stimulated dog thyrocytes in primary culture generated expression profiles of hundreds of genes simultaneously. Among the genes that displayed modified expression, we selected the transcription factor ID3, whose expression was increased by a cAMP-dependent stimulus. ID3 overexpression after TSH stimulation was first verified by Northern blotting analysis, and its mRNA regulation was then investigated in response to a variety of agents acting on thyrocyte proliferation and/or differentiation. We show that: (1) ID3 mRNA induction was stronger after stimulation of the cAMP cascade, but was not restricted to this signaling pathway, as phorbol myristate ester (TPA) and insulin also stimulated mRNA accumulation; (2) in contrast, powerful mitogens for thyroid cells, epidermal growth factor and hepatocyte growth factor, did not significantly modify ID3 mRNA levels; (3) ID3 protein levels closely parallelled mRNA levels, as revealed by immunofluorescence experiments showing a nuclear signal regulated by TSH; (4) in papillary thyroid carcinomas, ID3 mRNA was downregulated. Our results suggest that ID3 expression might be more related to the differentiating process induced by TSH than to the proliferative action of this hormone.

Evaluating The Ribosomal Internal Transcribed Spacer (ITS) as a Candidate Dinoflagellate Barcode Marker

DNA barcoding offers an efficient way to determine species identification and to measure biodiversity. For dinoflagellates, an ancient alveolate group of about 2000 described extant species, DNA barcoding studies have revealed large amounts of unrecognized species diversity, most of which is not represented in culture collections. To date, two mitochondrial gene markers, Cytochrome Oxidase I (COI) and Cytochrome b oxidase (COB), have been used to assess DNA barcoding in dinoflagellates, and both failed to amplify all taxa and suffered from low resolution. Nevertheless, both genes yielded many examples of morphospecies showing cryptic speciation and morphologically distinct named species being genetically similar, highlighting the need for a common marker. For example, a large number of cultured Symbiodinium strains have neither taxonomic identification, nor a common measure of diversity that can be used to compare this genus to other dinoflagellates.

Validation of the detection of Pseudo-nitzschia spp. using specific RNA probes tested in a microarray format: Calibration of signal based on variability of RNA content with environmental conditions.

Harmful algal blooms (HAB) occur worldwide and cause health problems and economic damage to fisheries and tourism. Monitoring for toxic algae is therefore essential but is based primarily on light microscopy, which is time consuming and can be limited by insufficient morphological characters such that more time is needed to examine critical features with electron microscopy. Monitoring with molecular tools is done in only a few places world-wide. EU FP7 MIDTAL (Microarray Detection of Toxic Algae) used SSU and LSU rRNA genes as targets on microarrays to identify toxic species. In order to comply with current monitoring requirements to report cell numbers as the relevant threshold measurement to trigger closure of fisheries, it was necessary to calibrate our microarray to convert the hybridisation signal obtained to cell numbers. Calibration curves for two species of Pseudo-nitzschia for use with the MIDTAL microarray are presented to obtain cell numbers following hybridisation. It complements work presented by Barra et al. (2012b. Environ. Sci. Pollut. Res. doi: 10.1007/s11356-012-1330-1v) for two other Pseudo-nitzschia spp., Dittami and Edvardsen (2012a. J. Phycol. 48, 1050) for Pseudochatonella, Blanco et al. (2013. Harmful Algae 24, 80) for Heterosigma, McCoy et al. (2013. FEMS. doi: 10.1111/1574-6941.12277) for Prymnesium spp., Karlodinium veneficum, and cf. Chatonella spp. and Taylor et al. (2014. Harmful Algae, in press) for Alexandrium.

Feeding selectivity of bivalve larvae on natural plankton assemblages in the Western English Channel

Meroplankton, including bivalve larvae, are an important and yet understudied component of coastal marine food webs. Understanding the baseline of meroplankton ecology is imperative to establish and predict their sensitivity to local and global marine stressors. Over an annual cycle (October 2009–September 2010), bivalve larvae were collected from the Western Channel Observatory time series station L4 (50°15.00′N, 4°13.02′W). The morphologically similar larvae were identified by analysis of the 18S nuclear small subunit ribosomal RNA gene, and a series of incubation experiments were conducted to determine larval ingestion rates on natural plankton assemblages. Complementary gut content analysis was performed using a PCR-based method for detecting prey DNA both from field-collected larvae and those from the feeding experiments. Molecular identification of bivalve larvae showed the community composition to change over the course of the sampling period with domination by Phaxas in winter and higher diversity in autumn. The larvae selected for nanoeukaryotes (2–20 µm) including coccolithophores (<20 µm) which together comprised >75 % of the bivalve larvae diet. Additionally, a small percentage of carbon ingested originated from heterotrophic ciliates (<30 µm). The molecular analysis of bivalve larvae gut content provided increased resolution of identification of prey consumed and demonstrated that the composition of prey consumed established through bottle incubations conferred with that established from in situ larvae. Despite changes in bivalve larvae community structure, clearance rates of each prey type did not change significantly over the course of the experiment, suggesting different bivalve larvae species may consume similar prey.