Whole exome sequencing in the rat.

BACKGROUND: The rat genome was sequenced in 2004 with the aim to improve human health altered by disease and environmental influences through gene discovery and animal model validation. Here, we report development and testing of a probe set for whole exome sequencing (WES) to detect sequence variants in exons and UTRs of the rat genome. Using an in-silico approach, we designed probes targeting the rat exome and compared captured mutations in cancer-related genes from four chemically induced rat tumor cell lines (C6, FAT7, DSL-6A/C1, NBTII) to validated cancer genes in the human database, Catalogue of Somatic Mutations in Cancer (COSMIC) as well as normal rat DNA. Paired, fresh frozen (FF) and formalin-fixed, paraffin-embedded (FFPE) liver tissue from naive rats were sequenced to confirm known dbSNP variants and identify any additional variants. RESULTS: Informatics analysis of available gene annotation from rat RGSC6.0/rn6 RefSeq and Ensembl transcripts provided 223,636 unique exons representing a total of 26,365 unique genes and untranslated regions. Using this annotation and the Rn6 reference genome, an in-silico probe design generated 826,878 probe sequences of which 94.2% were uniquely aligned to the rat genome without mismatches. Further informatics analysis revealed 25,249 genes (95.8%) covered by at least one probe and 23,603 genes (93.5%) had every exon covered by one or more probes. We report high performance metrics from exome sequencing of our probe set and Sanger validation of annotated, highly relevant, cancer gene mutations as cataloged in the human COSMIC database, in addition to several exonic variants in cancer-related genes. CONCLUSIONS: An in-silico probe set was designed to enrich the rat exome from isolated DNA. The platform was tested on rat tumor cell lines and normal FF and FFPE liver tissue. The method effectively captured target exome regions in the test DNA samples with exceptional sensitivity and specificity to obtain reliable sequencing data representing variants that are likely chemically induced somatic mutations. Genomic discovery conducted by means of high throughput WES queries should benefit investigators in discovering rat genomic variants in disease etiology and in furthering human translational research.

The RareCyte® platform for next-generation analysis of circulating tumor cells.

Circulating tumor cells (CTCs) can reliably be identified in cancer patients and are associated with clinical outcome. Next-generation "liquid biopsy" technologies will expand CTC diagnostic investigation to include phenotypic characterization and single-cell molecular analysis. We describe here a rare cell analysis platform designed to comprehensively collect and identify CTCs, enable multi-parameter assessment of individual CTCs, and retrieve single cells for molecular analysis. The platform has the following four integrated components: 1) density-based separation of the CTC-containing blood fraction and sample deposition onto microscope slides; 2) automated multiparameter fluorescence staining; 3) image scanning, analysis, and review; and 4) mechanical CTC retrieval. The open platform utilizes six fluorescence channels, of which four channels are used to identify CTC and two channels are available for investigational biomarkers; a prototype assay that allows three investigational biomarker channels has been developed. Single-cell retrieval from fixed slides is compatible with whole genome amplification methods for genomic analysis. © 2018 The Authors. Cytometry Part A published by Wiley Periodicals, Inc. on behalf of International Society for Advancement of Cytometry.

Let-7 family of microRNA is required for maturation and adult-like metabolism in stem cell-derived cardiomyocytes.

In metazoans, transition from fetal to adult heart is accompanied by a switch in energy metabolism-glycolysis to fatty acid oxidation. The molecular factors regulating this metabolic switch remain largely unexplored. We first demonstrate that the molecular signatures in 1-year (y) matured human embryonic stem cell-derived cardiomyocytes (hESC-CMs) are similar to those seen in in vivo-derived mature cardiac tissues, thus making them an excellent model to study human cardiac maturation. We further show that let-7 is the most highly up-regulated microRNA (miRNA) family during in vitro human cardiac maturation. Gain- and loss-of-function analyses of let-7g in hESC-CMs demonstrate it is both required and sufficient for maturation, but not for early differentiation of CMs. Overexpression of let-7 family members in hESC-CMs enhances cell size, sarcomere length, force of contraction, and respiratory capacity. Interestingly, large-scale expression data, target analysis, and metabolic flux assays suggest this let-7-driven CM maturation could be a result of down-regulation of the phosphoinositide 3 kinase (PI3K)/AKT protein kinase/insulin pathway and an up-regulation of fatty acid metabolism. These results indicate let-7 is an important mediator in augmenting metabolic energetics in maturing CMs. Promoting maturation of hESC-CMs with let-7 overexpression will be highly significant for basic and applied research.

A Distributed Network for Intensive Longitudinal Monitoring in Metastatic Triple-Negative Breast Cancer.

Accelerating cancer research is expected to require new types of clinical trials. This report describes the Intensive Trial of OMics in Cancer (ITOMIC) and a participant with triple-negative breast cancer metastatic to bone, who had markedly elevated circulating tumor cells (CTCs) that were monitored 48 times over 9 months. A total of 32 researchers from 14 institutions were engaged in the patient's evaluation; 20 researchers had no prior involvement in patient care and 18 were recruited specifically for this patient. Whole-exome sequencing of 3 bone marrow samples demonstrated a novel ROS1 variant that was estimated to be present in most or all tumor cells. After an initial response to cisplatin, a hypothesis of crizotinib sensitivity was disproven. Leukapheresis followed by partial CTC enrichment allowed for the development of a differential high-throughput drug screen and demonstrated sensitivity to investigational BH3-mimetic inhibitors of BCL-2 that could not be tested in the patient because requests to the pharmaceutical sponsors were denied. The number and size of CTC clusters correlated with clinical status and eventually death. Focusing the expertise of a distributed network of investigators on an intensively monitored patient with cancer can generate high-resolution views of the natural history of cancer and suggest new opportunities for therapy. Optimization requires access to investigational drugs.

Analysis of the DNA sequence and duplication history of human chromosome 15.

Here we present a finished sequence of human chromosome 15, together with a high-quality gene catalogue. As chromosome 15 is one of seven human chromosomes with a high rate of segmental duplication, we have carried out a detailed analysis of the duplication structure of the chromosome. Segmental duplications in chromosome 15 are largely clustered in two regions, on proximal and distal 15q; the proximal region is notable because recombination among the segmental duplications can result in deletions causing Prader-Willi and Angelman syndromes. Sequence analysis shows that the proximal and distal regions of 15q share extensive ancient similarity. Using a simple approach, we have been able to reconstruct many of the events by which the current duplication structure arose. We find that most of the intrachromosomal duplications seem to share a common ancestry. Finally, we demonstrate that some remaining gaps in the genome sequence are probably due to structural polymorphisms between haplotypes; this may explain a significant fraction of the gaps remaining in the human genome.

Amino acid primed mTOR activity is essential for heart regeneration.

Heart disease is the leading cause of death with no method to repair damaged myocardium due to the limited proliferative capacity of adult cardiomyocytes. Curiously, mouse neonates and zebrafish can regenerate their hearts via cardiomyocyte de-differentiation and proliferation. However, a molecular mechanism of why these cardiomyocytes can re-enter cell cycle is poorly understood. Here, we identify a unique metabolic state that primes adult zebrafish and neonatal mouse ventricular cardiomyocytes to proliferate. Zebrafish and neonatal mouse hearts display elevated glutamine levels, predisposing them to amino-acid-driven activation of TOR, and that TOR activation is required for zebrafish cardiomyocyte regeneration in vivo. Through a multi-omics approach with cellular validation we identify metabolic and mitochondrial changes during the first week of regeneration. These data suggest that regeneration of zebrafish myocardium is driven by metabolic remodeling and reveals a unique metabolic regulator, TOR-primed state, in which zebrafish and mammalian cardiomyocytes are regeneration competent.

The effect of smoking on MAOA promoter methylation in DNA prepared from lymphoblasts and whole blood.

Prior work using lymphoblast DNA prepared from 192 subjects from the Iowa Adoption Studies (IAS) demonstrated that decreased MAOA promoter methylation was associated with lifetime symptom count for nicotine dependence (ND) and provided suggestive evidence that the amount of methylation is genotype dependent. In the current investigation, we replicate and extend these prior findings in three ways using another 289 IAS subjects and the same methodologies. First, we show that methylation is dependent on current smoking status. Second, we introduce a factor analytic approach to DNA methylation, highlighting three distinct regions of the promoter that may function in somewhat different ways for males and females. Third, we directly compare the methylation signatures in DNA prepared from whole blood and lymphoblasts from a subset of these subjects and provide suggestive evidence favoring the use of lymphoblast DNA. We conclude that smoking reliably decreases MAOA methylation, but exact characterization of effects on level of methylation depend on genotype, smoking history, current smoking status, gender, and region of the promoter-associated CpG Island examined.

Author Correction: TFPa/HADHA is required for fatty acid beta-oxidation and cardiolipin re-modeling in human cardiomyocytes.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Metabolism as an early predictor of DPSCs aging.

Tissue resident adult stem cells are known to participate in tissue regeneration and repair that follows cell turnover, or injury. It has been well established that aging impedes the regeneration capabilities at the cellular level, but it is not clear if the different onset of stem cell aging between individuals can be predicted or prevented at an earlier stage. Here we studied the dental pulp stem cells (DPSCs), a population of adult stem cells that is known to participate in the repair of an injured tooth, and its properties can be affected by aging. The dental pulp from third molars of a diverse patient group were surgically extracted, generating cells that had a high percentage of mesenchymal stem cell markers CD29, CD44, CD146 and Stro1 and had the ability to differentiate into osteo/odontogenic and adipogenic lineages. Through RNA seq and qPCR analysis we identified homeobox protein, Barx1, as a marker for DPSCs. Furthermore, using high throughput transcriptomic and proteomic analysis we identified markers for DPSC populations with accelerated replicative senescence. In particular, we show that the transforming growth factor-beta (TGF-ß) pathway and the cytoskeletal proteins are upregulated in rapid aging DPSCs, indicating a loss of stem cell characteristics and spontaneous initiation of terminal differentiation. Importantly, using metabolic flux analysis, we identified a metabolic signature for the rapid aging DPSCs, prior to manifestation of senescence phenotypes. This metabolic signature therefore can be used to predict the onset of replicative senescence. Hence, the present study identifies Barx1 as a DPSCs marker and dissects the first predictive metabolic signature for DPSCs aging.

Discovery of high affinity, pan-allelic, and pan-mammalian reactive antibodies against the myeloid checkpoint receptor SIRPα.

Targeting the CD47-signal-regulatory protein α (SIRPα) pathway represents a novel therapeutic approach to enhance anti-cancer immunity by promoting both innate and adaptive immune responses. Unlike CD47, which is expressed ubiquitously, SIRPα expression is mainly restricted to myeloid cells and neurons. Therefore, compared to CD47-targeted therapies, targeting SIRPα may result in differential safety and efficacy profiles, potentially enabling lower effective doses and improved pharmacokinetics and pharmacodynamics. The development of effective SIRPα antagonists is restricted by polymorphisms within the CD47-binding domain of SIRPα, necessitating pan-allele reactive anti-SIRPα antibodies for therapeutic intervention in diverse patient populations. We immunized wild-type and human antibody transgenic chickens with a multi-allele and multi-species SIRPα regimen in order to discover pan-allelic and pan-mammalian reactive anti-SIRPα antibodies suitable for clinical translation. A total of 200 antibodies were isolated and screened for SIRPα reactivity from which approximately 70 antibodies with diverse SIRPα binding profiles, sequence families, and epitopes were selected for further characterization. A subset of anti-SIRPα antibodies bound to both human SIRPα v1 and v2 alleles with high affinity ranging from low nanomolar to picomolar, potently antagonized the CD47/SIRPα interaction, and potentiated macrophage-mediated antibody-dependent cellular phagocytosis in vitro. X-ray crystal structures of five anti-SIRPα antigen-binding fragments, each with unique epitopes, in complex with SIRPα (PDB codes 6NMV, 6NMU, 6NMT, 6NMS, and 6NMR) are reported. Furthermore, some of the anti-SIRPα antibodies cross-react with cynomolgus SIRPα and various mouse SIRPα alleles (BALB/c, NOD, BL/6), which can facilitate preclinical to clinical development. These properties provide an attractive rationale to advance the development of these anti-SIRPα antibodies as a novel therapy for advanced malignancies. Abbreviations: ADCC: antibody-dependent cellular cytotoxicity; ADCP: antibody-dependent cellular phagocytosis; CFSE: carboxyfluorescein succinimidyl ester; Fab: fragment antigen binding; Fc: fragment crystallizable; FcγR: Fcγ receptor; Ig: immunoglobulin; IND: investigational new drug; MDM⊘: monocyte-derived macrophage; NOD: non-obese diabetic; scFv: single chain fragment variable; SCID: severe combined immunodeficiency; SIRP: signal-regulatory protein.

TFPa/HADHA is required for fatty acid beta-oxidation and cardiolipin re-modeling in human cardiomyocytes.

Mitochondrial trifunctional protein deficiency, due to mutations in hydratase subunit A (HADHA), results in sudden infant death syndrome with no cure. To reveal the disease etiology, we generated stem cell-derived cardiomyocytes from HADHA-deficient hiPSCs and accelerated their maturation via an engineered microRNA maturation cocktail that upregulated the epigenetic regulator, HOPX.  Here we report, matured HADHA mutant cardiomyocytes treated with an endogenous mixture of fatty acids manifest the disease phenotype: defective calcium dynamics and repolarization kinetics which results in a pro-arrhythmic state. Single cell RNA-seq reveals a cardiomyocyte developmental intermediate, based on metabolic gene expression. This intermediate gives rise to mature-like cardiomyocytes in control cells but, mutant cells transition to a pathological state with reduced fatty acid beta-oxidation, reduced mitochondrial proton gradient, disrupted cristae structure and defective cardiolipin remodeling. This study reveals that HADHA (tri-functional protein alpha), a monolysocardiolipin acyltransferase-like enzyme, is required for fatty acid beta-oxidation and cardiolipin remodeling, essential for functional mitochondria in human cardiomyocytes.

Gene expression analysis of human otosclerotic stapedial footplates.

Otosclerosis is a complex disease that results in a common form of conductive hearing loss due to impaired mobility of the stapes. Stapedial motion becomes compromised secondary to invasion of otosclerotic foci into the stapedio-vestibular joint. Although environmental factors and genetic causes have been implicated in this process, the pathogenesis of otosclerosis remains poorly understood. To identify molecular contributors to otosclerosis we completed a microarray study of otosclerotic stapedial footplates. Stapes footplate samples from otosclerosis and control patients were used in the analysis. One-hundred-and-ten genes were found to be differentially expressed in otosclerosis samples. Ontological analysis of differentially expressed genes in otosclerosis provides evidence for the involvement of a number of pathways in the disease process that include interleukin signaling, inflammation and signal transduction, suggesting that aberrant regulation of these pathways leads to abnormal bone remodeling. Functional analyses of genes from this study will enhance our understanding of the pathogenesis of this disease.

Genome-wide analysis of epigenetic silencing identifies BEX1 and BEX2 as candidate tumor suppressor genes in malignant glioma.

Promoter hypermethylation and histone deacetylation are common epigenetic mechanisms implicated in the transcriptional silencing of tumor suppressor genes in human cancer. We treated two immortalized glioma cell lines, T98 and U87, and 10 patient-derived primary glioma cell lines with trichostatin A (TSA), a histone deacetylase inhibitor, or 5-aza-2'-deoxycytidine (5-AzaC), a DNA methyltransferase inhibitor, to comprehensively identify the cohort of genes reactivated through the pharmacologic reversal of these distinct but related epigenetic processes. Whole-genome microarray analysis identified genes induced by TSA (653) or 5-AzaC treatment (170). We selected a subset of reactivated genes that were markedly induced (greater than two-fold) after treatment with either TSA or 5-AzaC in a majority of glioma cell lines but not in cultured normal astrocytes. We then characterized the degree of promoter methylation and transcriptional silencing of selected genes in histologically confirmed human tumor and nontumor brain specimens. We identified two novel brain expressed genes, BEX1 and BEX2, which were silenced in all tumor specimens and exhibited extensive promoter hypermethylation. Viral-mediated reexpression of either BEX1 or BEX2 led to increased sensitivity to chemotherapy-induced apoptosis and potent tumor suppressor effects in vitro and in a xenograft mouse model. Using an integrated approach, we have established a novel platform for the genome-wide screening of epigenetically silenced genes in malignant glioma. This experimental paradigm provides a powerful new method for the identification of epigenetically silenced genes with potential function as tumor suppressors, biomarkers for disease diagnosis and detection, and therapeutically reversible modulators of critical regulatory pathways important in glioma pathogenesis.

The relationship of 5HTT (SLC6A4) methylation and genotype on mRNA expression and liability to major depression and alcohol dependence in subjects from the Iowa Adoption Studies.

Serotonin Transporter (5HTT or SLC6A4) mRNA transcription is regulated by both genetic and epigenetic mechanisms. Unfortunately, despite intense scrutiny, the exact identity and contribution of each of these regulatory mechanisms, and their relationship to behavioral illness remain unknown. This lack of knowledge is critical because alterations in SLC6A4 function are posited to be central to a wide variety of CNS disorders. In order to address this shortcoming, we quantified 5HTTLPR genotype, SLC6A4 mRNA production and CpG methylation using biomaterial from 192 lymphoblast cell lines derived from subjects who participated in the latest wave of the Iowa Adoption Studies. We then analyzed the resulting data with respect to clinical characteristics. We confirmed prior findings that the short (s) 5HTTLPR allele is associated with lower amounts of mRNA transcription, but there was no significant effect of the "Long G" allele on mRNA transcription. We also found that CpG methylation was higher (P < 0.0008) and mRNA production (P < 0.0001) was lower in females as compared to males. Those subjects with a lifetime history of Alcohol Dependence had higher levels of SLC6A4 mRNA. There was a trend for an association of increased overall methylation with lifetime history of major depression. Finally, we confirm our prior findings that the exact levels of 5HTT mRNA expression are dependent on how it is measured. We conclude that both genetic variation and epigenetic modifications contribute to the regulation of SLC6A4 function and that more in-depth studies of the molecular mechanisms controlling gene activity and the relationship of these mechanisms to behavioral illness are indicated.

MAOA methylation is associated with nicotine and alcohol dependence in women.

In recent years, the role of epigenetic phenomenon, such as methylation, in mediating vulnerability to behavioral illness has become increasingly appreciated. One prominent locus at which epigenetic phenomena are thought to be in play is the monoamine oxidase A (MAOA) locus. In order to examine the role of methylation at this locus, we performed quantitative methylation analysis across the promoter region of this gene in lymphoblast lines derived from 191 subjects participating in the Iowa Adoption Studies (IAS). We analyzed the resulting data with respect to genotype and lifetime symptom counts for the more common major behavioral disorders in the IAS, antisocial personality disorder (ASPD), and substance use disorders (alcohol (AD) and nicotine dependence (ND)). We found that methylation status was significantly associated with lifetime symptom counts for ND (P < 0.001) and AD (P < 0.008) in women, but not men. Furthermore, a trend was found for women homozygous for the 3,3 allele to have a higher degree of overall methylation than women homozygous for the 4,4 allele (P < 0.10). We conclude that methylation of MAOA may play a significant role in common psychiatric illness and that further examination of epigenetic processes at this locus is in order.

Dickkopf-1 is an epigenetically silenced candidate tumor suppressor gene in medulloblastoma.

Medulloblastoma is a heterogeneous pediatric brain tumor with significant therapy-related morbidity, its five-year survival rates ranging from 30% to 70%. Improvement in diagnosis and therapy requires better understanding of medulloblastoma pathology. We used whole-genome microarray analysis to identify putative tumor suppressor genes silenced by epigenetic mechanisms in medulloblastoma. This analysis yielded 714 up-regulated genes in immortalized medulloblastoma cell line D283 on treatment with histone deacetylase (HDAC) inhibitor trichostatin A (TSA). Dickkopf-1 (DKK1), a Wnt antagonist, was found to be up-regulated on HDAC inhibition. We examined DKK1 expression in primary medulloblastoma cells and patient samples by reverse transcriptase PCR and found it to be significantly down-regulated relative to normal cerebellum. Transfection of a DKK1 gene construct into D283 cell lines suppressed medulloblastoma tumor growth in colony focus assays by 60% (P < 0.001). In addition, adenoviral vector-mediated expression of DKK1 in medulloblastoma cells increased apoptosis fourfold (P < 0.001). These data reveal that inappropriate histone modifications might deregulate DKK1 expression in medulloblastoma tumorigenesis and block its tumor-suppressive activity.

Adaptation of a RAS pathway activation signature from FF to FFPE tissues in colorectal cancer.

BACKGROUND: The KRAS gene is mutated in about 40 % of colorectal cancer (CRC) cases, which has been clinically validated as a predictive mutational marker of intrinsic resistance to anti-EGFR inhibitor (EGFRi) therapy. Since nearly 60 % of patients with a wild type KRAS fail to respond to EGFRi combination therapies, there is a need to develop more reliable molecular signatures to better predict response. Here we address the challenge of adapting a gene expression signature predictive of RAS pathway activation, created using fresh frozen (FF) tissues, for use with more widely available formalin fixed paraffin-embedded (FFPE) tissues. METHODS: In this study, we evaluated the translation of an 18-gene RAS pathway signature score from FF to FFPE in 54 CRC cases, using a head-to-head comparison of five technology platforms. FFPE-based technologies included the Affymetrix GeneChip (Affy), NanoString nCounter™ (NanoS), Illumina whole genome RNASeq (RNA-Acc), Illumina targeted RNASeq (t-RNA), and Illumina stranded Total RNA-rRNA-depletion (rRNA). RESULTS: Using Affy_FF as the "gold" standard, initial analysis of the 18-gene RAS scores on all 54 samples shows varying pairwise Spearman correlations, with (1) Affy_FFPE (r = 0.233, p = 0.090); (2) NanoS_FFPE (r = 0.608, p < 0.0001); (3) RNA-Acc_FFPE (r = 0.175, p = 0.21); (4) t-RNA_FFPE (r = -0.237, p = 0.085); (5) and t-RNA (r = -0.012, p = 0.93). These results suggest that only NanoString has successful FF to FFPE translation. The subsequent removal of identified "problematic" samples (n = 15) and genes (n = 2) further improves the correlations of Affy_FF with three of the five technologies: Affy_FFPE (r = 0.672, p < 0.0001); NanoS_FFPE (r = 0.738, p < 0.0001); and RNA-Acc_FFPE (r = 0.483, p = 0.002). CONCLUSIONS: Of the five technology platforms tested, NanoString technology provides a more faithful translation of the RAS pathway gene expression signature from FF to FFPE than the Affymetrix GeneChip and multiple RNASeq technologies. Moreover, NanoString was the most forgiving technology in the analysis of samples with presumably poor RNA quality. Using this approach, the RAS signature score may now be reasonably applied to FFPE clinical samples.

A multigene mutation classification of 468 colorectal cancers reveals a prognostic role for APC.

Colorectal cancer (CRC) is a highly heterogeneous disease, for which prognosis has been relegated to clinicopathologic staging for decades. There is a need to stratify subpopulations of CRC on a molecular basis to better predict outcome and assign therapies. Here we report targeted exome-sequencing of 1,321 cancer-related genes on 468 tumour specimens, which identified a subset of 17 genes that best classify CRC, with APC playing a central role in predicting overall survival. APC may assume 0, 1 or 2 truncating mutations, each with a striking differential impact on survival. Tumours lacking any APC mutation carry a worse prognosis than single APC mutation tumours; however, two APC mutation tumours with mutant KRAS and TP53 confer the poorest survival among all the subgroups examined. Our study demonstrates a prognostic role for APC and suggests that sequencing of APC may have clinical utility in the routine staging and potential therapeutic assignment for CRC.

New perspectives on rickettsial evolution from new genome sequences of rickettsia, particularly R. canadensis, and Orientia tsutsugamushi.

The complete genome sequences available for eight species of Rickettsia and information for other near relatives in the Rickettsiales including Orientia and species of Anaplasmataceae are a rich resource for comparative analyses of the evolution of these obligate intracellular bacteria. Differences in these organisms have permitted them to colonize varied intracellular compartments, arthropod vectors, and vertebrate reservoirs in both pathogenic and symbiotic relationships. We summarize some comparative aspects of the genomes of these organisms, paying particular attention to the recently completed sequence for R. canadensis McKiel strain and an estimated two-thirds of the genome sequence for a Thailand patient isolate of Orientia tsutsugamushi. The Rickettsia genomes exhibit a high degree of synteny punctuated by distinctive chromosome inversions and consistent phylogenetic relationships regardless of whether protein coding sequences or RNA genes, concatenated open reading frames or gene regions, or whole genomes are used to construct phylogenetic trees. The aggregate characteristics (number, length, composition, repeat identity) of tandem repeat sequences of Rickettsia, which often exhibit recent and rapid divergence between closely related strains and species of bacteria, are also very conserved in Rickettsia but differed significantly in Orientia. O. tsutsugamushi shared no significant synteny to species of Rickettsia or Anaplasmataceae, supporting its placement in a unique genus. Like Rickettsia felis, Orientia has many transposases and ankyrin and tetratricopeptide repeat domains. Orientia shares the important ATP/ADP translocase and proline-betaine transporter multigene families with Rickettsia, but has more gene families that may be involved in regulatory and transporter responses to environmental stimuli.