Self-reported Sleep Problems Related to Amyloid Deposition in Cortical Regions with High HOMER1 Gene Expression

Sleep problems are related to the elevated levels of the Alzheimer's disease (AD) biomarker β-amyloid (Aβ). Hypotheses about the causes of this relationship can be generated from molecular markers of sleep problems identified in rodents. A major marker of sleep deprivation is Homer1a, a neural protein coded by the HOMER1 gene, which has also been implicated in brain Aβ accumulation. Here, we tested whether the relationship between cortical Aβ accumulation and self-reported sleep quality, as well as changes in sleep quality over 3 years, was stronger in cortical regions with high HOMER1 mRNA expression levels. In a sample of 154 cognitively healthy older adults, Aβ correlated with poorer sleep quality cross-sectionally and longitudinally (n = 62), but more strongly in the younger than in older individuals. Effects were mainly found in regions with high expression of HOMER1. The anatomical distribution of the sleep-Aβ relationship followed closely the Aβ accumulation pattern in 69 patients with mild cognitive impairment or AD. Thus, the results indicate that the relationship between sleep problems and Aβ accumulation may involve Homer1 activity in the cortical regions, where harbor Aβ deposits in AD. The findings may advance our understanding of the relationship between sleep problems and AD risk.

EOMES and IL-10 regulate antitumor activity of T regulatory type 1 CD4+ T cells in chronic lymphocytic leukemia

The transcription factor eomesodermin (EOMES) promotes interleukin (IL)-10 expression in CD4+ T cells, which has been linked to immunosuppressive and cytotoxic activities. We detected cytotoxic, programmed cell death protein-1 (PD-1) and EOMES co-expressing CD4+ T cells in lymph nodes (LNs) of patients with chronic lymphocytic leukemia (CLL) or diffuse large B-cell lymphoma. Transcriptome and flow cytometry analyses revealed that EOMES does not only drive IL-10 expression, but rather controls a unique transcriptional signature in CD4+ T cells, that is enriched in genes typical for T regulatory type 1 (TR1) cells. The TR1 cell identity of these CD4+ T cells was supported by their expression of interferon gamma and IL-10, as well as inhibitory receptors including PD-1. TR1 cells with cytotoxic capacity accumulate also in Eµ-TCL1 mice that develop CLL-like disease. Whereas wild-type CD4+ T cells control TCL1 leukemia development after adoptive transfer in leukopenic Rag2-/- mice, EOMES-deficient CD4+ T cells failed to do so. We further show that TR1 cell-mediated control of TCL1 leukemia requires IL-10 receptor (IL-10R) signaling, as Il10rb-deficient CD4+ T cells showed impaired antileukemia activity. Altogether, our data demonstrate that EOMES is indispensable for the development of IL-10-expressing, cytotoxic TR1 cells, which accumulate in LNs of CLL patients and control TCL1 leukemia in mice in an IL-10R-dependent manner.

Identifying genetic modifiers of age-associated penetrance in X-linked dystonia-parkinsonism

X-linked dystonia-parkinsonism is a neurodegenerative disorder caused by a founder retrotransposon insertion, in which a polymorphic hexanucleotide repeat accounts for ~50% of age at onset variability. Employing a genome-wide association study to identify additional factors modifying age at onset, we establish that three independent loci are significantly associated with age at onset (p < 5 × 10⁻⁸). The lead single nucleotide polymorphisms collectively account for 25.6% of the remaining variance not explained by the hexanucleotide repeat and 13.0% of the overall variance in age at onset in X-linked dystonia-parkinsonism with the protective alleles delaying disease onset by seven years. These regions harbor or lie adjacent to MSH3 and PMS2, the genes that were recently implicated in modifying age at onset in Huntington’s disease, likely through a common pathway influencing repeat instability. Our work indicates the existence of three modifiers of age at onset in X-linked dystonia-parkinsonism that likely affect the DNA mismatch repair pathway.

Age at onset of X-linked dystonia-parkinsonism is 50% explained by the length of a repeat in an SVA insert. The authors perform a GWAS for genetic modifiers and discover three more loci, accounting for another 13% of variability in age at onset with the protective alleles delaying onset by seven years.

Pathophysiological subtypes of Alzheimer's disease based on cerebrospinal fluid proteomics

Alzheimer's disease is biologically heterogeneous, and detailed understanding of the processes involved in patients is critical for development of treatments. CSF contains hundreds of proteins, with concentrations reflecting ongoing (patho)physiological processes. This provides the opportunity to study many biological processes at the same time in patients. We studied whether Alzheimer's disease biological subtypes can be detected in CSF proteomics using the dual clustering technique non-negative matrix factorization. In two independent cohorts (EMIF-AD MBD and ADNI) we found that 705 (77% of 911 tested) proteins differed between Alzheimer's disease (defined as having abnormal amyloid, n = 425) and controls (defined as having normal CSF amyloid and tau and normal cognition, n = 127). Using these proteins for data-driven clustering, we identified three robust pathophysiological Alzheimer's disease subtypes within each cohort showing (i) hyperplasticity and increased BACE1 levels; (ii) innate immune activation; and (iii) blood-brain barrier dysfunction with low BACE1 levels. In both cohorts, the majority of individuals were labelled as having subtype 1 (80, 36% in EMIF-AD MBD; 117, 59% in ADNI), 71 (32%) in EMIF-AD MBD and 41 (21%) in ADNI were labelled as subtype 2, and 72 (32%) in EMIF-AD MBD and 39 (20%) individuals in ADNI were labelled as subtype 3. Genetic analyses showed that all subtypes had an excess of genetic risk for Alzheimer's disease (all P > 0.01). Additional pathological comparisons that were available for a subset in ADNI suggested that subtypes showed similar severity of Alzheimer's disease pathology, and did not differ in the frequencies of co-pathologies, providing further support that found subtypes truly reflect Alzheimer's disease heterogeneity. Compared to controls, all non-demented Alzheimer's disease individuals had increased risk of showing clinical progression (all P < 0.01). Compared to subtype 1, subtype 2 showed faster clinical progression after correcting for age, sex, level of education and tau levels (hazard ratio = 2.5; 95% confidence interval = 1.2, 5.1; P = 0.01), and subtype 3 at trend level (hazard ratio = 2.1; 95% confidence interval = 1.0, 4.4; P = 0.06). Together, these results demonstrate the value of CSF proteomics in studying the biological heterogeneity in Alzheimer's disease patients, and suggest that subtypes may require tailored therapy.

Genome-wide association study of Alzheimer's disease CSF biomarkers in the EMIF-AD Multimodal Biomarker Discovery dataset

Alzheimer's disease (AD) is the most prevalent neurodegenerative disorder and the most common form of dementia in the elderly. Susceptibility to AD is considerably determined by genetic factors which hitherto were primarily identified using case-control designs. Elucidating the genetic architecture of additional AD-related phenotypic traits, ideally those linked to the underlying disease process, holds great promise in gaining deeper insights into the genetic basis of AD and in developing better clinical prediction models. To this end, we generated genome-wide single-nucleotide polymorphism (SNP) genotyping data in 931 participants of the European Medical Information Framework Alzheimer's Disease Multimodal Biomarker Discovery (EMIF-AD MBD) sample to search for novel genetic determinants of AD biomarker variability. Specifically, we performed genome-wide association study (GWAS) analyses on 16 traits, including 14 measures derived from quantifications of five separate amyloid-beta (Aβ) and tau-protein species in the cerebrospinal fluid (CSF). In addition to confirming the well-established effects of apolipoprotein E (APOE) on diagnostic outcome and phenotypes related to Aβ42, we detected novel potential signals in the zinc finger homeobox 3 (ZFHX3) for CSF-Aβ38 and CSF-Aβ40 levels, and confirmed the previously described sex-specific association between SNPs in geminin coiled-coil domain containing (GMNC) and CSF-tau. Utilizing the results from independent case-control AD GWAS to construct polygenic risk scores (PRS) revealed that AD risk variants only explain a small fraction of CSF biomarker variability. In conclusion, our study represents a detailed first account of GWAS analyses on CSF-Aβ and -tau-related traits in the EMIF-AD MBD dataset. In subsequent work, we will utilize the genomics data generated here in GWAS of other AD-relevant clinical outcomes ascertained in this unique dataset.

Genetic risk for Alzheimer disease predicts hippocampal volume through the human lifespan

Objective

To test the hypothesis that genetic risk for Alzheimer disease (AD) may represent a stable influence on the brain from early in life, rather than being primarily age dependent, we investigated in a lifespan sample of 1,181 persons with a total of 2,690 brain scans, whether higher polygenic risk score (PGS) for AD and presence of APOE ε4 was associated with lower hippocampal volumes to begin with, as an offset effect, or possibly faster decline in older age.

Methods

Using general additive mixed models, we assessed the relations of PGS for AD, including variants in APOE with hippocampal volume and its change in a cognitively healthy longitudinal lifespan sample (age range: 4–95 years, mean visit age 39.7 years, SD 26.9 years), followed for up to 11 years.

Results

AD-PGS and APOE ε4 in isolation showed a significant negative effect on hippocampal volume. The effect of a 1 sample SD increase in AD-PGS on hippocampal volume was estimated to –36.4 mm³ (confidence interval [CI]: –71.8, –1.04) and the effect of carrying ε4 allele(s) –107.0 mm³ (CI: –182.0, –31.5). Offset effects of AD-PGS and APOE ε4 were present in hippocampal development, and interactions between age and genetic risk on volume change were not consistently observed.

Conclusions

Endophenotypic manifestation of polygenic risk for AD may be seen across the lifespan in cognitively healthy persons, not being confined to clinical populations or older age. This emphasizes that a broader population and age range may be relevant targets for attempts to prevent AD.

snakePipes: facilitating flexible, scalable and integrative epigenomic analysis

SUMMARY

Due to the rapidly increasing scale and diversity of epigenomic data, modular and scalable analysis workflows are of wide interest. Here we present snakePipes, a workflow package for processing and downstream analysis of data from common epigenomic assays: ChIP-seq, RNA-seq, Bisulfite-seq, ATAC-seq, Hi-C and single-cell RNA-seq. snakePipes enables users to assemble variants of each workflow and to easily install and upgrade the underlying tools, via its simple command-line wrappers and yaml files.

AVAILABILITY AND IMPLEMENTATION

snakePipes can be installed via conda: `conda install -c mpi-ie -c bioconda -c conda-forge snakePipes'. Source code (https://github.com/maxplanck-ie/snakepipes) and documentation (https://snakepipes.readthedocs.io/en/latest/) are available online.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Dickkopf-1 Overexpression in vitro Nominates Candidate Blood Biomarkers Relating to Alzheimer's Disease Pathology

BACKGROUND

Previous studies suggest that Dickkopf-1 (DKK1), an inhibitor of Wnt signaling, plays a role in amyloid-induced toxicity and hence Alzheimer's disease (AD). However, the effect of DKK1 expression on protein expression, and whether such proteins are altered in disease, is unknown.

OBJECTIVE

We aim to test whether DKK1 induced protein signature obtained in vitro were associated with markers of AD pathology as used in the amyloid/tau/neurodegeneration (ATN) framework as well as with clinical outcomes.

METHODS

We first overexpressed DKK1 in HEK293A cells and quantified 1,128 proteins in cell lysates using aptamer capture arrays (SomaScan) to obtain a protein signature induced by DKK1. We then used the same assay to measure the DKK1-signature proteins in human plasma in two large cohorts, EMIF (n = 785) and ANM (n = 677).

RESULTS

We identified a 100-protein signature induced by DKK1 in vitro. Subsets of proteins, along with age and apolipoprotein E ɛ4 genotype distinguished amyloid pathology (A + T-N-, A+T+N-, A+T-N+, and A+T+N+) from no AD pathology (A-T-N-) with an area under the curve of 0.72, 0.81, 0.88, and 0.85, respectively. Furthermore, we found that some signature proteins (e.g., Complement C3 and albumin) were associated with cognitive score and AD diagnosis in both cohorts.

CONCLUSIONS

Our results add further evidence for a role of DKK regulation of Wnt signaling in AD and suggest that DKK1 induced signature proteins obtained in vitro could reflect theATNframework as well as predict disease severity and progression in vivo.

Discovery and validation of plasma proteomic biomarkers relating to brain amyloid burden by SOMAscan assay

Introduction

Plasma proteins have been widely studied as candidate biomarkers to predict brain amyloid deposition to increase recruitment efficiency in secondary prevention clinical trials for Alzheimer's disease. Most such biomarker studies are targeted to specific proteins or are biased toward high abundant proteins.

Methods

4001 plasma proteins were measured in two groups of participants (discovery group = 516, replication group = 365) selected from the European Medical Information Framework for Alzheimer's disease Multimodal Biomarker Discovery study, all of whom had measures of amyloid.

Results

A panel of proteins (n = 44), along with age and apolipoprotein E (APOE) ε4, predicted brain amyloid deposition with good performance in both the discovery group (area under the curve = 0.78) and the replication group (area under the curve = 0.68). Furthermore, a causal relationship between amyloid and tau was confirmed by Mendelian randomization.

Discussion

The results suggest that high-dimensional plasma protein testing could be a useful and reproducible approach for measuring brain amyloid deposition.

DOT1L promotes progenitor proliferation and primes neuronal layer identity in the developing cerebral cortex

Cortical development is controlled by transcriptional programs, which are orchestrated by transcription factors. Yet, stable inheritance of spatio-temporal activity of factors influencing cell fate and localization in different layers is only partly understood. Here we find that deletion of Dot1l in the murine telencephalon leads to cortical layering defects, indicating DOT1L activity and chromatin methylation at H3K79 impact on the cell cycle, and influence transcriptional programs conferring upper layer identity in early progenitors. Specifically, DOT1L prevents premature differentiation by increasing expression of genes that regulate asymmetric cell division (Vangl2, Cenpj). Loss of DOT1L results in reduced numbers of progenitors expressing genes including SoxB1 gene family members. Loss of DOT1L also leads to altered cortical distribution of deep layer neurons that express either TBR1, CTIP2 or SOX5, and less activation of transcriptional programs that are characteristic for upper layer neurons (Satb2, Pou3f3, Cux2, SoxC family members). Data from three different mouse models suggest that DOT1L balances transcriptional programs necessary for proper neuronal composition and distribution in the six cortical layers. Furthermore, because loss of DOT1L in the pre-neurogenic phase of development impairs specifically generation of SATB2-expressing upper layer neurons, our data suggest that DOT1L primes upper layer identity in cortical progenitors.

Alzheimer's disease risk SNPs show no strong effect on miRNA expression in human lymphoblastoid cell lines

The role of microRNAs (miRNAs) in the pathogenesis of Alzheimer's disease (AD) is currently extensively investigated. In this study, we assessed the potential impact of AD genetic risk variants on miRNA expression by performing large-scale bioinformatic data integration. Our analysis was based on genetic variants from 3 AD genome-wide association studies (GWASs). Association with miRNA expression was tested by expression quantitative trait locus analysis using next-generation miRNA sequencing data generated in lymphoblastoid cell lines. Although, overall, we did not identify a strong effect of AD GWAS variants on miRNA expression in this cell type, we highlight 2 notable outliers, that is, miR-29c-5p and miR-6840-5p. MiR-29c-5p was recently reported to be involved in the regulation of BACE1 and SORL1 expression. In conclusion, despite 2 exceptions, our large-scale assessment provides only limited support for the hypothesis that AD GWAS variants act as miRNA expression quantitative trait loci.

Critical evaluation of the DNA-methylation markers ABCG1 and SREBF1 for Type 2 diabetes stratification

Aim: Validation of epigenome-wide association studies is sparse. Therefore, we evaluated the methylation markers cg06500161 (ABCG1) and cg11024682 (SREBF1) as classifiers for diabetes stratification. Patients & methods: DNA methylation was measured in blood (n = 167), liver (n = 99) and visceral adipose tissue (n = 99) of nondiabetic or Type 2 diabetic subjects by bisulfite pyrosequencing. Results: DNA methylation at cg11024682 in blood and liver correlated with BMI. Methylation at cg06500161 was influenced by the adjacent SNP rs9982016. Insulin-resistant and sensitive subjects could be stratified by DNA methylation status in blood or visceral adipose tissue. Conclusion: DNA methylation at both loci in blood presents a promising approach for risk group stratification and could be valuable for personalized Type 2 diabetes risk prediction in the future.

Differential Methylation of H3K79 Reveals DOT1L Target Genes and Function in the Cerebellum In Vivo

The disruptor of telomeric silencing 1-like (DOT1L) mediates methylation of histone H3 at position lysine 79 (H3K79). Conditional knockout of Dot1l in mouse cerebellar granule cells (Dot1l-cKO^Atoh1) led to a smaller external granular layer with fewer precursors of granule neurons. Dot1l-cKO^Atoh1 mice had impaired proliferation and differentiation of granular progenitors, which resulted in a smaller cerebellum. Mutant mice showed mild ataxia in motor behavior tests. In contrast, Purkinje cell-specific conditional knockout mice showed no obvious phenotype. Genome-wide transcription analysis of Dot1l-cKO^Atoh1 cerebella using microarrays revealed changes in genes that function in cell cycle, cell migration, axon guidance, and metabolism. To identify direct DOT1L target genes, we used genome-wide profiling of H3K79me2 and transcriptional analysis. Analysis of differentially methylated regions (DR) and differentially expressed genes (DE) revealed in total 12 putative DOT1L target genes in Dot1l-cKO^Atoh1 affecting signaling (Tnfaip8l3, B3galt5), transcription (Otx1), cell migration and axon guidance (Sema4a, Sema5a, Robo1), cholesterol and lipid metabolism (Lss, Cyp51), cell cycle (Cdkn1a), calcium-dependent cell-adhesion or exocytosis (Pcdh17, Cadps2), and unknown function (Fam174b). Dysregulated expression of these target genes might be implicated in the ataxia phenotype observed in Dot1l-cKO^Atoh1.

deepTools2: a next generation web server for deep-sequencing data analysis

We present an update to our Galaxy-based web server for processing and visualizing deeply sequenced data. Its core tool set, deepTools, allows users to perform complete bioinformatic workflows ranging from quality controls and normalizations of aligned reads to integrative analyses, including clustering and visualization approaches. Since we first described our deepTools Galaxy server in 2014, we have implemented new solutions for many requests from the community and our users. Here, we introduce significant enhancements and new tools to further improve data visualization and interpretation. deepTools continue to be open to all users and freely available as a web service at deeptools.ie-freiburg.mpg.de The new deepTools2 suite can be easily deployed within any Galaxy framework via the toolshed repository, and we also provide source code for command line usage under Linux and Mac OS X. A public and documented API for access to deepTools functionality is also available.

High-resolution SAR11 ecotype dynamics at the Bermuda Atlantic Time-series Study site by phylogenetic placement of pyrosequences

Advances in next-generation sequencing technologies are providing longer nucleotide sequence reads that contain more information about phylogenetic relationships. We sought to use this information to understand the evolution and ecology of bacterioplankton at our long-term study site in the Western Sargasso Sea. A bioinformatics pipeline called PhyloAssigner was developed to align pyrosequencing reads to a reference multiple sequence alignment of 16S ribosomal RNA (rRNA) genes and assign them phylogenetic positions in a reference tree using a maximum likelihood algorithm. Here, we used this pipeline to investigate the ecologically important SAR11 clade of Alphaproteobacteria. A combined set of 2.7 million pyrosequencing reads from the 16S rRNA V1-V2 regions, representing 9 years at the Bermuda Atlantic Time-series Study (BATS) site, was quality checked and parsed into a comprehensive bacterial tree, yielding 929 036 Alphaproteobacteria reads. Phylogenetic structure within the SAR11 clade was linked to seasonally recurring spatiotemporal patterns. This analysis resolved four new SAR11 ecotypes in addition to five others that had been described previously at BATS. The data support a conclusion reached previously that the SAR11 clade diversified by subdivision of niche space in the ocean water column, but the new data reveal a more complex pattern in which deep branches of the clade diversified repeatedly across depth strata and seasonal regimes. The new data also revealed the presence of an unrecognized clade of Alphaproteobacteria, here named SMA-1 (Sargasso Mesopelagic Alphaproteobacteria, group 1), in the upper mesopelagic zone. The high-resolution phylogenetic analyses performed herein highlight significant, previously unknown, patterns of evolutionary diversification, within perhaps the most widely distributed heterotrophic marine bacterial clade, and strongly links to ecosystem regimes.

Multiple rearrangements in mitochondrial genomes of Isopoda and phylogenetic implications

In this study, we analyse the evolutionary dynamics and phylogenetic implications of gene order rearrangements in five newly sequenced mitochondrial (mt) genomes and four published mt genomes of isopod crustaceans. The sequence coverage is nearly complete for four of the five newly sequenced species, with only the control region and some tRNA genes missing, while in Janira maculosa only two thirds of the genome could be determined. Mitochondrial gene order in isopods seems to be more plastic than that in other crustacean lineages, making all nine known mt gene orders different. Especially the asellote Janira is characterized by many autapomorphies. The following inferred ancestral isopod mt gene order exists slightly modified in modern isopods: nad1, tnrL1, rrnS, control region, trnS1, cob, trnT, nad5, trnF. We consider the inferred gene translocation events leading to gene rearrangements as valuable characters in phylogenetic analyses. In this first study covering major isopod lineages, potential apomorphies were identified, e.g., a shared relative position of trnR in Valvifera. We also report one of the first findings of homoplasy in mitochondrial gene order, namely a shared relative position of trnV in unrelated isopod lineages. In addition to increased taxon sampling secondary structure, modification in tRNAs and GC-skew inversion may be potentially fruitful subjects for future mt genome studies in a phylogenetic context.

The Australian fresh water isopod (Phreatoicidea: Isopoda) allows insights into the early mitogenomic evolution of isopods

The complete mitochondrial (mt) genome sequence of the Australian fresh water isopod Eophreatoicus sp.-14 has been determined. The new species is a member of the taxon Phreatoicidea, a clade of particular interest, as it is often regarded as the sister group to all other Isopoda. Although the overall genome organization of Eophreatoicus sp.-14 conforms to the typical state of Metazoa--it is a circular ring of DNA hosting the usual 37 genes and one major non-coding region--it bears a number of derived characters that fall within the scope of "genome morphology". Earlier studies have indicated that the isopod mitochondrial gene order is not as conserved as that of other crustaceans. Indeed, the mt genome of Eophreatoicus sp.-14 shows an inversion of seven genes (including cox1), which is as far as we know unique. Even more interesting is the derived arrangement of nad1, trnL(CUN), rrnS, control region, cob, trnT, nad5 and trnF that is shared by nearly all available isopod mt genomes. A striking feature is the close proximity of the rearranged genes to the mt control region. Inferable gene translocation events are, however, more suitable to trace the evolution of mt genomes. Genes like nad1/trnL(CUN) and nad5/trnF, which retained their adjacent position after being rearranged, were most likely translocated together. A very good example for the need to understand the mechanisms of translocations is the remolding of trnL(UUR) to trnL(CUN). Both tRNA genes are adjacent and have a high sequence similarity, probably the result of a gene duplication and subsequent anticodon mutation. Modified secondary structures were found in three tRNAs of Eophreatoicus sp.-14, which are all characterized by the loss of the DHU-arm. This is common to crustaceans for tRNA Serine(AGY), while the arm-loss in tRNA Cysteine within Malacostraca is only shared by other isopods. Modification of the third tRNA, Isoleucine, is not known from any other related species. Nucleotide frequencies of genes have been found to be indirectly correlated to the orientation of the mitochondrial replication process. In Eophreatoicus sp.-14 and in other Isopoda the associated nucleotide bias is inversed to the state of other Malacostraca. This is a strong indication for an inversion of the control region that most likely evolved in the isopod ancestor.

The complete mitochondrial genome of the common sea slater, Ligia oceanica (Crustacea, Isopoda) bears a novel gene order and unusual control region features

Background

Sequence data and other characters from mitochondrial genomes (gene translocations, secondary structure of RNA molecules) are useful in phylogenetic studies among metazoan animals from population to phylum level. Moreover, the comparison of complete mitochondrial sequences gives valuable information about the evolution of small genomes, e.g. about different mechanisms of gene translocation, gene duplication and gene loss, or concerning nucleotide frequency biases.

The Peracarida (gammarids, isopods, etc.) comprise about 21,000 species of crustaceans, living in many environments from deep sea floor to arid terrestrial habitats. Ligia oceanica is a terrestrial isopod living at rocky seashores of the european North Sea and Atlantic coastlines.

Results

The study reveals the first complete mitochondrial DNA sequence from a peracarid crustacean. The mitochondrial genome of Ligia oceanica is a circular double-stranded DNA molecule, with a size of 15,289 bp. It shows several changes in mitochondrial gene order compared to other crustacean species. An overview about mitochondrial gene order of all crustacean taxa yet sequenced is also presented. The largest non-coding part (the putative mitochondrial control region) of the mitochondrial genome of Ligia oceanica is unexpectedly not AT-rich compared to the remainder of the genome. It bears two repeat regions (4× 10 bp and 3× 64 bp), and a GC-rich hairpin-like secondary structure. Some of the transfer RNAs show secondary structures which derive from the usual cloverleaf pattern. While some tRNA genes are putative targets for RNA editing, trnR could not be localized at all.

Conclusion

Gene order is not conserved among Peracarida, not even among isopods. The two isopod species Ligia oceanica and Idotea baltica show a similarly derived gene order, compared to the arthropod ground pattern and to the amphipod Parhyale hawaiiensis, suggesting that most of the translocation events were already present the last common ancestor of these isopods. Beyond that, the positions of three tRNA genes differ in the two isopod species. Strand bias in nucleotide frequency is reversed in both isopod species compared to other Malacostraca. This is probably due to a reversal of the replication origin, which is further supported by the fact that the hairpin structure typically found in the control region shows a reversed orientation in the isopod species, compared to other crustaceans.