The whale shark genome reveals how genomic and physiological properties scale with body size.

The endangered whale shark (Rhincodon typus) is the largest fish on Earth and a long-lived member of the ancient Elasmobranchii clade. To characterize the relationship between genome features and biological traits, we sequenced and assembled the genome of the whale shark and compared its genomic and physiological features to those of 83 animals and yeast. We examined the scaling relationships between body size, temperature, metabolic rates, and genomic features and found both general correlations across the animal kingdom and features specific to the whale shark genome. Among animals, increased lifespan is positively correlated to body size and metabolic rate. Several genomic traits also significantly correlated with body size, including intron and gene length. Our large-scale comparative genomic analysis uncovered general features of metazoan genome architecture: Guanine and cytosine (GC) content and codon adaptation index are negatively correlated, and neural connectivity genes are longer than average genes in most genomes. Focusing on the whale shark genome, we identified multiple features that significantly correlate with lifespan. Among these were very long gene length, due to introns being highly enriched in repetitive elements such as CR1-like long interspersed nuclear elements, and considerably longer neural genes of several types, including connectivity, activity, and neurodegeneration genes. The whale shark genome also has the second slowest evolutionary rate observed in vertebrates to date. Our comparative genomics approach uncovered multiple genetic features associated with body size, metabolic rate, and lifespan and showed that the whale shark is a promising model for studies of neural architecture and lifespan.

Genome-wide prediction of synthetic rescue mediators of resistance to targeted and immunotherapy.

Most patients with advanced cancer eventually acquire resistance to targeted therapies, spurring extensive efforts to identify molecular events mediating therapy resistance. Many of these events involve synthetic rescue (SR) interactions, where the reduction in cancer cell viability caused by targeted gene inactivation is rescued by an adaptive alteration of another gene (the rescuer). Here, we perform a genome-wide in silico prediction of SR rescuer genes by analyzing tumor transcriptomics and survival data of 10,000 TCGA cancer patients. Predicted SR interactions are validated in new experimental screens. We show that SR interactions can successfully predict cancer patients' response and emerging resistance. Inhibiting predicted rescuer genes sensitizes resistant cancer cells to therapies synergistically, providing initial leads for developing combinatorial approaches to overcome resistance proactively. Finally, we show that the SR analysis of melanoma patients successfully identifies known mediators of resistance to immunotherapy and predicts novel rescuers.

Staufen1-mediated mRNA decay functions in adipogenesis.

The double-stranded RNA binding protein Staufen1 (Stau1) is involved in diverse gene expression pathways. For Stau1-mediated mRNA decay (SMD) in mammals, Stau1 binds to the 3' untranslated region of target mRNA and recruits Upf1 to elicit rapid mRNA degradation. However, the events downstream of Upf1 recruitment and the biological importance of SMD remain unclear. Here we show that SMD involves PNRC2, decapping activity, and 5'-to-3' exonucleolytic activity. In particular, Upf1 serves as an adaptor protein for the association of PNRC2 and Stau1. During adipogenesis, Stau1 and PNRC2 increase in abundance, Upf1 becomes hyperphosphorylated, and consequently SMD efficiency is enhanced. Intriguingly, downregulation of SMD components attenuates adipogenesis in a way that is rescued by downregulation of an antiadipogenic factor, Krüppel-like factor 2 (KLF2), the mRNA of which is identified as a substrate of SMD. Our data thus identify a biological role for SMD in adipogenesis.

The genome of the giant Nomura's jellyfish sheds light on the early evolution of active predation.

BACKGROUND: Unique among cnidarians, jellyfish have remarkable morphological and biochemical innovations that allow them to actively hunt in the water column and were some of the first animals to become free-swimming. The class Scyphozoa, or true jellyfish, are characterized by a predominant medusa life-stage consisting of a bell and venomous tentacles used for hunting and defense, as well as using pulsed jet propulsion for mobility. Here, we present the genome of the giant Nomura's jellyfish (Nemopilema nomurai) to understand the genetic basis of these key innovations. RESULTS: We sequenced the genome and transcriptomes of the bell and tentacles of the giant Nomura's jellyfish as well as transcriptomes across tissues and developmental stages of the Sanderia malayensis jellyfish. Analyses of the Nemopilema and other cnidarian genomes revealed adaptations associated with swimming, marked by codon bias in muscle contraction and expansion of neurotransmitter genes, along with expanded Myosin type II family and venom domains, possibly contributing to jellyfish mobility and active predation. We also identified gene family expansions of Wnt and posterior Hox genes and discovered the important role of retinoic acid signaling in this ancient lineage of metazoans, which together may be related to the unique jellyfish body plan (medusa formation). CONCLUSIONS: Taken together, the Nemopilema jellyfish genome and transcriptomes genetically confirm their unique morphological and physiological traits, which may have contributed to the success of jellyfish as early multi-cellular predators.

SMG5-PNRC2 is functionally dominant compared with SMG5-SMG7 in mammalian nonsense-mediated mRNA decay.

In mammals, nonsense-mediated mRNA decay (NMD) functions in post-transcriptional gene regulation as well as mRNA surveillance. A key NMD factor, Upf1, becomes hyperphosphorylated by SMG1 kinase during the recognition of NMD substrates. Hyperphosphorylated Upf1 interacts with several factors including SMG5, SMG6, SMG7 and PNRC2 to trigger rapid mRNA degradation. However, the possible cross-talk among these factors and their selective use during NMD remain unknown. Here, we show that PNRC2 is preferentially complexed with SMG5, but not with SMG6 or SMG7, and that downregulation of PNRC2 abolishes the interaction between SMG5 and Dcp1a, a component of the decapping complex. In addition, tethering experiments reveal the function of Upf1, SMG5 and PNRC2 at the same step of NMD and the requirement of SMG6 for Upf1 for efficient mRNA degradation. Intriguingly, microarray results reveal the significant overlap of SMG5-dependent NMD substrates more with PNRC2-dependent NMD substrates than with SMG7-dependent NMD substrates, suggesting the functional dominance of SMG5-PNRC2, rather than SMG5-SMG7, under normal conditions. The results provide evidence that, to some extent, endogenous NMD substrates have their own binding preference for Upf1-interacting adaptors or effectors.

Conservation in first introns is positively associated with the number of exons within genes and the presence of regulatory epigenetic signals.

BACKGROUND: Genomes of higher eukaryotes have surprisingly long first introns and in some cases, the first introns have been shown to have higher conservation relative to other introns. However, the functional relevance of conserved regions in the first introns is poorly understood. Leveraging the recent ENCODE data, here we assess potential regulatory roles of conserved regions in the first intron of human genes. RESULTS: We first show that relative to other downstream introns, the first introns are enriched for blocks of highly conserved sequences. We also found that the first introns are enriched for several chromatin marks indicative of active regulatory regions and this enrichment of regulatory marks is correlated with enrichment of conserved blocks in the first intron; the enrichments of conservation and regulatory marks in first intron are not entirely explained by a general, albeit variable, bias for certain marks toward the 5' end of introns. Interestingly, conservation as well as proportions of active regulatory chromatin marks in the first intron of a gene correlates positively with the numbers of exons in the gene but the correlation is significantly weakened in second introns and negligible beyond the second intron. The first intron conservation is also positively correlated with the gene's expression level in several human tissues. Finally, a gene-wise analysis shows significant enrichments of active chromatin marks in conserved regions of first introns, relative to the conserved regions in other introns of the same gene. CONCLUSIONS: Taken together, our analyses strongly suggest that first introns are enriched for active transcriptional regulatory signals under purifying selection.

Interleukin-1 promotes coagulation, which is necessary for protective immunity in the lung against Streptococcus pneumoniae infection.

Interleukin (IL)-1 is a well-known cytokine for the initiation of innate immunity in bacterial infection. However, the underlying mechanism of IL-1 on the respiratory infection is not fully elucidated. We studied how IL-1 contributes to the host defense against Streptococcus pneumoniae. IL-1R(-/-) mice showed high mortality, local cytokine storm, and substantial infiltrates in the lower respiratory tract after intratracheal challenge with S. pneumoniae. The IL-1-deficient condition did not suppress the propagation of bacteria in the lung, although the recruitment and the bacteria-killing ability of neutrophils (CD11b(+)Ly6C(+)Ly6G(+)) were not defective compared with wild-type mice. Unexpectedly, we found that the transcription of fibrinogen alpha and gamma genes were highly activated in the lungs of wild-type mice after the infection, whereas no significant changes were found in IL-1R(-/-) mice. Of note, synthesis of fibrinogen was dependent on the IL-1-IL-6-Stat3 cascade. Treatment with recombinant fibrinogen improved survival and bacterial propagation in the IL-1R(-/-) mice and blockade of the coagulation increased the susceptibility of wild-type mice to pneumococcal pneumonia. Our findings suggest that IL-1 signaling leads to the synthesis of fibrinogen in the lung after pneumococcus infection and is followed by coagulation, which contributes to the control of bacterial infection in the pulmonary tract.

Expression breadth and expression abundance behave differently in correlations with evolutionary rates.

BACKGROUND: One of the main objectives of the molecular evolution and evolutionary systems biology field is to reveal the underlying principles that dictate protein evolutionary rates. Several studies argue that expression abundance is the most critical component in determining the rate of evolution, especially in unicellular organisms. However, the expression breadth also needs to be considered for multicellular organisms. RESULTS: In the present paper, we analyzed the relationship between the two expression variables and rates using two different genome-scale expression datasets, microarrays and ESTs. A significant positive correlation between the expression abundance (EA) and expression breadth (EB) was revealed by Kendall's rank correlation tests. A novel random shuffling approach was applied for EA and EB to compare the correlation coefficients obtained from real data sets to those estimated based on random chance. A novel method called a Fixed Group Analysis (FGA) was designed and applied to investigate the correlations between expression variables and rates when one of the two expression variables was evenly fixed. CONCLUSIONS: In conclusion, all of these analyses and tests consistently showed that the breadth rather than the abundance of gene expression is tightly linked with the evolutionary rate in multicellular organisms.

Korean Genome Project: 1094 Korean personal genomes with clinical information.

We present the initial phase of the Korean Genome Project (Korea1K), including 1094 whole genomes (sequenced at an average depth of 31×), along with data of 79 quantitative clinical traits. We identified 39 million single-nucleotide variants and indels of which half were singleton or doubleton and detected Korean-specific patterns based on several types of genomic variations. A genome-wide association study illustrated the power of whole-genome sequences for analyzing clinical traits, identifying nine more significant candidate alleles than previously reported from the same linkage disequilibrium blocks. Also, Korea1K, as a reference, showed better imputation accuracy for Koreans than the 1KGP panel. As proof of utility, germline variants in cancer samples could be filtered out more effectively when the Korea1K variome was used as a panel of normals compared to non-Korean variome sets. Overall, this study shows that Korea1K can be a useful genotypic and phenotypic resource for clinical and ethnogenetic studies.

Analysis of human disease genes in the context of gene essentiality.

The characteristics of human disease genes were investigated through a comparative analysis with mouse mutant phenotype data. Mouse orthologs with mutations that resulted in discernible phenotypes were separated from mutations with no phenotypic defect, listing 'phenotype' and 'no phenotype' genes. First, we showed that phenotype genes are more likely to be disease genes compared to no phenotype genes. Phenotype genes were further divided into 'embryonic lethal', 'postnatal lethal', and 'non-lethal phenotype' groups. Interestingly, embryonic lethal genes, the most essential genes in mouse, were less likely to be disease genes than postnatal lethal genes. These findings indicate that some extremely essential genes are less likely to be disease genes, although human disease genes tend to display characteristics of essential genes. We also showed that, in lethal groups, non-disease genes tend to evolve slower than disease genes indicating a strong purifying selection on non-disease genes in this group. In addition, phenotype and no phenotype groups showed differing types of disease mutations. Disease genes in the no phenotype group displayed a higher frequency of regulatory mutations while those in the phenotype group had more frequent coding mutations, indicating that the types of disease mutations vary depending on gene essentiality. Furthermore, missense disease mutations in no phenotype genes were found to be more radical amino acid substitutions than those in phenotype genes.

Comparison of exercise electrocardiography and magnetocardiography for detection of coronary artery disease using ST-segment fluctuation score.

BACKGROUND: Exercise electrocardiography (ECG) is frequently used as a diagnostic measure in patients with suspected coronary artery disease (CAD). However, it has low sensitivity for the detection of CAD. Magnetocardiography (MCG) has been proposed as an alternative tool to accurately diagnose CAD. OBJECTIVE: To date, a direct comparison of MCG to ECG has not been performed. This study sought to compare them for predicting the presence of significantly obstructive CAD. METHODS: The patients with chest pain or other symptoms suggestive of CAD were enrolled in the analysis. All the patients underwent a clinical evaluation, exercise ECG, MCG test, and coronary angiography (CA). CAD was defined as stenosis ≥70% in at least one major coronary artery on quantitative analysis of CA. RESULTS: We prospectively enrolled 202 consecutive patients who suggested CAD. The prevalence of CAD on CA was 39.1%. Sensitivity and accuracy for CAD diagnosis was higher for MCG compared with exercise ECG (sensitivities 68.4% and 40.5%, p <0.001, specificities 95.1% and 91.1%, p = 0.267, and accuracies 84.7% and 71.3%, p <0.001, respectively). There was no incremental diagnostic value of combined MCG and ECG to detect coronary artery disease (p = 0.357). CONCLUSIONS: For the patients with intermediate to high risk of CAD, MCG exercise test provides better diagnostic accuracy for the detection of relevant obstruction of the epicardial coronaries when directly compared to exercise ECG.

The Galleria mellonella Hologenome Supports Microbiota-Independent Metabolism of Long-Chain Hydrocarbon Beeswax.

The greater wax moth, Galleria mellonella, degrades wax and plastic molecules. Despite much interest, the genetic basis of these hallmark traits remains poorly understood. Herein, we assembled high-quality genome and transcriptome data from G. mellonella to investigate long-chain hydrocarbon wax metabolism strategies. Specific carboxylesterase and lipase and fatty-acid-metabolism-related enzymes in the G. mellonella genome are transcriptionally regulated during feeding on beeswax. Strikingly, G. mellonella lacking intestinal microbiota successfully decomposes long-chain fatty acids following wax metabolism, although the intestinal microbiome performs a supplementary role in short-chain fatty acid degradation. Notably, final wax derivatives were detected by gas chromatography even in the absence of gut microbiota. Our findings provide insight into wax moth adaptation and may assist in the development of unique wax-degradation strategies with a similar metabolic approach for a plastic molecule polyethylene biodegradation using organisms without intestinal microbiota.

Brca2 abrogation engages with the alternative lengthening of telomeres via break-induced replication.

A subset of cancer cells maintains their telomeres without telomerase through the recombination-based alternative lengthening of telomeres (ALT) pathway. Currently, it is not yet clear in what context ALT is induced and how the pathway choice is made. Here, we show that abrogation of Brca2 reinforces break-induced replication (BIR) and engages with ALT pathway. Brca2 depletion in telomerase-null mouse cells alleviated the growth defect, accompanied by telomere elongation, suggesting the induction of ALT. We also found that Brca2-depleted telomerase-null cells exhibited dynamic clustering of telomeres from G2 phase in Promyelocytic Nuclear (PML) bodies. For Brca2-deficient ALT induction, Rad51 was dispensable but Mre11 and Rad52 were required. Congruently, conservative telomeric DNA synthesis was apparent in mitosis, indicating that the absence of Brca2 directed towards Rad52-mediated BIR. Collectively, we propose that Brca2 abrogation can instigate ALT tumourigenesis through the induction of BIR. This study implies that inhibitors of BIR may be useful for BRCA2-associated ALT-type cancers. Assessing ALT features may be considered for the tailored therapy of BRCA2-associated cancers.

Depression and suicide risk prediction models using blood-derived multi-omics data.

More than 300 million people worldwide experience depression; annually, ~800,000 people die by suicide. Unfortunately, conventional interview-based diagnosis is insufficient to accurately predict a psychiatric status. We developed machine learning models to predict depression and suicide risk using blood methylome and transcriptome data from 56 suicide attempters (SAs), 39 patients with major depressive disorder (MDD), and 87 healthy controls. Our random forest classifiers showed accuracies of 92.6% in distinguishing SAs from MDD patients, 87.3% in distinguishing MDD patients from controls, and 86.7% in distinguishing SAs from controls. We also developed regression models for predicting psychiatric scales with R2 values of 0.961 and 0.943 for Hamilton Rating Scale for Depression-17 and Scale for Suicide Ideation, respectively. Multi-omics data were used to construct psychiatric status prediction models for improved mental health treatment.

The Draft Genome of an Octocoral, Dendronephthya gigantea.

Coral reefs composed of stony corals are threatened by global marine environmental changes. However, soft coral communities of octocorallian species, appear more resilient. The genomes of several cnidarians species have been published, including from stony corals, sea anemones, and hydra. To fill the phylogenetic gap for octocoral species of cnidarians, we sequenced the octocoral, Dendronephthya gigantea, a nonsymbiotic soft coral, commonly known as the carnation coral. The D. gigantea genome size is ∼276 Mb. A high-quality genome assembly was constructed from PacBio long reads (29.85 Gb with 108× coverage) and Illumina short paired-end reads (35.54 Gb with 128× coverage) resulting in the highest N50 value (1.4 Mb) reported thus far among cnidarian genomes. About 12% of the genome is repetitive elements and contained 28,879 predicted protein-coding genes. This gene set is composed of 94% complete BUSCO ortholog benchmark genes, which is the second highest value among the cnidarians, indicating high quality. Based on molecular phylogenetic analysis, octocoral and hexacoral divergence times were estimated at 544 MYA. There is a clear difference in Hox gene composition between these species: unlike hexacorals, the Antp superclass Evx gene was absent in D. gigantea. Here, we present the first genome assembly of a nonsymbiotic octocoral, D. gigantea to aid in the comparative genomic analysis of cnidarians, including stony and soft corals, both symbiotic and nonsymbiotic. The D. gigantea genome may also provide clues to mechanisms of differential coping between the soft and stony corals in response to scenarios of global warming.

Raptor genomes reveal evolutionary signatures of predatory and nocturnal lifestyles.

BACKGROUND: Birds of prey (raptors) are dominant apex predators in terrestrial communities, with hawks (Accipitriformes) and falcons (Falconiformes) hunting by day and owls (Strigiformes) hunting by night. RESULTS: Here, we report new genomes and transcriptomes for 20 species of birds, including 16 species of birds of prey, and high-quality reference genomes for the Eurasian eagle-owl (Bubo bubo), oriental scops owl (Otus sunia), eastern buzzard (Buteo japonicus), and common kestrel (Falco tinnunculus). Our extensive genomic analysis and comparisons with non-raptor genomes identify common molecular signatures that underpin anatomical structure and sensory, muscle, circulatory, and respiratory systems related to a predatory lifestyle. Compared with diurnal birds, owls exhibit striking adaptations to the nocturnal environment, including functional trade-offs in the sensory systems, such as loss of color vision genes and selection for enhancement of nocturnal vision and other sensory systems that are convergent with other nocturnal avian orders. Additionally, we find that a suite of genes associated with vision and circadian rhythm are differentially expressed in blood tissue between nocturnal and diurnal raptors, possibly indicating adaptive expression change during the transition to nocturnality. CONCLUSIONS: Overall, raptor genomes show genomic signatures associated with the origin and maintenance of several specialized physiological and morphological features essential to be apex predators.

Magnetocardiography scoring system to predict the presence of obstructive coronary artery disease.

BACKGROUND: Magnetocardiography (MCG) has been proposed as a non-invasive and functional technique with high accuracy for diagnosis of myocardial ischemia. OBJECTIVE: This study sought to develop a novel scoring system of MCG for predicting the presence of significant obstructive coronary artery disease (CAD). METHODS: In a training set of 108 subjects, predictors of ≥70% stenosis in at least one major coronary vessel were prospectively identified from MCG variables. The final model was then retrospectively validated in a separate set of 45 subjects. RESULTS: In the multivariable logistic regression, among those in the training set, elevated scores were predictive of ≥70% stenosis in all subjects (OR: 40.85; 95% CI: 6.28-265.90; p < 0.001). In the validation set, the score had an area under the receiver-operating characteristic curve of 0.91 (p < 0.001) for ≥70% stenosis. At an optimal cutoff, the score had 89% sensitivity, 77% specificity, 74% positive predictive value (PPV), 91% negative predictive value (NPV), and 82% accuracy for ≥70% stenosis. Partitioning the score into three levels of predicted risk, 91% of subjects could be identified or excluding CAD (81% PPV and 84% NPV). CONCLUSION: We described an MCG score with high accuracy for predicting the presence of anatomically significant CAD.

Putative functional genes in idiopathic dilated cardiomyopathy.

Idiopathic dilated cardiomyopathy (DCM) is a complex disorder with a genetic and an environmental component involving multiple genes, many of which are yet to be discovered. We integrate genetic, epigenetic, transcriptomic, phenotypic, and evolutionary features into a method - Hridaya, to infer putative functional genes underlying DCM in a genome-wide fashion, using 213 human heart genomes and transcriptomes. Many genes identified by Hridaya are experimentally shown to cause cardiac complications. We validate the top predicted genes, via five different genome-wide analyses: First, the predicted genes are associated with cardiovascular functions. Second, their knockdowns in mice induce cardiac abnormalities. Third, their inhibition by drugs cause cardiac side effects in human. Fourth, they tend to have differential exon usage between DCM and normal samples. Fifth, analyzing 213 individual genotypes, we show that regulatory polymorphisms of the predicted genes are associated with elevated risk of cardiomyopathy. The stratification of DCM patients based on cardiac expression of the functional genes reveals two subgroups differing in key cardiac phenotypes. Integrating predicted functional genes with cardiomyocyte drug treatment experiments reveals novel potential drug targets. We provide a list of investigational drugs that target the newly identified functional genes that may lead to cardiac side effects.

Harnessing synthetic lethality to predict the response to cancer treatment.

While synthetic lethality (SL) holds promise in developing effective cancer therapies, SL candidates found via experimental screens often have limited translational value. Here we present a data-driven approach, ISLE (identification of clinically relevant synthetic lethality), that mines TCGA cohort to identify the most likely clinically relevant SL interactions (cSLi) from a given candidate set of lab-screened SLi. We first validate ISLE via a benchmark of large-scale drug response screens and by predicting drug efficacy in mouse xenograft models. We then experimentally test a select set of predicted cSLi via new screening experiments, validating their predicted context-specific sensitivity in hypoxic vs normoxic conditions and demonstrating cSLi's utility in predicting synergistic drug combinations. We show that cSLi can successfully predict patients' drug treatment response and provide patient stratification signatures. ISLE thus complements existing actionable mutation-based methods for precision cancer therapy, offering an opportunity to expand its scope to the whole genome.

Osteogenesis for postoperative temporal bone defects using human ear adipose-derived stromal cells and tissue engineering: An animal model study.

Mastoidectomy, the removal of infected mastoid bones, is a common surgical procedure for the treatment of chronic otitis media. Persistent and recurrent otorrhea and accumulation of keratin debris following open cavity mastoidectomy are still bothersome issues for both patients and otologists. In this study, we used human ear adipose-derived stromal cells (hEASCs) in combination with polycaprolactone (PCL) scaffolds and osteogenic differentiation medium (ODM) to regenerate temporal bone defects. The hEASCs showed stem cell phenotypes, and these characteristics were maintained up to passage 5. Mastoid bulla and cranial bone defects were induced in Sprague-Dawley rats using AgNO3 and burr hole drilling, respectively, and the rats were then divided into five groups: (1) control, (2) hEASCs, (3) hEASCs + ODM, (4) hEASCs + PCL scaffolds, and (5) hEASCs + PCL scaffolds + ODM. Osteogenesis was evaluated by micro-computed tomography and histology. Compared with the control group, the groups transplanted with hEASCs and PCL scaffolds had significantly higher bone formation along the periphery of the mastoid bulla area. Moreover, ODM synergistically enhanced bone formation in mastoid bulla defects. Our results suggest that combining hEASCs with PCL scaffolds represents a promising method for anatomical and functional reconstruction of postoperative temporal bone defects following mastoidectomy. © 2017 Wiley Periodicals, Inc. J Biomed Mater Res Part A: 105A: 3493-3501, 2017.