Identification and characterization of structural variants related to meat quality in pigs using chromosome-level genome assemblies

BACKGROUND

Many studies have been performed to identify various genomic loci and genes associated with the meat quality in pigs. However, the full genetic architecture of the trait still remains unclear in part because of the lack of accurate identification of related structural variations (SVs) which resulted from the shortage of target breeds, the limitations of sequencing data, and the incompleteness of genome assemblies. The recent generation of a new pig breed with superior meat quality, called Nanchukmacdon, and its chromosome-level genome assembly (the NCMD assembly) has provided new opportunities.

RESULTS

By applying assembly-based SV calling approaches to various genome assemblies of pigs including Nanchukmacdon, the impact of SVs on meat quality was investigated. Especially, by checking the commonality of SVs with other pig breeds, a total of 13,819 Nanchukmacdon-specific SVs (NSVs) were identified, which have a potential effect on the unique meat quality of Nanchukmacdon. The regulatory potentials of NSVs for the expression of nearby genes were further examined using transcriptome- and epigenome-based analyses in different tissues.

CONCLUSIONS

Whole-genome comparisons based on chromosome-level genome assemblies have led to the discovery of SVs affecting meat quality in pigs, and their regulatory potentials were analyzed. The identified NSVs will provide new insights regarding genetic architectures underlying the meat quality in pigs. Finally, this study confirms the utility of chromosome-level genome assemblies and multi-omics analysis to enhance the understanding of unique phenotypes.

A chromosome-level genome assembly of the Korean crossbred pig Nanchukmacdon (Sus scrofa)

As plentiful high-quality genome assemblies have been accumulated, reference-guided genome assembly can be a good approach to reconstruct a high-quality assembly. Here, we present a chromosome-level genome assembly of the Korean crossbred pig called Nanchukmacdon (the NCMD assembly) using the reference-guided assembly approach with short and long reads. The NCMD assembly contains 20 chromosome-level scaffolds with a total size of 2.38 Gbp (N50: 138.77 Mbp). Its BUSCO score is 93.1%, which is comparable to the pig reference assembly, and a total of 20,588 protein-coding genes, 8,651 non-coding genes, and 996.14 Mbp of repetitive elements are annotated. The NCMD assembly was also used to close many gaps in the pig reference assembly. This NCMD assembly and annotation provide foundational resources for the genomic analyses of pig and related species.

Brain structure and allelic associations in Alzheimer's disease

BACKGROUND

Alzheimer's disease (AD), the most prevalent form of dementia, affects 6.5 million Americans and over 50 million people globally. Clinical, genetic, and phenotypic studies of dementia provide some insights of the observed progressive neurodegenerative processes, however, the mechanisms underlying AD onset remain enigmatic.

AIMS

This paper examines late-onset dementia-related cognitive impairment utilizing neuroimaging-genetics biomarker associations.

MATERIALS AND METHODS

The participants, ages 65-85, included 266 healthy controls (HC), 572 volunteers with mild cognitive impairment (MCI), and 188 Alzheimer's disease (AD) patients. Genotype dosage data for AD-associated single nucleotide polymorphisms (SNPs) were extracted from the imputed ADNI genetics archive using sample-major additive coding. Such 29 SNPs were selected, representing a subset of independent SNPs reported to be highly associated with AD in a recent AD meta-GWAS study by Jansen and colleagues.

RESULTS

We identified the significant correlations between the 29 genomic markers (GMs) and the 200 neuroimaging markers (NIMs). The odds ratios and relative risks for AD and MCI (relative to HC) were predicted using multinomial linear models.

DISCUSSION

In the HC and MCI cohorts, mainly cortical thickness measures were associated with GMs, whereas the AD cohort exhibited different GM-NIM relations. Network patterns within the HC and AD groups were distinct in cortical thickness, volume, and proportion of White to Gray Matter (pct), but not in the MCI cohort. Multinomial linear models of clinical diagnosis showed precisely the specific NIMs and GMs that were most impactful in discriminating between AD and HC, and between MCI and HC.

CONCLUSION

This study suggests that advanced analytics provide mechanisms for exploring the interrelations between morphometric indicators and GMs. The findings may facilitate further clinical investigations of phenotypic associations that support deep systematic understanding of AD pathogenesis.

Reference-based read clustering improves the de novo genome assembly of microbial strains

Constructing accurate microbial genome assemblies is necessary to understand genetic diversity in microbial genomes and its functional consequences. However, it still remains as a challenging task especially when only short-read sequencing technologies are used. Here, we present a new read-clustering algorithm, called RBRC, for improving de novo microbial genome assembly, by accurately estimating read proximity using multiple reference genomes. The performance of RBRC was confirmed by simulation-based evaluation in terms of assembly contiguity and the number of misassemblies, and was successfully applied to existing fungal and bacterial genomes by improving the quality of the assemblies without using additional sequencing data. RBRC is a very useful read-clustering algorithm that can be used (i) for generating high-quality genome assemblies of microbial strains when genome assemblies of related strains are available, and (ii) for upgrading existing microbial genome assemblies when the generation of additional sequencing data, such as long reads, is difficult.

Biological aging in maltreated children followed up into middle adulthood

BACKGROUND

Childhood adversity has been linked to many indicators of shorter healthy lifespan, including earlier onset of disease and disability as well as early mortality. These observations suggest the hypothesis that childhood maltreatment may accelerate aging.

OBJECTIVE

To characterize the relationship between childhood maltreatment and accelerated biological aging in a prospective cohort of 357 individuals with documented cases of childhood maltreatment and 250 controls matched on demographic and socioeconomic factors.

METHODS

Cases were drawn from juvenile and adult court records from the years 1967 through 1971 in a large Midwest metropolitan geographic area. Cases were defined as having court-substantiated cases of childhood physical or sexual abuse, or neglect occurring at age 11 or younger. Controls were selected from the same schools and hospitals of birth and matched on age, sex, race, and approximate socioeconomic status. We compared biological aging in these two groups using two blood-chemistry algorithms, the Klemera-Doubal method Biological Age (KDM BA) and the PhenoAge. Algorithms were developed and validated in data from the National Health and Nutrition Examination Surveys (NHANES) using published methods and publicly available software.

RESULTS

Participants (55% women, 49% non-White) had mean age of 41 years (SD=4). Those with court substantiated childhood maltreatment history exhibited more advanced biological aging as compared with matched controls, although this difference was statistically different for only the KDM BA measure (KDM BA Cohen's D=0.20, 95% CI=[0.03,0.36], p = 0.02; PhenoAge Cohen's D=0.09 95% CI=[-0.08,0.25], p = 0.296). In subgroup analyses, maltreatment effect sizes were larger for women as compared to men and for White participants as compared to non-White participants, although these differences were not statistically significant at the α= 0.05 level.

CONCLUSIONS AND RELEVANCE

As of midlife, effects of childhood maltreatment on biological aging are small in magnitude but discernible. Interventions to treat psychological and behavioral sequelae of exposure to childhood maltreatment, including in midlife adults, have potential to protect survivors from excess burden of disease, disability, and mortality in later life.

Generation and application of pseudo-long reads for metagenome assembly

BACKGROUND

Metagenomic assembly using high-throughput sequencing data is a powerful method to construct microbial genomes in environmental samples without cultivation. However, metagenomic assembly, especially when only short reads are available, is a complex and challenging task because mixed genomes of multiple microorganisms constitute the metagenome. Although long read sequencing technologies have been developed and have begun to be used for metagenomic assembly, many metagenomic studies have been performed based on short reads because the generation of long reads requires higher sequencing cost than short reads.

RESULTS

In this study, we present a new method called PLR-GEN. It creates pseudo-long reads from metagenomic short reads based on given reference genome sequences by considering small sequence variations existing in individual genomes of the same or different species. When applied to a mock community data set in the Human Microbiome Project, PLR-GEN dramatically extended short reads in length of 101 bp to pseudo-long reads with N50 of 33 Kbp and 0.4% error rate. The use of these pseudo-long reads generated by PLR-GEN resulted in an obvious improvement of metagenomic assembly in terms of the number of sequences, assembly contiguity, and prediction of species and genes.

CONCLUSIONS

PLR-GEN can be used to generate artificial long read sequences without spending extra sequencing cost, thus aiding various studies using metagenomes.

DLEB: a web application for building deep learning models in biological research

Deep learning has been applied for solving many biological problems, and it has shown outstanding performance. Applying deep learning in research requires knowledge of deep learning theories and programming skills, but researchers have developed diverse deep learning platforms to allow users to build deep learning models without programming. Despite these efforts, it is still difficult for biologists to use deep learning because of limitations of the existing platforms. Therefore, a new platform is necessary that can solve these challenges for biologists. To alleviate this situation, we developed a user-friendly and easy-to-use web application called DLEB (Deep Learning Editor for Biologists) that allows for building deep learning models specialized for biologists. DLEB helps researchers (i) design deep learning models easily and (ii) generate corresponding Python code to run directly in their machines. DLEB provides other useful features for biologists, such as recommending deep learning models for specific learning tasks and data, pre-processing of input biological data, and availability of various template models and example biological datasets for model training. DLEB can serve as a highly valuable platform for easily applying deep learning to solve many important biological problems. DLEB is freely available at http://dleb.konkuk.ac.kr/.

Long-term natural history of transient constrictive pericarditis

Background

In recent years, there has been increasing recognition of a potentially-reversible, transient/ subacute form of constrictive pericarditis (CP). To date, studies have been small with lack of long-term longitudinal follow-up.

Purpose

We aimed to elucidate the causes and natural history of subacute CP.

Methods

Patients were included if (1) they had a diagnosis of CP, (2) had cardiac magnetic resonance (CMR) within 12 months of symptom onset with evidence of pericardial delayed enhancement/ inflammation (Figure) (3) received anti-inflammatory medications.

Results

A total of 78 individuals were included, comprising 61 men (78%) with a mean age of 59±14 years. Causes of subacute CP included idiopathic/ viral pericarditis (58%), post-pericardiotomy (29%), autoimmune (6%), radiation therapy (3%) and others (4%). After median follow-up of 4.4 years, 31 (40%) required pericardiectomy. There were no deaths.

Patients who underwent pericardiectomy had longer duration of symptoms at presentation [6 (4–9) vs. 3 (2–5) months, P<0.01], were more likely to be on diuretic therapy (87 vs. 45%, P<0.001), had lower ultra-sensitive C-reactive protein [4.4 (2.6–13.1) vs. 11.95 (1.8–61.55) mg/dl, P<0.001] and lower erythrocyte sedimentation rate [5 (2–10) vs. 25 (6–43 mm/hr), P=0.031] compared to those who were managed medically. There were no other significant differences in clinical characteristics or baseline anti-inflammatory therapy. The presence of elevated inflammatory markers (HR: 0.18; 95% CI: 0.06–0.58, P<0.01) was an independent predictor of freedom from pericardiectomy after adjustment for relevant clinical and imaging parameters.

Conclusions

We present the original observations of the largest cohort of patients with transient CP to date and demonstrate that increased inflammatory markers were independently associated with long-term freedom from pericardiectomy. Our results suggest that a trial of anti-inflammatory therapy in the setting of elevated inflammatory markers may be appropriate prior to referral for surgery given the possible reversibility.

Pericardial delayed enhancement

Funding Acknowledgement

Type of funding source: None

0558 Sleep Disordered Breathing and Right Ventricular Electrocardiographic and Functional Characteristics in Group 1 Pulmonary Arterial Hypertension

Introduction

Right ventricular (RV) electrophysiologic and functional alterations related to sleep disordered breathing (SDB) in pulmonary arterial hypertension (PAH) are not well understood. We hypothesize an association between SDB and RV electrophysiological/functional measures in World Symposium of Pulmonary Hypertension (WSPH) Group 1 PAH.

Methods

The NHLBI multicenter PVDOMICS study (NCT02980887) enrolls patients with PAH undergoing a battery of assessments including home sleep apnea testing(NOX-T3, Carefusion®) or with historical sleep study data. Logistic(OR,95%CI) and linear(beta coefficients,95%CI) regression models adjusted for age, sex, race, body mass index (BMI, kg/m2), PAH medications, supplemental oxygen(O2), positive airway pressure(PAP) were used to assess associations of SDB(apnea hypopnea index,(AHI), ≥3% desaturations(hypopnea), percentage recording time with SaO2<90% (TRT<90%) with electrocardiographic measures: RV hypertrophy(RVH), right bundle branch block(RBBB), and right axis deviation(RAD), echocardiographically-derived RV systolic pressure(RVSP) and RV ejection fraction(RVEF) from cardiac MRI. Analyses were performed based on an overall significance level of 0.05, using SAS software (version 9.4, Cary, NC).

Results

The analysis consisted of 182 PAH participants with age: 52.5±13.9 years, 71.4% female, 88.9% Caucasian, BMI:30.3±7.8 kg/m2, RVEF: 37.3±11.6, and RVSP: 67.0±23.4. None of the electrocardiographic measures were associated with AHI and only RVH was significantly associated with TRT<90% (1.25:1.09,1.43),p=0.001. Although AHI was not associated with RVSP, a 10% increase in TRT<90% was associated with a 2.60mmHg increase in RVSP (2.60:1.44,3.76),p<0.001. Each 10-unit increase in AHI was associated with a 2.72% reduction of RVEF (-2.72:-4.89,-0.56),p=0.014, and each 10-unit increase in TRT<90% was associated with a 0.72% reduction of RVEF (-0.72:-1.38,-0.06),p=0.033.

Conclusion

We identify nocturnal hypoxia as a predictor of RV electrophysiological and functional alterations even after consideration of confounding factors. SDB as determined by AHI was also more so associated with reduced RVEF than hypoxia. Future mechanistic studies should focus on further elucidation of SDB and nocturnal hypoxia on pathogenesis of RV dysfunction in PAH.

Support

U01HL125218/U01HL125205/U01HL125212/U01HL125208/U01HL125175/U01HL125215, U01HL125177/Pulmonary Hypertension Association

TAMA: improved metagenomic sequence classification through meta-analysis

BACKGROUND

Microorganisms are important occupants of many different environments. Identifying the composition of microbes and estimating their abundance promote understanding of interactions of microbes in environmental samples. To understand their environments more deeply, the composition of microorganisms in environmental samples has been studied using metagenomes, which are the collections of genomes of the microorganisms. Although many tools have been developed for taxonomy analysis based on different algorithms, variability of analysis outputs of existing tools from the same input metagenome datasets is the main obstacle for many researchers in this field.

RESULTS

Here, we present a novel meta-analysis tool for metagenome taxonomy analysis, called TAMA, by intelligently integrating outputs from three different taxonomy analysis tools. Using an integrated reference database, TAMA performs taxonomy assignment for input metagenome reads based on a meta-score by integrating scores of taxonomy assignment from different taxonomy classification tools. TAMA outperformed existing tools when evaluated using various benchmark datasets. It was also successfully applied to obtain relative species abundance profiles and difference in composition of microorganisms in two types of cheese metagenome and human gut metagenome.

CONCLUSION

TAMA can be easily installed and used for metagenome read classification and the prediction of relative species abundance from multiple numbers and types of metagenome read samples. TAMA can be used to more accurately uncover the composition of microorganisms in metagenome samples collected from various environments, especially when the use of a single taxonomy analysis tool is unreliable. TAMA is an open source tool, and can be downloaded at https://github.com/jkimlab/TAMA.

Using INTERSPIA to Explore the Dynamics of Protein-Protein Interactions Among Multiple Species

INTER-Species Protein Interaction Analysis (INTERSPIA) is a web application for identifying diverse patterns of protein-protein interactions (PPIs) in different species. Given a set of proteins of interest to the user, INTERSPIA first discovers additional proteins that are functionally associated with the input proteins as well as different or common patterns of PPIs among the proteins in multiple species through a server-side pipeline. Second, it visualizes the dynamics of PPIs in multiple species via an easy-to-use web interface. This article contains a basic protocol describing how to visualize diverse patterns of PPIs of input proteins in multiple species, and how to use them for functional analysis in the web interface. INTERSPIA is freely available at http://bioinfo.konkuk.ac.kr/INTERSPIA/. © 2019 by John Wiley & Sons, Inc. Basic Protocol: Running INTERSPIA using a list of input proteins.

Integrative Meta-Assembly Pipeline (IMAP): Chromosome-level genome assembler combining multiple de novo assemblies

Background

Genomic data have become major resources to understand complex mechanisms at fine-scale temporal and spatial resolution in functional and evolutionary genetic studies, including human diseases, such as cancers. Recently, a large number of whole genomes of evolving populations of yeast (Saccharomyces cerevisiae W303 strain) were sequenced in a time-dependent manner to identify temporal evolutionary patterns. For this type of study, a chromosome-level sequence assembly of the strain or population at time zero is required to compare with the genomes derived later. However, there is no fully automated computational approach in experimental evolution studies to establish the chromosome-level genome assembly using unique features of sequencing data.

Methods and results

In this study, we developed a new software pipeline, the integrative meta-assembly pipeline (IMAP), to build chromosome-level genome sequence assemblies by generating and combining multiple initial assemblies using three de novo assemblers from short-read sequencing data. We significantly improved the continuity and accuracy of the genome assembly using a large collection of sequencing data and hybrid assembly approaches. We validated our pipeline by generating chromosome-level assemblies of yeast strains W303 and SK1, and compared our results with assemblies built using long-read sequencing and various assembly evaluation metrics. We also constructed chromosome-level sequence assemblies of S. cerevisiae strain Sigma1278b, and three commonly used fungal strains: Aspergillus nidulans A713, Neurospora crassa 73, and Thielavia terrestris CBS 492.74, for which long-read sequencing data are not yet available. Finally, we examined the effect of IMAP parameters, such as reference and resolution, on the quality of the final assembly of the yeast strains W303 and SK1.

Conclusions

We developed a cost-effective pipeline to generate chromosome-level sequence assemblies using only short-read sequencing data. Our pipeline combines the strengths of reference-guided and meta-assembly approaches. Our pipeline is available online at http://github.com/jkimlab/IMAP including a Docker image, as well as a Perl script, to help users install the IMAP package, including several prerequisite programs. Users can use IMAP to easily build the chromosome-level assembly for the genome of their interest.

INTERSPIA: a web application for exploring the dynamics of protein-protein interactions among multiple species

Proteins perform biological functions through cascading interactions with each other by forming protein complexes. As a result, interactions among proteins, called protein-protein interactions (PPIs) are not completely free from selection constraint during evolution. Therefore, the identification and analysis of PPI changes during evolution can give us new insight into the evolution of functions. Although many algorithms, databases and websites have been developed to help the study of PPIs, most of them are limited to visualize the structure and features of PPIs in a chosen single species with limited functions in the visualization perspective. This leads to difficulties in the identification of different patterns of PPIs in different species and their functional consequences. To resolve these issues, we developed a web application, called INTER-Species Protein Interaction Analysis (INTERSPIA). Given a set of proteins of user's interest, INTERSPIA first discovers additional proteins that are functionally associated with the input proteins and searches for different patterns of PPIs in multiple species through a server-side pipeline, and second visualizes the dynamics of PPIs in multiple species using an easy-to-use web interface. INTERSPIA is freely available at http://bioinfo.konkuk.ac.kr/INTERSPIA/.

GMASS: a novel measure for genome assembly structural similarity

BACKGROUND

Thanks to the recent advancements in next-generation sequencing (NGS) technologies, large amount of genomic data, which are short DNA sequences known as reads, has been accumulating. Diverse assemblers have been developed to generate high quality de novo assemblies using the NGS reads, but their output is very different because of algorithmic differences. However, there are not properly structured measures to show the similarity or difference in assemblies.

RESULTS

We developed a new measure, called the GMASS score, for comparing two genome assemblies in terms of their structure. The GMASS score was developed based on the distribution pattern of the number and coverage of similar regions between a pair of assemblies. The new measure was able to show structural similarity between assemblies when evaluated by simulated assembly datasets. The application of the GMASS score to compare assemblies in recently published benchmark datasets showed the divergent performance of current assemblers as well as its ability to compare assemblies.

CONCLUSION

The GMASS score is a novel measure for representing structural similarity between two assemblies. It will contribute to the understanding of assembly output and developing de novo assemblers.