Scalable Parallel Algorithm for Graph Neural Network Interatomic Potentials in Molecular Dynamics Simulations.

Message-passing graph neural network interatomic potentials (GNN-IPs), particularly those with equivariant representations such as NequIP, are attracting significant attention due to their data efficiency and high accuracy. However, parallelizing GNN-IPs poses challenges because multiple message-passing layers complicate data communication within the spatial decomposition method, which is preferred by many molecular dynamics (MD) packages. In this article, we propose an efficient parallelization scheme compatible with GNN-IPs and develop a package, SevenNet (Scalable EquiVariance-Enabled Neural NETwork), based on the NequIP architecture. For MD simulations, SevenNet interfaces with the LAMMPS package. Through benchmark tests on a 32-GPU cluster with examples of SiO2, SevenNet achieves over 80% parallel efficiency in weak-scaling scenarios and exhibits nearly ideal strong-scaling performance as long as GPUs are fully utilized. However, the strong-scaling performance significantly declines with suboptimal GPU utilization, particularly affecting parallel efficiency in cases involving lightweight models or simulations with small numbers of atoms. We also pretrain SevenNet with a vast data set from the Materials Project (dubbed "SevenNet-0") and assess its performance on generating amorphous Si3N4 containing more than 100,000 atoms. By developing scalable GNN-IPs, this work aims to bridge the gap between advanced machine-learning models and large-scale MD simulations, offering researchers a powerful tool to explore complex material systems with high accuracy and efficiency.

Stability and Equilibrium Structures of Unknown Ternary Metal Oxides Explored by Machine-Learned Potentials.

Ternary metal oxides are crucial components in a wide range of applications and have been extensively cataloged in experimental materials databases. However, there still exist cation combinations with unknown stability and structures of their compounds in oxide forms. In this study, we employ extensive crystal structure prediction methods, accelerated by machine-learned potentials, to investigate these untapped chemical spaces. We examine 181 ternary metal oxide systems, encompassing most cations except for partially filled 3d or f shells, and determine their lowest-energy crystal structures with representative stoichiometry derived from prevalent oxidation states or recommender systems. Consequently, we discover 45 ternary oxide systems containing stable compounds against decomposition into binary or elemental phases, the majority of which incorporate noble metals. Comparisons with other theoretical databases highlight the strengths and limitations of informatics-based material searches. With a relatively modest computational resource requirement, we contend that heuristic-based structure searches, as demonstrated in this study, offer a promising approach for future materials discovery endeavors.

Introduction of the Korea BioData Station (K-BDS) for sharing biological data.

A wave of new technologies has created opportunities for the cost-effective generation of high-throughput profiles of biological systems, foreshadowing a "data-driven science" era. The large variety of data available from biological research is also a rich resource that can be used for innovative endeavors. However, we are facing considerable challenges in big data deposition, integration, and translation due to the complexity of biological data and its production at unprecedented exponential rates. To address these problems, in 2020, the Korean government officially announced a national strategy to collect and manage the biological data produced through national R&D fund allocations and provide the collected data to researchers. To this end, the Korea Bioinformation Center (KOBIC) developed a new biological data repository, the Korea BioData Station (K-BDS), for sharing data from individual researchers and research programs to create a data-driven biological study environment. The K-BDS is dedicated to providing free open access to a suite of featured data resources in support of worldwide activities in both academia and industry.

Transcriptome-metabolome-wide association study (TMWAS) in rats revealed a potential carcinogenic effect of DEHP in thyroid associated with eicosanoids.

The incidence of thyroid cancer (TC) has increased considerably in the last few decades. Environmental factors, including plasticizers, are recognized as potential risks leading to thyroid cancer in humans. In this study, we used a transcriptome-metabolome-wide association study to find the unidentified carcinogenic mechanism of di-2-ethylhexyl phthalate (DEHP) in thyroid and biomarkers for non-invasive diagnosis. Rats were treated with different doses of DEHP (0, 0.3, 3, 30, 150 mg DEHP/kg bw/day) for 13 weeks. Then, the thyroids were processed for Ki67 staining and RNA-seq. Also, 17-h urine samples were collected for high-resolution metabolomics analysis. After a high dose of DEHP exposure, the terminal body weights and the thyroid and parathyroid glands weights were not altered. However, the liver weights and numbers of Ki67-positive cells were increased. Further, multivariate statistical analysis revealed that metabolic shifts were considerably altered above 30 mg DEHP/kg bw/day. In RNA-seq analysis, some cancer-related genes were altered, including 18 upregulated and 9 downregulated transcripts. These cancer transcripts and whole metabolome data were integrated to uncover thyroid cancer-related metabolic pathways, which revealed that cancer-related transcripts had a network structure linked to eicosanoids such as leukotriene D4 and prostaglandin. In brief, our study demonstrated that DEHP can induce thyroid hyperplasia through the eicosanoid-associated pathway, providing further insight into the mechanism of DEHP-associated thyroid cancer.

IL-10 Protects Mice From the Lung Infection of Acinetobacter baumannii and Contributes to Bacterial Clearance by Regulating STAT3-Mediated MARCO Expression in Macrophages.

Interleukin-10 plays important, yet contrasting, roles in host protection against bacterial infections and in the septic response. To determine the role of IL-10 in the host defense against Acinetobacter baumannii infection, wild-type (WT) and IL-10-deficient mice were infected intranasally with the bacteria. IL-10-deficient mice exhibited increased mortality, severe pathology, and excess production of proinflammatory cytokines and chemokines in the lungs, and increased bacterial burdens in bronchoalveolar lavage (BAL) fluids and lung homogenates after A. baumannii infection, compared to WT mice. Intranasal administration of recombinant IL-10 rescued mice from the lethality of the bacterial infection by promoting bacterial clearance and reducing production of cytokines and chemokines in the lungs. In vitro experiments revealed that IL-10 enhanced phagocytosis and bacterial killing by macrophages by upregulating the macrophage receptor with collagenous structure (MARCO). In addition, A. baumannii-induced activation of STAT3 was impaired in IL-10-deficient macrophages, which was essential for expression of MARCO. Intranasal adoptive transfer of WT macrophages resulted in significant increases in mice survival and bacterial clearance in IL-10-deficient mice infected with A. baumannii. Our results show that IL-10 played an important role in the host defense against pulmonary infection of A. baumannii by promoting the antibacterial function of macrophages by regulating MARCO expression through the STAT3-mediated pathway.

Novel genes exhibiting DNA methylation alterations in Korean patients with chronic lymphocytic leukaemia: a methyl-CpG-binding domain sequencing study.

Chronic lymphocytic leukaemia (CLL) exhibits differences between Asians and Caucasians in terms of incidence rate, age at onset, immunophenotype, and genetic profile. We performed genome-wide methylation profiling of CLL in an Asian cohort for the first time. Eight Korean patients without somatic immunoglobulin heavy chain gene hypermutations underwent methyl-CpG-binding domain sequencing (MBD-seq), as did five control subjects. Gene Ontology, pathway analysis, and network-based prioritization of differentially methylated genes were also performed. More regions were hypomethylated (2,062 windows) than were hypermethylated (777 windows). Promoters contained the highest proportion of differentially methylated regions (0.08%), while distal intergenic and intron regions contained the largest number of differentially methylated regions. Protein-coding genes were the most abundant, followed by long noncoding and short noncoding genes. The most significantly over-represented signalling pathways in the differentially methylated gene list included immune/cancer-related pathways and B-cell receptor signalling. Among the top 10 hub genes identified via network-based prioritization, four (UBC, GRB2, CREBBP, and GAB2) had no known relevance to CLL, while the other six (STAT3, PTPN6, SYK, STAT5B, XPO1, and ABL1) have previously been linked to CLL in Caucasians. As such, our analysis identified four novel candidate genes of potential significance to Asian patients with CLL.

Characterization of Distinct CyanoHABs-Related Modules in Microbial Recurrent Association Network.

To elucidate the interspecies connectivity between cyanobacteria and other bacteria (non-cyanobacteria) during cyanobacterial harmful algal blooms (cyanoHABs), samples were collected from the Nakdong River, Korea, from June 2016 to August 2017, and microbial recurrent association network (MRAN) analysis was performed to overcome the limitations of conventional network analysis. Microcystis blooms were tightly linked with Pseudanabaena in summer and were accompanied by significant changes in the non-cyanobacterial community composition (nCCC) compared to non-bloom period. Riverine bacterial communities could be clearly separated into modules that were involved in the formation, maintenance, and decomposition of cyanoHABs. Roseomonas and Herbaspirillum were directly linked with major cyanobacteria and assigned to connector and module hub in cyanoHABs-related modules, respectively. The functional profiles of the cyanoHABs-related modules suggested that nitrate reduction, aerobic ammonia oxidation, fermentation, and hydrocarbon degradation could be increased during the Microcystis bloom periods. In conclusion, MRAN analysis revealed that specific bacteria belonging to cyanoHABs-related module, including connectors and module hubs, appeared to contribute to the development and collapse of cyanoHABs. Therefore, to understand cyanoHABs, a modular microbial perspective may be more helpful than a single bacterial species perspective.

The water depth-dependent co-occurrence patterns of marine bacteria in shallow and dynamic Southern Coast, Korea.

To investigate the interactions between bacterial species in relation to the biotic and abiotic environmental fluctuations, free-living (FL), nanoparticle-associated (NP), and microparticle-associated (MP) bacterial community compositions (BCCs) were analyzed. A total of 267 samples were collected from July to December 2016 in the dynamic and shallow southern coastal water of Korea. The variations in BCC mostly depended on planktonic size fraction. Network analysis revealed water depth-dependent co-occurrence patterns of coastal bacterial communities. Higher interspecies connectivity was observed within FL bacteria than NP/MP bacteria, suggesting that FL bacteria with a streamlined genome may need other bacterial metabolites for survival, while the NP/MP copiotrophs may have the self-supporting capacity to produce the vital nutrients. The analysis of topological roles of individual OTUs in the network revealed that several groups of metabolically versatile bacteria (the marine Roseobacters, Flavobacteriales, Desulfobacterales, and SAR406 clade) acted as module hubs in different water depth. In conclusion, interspecies interactions dominated in FL bacteria, compared to NP and MP bacteria; modular structures of bacterial communities and keystone species strongly depended on the water depth-derived environmental factors. Furthermore, the multifunctional, versatile FL bacteria could play pivotal roles in dynamic shallow coastal ecosystems.

Low Dose Exposure to Di-2-Ethylhexylphthalate in Juvenile Rats Alters the Expression of Genes Related with Thyroid Hormone Regulation.

Phthalates widely used in the manufacture of plastics have deeply penetrated into our everyday lives. Recently, a concern over the toxicity of phthalates on thyroid, has been raised but in most of cases, the doses employed were unrealistically high. To investigate the effects of phthalates on thyroid, we investigated the effects of the repeated oral exposure to low to high doses (0.3, 3, 30 and 150 mg/kg) di-2-ethylhexylphthalate (DEHP) from weaning to maturity for 90 days in juvenile rats on the thyroid. The histological examination revealed that DEHP significantly induced hyperplasia in the thyroid from the doses of 30 mg/kg, which was confirmed with Ki67 staining. In line with this finding, increased mRNA expression of thyrotropin releasing hormone (Trh) was observed in the thyroid of female at 0.3 mg/kg and 150 mg/kg as determined by RNAseq analysis. Moreover, significantly increased expression of parathyroid hormone (Pth) in the female at 0.3 mg/kg, and thyroglobulin (Tg) and thyroid hormone responsive (Thrsp) in the male at 0.3 mg/kg were noted in the blood, of which changes were substantially attenuated at 150 m/kg, alluding the meaningful effects of low dose DEHP on the thyroid hormone regulation. Urinary excretion of mono-2-ethylhexyl-phthalate (MEHP), a major metabolite of DEHP was determined to be 4.10 and 12.26 ppb in male, 6.65 and 324 ppb in female at 0.3 and 30 mg/kg DEHP, respectively, which fell within reported human urine levels. Collectively, these results suggest a potential adverse effects of low dose phthalates on the thyroid.

On-Chip Fluorescence Switching System for Constructing a Rewritable Random Access Data Storage Device.

We report the development of on-chip fluorescence switching system based on DNA strand displacement and DNA hybridization for the construction of a rewritable and randomly accessible data storage device. In this study, the feasibility and potential effectiveness of our proposed system was evaluated with a series of wet experiments involving 40 bits (5 bytes) of data encoding a 5-charactered text (KRIBB). Also, a flexible data rewriting function was achieved by converting fluorescence signals between "ON" and "OFF" through DNA strand displacement and hybridization events. In addition, the proposed system was successfully validated on a microfluidic chip which could further facilitate the encoding and decoding process of data. To the best of our knowledge, this is the first report on the use of DNA hybridization and DNA strand displacement in the field of data storage devices. Taken together, our results demonstrated that DNA-based fluorescence switching could be applicable to construct a rewritable and randomly accessible data storage device through controllable DNA manipulations.

Long-Term Stability and Integrity of Plasmid-Based DNA Data Storage.

Validation of long-term DNA stability and integrity are essential for the use of DNA in data storage applications. Because of this, we evaluated the plasmid-based DNA data storage in a manner that preserves DNA stability and integrity. A document consisting of 2046 words was encoded with DNA sequences using Perl script, and the encoded DNA sequences were synthesized for information storage. The DNA comprised a total of 22 chemically synthesized DNA fragments with 400 nucleotides each, which were incorporated into a plasmid vector. A long-term DNA stability study demonstrated that 3-year stored plasmid containing text information showed DNA stability at controlled conditions of -20 °C. The plasmid DNA under accelerated aging conditions (AAC) up to 65 °C for 20 days, which corresponds to approximately 20 years of storage at -20 °C, also exhibited no significant differences in DNA stability compared to newly produced plasmid. Also, the 3-year old plasmid stored at -20 °C and the AAC-tested plasmid stored up to 65 °C for 20 days had functional integrity and nucleotide integrity comparable to control sample, thereby allowing for retrieval of the original error-free text data. Finally, the nucleotides were sequenced, and then decoded to retrieve the original data, thereby allowing us to read the text with 100% accuracy, and amplify the DNA with a simple and quick bacterial transformation. To the best of our knowledge, this is the first report on examining the long-term stability and integrity of plasmid-based DNA data storage. Taken together, our results indicate that plasmid DNA data storage can be useful for long-term archival storage to recover the source text in a reproducible and accountable manner.

Idiopathic hypereosinophilia is clonal disorder? Clonality identified by targeted sequencing.

Idiopathic hypereosinophilia (IHE)/idiopathic hypereosinophilic syndrome (IHES) has been defined by a persistent elevation of the blood eosinophil count exceeding 1.5×103/µL, without evidence of reactive or clonal causes. While T-cell clonality assessment has been recommended for unexplained hypereosinophilia, this approach is not often applied to routine practice in the clinic. We hypothesized that the clonality would exist in a subset of IHE/IHES patients. We aimed to investigate the candidate mutations and T-cell clonality in IHE/IHES and to explore the role of mutations in eosinophil proliferation. We performed targeted capture sequencing for 88 genes using next-generation sequencing, T-cell receptor (TCR) gene rearrangement assays, and pathway network analysis in relation to eosinophil proliferation. By targeted sequencing, 140 variants in 59 genes were identified. Sixteen out of 30 patients (53.3%) harbored at least one candidate mutation. The most frequently affected genes were NOTCH1 (26.7%), SCRIB and STAG2 (16.7%), and SH2B3 (13.3%). Network analysis revealed that our 21 candidate genes (BIRC3, BRD4, CSF3R, DNMT3A, EGR2, EZH2, FAT4, FLT3, GATA2, IKZF, JAK2, MAPK1, MPL, NF1, NOTCH1, PTEN, RB1, RUNX1, TET2, TP53 and WT1) are functionally linked to the eosinophilopoietic pathway. Among the 21 candidate genes, five genes (MAPK1, RUNX1, GATA2, NOTCH1 and TP53) with the highest number of linkages were considered major genes. A TCR assay revealed that four patients (13.3%) had a clonal TCR rearrangement. NOTCH1 was the most frequently mutated gene and was shown to be a common node for eosinophilopoiesis in our network analysis, while the possibility of hidden T cell malignancy was indwelling in the presence of NOTCH1 mutation, though not revealed by aberrant T cell study. Collectively, these results provide new evidence that mutations affecting eosinophilopoiesis underlie a subset of IHE/IHES, and the candidate genes are inferred to act their potential roles in the eosinophilopoietic pathway.

Dysregulated expression of proteins associated with ER stress, autophagy and apoptosis in tissues from nonalcoholic fatty liver disease.

Nonalcoholic fatty liver disease (NAFLD) is categorized into nonalcoholic fatty liver (NAFL) and nonalcoholic steatohepatitis (NASH) and has emerged as a risk factor for more critical clinical conditions. However, the underlying mechanisms of NAFLD pathogenesis are not fully understood. In this study, expression of proteins associated with endoplasmic reticulum (ER) stress, apoptosis and autophagy were analyzed in normal, NAFL and NASH human livers by western blotting. Levels of some ER stress-transducing transcription factors, including cleaved activating transcription factor 6, were higher in NASH than in the normal tissues. However, the expression of a majority of the ER chaperones and foldases analyzed, including glucose-regulated protein 78 and ER protein 44, was lower in NASH than in the normal tissues. Levels of apoptosis markers, such as cleaved poly (ADP-ribose) polymerase, were also lower in NASH tissues, in which expression of some B-cell lymphoma-2 family proteins was up- or down-regulated compared to the normal tissues. The level of the autophagy substrate p62 was not different in NASH and normal tissues, although some autophagy regulators were up- or down-regulated in the NASH tissues compared to the normal tissues. Levels of most of the proteins analyzed in NAFL tissues were either similar to those in one of the other two types, NASH and normal, or were somewhere in between. Together, these findings suggest that regulation of certain important tissues processes involved in protein quality control and cell survival were broadly compromised in the NAFLD tissues.

Expression of microRNAs in Horse Plasma and Their Characteristic Nucleotide Composition.

MicroRNAs (miRNAs) in blood plasma are stable under high levels of ribonuclease activity and could function in tissue-to-tissue communication, suggesting that they may have distinctive structural characteristics compared with non-circulating miRNAs. In this study, the expression of miRNAs in horse plasma and their characteristic nucleotide composition were examined and compared with non-plasma miRNAs. Highly expressed plasma miRNA species were not part of the abundant group of miRNAs in non-plasma tissues, except for the eca-let-7 family. eca-miR-486-5p, -92a, and -21 were among the most abundant plasma miRNAs, and their human orthologs also belong to the most abundant group of miRNAs in human plasma. Uracil and guanine were the most common nucleotides of both plasma and non-plasma miRNAs. Cytosine was the least common in plasma and non-plasma miRNAs, although levels were higher in plasma miRNAs. Plasma miRNAs also showed higher expression levels of miRNAs containing adenine and cytosine repeats, compared with non-plasma miRNAs. These observations indicate that miRNAs in the plasma have a unique nucleotide composition.

Whole genome sequencing of an ethnic Pathan (Pakhtun) from the north-west of Pakistan.

BACKGROUND: Pakistan covers a key geographic area in human history, being both part of the Indus River region that acted as one of the cradles of civilization and as a link between Western Eurasia and Eastern Asia. This region is inhabited by a number of distinct ethnic groups, the largest being the Punjabi, Pathan (Pakhtuns), Sindhi, and Baloch. RESULTS: We analyzed the first ethnic male Pathan genome by sequencing it to 29.7-fold coverage using the Illumina HiSeq2000 platform. A total of 3.8 million single nucleotide variations (SNVs) and 0.5 million small indels were identified by comparing with the human reference genome. Among the SNVs, 129,441 were novel, and 10,315 nonsynonymous SNVs were found in 5,344 genes. SNVs were annotated for health consequences and high risk diseases, as well as possible influences on drug efficacy. We confirmed that the Pathan genome presented here is representative of this ethnic group by comparing it to a panel of Central Asians from the HGDP-CEPH panels typed for ~650 k SNPs. The mtDNA (H2) and Y haplogroup (L1) of this individual were also typical of his geographic region of origin. Finally, we reconstruct the demographic history by PSMC, which highlights a recent increase in effective population size compatible with admixture between European and Asian lineages expected in this geographic region. CONCLUSIONS: We present a whole-genome sequence and analyses of an ethnic Pathan from the north-west province of Pakistan. It is a useful resource to understand genetic variation and human migration across the whole Asian continent.

Increased expression of interferon signaling genes in the bone marrow microenvironment of myelodysplastic syndromes.

INTRODUCTION: The bone marrow (BM) microenvironment plays an important role in the pathogenesis of myelodysplastic syndromes (MDS) through a reciprocal interaction with resident BM hematopoietic cells. We investigated the differences between BM mesenchymal stromal cells (MSCs) in MDS and normal individuals and identified genes involved in such differences. MATERIALS AND METHODS: BM-derived MSCs from 7 MDS patients (3 RCMD, 3 RAEB-1, and 1 RAEB-2) and 7 controls were cultured. Global gene expression was analyzed using a microarray. RESULT: We found 314 differentially expressed genes (DEGs) in RCMD vs. control, 68 in RAEB vs. control, and 51 in RAEB vs. RCMD. All comparisons were clearly separated from one another by hierarchical clustering. The overall similarity between differential expression signatures from the RCMD vs. control comparison and the RAEB vs. control comparison was highly significant (p = 0), which indicates a common transcriptomic response in these two MDS subtypes. RCMD and RAEB simultaneously showed an up-regulation of interferon alpha/beta signaling and the ISG15 antiviral mechanism, and a significant fraction of the RAEB vs. control DEGs were also putative targets of transcription factors IRF and ICSBP. Pathways that involved RNA polymerases I and III and mitochondrial transcription were down-regulated in RAEB compared to RCMD. CONCLUSION: Gene expression in the MDS BM microenvironment was different from that in normal BM and exhibited altered expression according to disease progression. The present study provides genetic evidence that inflammation and immune dysregulation responses that involve the interferon signaling pathway in the BM microenvironment are associated with MDS pathogenesis, which suggests BM MSCs as a possible therapeutic target in MDS.

Short-term differential adaptation to anaerobic stress via genomic mutations by Escherichia coli strains K-12 and B lacking alcohol dehydrogenase.

Microbial adaptations often occur via genomic mutations under adverse environmental conditions. This study used Escherichia coli ΔadhE cells as a model system to investigate adaptation to anaerobic conditions, which we then compared with the adaptive mechanisms of two closely related E. coli strains, K-12 and B. In contrast to K-12 ΔadhE cells, the E. coli B ΔadhE cells exhibited significantly delayed adaptive growth under anaerobic conditions. Adaptation by the K-12 and B strains mainly employed anaerobic lactate fermentation to restore cellular growth. Several mutations were identified in the pta or pflB genes of adapted K-12 cells, but mostly in the pta gene of the B strains. However, the types of mutation in the adapted K-12 and B strains were similar. Cellular viability was affected directly by severe redox imbalance in B ΔadhE cells, which also impaired their ability to adapt to anaerobic conditions. This study demonstrates that closely related microorganisms may undergo different adaptations under the same set of adverse conditions, which might be associated with the specific metabolic characteristics of each strain. This study provides new insights into short-term microbial adaptation to stressful conditions, which may reflect dynamic microbial population changes in nature.

Variation block-based genomics method for crop plants.

BACKGROUND: In contrast with wild species, cultivated crop genomes consist of reshuffled recombination blocks, which occurred by crossing and selection processes. Accordingly, recombination block-based genomics analysis can be an effective approach for the screening of target loci for agricultural traits. RESULTS: We propose the variation block method, which is a three-step process for recombination block detection and comparison. The first step is to detect variations by comparing the short-read DNA sequences of the cultivar to the reference genome of the target crop. Next, sequence blocks with variation patterns are examined and defined. The boundaries between the variation-containing sequence blocks are regarded as recombination sites. All the assumed recombination sites in the cultivar set are used to split the genomes, and the resulting sequence regions are termed variation blocks. Finally, the genomes are compared using the variation blocks. The variation block method identified recurring recombination blocks accurately and successfully represented block-level diversities in the publicly available genomes of 31 soybean and 23 rice accessions. The practicality of this approach was demonstrated by the identification of a putative locus determining soybean hilum color. CONCLUSIONS: We suggest that the variation block method is an efficient genomics method for the recombination block-level comparison of crop genomes. We expect that this method will facilitate the development of crop genomics by bringing genomics technologies to the field of crop breeding.

The tiger genome and comparative analysis with lion and snow leopard genomes.

Tigers and their close relatives (Panthera) are some of the world's most endangered species. Here we report the de novo assembly of an Amur tiger whole-genome sequence as well as the genomic sequences of a white Bengal tiger, African lion, white African lion and snow leopard. Through comparative genetic analyses of these genomes, we find genetic signatures that may reflect molecular adaptations consistent with the big cats' hypercarnivorous diet and muscle strength. We report a snow leopard-specific genetic determinant in EGLN1 (Met39>Lys39), which is likely to be associated with adaptation to high altitude. We also detect a TYR260G>A mutation likely responsible for the white lion coat colour. Tiger and cat genomes show similar repeat composition and an appreciably conserved synteny. Genomic data from the five big cats provide an invaluable resource for resolving easily identifiable phenotypes evident in very close, but distinct, species.