IGHV allele similarity clustering improves genotype inference from adaptive immune receptor repertoire sequencing data.

In adaptive immune receptor repertoire analysis, determining the germline variable (V) allele associated with each T- and B-cell receptor sequence is a crucial step. This process is highly impacted by allele annotations. Aligning sequences, assigning them to specific germline alleles, and inferring individual genotypes are challenging when the repertoire is highly mutated, or sequence reads do not cover the whole V region. Here, we propose an alternative naming scheme for the V alleles, as well as a novel method to infer individual genotypes. We demonstrate the strengths of the two by comparing their outcomes to other genotype inference methods. We validate the genotype approach with independent genomic long-read data. The naming scheme is compatible with current annotation tools and pipelines. Analysis results can be converted from the proposed naming scheme to the nomenclature determined by the International Union of Immunological Societies (IUIS). Both the naming scheme and the genotype procedure are implemented in a freely available R package (PIgLET https://bitbucket.org/yaarilab/piglet). To allow researchers to further explore the approach on real data and to adapt it for their uses, we also created an interactive website (https://yaarilab.github.io/IGHV_reference_book).

Integrated Proteomics Analysis of Baseline Protein Expression in Pig Tissues.

The availability of an increasingly large amount of public proteomics data sets presents an opportunity for performing combined analyses to generate comprehensive organism-wide protein expression maps across different organisms and biological conditions. Sus scrofa, a domestic pig, is a model organism relevant for food production and for human biomedical research. Here, we reanalyzed 14 public proteomics data sets from the PRIDE database coming from pig tissues to assess baseline (without any biological perturbation) protein abundance in 14 organs, encompassing a total of 20 healthy tissues from 128 samples. The analysis involved the quantification of protein abundance in 599 mass spectrometry runs. We compared protein expression patterns among different pig organs and examined the distribution of proteins across these organs. Then, we studied how protein abundances were compared across different data sets and studied the tissue specificity of the detected proteins. Of particular interest, we conducted a comparative analysis of protein expression between pig and human tissues, revealing a high degree of correlation in protein expression among orthologs, particularly in brain, kidney, heart, and liver samples. We have integrated the protein expression results into the Expression Atlas resource for easy access and visualization of the protein expression data individually or alongside gene expression data.

Change in cyclone disaster vulnerability and response in coastal Bangladesh.

The number of deaths owing to tropical cyclones in Bangladesh has significantly reduced. Category 4 Cyclone Gorky in 1991 and Sidr in 2007 caused 147,000 and 4,500 deaths respectively, whereas Category 1 Cyclone Mora in 2017 resulted in six. Face-to-face interviews with 362 residents, participant observation, and focus-group discussions answer a research question about how change in coastal areas has contributed to this outcome. The study considered institutional approaches of disaster risk management through legal frameworks, administrative arrangements, cyclone preparedness activities, cyclone detection and early warning dissemination, construction of shelter centres, strengthening of various types of coastal embankments, paved roads, and pre-cyclone evacuation. The findings indicate significant improvement in house structures and design, income levels and diversification, education, awareness, individual capacity, poverty reduction, and lowering dependency on agriculture-based earning. Furthermore, the availability of mobile telephones, radio, television, and social media platforms enhanced social connectivity and greater gender equality and empowerment helped to facilitate disaster preparedness, evacuation, and response.

Assessing Multiple Evidence Streams to Decide on Confidence for Identification of Post-Translational Modifications, within and Across Data Sets.

Phosphorylation is a post-translational modification of great interest to researchers due to its relevance in many biological processes. LC-MS/MS techniques have enabled high-throughput data acquisition, with studies claiming identification and localization of thousands of phosphosites. The identification and localization of phosphosites emerge from different analytical pipelines and scoring algorithms, with uncertainty embedded throughout the pipeline. For many pipelines and algorithms, arbitrary thresholding is used, but little is known about the actual global false localization rate in these studies. Recently, it has been suggested to use decoy amino acids to estimate global false localization rates of phosphosites, among the peptide-spectrum matches reported. Here, we describe a simple pipeline aiming to maximize the information extracted from these studies by objectively collapsing from peptide-spectrum match to the peptidoform-site level, as well as combining findings from multiple studies while maintaining track of false localization rates. We show that the approach is more effective than current processes that use a simpler mechanism for handling phosphosite identification redundancy within and across studies. In our case study using eight rice phosphoproteomics data sets, 6368 unique sites were confidently identified using our decoy approach compared to 4687 using traditional thresholding in which false localization rates are unknown.

Integrated View of Baseline Protein Expression in Human Tissues.

The availability of proteomics datasets in the public domain, and in the PRIDE database, in particular, has increased dramatically in recent years. This unprecedented large-scale availability of data provides an opportunity for combined analyses of datasets to get organism-wide protein abundance data in a consistent manner. We have reanalyzed 24 public proteomics datasets from healthy human individuals to assess baseline protein abundance in 31 organs. We defined tissue as a distinct functional or structural region within an organ. Overall, the aggregated dataset contains 67 healthy tissues, corresponding to 3,119 mass spectrometry runs covering 498 samples from 489 individuals. We compared protein abundances between different organs and studied the distribution of proteins across these organs. We also compared the results with data generated in analogous studies. Additionally, we performed gene ontology and pathway-enrichment analyses to identify organ-specific enriched biological processes and pathways. As a key point, we have integrated the protein abundance results into the resource Expression Atlas, where they can be accessed and visualized either individually or together with gene expression data coming from transcriptomics datasets. We believe this is a good mechanism to make proteomics data more accessible for life scientists.

Predicting Pancreatic Cancer in the UK Biobank Cohort Using Polygenic Risk Scores and Diabetes Mellitus.

BACKGROUND & AIMS: Diabetes mellitus (DM) is known to be associated with pancreatic ductal adenocarcinoma (PDAC), particularly new-onset DM (NODM). Others have developed polygenic risk scores (PRS) associated with PDAC risk. We aimed to compare the performance of these PRS in an independent cohort to determine if they can discriminate between NODM and long-standing DM patients with PDAC. METHODS: Cases (1042) and matched cancer-free controls (10,420) were drawn from the UK Biobank. Five PRS models were calculated using single nucleotide polymorphisms (SNPs) from previous studies (Nakatochi, Galeotti, Molina, Jia, and Rashkin) and a combination of these. Regression models were used to assess the association between PDAC and PRS adjusted for ancestry, smoking, DM, waist circumference, and family history of digestive cancer. Receiver operator characteristic curves and area under the curve metrics (AUC) were used to assess the performance of each PRS for classifying PDAC risk. RESULTS: The combined PRS model achieved the highest AUC (0.605), and significantly improved a clinical risk model in this cohort (AUC = 0.83; P = .0002). Individuals within the fifth quintile have a 2.74-fold increased risk of developing PDAC vs those in the first quintile (P < .001), and have a 3.05-fold increased risk of developing PDAC if they have DM vs those without DM (P < .001). The positive predictive value was 11.9% in participants without DM, 23.9% with long-standing DM, and 86.7% with NODM. CONCLUSIONS: The PDAC-related common genetic variants are more strongly associated with DM. This PRS has the potential for targeting individuals with NODM for PDAC secondary screening measures.

Essentiality-specific pathogenicity prioritization gene score to improve filtering of disease sequence data.

The causal genetic variants underlying more than 50% of single gene (monogenic) disorders are yet to be discovered. Many patients with conditions likely to have a monogenic basis do not receive a confirmed molecular diagnosis which has potential impacts on clinical management. We have developed a gene-specific score, essentiality-specific pathogenicity prioritization (ESPP), to guide the recognition of genes likely to underlie monogenic disease variation to assist in filtering of genome sequence data. When a patient genome is sequenced, there are frequently several plausibly pathogenic variants identified in different genes. Recognition of the single gene most likely to include pathogenic variation can guide the identification of a causal variant. The ESPP score integrates gene-level scores which are broadly related to gene essentiality. Previous work towards the recognition of monogenic disease genes proposed a model with increasing gene essentiality from 'non-essential' to 'essential' genes (for which pathogenic variation may be incompatible with survival) with genes liable to contain disease variation positioned between these two extremes. We demonstrate that the ESPP score is useful for recognizing genes with high potential for pathogenic disease-related variation. Genes classed as essential have particularly high scores, as do genes recently recognized as strong candidates for developmental disorders. Through the integration of individual gene-specific scores, which have different properties and assumptions, we demonstrate the utility of an essentiality-based gene score to improve sequence genome filtering.

Visual comet scoring revisited: a guide to scoring comet assay slides and obtaining reliable results.

Measurement of DNA migration in the comet assay can be done by image analysis or visual scoring. The latter accounts for 20%-25% of the published comet assay results. Here we assess the intra- and inter-investigator variability in visual scoring of comets. We include three training sets of comet images, which can be used as reference for researchers who wish to use visual scoring of comets. Investigators in 11 different laboratories scored the comet images using a five-class scoring system. There is inter-investigator variation in the three training sets of comets (i.e. coefficient of variation (CV) = 9.7%, 19.8%, and 15.2% in training sets I-III, respectively). However, there is also a positive correlation of inter-investigator scoring in the three training sets (r = 0.60). Overall, 36% of the variation is attributed to inter-investigator variation and 64% stems from intra-investigator variation in scoring between comets (i.e. the comets in training sets I-III look slightly different and this gives rise to heterogeneity in scoring). Intra-investigator variation in scoring was also assessed by repeated analysis of the training sets by the same investigator. There was larger variation when the training sets were scored over a period of six months (CV = 5.9%-9.6%) as compared to 1 week (CV = 1.3%-6.1%). A subsequent study revealed a high inter-investigator variation when premade slides, prepared in a central laboratory, were stained and scored by investigators in different laboratories (CV = 105% and 18%-20% in premade slides with comets from unexposed and hydrogen peroxide-exposed cells, respectively). The results indicate that further standardization of visual scoring is desirable. Nevertheless, the analysis demonstrates that visual scoring is a reliable way of analysing DNA migration in comets.

Inter-laboratory variation in measurement of DNA damage by the alkaline comet assay in the hCOMET ring trial.

The comet assay is a simple and versatile method for measurement of DNA damage in eukaryotic cells. More specifically, the assay detects DNA migration from agarose gel-embedded nucleoids, which depends on assay conditions and the level of DNA damage. Certain steps in the comet assay procedure have substantial impact on the magnitude of DNA migration (e.g. electric potential and time of electrophoresis). Inter-laboratory variation in DNA migration levels occurs because there is no agreement on optimal assay conditions or suitable assay controls. The purpose of the hCOMET ring trial was to test potassium bromate (KBrO3) as a positive control for the formamidopyrimidine DNA glycosylase (Fpg)-modified comet assay. To this end, participating laboratories used semi-standardized protocols for cell culture (i.e. cell culture, KBrO3 exposure, and cryopreservation of cells) and comet assay procedures, whereas the data acquisition was not standardized (i.e. staining of comets and image analysis). Segregation of the total variation into partial standard deviation (SD) in % Tail DNA units indicates the importance of cell culture procedures (SD = 10.9), comet assay procedures (SD = 12.3), staining (SD = 7.9) and image analysis (SD = 0.5) on the overall inter-laboratory variation of DNA migration (SD = 18.2). Future studies should assess sources of variation in each of these steps. On the positive side, the hCOMET ring trial demonstrates that KBrO3 is a robust positive control for the Fpg-modified comet assay. In conclusion, the hCOMET ring trial has demonstrated a high reproducibility of detecting genotoxic effects by the comet assay, but inter-laboratory variation of DNA migration levels is a concern.

DNA strand break levels in cryopreserved mononuclear blood cell lines measured by the alkaline comet assay: results from the hCOMET ring trial.

The comet assay is widely used in biomonitoring studies for the analysis of DNA damage in leukocytes and peripheral blood mononuclear cells. Rather than processing blood samples directly, it can be desirable to cryopreserve whole blood or isolated cells for later analysis by the comet assay. However, this creates concern about artificial accumulation of DNA damage during cryopreservation. In this study, 10 laboratories used standardized cryopreservation and thawing procedures of monocytic (THP-1) or lymphocytic (TK6) cells. Samples were cryopreserved in small aliquots in 50% foetal bovine serum, 40% cell culture medium, and 10% dimethyl sulphoxide. Subsequently, cryopreserved samples were analysed by the standard comet assay on three occasions over a 3-year period. Levels of DNA strand breaks in THP-1 cells were increased (four laboratories), unaltered (four laboratories), or decreased (two laboratories) by long-term storage. Pooled analysis indicates only a modest positive association between storage time and levels of DNA strand breaks in THP-1 cells (0.37% Tail DNA per year, 95% confidence interval: -0.05, 0.78). In contrast, DNA strand break levels were not increased by cryopreservation in TK6 cells. There was inter-laboratory variation in levels of DNA strand breaks in THP-1 cells (SD = 3.7% Tail DNA) and TK6 reference sample cells (SD = 9.4% Tail DNA), whereas the intra-laboratory residual variation was substantially smaller (i.e. SD = 0.4%-2.2% Tail DNA in laboratories with the smallest and largest variation). In conclusion, the study shows that accumulation of DNA strand breaks in cryopreserved mononuclear blood cell lines is not a matter of concern.

Long-term cryopreservation of potassium bromate positive assay controls for measurement of oxidatively damaged DNA by the Fpg-modified comet assay: results from the hCOMET ring trial.

The formamidopyrimidine DNA glycosylase (Fpg)-modified comet assay is widely used for the measurement of oxidatively generated damage to DNA. However, there has not been a recommended long-term positive control for this version of the comet assay. We have investigated potassium bromate as a positive control for the Fpg-modified comet assay because it generates many Fpg-sensitive sites with a little concurrent generation of DNA strand breaks. Eight laboratories used the same procedure for the treatment of monocytic THP-1 cells with potassium bromate (0, 0.5, 1.5, and 4.5 mM) and subsequent cryopreservation in a freezing medium consisting of 50% foetal bovine serum, 40% RPMI-1640 medium, and 10% dimethyl sulphoxide. The samples were analysed by the Fpg-modified comet assay three times over a 3-year period. All laboratories obtained a positive concentration-response relationship in cryopreserved samples (linear regression coefficients ranging from 0.79 to 0.99). However, there was a wide difference in the levels of Fpg-sensitive sites between the laboratory with the lowest (4.2% Tail DNA) and highest (74% Tail DNA) values in THP-1 cells after exposure to 4.5 mM KBrO3. In an attempt to assess sources of inter-laboratory variation in Fpg-sensitive sites, comet images from one experiment in each laboratory were forwarded to a central laboratory for visual scoring. There was high consistency between measurements of %Tail DNA values in each laboratory and the visual score of the same comets done in the central laboratory (r = 0.98, P < 0.001, linear regression). In conclusion, the results show that potassium bromate is a suitable positive comet assay control.

Integrated view and comparative analysis of baseline protein expression in mouse and rat tissues.

The increasingly large amount of proteomics data in the public domain enables, among other applications, the combined analyses of datasets to create comparative protein expression maps covering different organisms and different biological conditions. Here we have reanalysed public proteomics datasets from mouse and rat tissues (14 and 9 datasets, respectively), to assess baseline protein abundance. Overall, the aggregated dataset contained 23 individual datasets, including a total of 211 samples coming from 34 different tissues across 14 organs, comprising 9 mouse and 3 rat strains, respectively. In all cases, we studied the distribution of canonical proteins between the different organs. The number of canonical proteins per dataset ranged from 273 (tendon) and 9,715 (liver) in mouse, and from 101 (tendon) and 6,130 (kidney) in rat. Then, we studied how protein abundances compared across different datasets and organs for both species. As a key point we carried out a comparative analysis of protein expression between mouse, rat and human tissues. We observed a high level of correlation of protein expression among orthologs between all three species in brain, kidney, heart and liver samples, whereas the correlation of protein expression was generally slightly lower between organs within the same species. Protein expression results have been integrated into the resource Expression Atlas for widespread dissemination.

Assay conditions for estimating differences in base excision repair activity with Fpg-modified comet assay.

DNA repair is an essential agent in cancer development, progression, prognosis, and response to therapy. We have adapted a cellular repair assay based on the formamidopyrimidine DNA glycosylase (Fpg)-modified comet assay to assess DNA repair kinetics. The removal of oxidized nucleobases over time (0-480 min) was analyzed in peripheral blood mononuclear cells (PBMCs) and 8 cell lines. DNA damage was induced by exposure to either Ro19-8022 plus visible light or potassium bromate (KBrO3). The initial amount of damage induced by Ro 19-8022 plus light varied between cell lines, and this was apparently associated with the rate of repair. However, the amount of DNA damage induced by KBrO3 varied less between cell types, so we used this agent to study the kinetics of DNA repair. We found an early phase of ca. 60 min with fast removal of Fpg-sensitive sites, followed by slower removal over the following 7 h. In conclusion, adjusting the initial damage at T0 to an equal level can be achieved by the use of KBrO3, which allows for accurate analysis of subsequent cellular DNA repair kinetics in the first hour after exposure.

Comparison of comet-based approaches to assess base excision repair.

DNA repair plays an essential role in maintaining genomic stability, and can be assessed by various comet assay-based approaches, including the cellular repair assay and the in vitro repair assay. In the cellular repair assay, cells are challenged with a DNA-damaging compound and DNA damage removal over time is assessed. In the in vitro repair assay, an early step in the repair process is assessed as the ability of a cellular extract to recognize and incise damaged DNA in substrate nucleoids from cells treated with a DNA-damaging compound. Our direct comparison of both assays in eight cell lines and human peripheral blood lymphocytes indicated no significant relationship between these DNA repair assays (R2 = 0.084, P = 0.52). The DNA incision activity of test cells measured with the in vitro repair assay correlated with the background level of DNA damage in the untreated test cells (R2 = 0.621, P = 0.012). When extracts were prepared from cells exposed to DNA-damaging agents (10 mM KBrO3 or 1 µM Ro 19-8022 plus light), the incision activity was significantly increased, which is in line with the notion that base excision repair is inducible. The data presented suggest that the two assays do not measure the same endpoint of DNA repair and should be considered as complementary.

Abdominal lymphocele following multi-level anterior lumbar interbody fusion (ALIF) managed with a laparoscopic peritoneal window: case report and review of the literature.

PURPOSE: Lymphocele formation following anterior lumbar interbody fusion (ALIF) is not common, but it can pose diagnostic and treatment challenges. The purpose of this case is to report for the first time the treatment of a postoperative lymphocele following a multi-level ALIF using a peritoneal window made through a minimally invasive laparoscopic approach. METHODS: Case report. RESULTS: A 74-year-old male with a history of prostatectomy and pelvic radiation underwent a staged L3-S1 ALIF (left paramedian approach) and T10-pelvis posterior instrumented with L1-5 decompression/posterior column osteotomies for degenerative scoliosis and neurogenic claudication. Three weeks after surgery, swelling of the left abdomen and entire left leg was reported. Computed tomography of the abdomen/pelvis demonstrated a large (19.2 × 12.0 × 15.4 cm) retroperitoneal fluid collection with compression of the left ureter and left common iliac vein. Fluid analysis (80% lymphocytes) was consistent with a lymphocele. Percutaneous drainage for 4 days was ineffective at clearing the lymphocele. For more definitive management, the patient underwent an uncomplicated laparoscopic creation of a peritoneal window to allow passive drainage of lymphatic fluid into the abdomen. Three years after surgery, he had no back or leg pain, had achieved spinal union, and had no abdominal swelling or left leg swelling. Advanced imaging also confirmed resolution of the lymphocele. CONCLUSIONS: In this case report, creation of a peritoneal window minimally invasively via a laparoscope allowing passive drainage of lymphatic fluid into the abdomen was safe and effective for management of an abdominal lymphocele following a multi-level ALIF.

Factors associated with success of ream-and-run arthroplasty at a minimum of 5 years.

BACKGROUND: Ream-and-run arthroplasty offers improvements in shoulder pain and function for patients with primary glenohumeral arthritis who wish to avoid limitations associated with a polyethylene glenoid component. Longer-term clinical outcome assessments of the ream-and-run procedure are sparse in the literature. This study aimed to present minimum 5-year functional outcomes in a large cohort following ream-and-run arthroplasty to determine factors associated with clinical success and reoperation. METHODS: Patients who underwent ream-and-run surgery with a minimum of 5 years of follow-up were identified through a retrospective review of a prospectively maintained database from a single academic institution. To assess clinical outcomes, the Simple Shoulder Test (SST) was administered and assessed to determine achievement of the minimal clinically important difference, as well as the necessity for open revision surgery. Factors associated with P < .1 on univariate analysis were included in multivariate analysis. RESULTS: Of 228 patients who consented to undergo long-term follow-up, 201 (88%) were included in our analysis (mean follow-up period, 7.6 ± 2.1 years). The average age was 59.4 ± 8.8 years, and 93% of patients were men. The majority of patients had a diagnosis of osteoarthritis (79%) or capsulorrhaphy arthropathy (10%). SST scores improved from a mean of 4.9 ± 2.5 preoperatively to a mean of 10.2 ± 2.6 at latest follow-up. The minimal clinically important difference in the SST score of 2.6 was reached by 165 patients (82%). Male sex (P = .020), nondiabetic status (P = .080), and lower preoperative SST score (P < .001) were included in a multivariate analysis. Male sex (P = .010) and lower preoperative SST score (P < .001) were associated with clinically important improvements in SST scores on multivariate analysis. Open revision surgery was required in 22 patients (11%). Younger age (P < .001), female sex (P = .055), and higher preoperative pain score (P = .023) were included in a multivariate analysis. Only younger age was predictive of open revision surgery (P = .003). CONCLUSIONS: Ream-and-run arthroplasty can provide significant and clinically important improvements in clinical outcomes at minimum 5-year follow-up. Successful clinical outcomes were significantly associated with male sex and lower preoperative SST scores. Reoperation was more common in younger patients.

Immunoglobulin genes, reproductive isolation and vertebrate speciation.

Reproductive isolation drives the formation of new species, and many genes contribute to this through Dobzhansky-Muller incompatibilities (DMIs). These incompatibilities occur when gene divergence affects loci encoding interacting products such as receptors and their ligands. We suggest here that the nature of vertebrate immunoglobulin (IG) genes must make them prone to DMIs. The genes of these complex loci form functional genes through the process of recombination, giving rise to a repertoire of heterodimeric receptors of incredible diversity. This repertoire, within individuals and within species, must defend against pathogens but must also avoid pathogenic self-reactivity. We suggest that this avoidance of autoimmunity is only achieved through a coordination of evolution between heavy- and light-chain genes, and between these genes and the rest of the genome. Without coordinated evolution, the hybrid offspring of two diverging populations will carry a heavy burden of DMIs, resulting in a loss of fitness. Critical incompatibilities could manifest as incompatibilities between a mother and her divergent offspring. During fetal development, biochemical differences between the parents of hybrid offspring could make their offspring a target of the maternal immune system. This hypothesis was conceived in the light of recent insights into the population genetics of IG genes. This has suggested that antibody genes are probably as susceptible to evolutionary forces as other parts of the genome. Further repertoire studies in human and nonhuman species should now help determine whether antibody genes have been part of the evolutionary forces that drive the development of species.

OGRDB: a reference database of inferred immune receptor genes.

High-throughput sequencing of the adaptive immune receptor repertoire (AIRR-seq) is providing unprecedented insights into the immune response to disease and into the development of immune disorders. The accurate interpretation of AIRR-seq data depends on the existence of comprehensive germline gene reference sets. Current sets are known to be incomplete and unrepresentative of the degree of polymorphism and diversity in human and animal populations. A key issue is the complexity of the genomic regions in which they lie, which, because of the presence of multiple repeats, insertions and deletions, have not proved tractable with short-read whole genome sequencing. Recently, tools and methods for inferring such gene sequences from AIRR-seq datasets have become available, and a community approach has been developed for the expert review and publication of such inferences. Here, we present OGRDB, the Open Germline Receptor Database (https://ogrdb.airr-community.org), a public resource for the submission, review and publication of previously unknown receptor germline sequences together with supporting evidence.

VDJbase: an adaptive immune receptor genotype and haplotype database.

VDJbase is a publicly available database that offers easy searching of data describing the complete sets of gene sequences (genotypes and haplotypes) inferred from adaptive immune receptor repertoire sequencing datasets. VDJbase is designed to act as a resource that will allow the scientific community to explore the genetic variability of the immunoglobulin (Ig) and T cell receptor (TR) gene loci. It can also assist in the investigation of Ig- and TR-related genetic predispositions to diseases. Our database includes web-based query and online tools to assist in visualization and analysis of the genotype and haplotype data. It enables users to detect those alleles and genes that are significantly over-represented in a particular population, in terms of genotype, haplotype and gene expression. The database website can be freely accessed at https://www.vdjbase.org/, and no login is required. The data and code use creative common licenses and are freely downloadable from https://bitbucket.org/account/user/yaarilab/projects/GPHP.

Symbiotic organs shaped by distinct modes of genome evolution in cephalopods.

Microbes have been critical drivers of evolutionary innovation in animals. To understand the processes that influence the origin of specialized symbiotic organs, we report the sequencing and analysis of the genome of Euprymna scolopes, a model cephalopod with richly characterized host-microbe interactions. We identified large-scale genomic reorganization shared between E. scolopes and Octopus bimaculoides and posit that this reorganization has contributed to the evolution of cephalopod complexity. To reveal genomic signatures of host-symbiont interactions, we focused on two specialized organs of E. scolopes: the light organ, which harbors a monoculture of Vibrio fischeri, and the accessory nidamental gland (ANG), a reproductive organ containing a bacterial consortium. Our findings suggest that the two symbiotic organs within E. scolopes originated by different evolutionary mechanisms. Transcripts expressed in these microbe-associated tissues displayed their own unique signatures in both coding sequences and the surrounding regulatory regions. Compared with other tissues, the light organ showed an abundance of genes associated with immunity and mediating light, whereas the ANG was enriched in orphan genes known only from E. scolopes Together, these analyses provide evidence for different patterns of genomic evolution of symbiotic organs within a single host.