Discovering Frequent Poly-Regions in DNA Sequences

The problem of discovering frequent arrangements of regions of high occurrence of one or more items of a given alphabet in a sequence is studied, and two efficient approaches are proposed to solve it. The first approach is entropy-based and uses an existing recursive segmentation technique to split the input sequence into a set of homogeneous segments. The key idea of the second approach is to use a set of sliding windows over the sequence. Each sliding window keeps a set of statistics of a sequence segment that mainly includes the number of occurrences of each item in that segment. Combining these statistics efficiently yields the complete set of regions of high occurrence of the items of the given alphabet. After identifying these regions, the sequence is converted to a sequence of labeled intervals (each one corresponding to a region). An efficient algorithm for mining frequent arrangements of temporal intervals on a single sequence is applied on the converted sequence to discover frequently occurring arrangements of these regions. The proposed algorithms are tested on various DNA sequences producing results with significant biological meaning.

Compilation and analysis of integrated regional input-output tables for NUTS 2 regions in Ireland

In 1966, Roy Geary, Director of the ESRI, noted “the absence of any kind of import and export statistics for regions is a grave lacuna” and further noted that if regional analyses were to be developed then regional Input-Output Tables must be put on the “regular statistical assembly line”. Forty-five years later, the lacuna lamented by Geary still exists and remains the most significant challenge to the construction of regional Input-Output Tables in Ireland. The continued paucity of sufficient regional data to compile effective regional Supply and Use and Input-Output Tables has retarded the capacity to construct sound regional economic models and provide a robust evidence base with which to formulate and assess regional policy. This study makes a first step towards addressing this gap by presenting the first set of fully integrated, symmetric, Supply and Use and domestic Input-Output Tables compiled for the NUTS 2 regions in Ireland: The Border, Midland and Western region and the Southern & Eastern region. These tables are general purpose in nature and are consistent fully with the official national Supply & Use and Input-Output Tables, and the regional accounts. The tables are constructed using a survey-based or bottom-up approach rather than employing modelling techniques, yielding more robust and credible tables. These tables are used to present a descriptive statistical analysis of the two administrative NUTS 2 regions in Ireland, drawing particular attention to the underlying structural differences of regional trade balances and composition of Gross Value Added in those regions. By deriving regional employment multipliers, Domestic Demand Employment matrices are constructed to quantify and illustrate the supply chain impact on employment. In the final part of the study, the predictive capability of the Input-Output framework is tested over two time periods. For both periods, the static Leontief production function assumptions are relaxed to allow for labour productivity. Comparative results from this experiment are presented.

Identification of RNA editing in the human exome and development of DARNED DatabaseSF

RNA editing is a biological phenomena that alters nascent RNA transcripts by insertion, deletion and/or substitution of one or a few nucleotides. It is ubiquitous in all kingdoms of life and in viruses. The predominant editing event in organisms with a developed central nervous system is Adenosine to Inosine deamination. Inosine is recognized as Guanosine by the translational machinery and reverse-transcriptase. In primates, RNA editing occurs frequently in transcripts from repetitive regions of the genome. In humans, more than 500,000 editing instances have been identified, by applying computational pipelines on available ESTs and high-throughput sequencing data, and by using chemical methods. However, the functions of only a small number of cases have been studied thoroughly. RNA editing instances have been found to have roles in peptide variants synthesis by non-synonymous codon substitutions, transcript variants by alterations in splicing sites and gene silencing by miRNAs sequence modifications. We established the Database of RNA EDiting (DARNED) to accommo-date the reference genomic coordinates of substitution editing in human, mouse and fly transcripts from published literatures, with additional information on edited genomic coordinates collected from various databases e.g. UCSC, NCBI. DARNED contains mostly Adenosine to Inosine editing and allows searches based on genomic region, gene ID, and user provided sequence. The Database is accessible at http://darned.ucc.ie RNA editing instances in coding region are likely to result in recoding in protein synthesis. This encouraged me to focus my research on the occurrences of RNA editing specific CDS and non-Alu exonic regions. By applying various filters on discrepancies between available ESTs and their corresponding reference genomic sequences, putative RNA editing candidates were identified. High-throughput sequencing was used to validate these candidates. All predicted coordinates appeared to be either SNPs or unedited.

The biodiversity, control and genetic characterisation of the lactococcal 936 group phages

Phages belonging to the 936 group represent one of the most prevalent and frequently isolated phages in dairy fermentation processes using Lactococcus lactis as the primary starter culture. In recent years extensive research has been carried out to characterise this phage group at a genomic level in an effort to understand how the 936 group phages dominate this particular niche and cause regular problems during large scale milk fermentations. This thesis describes a large scale screening of industrial whey samples, leading to the isolation of forty three genetically different lactococcal phages. Using multiplex PCR, all phages were identified as members of the 936 group. The complete genome of thirty eight of these phages was determined using next generation sequencing technologies which identified several regions of divergence. These included the structural region surrounding the major tail protein, the replication region as well as the genes involved in phage DNA packing. For a number of phages the latter genomic region was found to harbour genes encoding putative orphan methyltransferases. Using small molecule real time (SMRT) sequencing and heterologous gene expression, the target motifs for several of these MTases were determined and subsequently shown to actively protect phage DNA from restriction endonuclease activity. Comparative analysis of the thirty eight phages with fifty two previously sequenced members of this group showed that the core genome consists of 28 genes, while the non-core genome was found to fluctuate irrespective of geographical location or time of isolation. This study highlights the continued need to perform large scale characterisation of the bacteriophage populations infecting industrial fermentation facilities in effort to further our understanding dairy phages and ways to control their proliferation.

Similar T-cell immune responses induced by group M consensus env immunogens with wild-type or minimum consensus variable regions.

Consensus HIV-1 genes can decrease the genetic distances between candidate immunogens and field virus strains. To ensure the functionality and optimal presentation of immunologic epitopes, we generated two group-M consensus env genes that contain variable regions either from a wild-type B/C recombinant virus isolate (CON6) or minimal consensus elements (CON-S) in the V1, V2, V4, and V5 regions. C57BL/6 and BALB/c mice were primed twice with CON6, CON-S, and subtype control (92UG37_A and HXB2/Bal_B) DNA and boosted with recombinant vaccinia virus (rVV). Mean antibody titers against 92UG37_A, 89.6_B, 96ZM651_C, CON6, and CON-S Env protein were determined. Both CON6 and CON-S induced higher mean antibody titers against several of the proteins, as compared with the subtype controls. However, no significant differences were found in mean antibody titers in animals immunized with CON6 or CON-S. Cellular immune responses were measured by using five complete Env overlapping peptide sets: subtype A (92UG37_A), subtype B (MN_B, 89.6_B and SF162_B), and subtype C (Chn19_C). The intensity of the induced cellular responses was measured by using pooled Env peptides; T-cell epitopes were identified by using matrix peptide pools and individual peptides. No significant differences in T-cell immune-response intensities were noted between CON6 and CON-S immunized BALB/c and C57BL/6 mice. In BALB/c mice, 10 and eight nonoverlapping T-cell epitopes were identified in CON6 and CON-S, whereas eight epitopes were identified in 92UG37_A and HXB2/BAL_B. In C57BL/6 mice, nine and six nonoverlapping T-cell epitopes were identified after immunization with CON6 and CON-S, respectively, whereas only four and three were identified in 92UG37_A and HXB2/BAL_B, respectively. When combined together from both mouse strains, 18 epitopes were identified. The group M artificial consensus env genes, CON6 and CON-S, were equally immunogenic in breadth and intensity for inducing humoral and cellular immune responses.

A longitudinal study of epigenetic variation in twins.

DNA methylation is a key epigenetic mechanism involved in the developmental regulation of gene expression. Alterations in DNA methylation are established contributors to inter-individual phenotypic variation and have been associated with disease susceptibility. The degree to which changes in loci-specific DNA methylation are under the influence of heritable and environmental factors is largely unknown. In this study, we quantitatively measured DNA methylation across the promoter regions of the dopamine receptor 4 gene (DRD4), the serotonin transporter gene (SLC6A4/SERT) and the X-linked monoamine oxidase A gene (MAOA) using DNA sampled at both ages 5 and 10 years in 46 MZ twin-pairs and 45 DZ twin-pairs (total n=182). Our data suggest that DNA methylation differences are apparent already in early childhood, even between genetically identical individuals, and that individual differences in methylation are not stable over time. Our longitudinal-developmental study suggests that environmental influences are important factors accounting for interindividual DNA methylation differences, and that these influences differ across the genome. The observation of dynamic changes in DNA methylation over time highlights the importance of longitudinal research designs for epigenetic research.

The genomic analysis of erythrocyte microRNA expression in sickle cell diseases.

BACKGROUND: Since mature erythrocytes are terminally differentiated cells without nuclei and organelles, it is commonly thought that they do not contain nucleic acids. In this study, we have re-examined this issue by analyzing the transcriptome of a purified population of human mature erythrocytes from individuals with normal hemoglobin (HbAA) and homozygous sickle cell disease (HbSS). METHODS AND FINDINGS: Using a combination of microarray analysis, real-time RT-PCR and Northern blots, we found that mature erythrocytes, while lacking ribosomal and large-sized RNAs, contain abundant and diverse microRNAs. MicroRNA expression of erythrocytes was different from that of reticulocytes and leukocytes, and contributed the majority of the microRNA expression in whole blood. When we used microRNA microarrays to analyze erythrocytes from HbAA and HbSS individuals, we noted a dramatic difference in their microRNA expression pattern. We found that miR-320 played an important role for the down-regulation of its target gene, CD71 during reticulocyte terminal differentiation. Further investigation revealed that poor expression of miR-320 in HbSS cells was associated with their defective downregulation CD71 during terminal differentiation. CONCLUSIONS: In summary, we have discovered significant microRNA expression in human mature erythrocytes, which is dramatically altered in HbSS erythrocytes and their defect in terminal differentiation. Thus, the global analysis of microRNA expression in circulating erythrocytes can provide mechanistic insights into the disease phenotypes of erythrocyte diseases.

Histone modifications within the human X centromere region.

Human centromeres are multi-megabase regions of highly ordered arrays of alpha satellite DNA that are separated from chromosome arms by unordered alpha satellite monomers and other repetitive elements. Complexities in assembling such large repetitive regions have limited detailed studies of centromeric chromatin organization. However, a genomic map of the human X centromere has provided new opportunities to explore genomic architecture of a complex locus. We used ChIP to examine the distribution of modified histones within centromere regions of multiple X chromosomes. Methylation of H3 at lysine 4 coincided with DXZ1 higher order alpha satellite, the site of CENP-A localization. Heterochromatic histone modifications were distributed across the 400-500 kb pericentromeric regions. The large arrays of alpha satellite and gamma satellite DNA were enriched for both euchromatic and heterochromatic modifications, implying that some pericentromeric repeats have multiple chromatin characteristics. Partial truncation of the X centromere resulted in reduction in the size of the CENP-A/Cenp-A domain and increased heterochromatic modifications in the flanking pericentromere. Although the deletion removed approximately 1/3 of centromeric DNA, the ratio of CENP-A to alpha satellite array size was maintained in the same proportion, suggesting that a limited, but defined linear region of the centromeric DNA is necessary for kinetochore assembly. Our results indicate that the human X centromere contains multiple types of chromatin, is organized similarly to smaller eukaryotic centromeres, and responds to structural changes by expanding or contracting domains.

Evolution of the sex-related locus and genomic features shared in microsporidia and fungi.

BACKGROUND: Microsporidia are obligate intracellular, eukaryotic pathogens that infect a wide range of animals from nematodes to humans, and in some cases, protists. The preponderance of evidence as to the origin of the microsporidia reveals a close relationship with the fungi, either within the kingdom or as a sister group to it. Recent phylogenetic studies and gene order analysis suggest that microsporidia share a particularly close evolutionary relationship with the zygomycetes. METHODOLOGY/PRINCIPAL FINDINGS: Here we expanded this analysis and also examined a putative sex-locus for variability between microsporidian populations. Whole genome inspection reveals a unique syntenic gene pair (RPS9-RPL21) present in the vast majority of fungi and the microsporidians but not in other eukaryotic lineages. Two other unique gene fusions (glutamyl-prolyl tRNA synthetase and ubiquitin-ribosomal subunit S30) that are present in metazoans, choanoflagellates, and filasterean opisthokonts are unfused in the fungi and microsporidians. One locus previously found to be conserved in many microsporidian genomes is similar to the sex locus of zygomycetes in gene order and architecture. Both sex-related and sex loci harbor TPT, HMG, and RNA helicase genes forming a syntenic gene cluster. We sequenced and analyzed the sex-related locus in 11 different Encephalitozoon cuniculi isolates and the sibling species E. intestinalis (3 isolates) and E. hellem (1 isolate). There was no evidence for an idiomorphic sex-related locus in this Encephalitozoon species sample. According to sequence-based phylogenetic analyses, the TPT and RNA helicase genes flanking the HMG genes are paralogous rather than orthologous between zygomycetes and microsporidians. CONCLUSION/SIGNIFICANCE: The unique genomic hallmarks between microsporidia and fungi are independent of sequence based phylogenetic comparisons and further contribute to define the borders of the fungal kingdom and support the classification of microsporidia as unusual derived fungi. And the sex/sex-related loci appear to have been subject to frequent gene conversion and translocations in microsporidia and zygomycetes.

Cholera in Haiti and other Caribbean regions, 19th century.

Medical journals and other sources do not show evidence that cholera occurred in Haiti before 2010, despite the devastating effect of this disease in the Caribbean region in the 19th century. Cholera occurred in Cuba in 1833-1834; in Jamaica, Cuba, Puerto Rico, St. Thomas, St. Lucia, St. Kitts, Nevis, Trinidad, the Bahamas, St. Vincent, Granada, Anguilla, St. John, Tortola, the Turks and Caicos, the Grenadines (Carriacou and Petite Martinique), and possibly Antigua in 1850-1856; and in Guadeloupe, Cuba, St. Thomas, the Dominican Republic, Dominica, Martinique, and Marie Galante in 1865-1872. Conditions associated with slavery and colonial military control were absent in independent Haiti. Clustered populations, regular influx of new persons, and close quarters of barracks living contributed to spread of cholera in other Caribbean locations. We provide historical accounts of the presence and spread of cholera epidemics in Caribbean islands.

Geography and Genetics of Ecological Speciation in Pea Aphids

During ecological speciation, divergent natural selection drives evolution of ecological specialization and genetic differentiation of populations on alternate environments. Populations diverging onto the same alternate environments may be geographically widespread, so that divergence may occur at an array of locations simultaneously. Spatial variation in the process of divergence may produce a pattern of differences in divergence among locations called the Geographic Mosaic of Divergence. Diverging populations may vary in their degree of genetic differentiation and ecological specialization among locations. My dissertation examines the pattern and evolutionary processes of divergence in pea aphids (Acyrthosiphon pisum) on alfalfa (Medicago sativa) and clover (Trifolium pretense). In Chapter One, I examined differences among North American aphid populations in genetic differentiation at nuclear, sequence-based markers and in ecological specialization, measured as aphid fecundity on each host plant. In the East, aphids showed high host-plant associated ecological specialization and high genetic differentiation. In the West, aphids from clover were genetically indistinguishable from aphids on alfalfa, and aphids from clover were less specialized. Thus, the pattern of divergence differed among locations, suggesting a Geographic Mosaic of Divergence. In Chapter Two, I examined genomic heterogeneity in divergence in aphids on alfalfa and clover across North America using amplified fragment length polymorphisms (AFLPs). The degree of genetic differentiation varied greatly among markers, suggesting that divergent natural selection drives aphid divergence in all geographic locations. Three of the same genetic markers were identified as evolving under divergent selection in the eastern and western regions, and additional divergent markers were identified in the East. In Chapter Three, I investigated population structure of aphids in North America, France, and Sweden using AFLPs. Aphids on the same host plant were genetically similar across many parts of their range, so the evolution of host plant specialization does not appear to have occurred independently in every location. While aphids on alfalfa and clover were genetically differentiated in most locations, aphids from alfalfa and clover were genetically similar in both western North America and Sweden. High gene flow from alfalfa onto clover may constrain divergence in these locations.

Functional parcellation of attentional control regions of the brain.

Recently, a number of investigators have examined the neural loci of psychological processes enabling the control of visual spatial attention using cued-attention paradigms in combination with event-related functional magnetic resonance imaging. Findings from these studies have provided strong evidence for the involvement of a fronto-parietal network in attentional control. In the present study, we build upon this previous work to further investigate these attentional control systems. In particular, we employed additional controls for nonattentional sensory and interpretative aspects of cue processing to determine whether distinct regions in the fronto-parietal network are involved in different aspects of cue processing, such as cue-symbol interpretation and attentional orienting. In addition, we used shorter cue-target intervals that were closer to those used in the behavioral and event-related potential cueing literatures. Twenty participants performed a cued spatial attention task while brain activity was recorded with functional magnetic resonance imaging. We found functional specialization for different aspects of cue processing in the lateral and medial subregions of the frontal and parietal cortex. In particular, the medial subregions were more specific to the orienting of visual spatial attention, while the lateral subregions were associated with more general aspects of cue processing, such as cue-symbol interpretation. Additional cue-related effects included differential activations in midline frontal regions and pretarget enhancements in the thalamus and early visual cortical areas.

Expansion of the alpha 2-adrenergic receptor family: cloning and characterization of a human alpha 2-adrenergic receptor subtype, the gene for which is located on chromosome 2.

Pharmacologic, biochemical, and genetic analyses have demonstrated the existence of multiple alpha 2-adrenergic receptor (alpha 2AR) subtypes. We have cloned a human alpha 2AR by using the polymerase chain reaction with oligonucleotide primers homologous to conserved regions of the previously cloned alpha 2ARs, the genes for which are located on human chromosomes 4 (C4) and 10 (C10). The deduced amino acid sequence encodes a protein of 450 amino acids whose putative topology is similar to that of the family of guanine nucleotide-binding protein-coupled receptors, but whose structure most closely resembles that of the alpha 2ARs. Competition curve analysis of the binding properties of the receptor expressed in COS-7 cells with a variety of adrenergic ligands demonstrates a unique alpha 2AR pharmacology. Hybridization with somatic cell hybrids shows that the gene for this receptor is located on chromosome 2. Northern blot analysis of various rat tissues shows expression in liver and kidney. The unique pharmacology and tissue localization of this receptor suggest that this is an alpha 2AR subtype not previously identified by classical pharmacological or ligand binding approaches.

Regions of the alpha 1-adrenergic receptor involved in coupling to phosphatidylinositol hydrolysis and enhanced sensitivity of biological function.

Regions of the hamster alpha 1-adrenergic receptor (alpha 1 AR) that are important in GTP-binding protein (G protein)-mediated activation of phospholipase C were determined by studying the biological functions of mutant receptors constructed by recombinant DNA techniques. A chimeric receptor consisting of the beta 2-adrenergic receptor (beta 2AR) into which the putative third cytoplasmic loop of the alpha 1AR had been placed activated phosphatidylinositol metabolism as effectively as the native alpha 1AR, as did a truncated alpha 1AR lacking the last 47 residues in its cytoplasmic tail. Substitutions of beta 2AR amino acid sequence in the intermediate portions of the third cytoplasmic loop of the alpha 1AR or at the N-terminal portion of the cytoplasmic tail caused marked decreases in receptor coupling to phospholipase C. Conservative substitutions of two residues in the C terminus of the third cytoplasmic loop (Ala293----Leu, Lys290----His) increased the potency of agonists for stimulating phosphatidylinositol metabolism by up to 2 orders of magnitude. These data indicate (i) that the regions of the alpha 1AR that determine coupling to phosphatidylinositol metabolism are similar to those previously shown to be involved in coupling of beta 2AR to adenylate cyclase stimulation and (ii) that point mutations of a G-protein-coupled receptor can cause remarkable increases in sensitivity of biological response.

Reconstruction of gross avian genome structure, organization and evolution suggests that the chicken lineage most closely resembles the dinosaur avian ancestor.

BACKGROUND: The availability of multiple avian genome sequence assemblies greatly improves our ability to define overall genome organization and reconstruct evolutionary changes. In birds, this has previously been impeded by a near intractable karyotype and relied almost exclusively on comparative molecular cytogenetics of only the largest chromosomes. Here, novel whole genome sequence information from 21 avian genome sequences (most newly assembled) made available on an interactive browser (Evolution Highway) was analyzed. RESULTS: Focusing on the six best-assembled genomes allowed us to assemble a putative karyotype of the dinosaur ancestor for each chromosome. Reconstructing evolutionary events that led to each species' genome organization, we determined that the fastest rate of change occurred in the zebra finch and budgerigar, consistent with rapid speciation events in the Passeriformes and Psittaciformes. Intra- and interchromosomal changes were explained most parsimoniously by a series of inversions and translocations respectively, with breakpoint reuse being commonplace. Analyzing chicken and zebra finch, we found little evidence to support the hypothesis of an association of evolutionary breakpoint regions with recombination hotspots but some evidence to support the hypothesis that microchromosomes largely represent conserved blocks of synteny in the majority of the 21 species analyzed. All but one species showed the expected number of microchromosomal rearrangements predicted by the haploid chromosome count. Ostrich, however, appeared to retain an overall karyotype structure of 2n=80 despite undergoing a large number (26) of hitherto un-described interchromosomal changes. CONCLUSIONS: Results suggest that mechanisms exist to preserve a static overall avian karyotype/genomic structure, including the microchromosomes, with widespread interchromosomal change occurring rarely (e.g., in ostrich and budgerigar lineages). Of the species analyzed, the chicken lineage appeared to have undergone the fewest changes compared to the dinosaur ancestor.