Using iterative random forest to find geospatial environmental and Sociodemographic predictors of suicide attempts

Introduction

Despite a recent global decrease in suicide rates, death by suicide has increased in the United States. It is therefore imperative to identify the risk factors associated with suicide attempts to combat this growing epidemic. In this study, we aim to identify potential risk factors of suicide attempt using geospatial features in an Artificial intelligence framework.

Methods

We use iterative Random Forest, an explainable artificial intelligence method, to predict suicide attempts using data from the Million Veteran Program. This cohort incorporated 405,540 patients with 391,409 controls and 14,131 attempts. Our predictive model incorporates multiple climatic features at ZIP-code-level geospatial resolution. We additionally consider demographic features from the American Community Survey as well as the number of firearms and alcohol vendors per 10,000 people to assess the contributions of proximal environment, access to means, and restraint decrease to suicide attempts. In total 1,784 features were included in the predictive model.

Results

Our results show that geographic areas with higher concentrations of married males living with spouses are predictive of lower rates of suicide attempts, whereas geographic areas where males are more likely to live alone and to rent housing are predictive of higher rates of suicide attempts. We also identified climatic features that were associated with suicide attempt risk by age group. Additionally, we observed that firearms and alcohol vendors were associated with increased risk for suicide attempts irrespective of the age group examined, but that their effects were small in comparison to the top features.

Discussion

Taken together, our findings highlight the importance of social determinants and environmental factors in understanding suicide risk among veterans.

Identification of Novel, Replicable Genetic Risk Loci for Suicidal Thoughts and Behaviors Among US Military Veterans

Importance

Suicide is a leading cause of death; however, the molecular genetic basis of suicidal thoughts and behaviors (SITB) remains unknown.

Objective

To identify novel, replicable genomic risk loci for SITB.

Design, Setting, and Participants

This genome-wide association study included 633 778 US military veterans with and without SITB, as identified through electronic health records. GWAS was performed separately by ancestry, controlling for sex, age, and genetic substructure. Cross-ancestry risk loci were identified through meta-analysis. Study enrollment began in 2011 and is ongoing. Data were analyzed from November 2021 to August 2022.

Main Outcome and Measures

SITB.

Results

A total of 633 778 US military veterans were included in the analysis (57 152 [9%] female; 121 118 [19.1%] African ancestry, 8285 [1.3%] Asian ancestry, 452 767 [71.4%] European ancestry, and 51 608 [8.1%] Hispanic ancestry), including 121 211 individuals with SITB (19.1%). Meta-analysis identified more than 200 GWS (P < 5 × 10-8) cross-ancestry risk single-nucleotide variants for SITB concentrated in 7 regions on chromosomes 2, 6, 9, 11, 14, 16, and 18. Top single-nucleotide variants were largely intronic in nature; 5 were independently replicated in ISGC, including rs6557168 in ESR1, rs12808482 in DRD2, rs77641763 in EXD3, rs10671545 in DCC, and rs36006172 in TRAF3. Associations for FBXL19 and AC018880.2 were not replicated. Gene-based analyses implicated 24 additional GWS cross-ancestry risk genes, including FURIN, TSNARE1, and the NCAM1-TTC12-ANKK1-DRD2 gene cluster. Cross-ancestry enrichment analyses revealed significant enrichment for expression in brain and pituitary tissue, synapse and ubiquitination processes, amphetamine addiction, parathyroid hormone synthesis, axon guidance, and dopaminergic pathways. Seven other unique European ancestry-specific GWS loci were identified, 2 of which (POM121L2 and METTL15/LINC02758) were replicated. Two additional GWS ancestry-specific loci were identified within the African ancestry (PET112/GATB) and Hispanic ancestry (intergenic locus on chromosome 4) subsets, both of which were replicated. No GWS loci were identified within the Asian ancestry subset; however, significant enrichment was observed for axon guidance, cyclic adenosine monophosphate signaling, focal adhesion, glutamatergic synapse, and oxytocin signaling pathways across all ancestries. Within the European ancestry subset, genetic correlations (r > 0.75) were observed between the SITB phenotype and a suicide attempt-only phenotype, depression, and posttraumatic stress disorder. Additionally, polygenic risk score analyses revealed that the Million Veteran Program polygenic risk score had nominally significant main effects in 2 independent samples of veterans of European and African ancestry.

Conclusions and Relevance

The findings of this analysis may advance understanding of the molecular genetic basis of SITB and provide evidence for ESR1, DRD2, TRAF3, and DCC as cross-ancestry candidate risk genes. More work is needed to replicate these findings and to determine if and how these genes might impact clinical care.

Structural variants identified using non-Mendelian inheritance patterns advance the mechanistic understanding of autism spectrum disorder

The heritability of autism spectrum disorder (ASD), based on 680,000 families and five countries, is estimated to be nearly 80%, yet heritability reported from SNP-based studies are consistently lower, and few significant loci have been identified with genome-wide association studies. This gap in genomic information may reside in rare variants, interaction among variants (epistasis), or cryptic structural variation (SV) and may provide mechanisms that underlie ASD. Here we use a method to identify potential SVs based on non-Mendelian inheritance patterns in pedigrees using parent-child genotypes from ASD families and demonstrate that they are enriched in ASD-risk genes. Most are in non-coding genic space and are over-represented in expression quantitative trait loci, suggesting that they affect gene regulation, which we confirm with their overlap of differentially expressed genes in postmortem brain tissue of ASD individuals. We then identify an SV in the GRIK2 gene that alters RNA splicing and a regulatory region of the ACMSD gene in the kynurenine pathway as significantly associated with a non-verbal ASD phenotype, supporting our hypothesis that these currently excluded loci can provide a clearer mechanistic understanding of ASD. Finally, we use an explainable artificial intelligence approach to define subgroups demonstrating their use in the context of precision medicine.

Structural and functional characterization of NEMO cleavage by SARS-CoV-2 3CLpro

In addition to its essential role in viral polyprotein processing, the SARS-CoV-2 3C-like protease (3CLpro) can cleave human immune signaling proteins, like NF-κB Essential Modulator (NEMO) and deregulate the host immune response. Here, in vitro assays show that SARS-CoV-2 3CLpro cleaves NEMO with fine-tuned efficiency. Analysis of the 2.50 Å resolution crystal structure of 3CLpro C145S bound to NEMO226-234 reveals subsites that tolerate a range of viral and host substrates through main chain hydrogen bonds while also enforcing specificity using side chain hydrogen bonds and hydrophobic contacts. Machine learning- and physics-based computational methods predict that variation in key binding residues of 3CLpro-NEMO helps explain the high fitness of SARS-CoV-2 in humans. We posit that cleavage of NEMO is an important piece of information to be accounted for, in the pathology of COVID-19.

NF-κB perturbation reveals unique immunomodulatory functions in Prx1+ fibroblasts that promote development of atopic dermatitis

Skin is composed of diverse cell populations that cooperatively maintain homeostasis. Up-regulation of the nuclear factor κB (NF-κB) pathway may lead to the development of chronic inflammatory disorders of the skin, but its role during the early events remains unclear. Through analysis of single-cell RNA sequencing data via iterative random forest leave one out prediction, an explainable artificial intelligence method, we identified an immunoregulatory role for a unique paired related homeobox-1 (Prx1)⁺ fibroblast subpopulation. Disruption of Ikkb-NF-κB under homeostatic conditions in these fibroblasts paradoxically induced skin inflammation due to the overexpression of C-C motif chemokine ligand 11 (CCL11; or eotaxin-1) characterized by eosinophil infiltration and a subsequent T_H2 immune response. Because the inflammatory phenotype resembled that seen in human atopic dermatitis (AD), we examined human AD skin samples and found that human AD fibroblasts also overexpressed CCL11 and that perturbation of Ikkb-NF-κB in primary human dermal fibroblasts up-regulated CCL11. Monoclonal antibody treatment against CCL11 was effective in reducing the eosinophilia and T_H2 inflammation in a mouse model. Together, the murine model and human AD specimens point to dysregulated Prx1⁺ fibroblasts as a previously unrecognized etiologic factor that may contribute to the pathogenesis of AD and suggest that targeting CCL11 may be a way to treat AD-like skin lesions.

Response to comment on 'SARS-CoV-2 suppresses anticoagulant and fibrinolytic gene expression in the lung'

Early in the SARS-CoV-2 pandemic, we compared transcriptome data from hospitalized COVID-19 patients and control patients without COVID-19. We found changes in procoagulant and fibrinolytic gene expression in the lungs of COVID-19 patients (Mast et al., 2021). These findings have been challenged based on issues with the samples (Fitzgerald and Jamieson, 2022). We have revisited our previous analyses in the light of this challenge and find that these new analyses support our original conclusions.

Antiviral Strategies Against SARS-CoV-2: A Systems Biology Approach

The unprecedented scientific achievements in combating the COVID-19 pandemic reflect a global response informed by unprecedented access to data. We now have the ability to rapidly generate a diversity of information on an emerging pathogen and, by using high-performance computing and a systems biology approach, we can mine this wealth of information to understand the complexities of viral pathogenesis and contagion like never before. These efforts will aid in the development of vaccines, antiviral medications, and inform policymakers and clinicians. Here we detail computational protocols developed as SARS-CoV-2 began to spread across the globe. They include pathogen detection, comparative structural proteomics, evolutionary adaptation analysis via network and artificial intelligence methodologies, and multiomic integration. These protocols constitute a core framework on which to build a systems-level infrastructure that can be quickly brought to bear on future pathogens before they evolve into pandemic proportions.

Structural and functional characterization of NEMO cleavage by SARS-CoV-2 3CLpro.. [PMID: 34816264 PMCID: PMC8609902 DOI: 10.1101/2021.11.11.468228] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

In addition to its essential role in viral polyprotein processing, the SARS-CoV-2 3C-like (3CLpro) protease can cleave human immune signaling proteins, like NF-κB Essential Modulator (NEMO) and deregulate the host immune response. Here, in vitro assays show that SARS-CoV-2 3CLpro cleaves NEMO with fine-tuned efficiency. Analysis of the 2.14 Å resolution crystal structure of 3CLpro C145S bound to NEMO_226–235 reveals subsites that tolerate a range of viral and host substrates through main chain hydrogen bonds while also enforcing specificity using side chain hydrogen bonds and hydrophobic contacts. Machine learning- and physics-based computational methods predict that variation in key binding residues of 3CLpro-NEMO helps explain the high fitness of SARS-CoV-2 in humans. We posit that cleavage of NEMO is an important piece of information to be accounted for in the pathology of COVID-19.

SARS-CoV-2 suppresses anticoagulant and fibrinolytic gene expression in the lung

Extensive fibrin deposition in the lungs and altered levels of circulating blood coagulation proteins in COVID-19 patients imply local derangement of pathways that limit fibrin formation and/or promote its clearance. We examined transcriptional profiles of bronchoalveolar lavage fluid (BALF) samples to identify molecular mechanisms underlying these coagulopathies. mRNA levels for regulators of the kallikrein-kinin (C1-inhibitor), coagulation (thrombomodulin, endothelial protein C receptor), and fibrinolytic (urokinase and urokinase receptor) pathways were significantly reduced in COVID-19 patients. While transcripts for several coagulation proteins were increased, those encoding tissue factor, the protein that initiates coagulation and whose expression is frequently increased in inflammatory disorders, were not increased in BALF from COVID-19 patients. Our analysis implicates enhanced propagation of coagulation and decreased fibrinolysis as drivers of the coagulopathy in the lungs of COVID-19 patients.

Potential Pathogenicity Determinants Identified from Structural Proteomics of SARS-CoV and SARS-CoV-2

Despite SARS-CoV and SARS-CoV-2 being equipped with highly similar protein arsenals, the corresponding zoonoses have spread among humans at extremely different rates. The specific characteristics of these viruses that led to such distinct outcomes remain unclear. Here, we apply proteome-wide comparative structural analysis aiming to identify the unique molecular elements in the SARS-CoV-2 proteome that may explain the differing consequences. By combining protein modeling and molecular dynamics simulations, we suggest nonconservative substitutions in functional regions of the spike glycoprotein (S), nsp1, and nsp3 that are contributing to differences in virulence. Particularly, we explain why the substitutions at the receptor-binding domain of S affect the structure–dynamics behavior in complexes with putative host receptors. Conservation of functional protein regions within the two taxa is also noteworthy. We suggest that the highly conserved main protease, nsp5, of SARS-CoV and SARS-CoV-2 is part of their mechanism of circumventing the host interferon antiviral response. Overall, most substitutions occur on the protein surfaces and may be modulating their antigenic properties and interactions with other macromolecules. Our results imply that the striking difference in the pervasiveness of SARS-CoV-2 and SARS-CoV among humans seems to significantly derive from molecular features that modulate the efficiency of viral particles in entering the host cells and blocking the host immune response.

Potentially adaptive SARS-CoV-2 mutations discovered with novel spatiotemporal and explainable AI models

BACKGROUND

A mechanistic understanding of the spread of SARS-CoV-2 and diligent tracking of ongoing mutagenesis are of key importance to plan robust strategies for confining its transmission. Large numbers of available sequences and their dates of transmission provide an unprecedented opportunity to analyze evolutionary adaptation in novel ways. Addition of high-resolution structural information can reveal the functional basis of these processes at the molecular level. Integrated systems biology-directed analyses of these data layers afford valuable insights to build a global understanding of the COVID-19 pandemic.

RESULTS

Here we identify globally distributed haplotypes from 15,789 SARS-CoV-2 genomes and model their success based on their duration, dispersal, and frequency in the host population. Our models identify mutations that are likely compensatory adaptive changes that allowed for rapid expansion of the virus. Functional predictions from structural analyses indicate that, contrary to previous reports, the Asp614Gly mutation in the spike glycoprotein (S) likely reduced transmission and the subsequent Pro323Leu mutation in the RNA-dependent RNA polymerase led to the precipitous spread of the virus. Our model also suggests that two mutations in the nsp13 helicase allowed for the adaptation of the virus to the Pacific Northwest of the USA. Finally, our explainable artificial intelligence algorithm identified a mutational hotspot in the sequence of S that also displays a signature of positive selection and may have implications for tissue or cell-specific expression of the virus.

CONCLUSIONS

These results provide valuable insights for the development of drugs and surveillance strategies to combat the current and future pandemics.

A mechanistic model and therapeutic interventions for COVID-19 involving a RAS-mediated bradykinin storm

Neither the disease mechanism nor treatments for COVID-19 are currently known. Here, we present a novel molecular mechanism for COVID-19 that provides therapeutic intervention points that can be addressed with existing FDA-approved pharmaceuticals. The entry point for the virus is ACE2, which is a component of the counteracting hypotensive axis of RAS. Bradykinin is a potent part of the vasopressor system that induces hypotension and vasodilation and is degraded by ACE and enhanced by the angiotensin1-9 produced by ACE2. Here, we perform a new analysis on gene expression data from cells in bronchoalveolar lavage fluid (BALF) from COVID-19 patients that were used to sequence the virus. Comparison with BALF from controls identifies a critical imbalance in RAS represented by decreased expression of ACE in combination with increases in ACE2, renin, angiotensin, key RAS receptors, kinogen and many kallikrein enzymes that activate it, and both bradykinin receptors. This very atypical pattern of the RAS is predicted to elevate bradykinin levels in multiple tissues and systems that will likely cause increases in vascular dilation, vascular permeability and hypotension. These bradykinin-driven outcomes explain many of the symptoms being observed in COVID-19.

Mitochondrial variation in small brown planthoppers linked to multiple traits and probably reflecting a complex evolutionary trajectory

While it has been proposed in several taxa that the mitochondrial genome is associated with adaptive evolution to different climatic conditions, making links between mitochondrial haplotypes and organismal phenotypes remains a challenge. Mitonuclear discordance occurs in the small brown planthopper (SBPH), Laodelphax striatellus, with one mitochondrial haplogroup (HGI) more common in the cold climate region of China relative to another form (HGII) despite strong nuclear gene flow, providing a promising model to investigate climatic adaptation of mitochondrial genomes. We hypothesized that cold adaptation through HGI may be involved, and considered mitogenome evolution, population genetic analyses, and bioassays to test this hypothesis. In contrast to our hypothesis, chill-coma recovery tests and population genetic tests of selection both pointed to HGII being involved in cold adaptation. Phylogenetic analyses revealed that HGII is nested within HGI, and has three nonsynonymous changes in ND2, ND5 and CYTB in comparison to HGI. These molecular changes likely increased mtDNA copy number, cold tolerance and fecundity of SBPH, particularly through a function-altering amino acid change involving M114T in ND2. Nuclear background also influenced fecundity and chill recovery (i.e., mitonuclear epistasis) and protein modelling indicates possible nuclear interactions for the two nonsynonymous changes in ND2 and CYTB. The high occurrence frequency of HGI in the cold climate region of China remains unexplained, but several possible reasons are discussed. Overall, our study points to a link between mtDNA variation and organismal-level evolution and suggests a possible role of mitonuclear interactions in maintaining mtDNA diversity.

Genotype to phenotype: Diet-by-mitochondrial DNA haplotype interactions drive metabolic flexibility and organismal fitness

Diet may be modified seasonally or by biogeographic, demographic or cultural shifts. It can differentially influence mitochondrial bioenergetics, retrograde signalling to the nuclear genome, and anterograde signalling to mitochondria. All these interactions have the potential to alter the frequencies of mtDNA haplotypes (mitotypes) in nature and may impact human health. In a model laboratory system, we fed four diets varying in Protein: Carbohydrate (P:C) ratio (1:2, 1:4, 1:8 and 1:16 P:C) to four homoplasmic Drosophila melanogaster mitotypes (nuclear genome standardised) and assayed their frequency in population cages. When fed a high protein 1:2 P:C diet, the frequency of flies harbouring Alstonville mtDNA increased. In contrast, when fed the high carbohydrate 1:16 P:C food the incidence of flies harbouring Dahomey mtDNA increased. This result, driven by differences in larval development, was generalisable to the replacement of the laboratory diet with fruits having high and low P:C ratios, perturbation of the nuclear genome and changes to the microbiome. Structural modelling and cellular assays suggested a V161L mutation in the ND4 subunit of complex I of Dahomey mtDNA was mildly deleterious, reduced mitochondrial functions, increased oxidative stress and resulted in an increase in larval development time on the 1:2 P:C diet. The 1:16 P:C diet triggered a cascade of changes in both mitotypes. In Dahomey larvae, increased feeding fuelled increased β-oxidation and the partial bypass of the complex I mutation. Conversely, Alstonville larvae upregulated genes involved with oxidative phosphorylation, increased glycogen metabolism and they were more physically active. We hypothesise that the increased physical activity diverted energy from growth and cell division and thereby slowed development. These data further question the use of mtDNA as an assumed neutral marker in evolutionary and population genetic studies. Moreover, if humans respond similarly, we posit that individuals with specific mtDNA variations may differentially metabolise carbohydrates, which has implications for a variety of diseases including cardiovascular disease, obesity, and perhaps Parkinson's Disease.

Sex-specific influences of mtDNA mitotype and diet on mitochondrial functions and physiological traits in Drosophila melanogaster

Here we determine the sex-specific influence of mtDNA type (mitotype) and diet on mitochondrial functions and physiology in two Drosophila melanogaster lines. In many species, males and females differ in aspects of their energy production. These sex-specific influences may be caused by differences in evolutionary history and physiological functions. We predicted the influence of mtDNA mutations should be stronger in males than females as a result of the organelle's maternal mode of inheritance in the majority of metazoans. In contrast, we predicted the influence of diet would be greater in females due to higher metabolic flexibility. We included four diets that differed in their protein: carbohydrate (P:C) ratios as they are the two-major energy-yielding macronutrients in the fly diet. We assayed four mitochondrial function traits (Complex I oxidative phosphorylation, reactive oxygen species production, superoxide dismutase activity, and mtDNA copy number) and four physiological traits (fecundity, longevity, lipid content, and starvation resistance). Traits were assayed at 11 d and 25 d of age. Consistent with predictions we observe that the mitotype influenced males more than females supporting the hypothesis of a sex-specific selective sieve in the mitochondrial genome caused by the maternal inheritance of mitochondria. Also, consistent with predictions, we found that the diet influenced females more than males.

Differential Expression of Genes that Control Respiration Contribute to Thermal Adaptation in Redband Trout (Oncorhynchus mykiss gairdneri)

Organisms can adapt to local environmental conditions as a plastic response or become adapted through natural selection on genetic variation. The ability to adapt to increased water temperatures will be of paramount importance for many fish species as the climate continues to warm and water resources become limited. Because increased water temperatures will reduce the dissolved oxygen available for fish, we hypothesized that adaptation to low oxygen environments would involve improved respiration through oxidative phosphorylation (OXPHOS). To test this hypothesis, we subjected individuals from two ecologically divergent populations of inland (redband) rainbow trout (Oncorhynchus mykiss gairdneri) with historically different temperature regimes (desert and montane) and their F1 progeny to diel cycles of temperature stress and then examined gene expression data for 80 nuclear- and mitochondrial-encoded OXPHOS subunits that participate in respiration. Of the 80 transcripts, 7 showed ≥ 2-fold difference in expression levels in gill tissue from desert fish under heat stress whereas the montane fish had none and the F1 only had one differentially expressed gene. A structural analysis of the proteins encoded by those genes suggests that the response could coordinate the formation of supercomplexes and oligomers. Supercomplexes may increase the efficiency of respiration because complexes I, III, and IV are brought into close proximity and oligomerization of complex V alters the macrostructure of mitochondria to improve respiration. Significant differences in gene expression patterns in response to heat stress in a common environment indicate that the response was not due to plasticity but had a genetic basis.

Evolution: are the monkeys' typewriters rigged?

Evolution is presumed to proceed by random mutations, which increase an individual's fitness. Increased fitness produces a higher survival rate for those individuals within populations and drives the variants to fixation over large timescales to produce new species. We recently identified positively selected sites in mitochondrial complex I in numerous, diverse taxa. In one taxon, a simple sequence repeat (SSR) encompassed the positively selected sites. We hypothesized a model in which: (i) slip-strand mis-pairing during replication due to the SSR increases the mutation rate at these sites, and (ii) a functional constraint at the protein level maintains the SSR and therefore a higher mutation rate at this site over large time scales to drive evolution. We tested this model by identifying SSRs in a mitochondrial-encoded protein in species from our previous work and determined that nearly all of the positively selected sites encompass an SSR. Furthermore, we show that our proposed model accounts for most of the mutations at neutral sites but it is probably the predominant mechanism at positively selected sites. This suggests that evolution does not proceed by simple random processes but is guided by physical properties of the DNA itself and functional constraint of the proteins encoded by the DNA.

Application of single nucleotide polymorphisms to non-model species: a technical review

Single nucleotide polymorphisms (SNPs) have gained wide use in humans and model species and are becoming the marker of choice for applications in other species. Technology that was developed for work in model species may provide useful tools for SNP discovery and genotyping in non-model organisms. However, SNP discovery can be expensive, labour intensive, and introduce ascertainment bias. In addition, the most efficient approaches to SNP discovery will depend on the research questions that the markers are to resolve as well as the focal species. We discuss advantages and disadvantages of several past and recent technologies for SNP discovery and genotyping and summarize a variety of SNP discovery and genotyping studies in ecology and evolution.

Recent physical connections may explain weak genetic structure in western Alaskan chum salmon (Oncorhynchus keta) populations

Low genetic divergence at neutral loci among populations is often the result of high levels of contemporary gene flow. Western Alaskan summer-run chum salmon (Oncorhynchus keta) populations demonstrate weak genetic structure, but invoking contemporary gene flow as the basis for the low divergence is problematic because salmon home to their natal streams and some of the populations are thousands of kilometers apart. We used genotypes from microsatellite and single nucleotide polymorphism loci to investigate alternative explanations for the current genetic structure of chum salmon populations from western Alaska. We also estimated current levels of gene flow among Kuskokwim River populations. Our results suggest that weak genetic structure is best explained by physical connections that occurred after the Holocene Thermal Maximum among the Yukon, Kuskokwim, and Nushagak drainages that allowed gene flow to occur among now distant populations.

Positive Darwinian selection in the piston that powers proton pumps in complex I of the mitochondria of Pacific salmon

The mechanism of oxidative phosphorylation is well understood, but evolution of the proteins involved is not. We combined phylogenetic, genomic, and structural biology analyses to examine the evolution of twelve mitochondrial encoded proteins of closely related, yet phenotypically diverse, Pacific salmon. Two separate analyses identified the same seven positively selected sites in ND5. A strong signal was also detected at three sites of ND2. An energetic coupling analysis revealed several structures in the ND5 protein that may have co-evolved with the selected sites. These data implicate Complex I, specifically the piston arm of ND5 where it connects the proton pumps, as important in the evolution of Pacific salmon. Lastly, the lineage to Chinook experienced rapid evolution at the piston arm.

Application of single nucleotide polymorphism markers to chum salmon Oncorhynchus keta: discovery, genotyping and linkage phase resolution

This study describes (1) the application of new methods to the discovery of informative single nucleotide polymorphism (SNP) markers in chum salmon Oncorhynchus keta, (2) a method to resolve the linkage phase of closely linked SNPs and (3) a method to inexpensively genotype them. Finally, it demonstrates that these SNPs provide information that discriminates among O. keta populations from different geographical regions of the northern Pacific Ocean. These informative markers can be used in conjunction with mixed-stock analysis to learn about the spatial and temporal marine distributions of O. keta and the factors that influence the distributions.

Single nucleotide polymorphisms in chum salmon (Oncorhynchus keta) mitochondrial DNA derived from restriction site haplotype information

Single nucleotide polymorphisms (SNPs) are useful genetic markers for the management and conservation of commercially important species such as salmon. Informative markers can be derived from data obtained for other purposes. We used restriction endonuclease data from earlier work to identify potentially useful restriction sites in chum salmon (Oncorhynchus keta). With the aid of a newly generated complete mitochondrial DNA sequence (accession number AP010773), we identified the SNP responsible for each restriction site variant, designed rapid genotyping assays, and surveyed the SNPs in more than 400 individuals. The restriction site analysis and the SNP genotyping assays were almost perfectly concordant. Some reasons for the non-concordance were identified and discussed.

ABCA1 is the cAMP-inducible apolipoprotein receptor that mediates cholesterol secretion from macrophages

Lipid-poor high density lipoprotein apolipoproteins remove cholesterol and phospholipids from cells by an active secretory pathway controlled by an ABC transporter called ABCA1. This pathway is induced by cholesterol and cAMP analogs in a cell-specific manner. Here we provide evidence that increased plasma membrane ABCA1 accounts for the enhanced apolipoprotein-mediated lipid secretion from macrophages induced by cAMP analogs. Treatment of RAW264 macrophages with 8-bromo-cAMP caused parallel increases in apoA-I-mediated cholesterol efflux, ABCA1 mRNA and protein levels, incorporation of ABCA1 into the plasma membrane, and binding of apoA-I to cell-surface ABCA1. All of these parameters declined to near base-line values within 6 h after removal of 8-bromo-cAMP, indicating that ABCA1 is highly unstable and is degraded rapidly in the absence of inducer. Thus, ABCA1 is likely to be the cAMP-inducible apolipoprotein receptor that promotes removal of cholesterol and phospholipids from macrophages.

Replication in restenotic atherectomy tissue

Previously, we demonstrated that replication in restenotic coronary atherectomy specimens was an infrequent and modest event. In general, this data was interpreted with caution, as immunocytochemistry for the proliferating cell nuclear antigen (PCNA) was used to subjectively assess proliferation and most of the tissue specimens were resected more than 3 months after the initial interventional procedure. The purpose of the present study was to use a more sensitive method of detecting replication, in situ hybridization for histone 3 (H3) mRNA, to determine the replication profile of human directional atherectomy specimens. Restenotic directional coronary atherectomy specimens from lesions that had undergone an interventional procedure within the preceding 3 months were studied. In addition, larger atherectomy specimens from peripheral arterial lesions were assessed to ensure that pockets of replication were not being overlooked in the smaller coronary specimens. We found evidence for replication in tissue resected from 2/17 coronary and 9/12 peripheral artery restenotic lesions. In contrast, 3/11 specimens resected from primary lesions of peripheral arteries also expressed H3 mRNA. We estimated that the maximum percentage of cells that were replicating in restenotic coronary, restenotic peripheral and primary peripheral artery tissue slides to be <0.5, < or =1.2 and <0.01%, respectively. Replication was found in tissue specimens resected both early and late after a previous interventional procedure. For specimens with >15 replicating cells per slide we found high levels of focal replication. Therefore, cell replication, as assessed by the expression of H3 mRNA, was infrequent in restenotic coronary artery specimens, whereas peripheral restenotic lesions had more frequent and higher levels of replication regardless of the interval from the previous interventional procedure. For all specimens the percentage of cells that were replicating was low, however focal areas with relatively high replication indices were presented. Although replication was more abundant in restenotic lesions it does not appear to be a dominant event in the pathophysiology of restenosis.

The Tangier disease gene product ABC1 controls the cellular apolipoprotein-mediated lipid removal pathway

The ABC1 transporter was identified as the defect in Tangier disease by a combined strategy of gene expression microarray analysis, genetic mapping, and biochemical studies. Patients with Tangier disease have a defect in cellular cholesterol removal, which results in near zero plasma levels of HDL and in massive tissue deposition of cholesteryl esters. Blocking the expression or activity of ABC1 reduces apolipoprotein-mediated lipid efflux from cultured cells, and increasing expression of ABC1 enhances it. ABC1 expression is induced by cholesterol loading and cAMP treatment and is reduced upon subsequent cholesterol removal by apolipoproteins. The protein is incorporated into the plasma membrane in proportion to its level of expression. Different mutations were detected in the ABC1 gene of 3 unrelated patients. Thus, ABC1 has the properties of a key protein in the cellular lipid removal pathway, as emphasized by the consequences of its defect in patients with Tangier disease.

Arterial expression of the plasminogen activator system early after cardiac transplantation

OBJECTIVES

Recent studies suggest that alterations in tissue thrombolysis as well as the inward migration of cells may be specific events that contribute to coronary artery narrowing after cardiac transplantation. Plasminogen activators and inhibitors play a central role in governing not only tissue thrombolysis, but also vascular cell migration. The purpose of this study was to examine arterial wall expression of the plasminogen activation system in coronary arteries during graft vascular disease initiation and progression.

METHODS

Using in situ hybridization and immunocytochemistry, the expression patterns of uPA and PAI-1 in coronary arteries from cardiac allografts were compared to those of young individuals without disease.

RESULTS

Both PAI-1 and uPA were over-expressed early after transplantation and as late as 27 months post grafting. Over-expression of these molecules preceded morphological evidence of graft vascular disease. Of special note was the adventitial expression of uPA and PAI-1 in microvessels and myofibroblasts. In contrast, the expression of uPA and PAI-1 in normal coronary arteries was confined to endothelial cells of the central lumen, as well as low levels of expression in intimal and medial smooth muscle cells.

CONCLUSIONS

Despite morphologic similarities between normal and transplant coronary arteries, differences were noted in the vascular expression pattern of uPA and PAI-1. The exact role of these molecules in graft vascular disease requires further study; however, it is intriguing to consider that a local imbalance in the plasminogen system may contribute to arterial wall thrombosis and/or excessive cell migration and the genesis of complex vascular lesions.

Beta ig-h3, a transforming growth factor-beta-inducible gene, is overexpressed in atherosclerotic and restenotic human vascular lesions

Transforming growth factor-beta (TGF-beta) plays an important role in vascular lesion formation and possibly the renarrowing process ("restenosis") that occurs after balloon angioplasty. Secreted in a latent form by most cells, TFG-beta requires enzymatic conversion before it is biologically active. TGF-beta-inducible gene h3 (beta ig-h3) is a novel molecule that is induced when cells are treated with TGF-beta1. This study examined the expression of beta ig-h3 in normal and diseased human vascular tissue. To determine the expression pattern of beta ig-h3 in human arteries, immunocytochemistry was performed on tissue sections from (1) normal internal mammary arteries, (2) the proximal left anterior descending coronary artery (with minimal intimal thickening) of 15 patients aged 18 to 40 years, (3) primary and restenotic coronary lesions from 7 patients, and (4) fresh directional atherectomy tissue from 11 patients. A polyclonal antibody consistently immunodetected beta ig-h3 protein in endothelial cells of all vascular tissue. In normal coronary arteries of young individuals, beta ig-h3 protein was absent from the intima and media but was found in the subendothelial smooth muscle cells of some arteries with modest intimal thickening. In diseased arteries beta ig-h3 protein was more abundant in the intima than the media. Restenotic coronary lesions tended to show higher levels of immunodetectable beta ig-h3 protein, especially in areas of dense fibrous connective tissue. Beta ig-h3 protein was immunodetected in the cytoplasm of plaque macrophages as well as smooth muscle and endothelial cells. By using in situ hybridization on fresh directional atherectomy specimens, we found beta ig-h3 mRNA to be overexpressed by plaque macrophages and smooth muscle cells. Nondiseased human internal mammary arteries also expressed beta ig-h3 mRNA in endothelial cells but not in the smooth muscle cells of the normal intima and media. These results document the expression of beta ig-h3 in diseased human arterial tissue and support the hypothesis that active TGF-beta plays a role in atherogenesis and restenosis.

Osteopontin is synthesized by macrophage, smooth muscle, and endothelial cells in primary and restenotic human coronary atherosclerotic plaques

How an atherosclerotic plaque evolves from minimal diffuse intimal hyperplasia to a critical lesion is not well understood. Cellular proliferation is a relatively infrequent and modest event in both primary and restenotic coronary atherectomy specimens, leading us to believe that other processes, such as the formation of extracellular matrix, cell migration, neovascularization, and calcification might be more important for lesion formation. The investigation of proteins that are overexpressed in plaque compared with the normal vessel wall may provide clues that will help determine which of these processes are key to lesion pathogenesis. One such molecule, osteopontin (OPN), is an arginine-glycine-aspartate-containing acidic phosphoprotein recently shown to be a novel component of human atherosclerotic plaques and selectively expressed in the rat neointima following balloon angioplasty. Using in situ hybridization and immunohistochemical methods, we demonstrate that in addition to macrophages, smooth muscle and endothelial cells synthesize OPN mRNA and protein in human coronary atherosclerotic plaque specimens obtained by directional atherectomy. In contrast, OPN mRNA and protein were not detected in nondiseased vessel walls. Furthermore, extracellular OPN protein collocalized with sites of early calcification in the plaque that were identified with a sensitive modification of the von Kossa staining technique. These findings, combined with studies showing that OPN has adhesive, chemotactic, and calcium-binding properties, suggest that OPN may contribute to cellular accumulations and dystrophic calcification in atherosclerotic plaques.

Angiogenesis in human coronary atherosclerotic plaques

Neovascularization in the walls of coronary arteries is associated with the presence of atherosclerotic plaque. The mechanisms responsible for the formation of these intraplaque microvessels are not understood. The purpose of this study is to examine the prevalence of endothelial cell replication in plaque microvessels. Two hundred and one primary and restenotic coronary atherectomy specimens were analyzed for the presence of microvessels and proliferation as reflected by positive immunolabeling for Ulex agglutinin and the proliferating cell nuclear antigen, respectively. In primary but not restenotic specimens, proliferation of any cell type was associated with the detection of microvessels on the same slide. However, intraplaque microvessels were more commonly found in restenotic compared to primary specimens (P = 0.004). Twelve highly vascularized specimens with evidence of replication were subjected to detailed histomorphological and quantitative image analyses. At 200 x, the most vascular optical field of each slide was identified and consistently included plaque macrophages. Total slide endothelial cell replication indices for these specimens varied, but in some instances were remarkably elevated (eg, 43.5%). The role of intraplaque angiogenesis may be analogous to that of tumor or wound angiogenesis and be important in development and progression of coronary artery lesions and restenosis.