Characterization and identification of hidden rare variants in the human genome

Background

By examining the genotype calls generated by the 1000 Genomes Project we discovered that the human reference genome GRCh37 contains almost 20,000 loci in which the reference allele has never been observed in healthy individuals and around 70,000 loci in which it has been observed only in the heterozygous state.

Results

We show that a large fraction of this rare reference allele (RRA) loci belongs to coding, functional and regulatory elements of the genome and could be linked to rare Mendelian disorders as well as cancer. We also demonstrate that classical germline and somatic variant calling tools are not capable to recognize the rare allele when present in these loci. To overcome such limitations, we developed a novel tool, named RAREVATOR, that is able to identify and call the rare allele in these genomic positions. By using a small cancer dataset we compared our tool with two state-of-the-art callers and we found that RAREVATOR identified more than 1,500 germline and 22 somatic RRA variants missed by the two methods and which belong to significantly mutated pathways.

Conclusions

These results show that, to date, the investigation of around 100,000 loci of the human genome has been missed by re-sequencing experiments based on the GRCh37 assembly and that our tool can fill the gap left by other methods. Moreover, the investigation of the latest version of the human reference genome, GRCh38, showed that although the GRC corrected almost all insertions and a small part of SNVs and deletions, a large number of functionally relevant RRAs still remain unchanged. For this reason, also future resequencing experiments, based on GRCh38, will benefit from RAREVATOR analysis results. RAREVATOR is freely available at http://sourceforge.net/projects/rarevator.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1481-9) contains supplementary material, which is available to authorized users.     

Synthesis,Structural Characterization,and Chiroptical Studies of Bidentate Salen‐Type Lanthanide (III) Complexes

The salen‐type ligand prepared with (R,R) diphenylethan‐1,2‐diamine and salicylaldehyde provides stable and inert complexes KLnL₂ upon simple reaction with lanthanide halides or pseudohalides LnX₃ (Ln = Tb³⁺‐Lu³⁺; X = Cl^? or TfO^?) of its potassium salt. All the complexes were completely characterized through nuclear magnetic resonance (NMR), electronic circular dichroism (ECD) in the UV and some (Er³⁺, Tm³⁺, Yb³⁺) also with Near‐IR ECD (NIR‐ECD) and luminescence (Tb³⁺, Tm³⁺). Careful analysis of the NMR shifts demonstrated that the complexes are isostructural in solution and afforded an accurate geometry. This was further confirmed by means of Density Functional Theory (DFT) optimization of the Lu³⁺ complex, and by comparing the ligand‐centered experimental and time‐dependent TD‐DFT computed UV‐ECD spectra. As final validation, we used the NIR‐ECD spectrum of the Yb³⁺ derivative calculated by means of Richardson's equations. The excellent match between calculated and experimental ECD spectra confirm the quality of the NMR structure.  Chirality 27:857–863, 2015. © 2015 Wiley Periodicals, Inc.     

Loss of lysosomal membrane protein NCU-G1 in mice results in spontaneous liver fibrosis with accumulation of lipofuscin and iron in Kupffer cells

Human kidney predominant protein, NCU-G1, is a highly conserved protein with an unknown biological function. Initially described as a nuclear protein, it was later shown to be a bona fide lysosomal integral membrane protein. To gain insight into the physiological function of NCU-G1, mice with no detectable expression of this gene were created using a gene-trap strategy, and Ncu-g1^gt/gt mice were successfully characterized. Lysosomal disorders are mainly caused by lack of or malfunctioning of proteins in the endosomal-lysosomal pathway. The clinical symptoms vary, but often include liver dysfunction. Persistent liver damage activates fibrogenesis and, if unremedied, eventually leads to liver fibrosis/cirrhosis and death. We demonstrate that the disruption of Ncu-g1 results in spontaneous liver fibrosis in mice as the predominant phenotype. Evidence for an increased rate of hepatic cell death, oxidative stress and active fibrogenesis were detected in Ncu-g1^gt/gt liver. In addition to collagen deposition, microscopic examination of liver sections revealed accumulation of autofluorescent lipofuscin and iron in Ncu-g1^gt/gt Kupffer cells. Because only a few transgenic mouse models have been identified with chronic liver injury and spontaneous liver fibrosis development, we propose that the Ncu-g1^gt/gt mouse could be a valuable new tool in the development of novel treatments for the attenuation of fibrosis due to chronic liver damage.KEY WORDS: NCU-G1, Lysosome, Fibrosis     

Maternal intake of vitamin E and birth defects,national birth defects prevention study, 1997 to 2005

BACKGROUND In a recent study, high maternal periconceptional intake of vitamin E was found to be associated with risk of congenital heart defects (CHDs). To explore this association further, we investigated the association between total daily vitamin E intake and selected birth defects. METHODS: We analyzed data from 4525 controls and 8665 cases from the 1997 to 2005 National Birth Defects Prevention Study. We categorized estimated periconceptional energy‐adjusted total daily vitamin E intake from diet and supplements into quartiles (referent, lowest quartile). Associations between quartiles of energy‐adjusted vitamin E intake and selected birth defects were adjusted for demographic, lifestyle, and nutritional factors. RESULTS: We observed a statistically significant association with the third quartile of vitamin E intake (odds ratio [OR], 1.17; 95% confidence interval [CI], 1.01–1.35) and all CHDs combined. Among CHD sub‐types, we observed associations with left ventricular outflow tract obstruction defects, and its sub‐type, coarctation of the aorta and the third quartile of vitamin E intake. Among defects other than CHDs, we observed associations between anorectal atresia and the third quartile of vitamin E intake (OR, 1.66; 95% CI, 1.01–2.72) and hypospadias and the fourth quartile of vitamin E intake (OR, 1.42; 95% CI, 1.09–1.87). CONCLUSION: Selected quartiles of energy‐adjusted estimated total daily vitamin E intake were associated with selected birth defects. However, because these few associations did not exhibit exposure‐response patterns consistent with increasing risk associated with increasing intake of vitamin E, further studies are warranted to corroborate our findings. Birth Defects Research (Part A), 100:647–657, 2014. © 2014 Wiley Periodicals, Inc.     

Two-Component System Cross-Regulation Integrates Bacillus anthracis Response to Heme and Cell Envelope Stress

Two-component signaling systems (TCSs) are one of the mechanisms that bacteria employ to sense and adapt to changes in the environment. A prototypical TCS functions as a phosphorelay from a membrane-bound sensor histidine kinase (HK) to a cytoplasmic response regulator (RR) that controls target gene expression. Despite significant homology in the signaling domains of HKs and RRs, TCSs are thought to typically function as linear systems with little to no cross-talk between non-cognate HK-RR pairs. Here we have identified several cell envelope acting compounds that stimulate a previously uncharacterized Bacillus anthracis TCS. Furthermore, this TCS cross-signals with the heme sensing TCS HssRS; therefore, we have named it HssRS interfacing TCS (HitRS). HssRS reciprocates cross-talk to HitRS, suggesting a link between heme toxicity and cell envelope stress. The signaling between HssRS and HitRS occurs in the parental B. anthracis strain; therefore, we classify HssRS-HitRS interactions as cross-regulation. Cross-talk between HssRS and HitRS occurs at both HK-RR and post-RR signaling junctions. Finally, HitRS also regulates a previously unstudied ABC transporter implicating this transporter in the response to cell envelope stress. This chemical biology approach to probing TCS signaling provides a new model for understanding how bacterial signaling networks are integrated to enable adaptation to complex environments such as those encountered during colonization of the vertebrate host.     

Rotavirus Increases Levels of Lipidated LC3 Supporting Accumulation of Infectious Progeny Virus without Inducing Autophagosome Formation

Replication of many RNA viruses benefits from subversion of the autophagic pathway through many different mechanisms. Rotavirus, the main etiologic agent of pediatric gastroenteritis worldwide, has been recently described to induce accumulation of autophagosomes as a mean for targeting viral proteins to the sites of viral replication. Here we show that the viral-induced increase of the lipidated form of LC3 does not correlate with an augmented formation of autophagosomes, as detected by immunofluorescence and electron microscopy. The LC3-II accumulation was found to be dependent on active rotavirus replication through the use of antigenically intact inactivated viral particles and of siRNAs targeting viral genes that are essential for viral replication. Silencing expression of LC3 or of Atg7, a protein involved in LC3 lipidation, resulted in a significant impairment of viral titers, indicating that these elements of the autophagic pathway are required at late stages of the viral cycle.     

The Personalized Advantage Index: Translating Research on Prediction into Individualized Treatment Recommendations. A Demonstration

Background

Advances in personalized medicine require the identification of variables that predict differential response to treatments as well as the development and refinement of methods to transform predictive information into actionable recommendations.

Objective

To illustrate and test a new method for integrating predictive information to aid in treatment selection, using data from a randomized treatment comparison.

Method

Data from a trial of antidepressant medications (N = 104) versus cognitive behavioral therapy (N = 50) for Major Depressive Disorder were used to produce predictions of post-treatment scores on the Hamilton Rating Scale for Depression (HRSD) in each of the two treatments for each of the 154 patients. The patient''s own data were not used in the models that yielded these predictions. Five pre-randomization variables that predicted differential response (marital status, employment status, life events, comorbid personality disorder, and prior medication trials) were included in regression models, permitting the calculation of each patient''s Personalized Advantage Index (PAI), in HRSD units.

Results

For 60% of the sample a clinically meaningful advantage (PAI≥3) was predicted for one of the treatments, relative to the other. When these patients were divided into those randomly assigned to their “Optimal” treatment versus those assigned to their “Non-optimal” treatment, outcomes in the former group were superior (d = 0.58, 95% CI .17—1.01).

Conclusions

This approach to treatment selection, implemented in the context of two equally effective treatments, yielded effects that, if obtained prospectively, would rival those routinely observed in comparisons of active versus control treatments.     

Assessment of Acropora palmata in the Mesoamerican Reef System

The once-dominant shallow reef-building coral Acropora palmata has suffered drastic geographical declines in the wider Caribbean from a disease epidemic that began in the late 1970s. At present there is a lack of quantitative data to determine whether this species is recovering over large spatial scales. Here, we use quantitative surveys conducted in 107 shallow-water reef sites between 2010 and 2012 to investigate the current distribution and abundance of A. palmata along the Mesoamerican Reef System (MRS). Using historical data we also explored how the distribution and abundance of this species has changed in the northern portion of the MRS between 1985 and 2010–2012. A. palmata was recorded in only a fifth of the surveyed reef sites in 2010–2012. In the majority of these reef sites the presence of A. palmata was patchy and rare. Only one site (Limones reef), in the northernmost portion of the MRS, presented considerably high A. palmata cover (mean: 34.7%, SD: 24.5%). At this site, the size-frequency distribution of A. palmata colonies was skewed towards small colony sizes; 84% of the colonies were healthy, however disease prevalence increased with colony size. A comparison with historical data showed that in the northern portion of the MRS, in 1985, A. palmata occurred in 74% of the 31 surveyed sites and had a mean cover of 7.7% (SD = 9.0), whereas in 2010–2012 this species was recorded in 48% of the sites with a mean cover of 2.9% (SD = 7.5). A. palmata populations along the MRS are failing to recover the distribution and abundance they had prior to the 1980s. Investigating the biological (e.g., population genetics) and environmental conditions (e.g., sources of stress) of the few standing reefs with relatively high A. palmata cover is crucial for the development of informed restoration models for this species.     

Exploiting the repertoire of CK2 inhibitors to target DYRK and PIM kinases

Advantage has been taken of the relative promiscuity of commonly used inhibitors of protein kinase CK2 to develop compounds that can be exploited for the selective inhibition of druggable kinases other than CK2 itself. Here we summarize data obtained by altering the scaffold of CK2 inhibitors to give rise to novel selective inhibitors of DYRK1A and to a powerful cell permeable dual inhibitor of PIM1 and CK2. In the former case one of the new compounds, C624 (naphto [1,2-b]benzofuran-5,9-diol) displays a potency comparable to that of the first-in-class DYRK1A inhibitor, harmine, lacking however the drawback of drastically inhibiting monoamine oxidase-A (MAO-A) as harmine does. On the other hand the promiscuous CK2 inhibitor 4,5,6,7-tetrabromo-1H-benzimidazole (TBI,TBBz) has been derivatized with a sugar moiety to generate a 1-(β-D-2′-deoxyribofuranosyl)-4,5,6,7-tetrabromo-1H-benzimidazole (TDB) compound which inhibits PIM1 and CK2 with comparably high efficacy (IC₅₀ values < 100 nM) and remarkable selectivity. TDB, unlike other dual PIM1/CK2 inhibitors described in the literature is readily cell permeable and displays a cytotoxic effect on cancer cells consistent with concomitant inhibition of both its onco-kinase targets. This article is part of a Special Issue entitled: Inhibitors of Protein Kinases (2012).     

Pharmacological and genetical evidence supporting nitric oxide requirement for 2,4-epibrassinolide regulation of root architecture in Arabidopsis thaliana

Brassinosteroids (BRs) regulate various physiological processes, such as tolerance to stresses and root growth. Recently, a connection was reported between BRs and nitric oxide (NO) in plant responses to abiotic stress. Here we present evidence supporting NO functions in BR signaling during root growth process. Arabidopsis seedlings treated with BR 24-epibrassinolide (BL) show increased lateral roots (LR) density, inhibition of primary root (PR) elongation and NO accumulation. Similar effects were observed adding the NO donor GSNO to BR-receptor mutant bri1-1. Furthermore, BL-induced responses in the root were abolished by the specific NO scavenger c-PTIO. The activities of nitrate reductase (NR) and nitric oxide synthase (NOS)-like, two NO generating enzymes were involved in BR signaling. These results demonstrate that BR increases the NO concentration in root cells, which is required for BR-induced changes in root architecture.