Differential Roles of Cell Death-inducing DNA Fragmentation Factor-α-like Effector (CIDE) Proteins in Promoting Lipid Droplet Fusion and Growth in Subpopulations of Hepatocytes.

Lipid droplets (LDs) are dynamic subcellular organelles whose growth is closely linked to obesity and hepatic steatosis. Cell death-inducing DNA fragmentation factor-α-like effector (CIDE) proteins, including Cidea, Cideb, and Cidec (also called Fsp27), play important roles in lipid metabolism. Cidea and Cidec are LD-associated proteins that promote atypical LD fusion in adipocytes. Here, we find that CIDE proteins are all localized to LD-LD contact sites (LDCSs) and promote lipid transfer, LD fusion, and growth in hepatocytes. We have identified two types of hepatocytes, one with small LDs (small LD-containing hepatocytes, SLHs) and one with large LDs (large LD-containing hepatocytes, LLHs) in the liver. Cideb is localized to LDCSs and promotes lipid exchange and LD fusion in both SLHs and LLHs, whereas Cidea and Cidec are specifically localized to the LDCSs and promote lipid exchange and LD fusion in LLHs. Cideb-deficient SLHs have reduced LD sizes and lower lipid exchange activities. Fasting dramatically induces the expression of Cidea/Cidec and increases the percentage of LLHs in the liver. The majority of the hepatocytes from the liver of obese mice are Cidea/Cidec-positive LLHs. Knocking down Cidea or Cidec significantly reduced lipid storage in the livers of obese animals. Our data reveal that CIDE proteins play differential roles in promoting LD fusion and lipid storage; Cideb promotes lipid storage under normal diet conditions, whereas Cidea and Cidec are responsible for liver steatosis under fasting and obese conditions.

A sequence-based evolutionary distance method for Phylogenetic analysis of highly divergent proteins.

Because of the limited effectiveness of prevailing phylogenetic methods when applied to highly divergent protein sequences, the phylogenetic analysis problem remains challenging. Here, we propose a sequence-based evolutionary distance algorithm termed sequence distance (SD), which innovatively incorporates site-to-site correlation within protein sequences into the distance estimation. In protein superfamilies, SD can effectively distinguish evolutionary relationships both within and between protein families, producing phylogenetic trees that closely align with those based on structural information, even with sequence identity less than 20%. SD is highly correlated with the similarity of the protein structure, and can calculate evolutionary distances for thousands of protein pairs within seconds using a single CPU, which is significantly faster than most protein structure prediction methods that demand high computational resources and long run times. The development of SD will significantly advance phylogenetics, providing researchers with a more accurate and reliable tool for exploring evolutionary relationships.

Poor outcome of pediatric B-cell acute lymphoblastic leukemia associated with high level of CRLF2 gene expression in distinct molecular subtypes.

Background: Overexpression of the cytokine receptor-like factor 2 (CRLF2) gene is the most common feature in the Philadelphia chromosome (Ph)-like subtype of B-cell acute lymphoblastic leukemia (B-ALL). However, the predictive value of CRLF2 overexpression for the prognosis of pediatric B-ALL patients remain controversial. The molecular mechanisms that upregulate CRLF2 expression level in patients has not been fully elucidated. Methods: In this study, the prognostic impact of CRLF2 expression level on molecular types of B-ALL in pediatric patients from Zhujiang Hospital (n = 111) was retrospectively analyzed. Youden index analysis was used to categorize CRLF2 expression into 3 groups, and these categories more precisely described the differences in the prognosis of patients with varying expression levels of CRLF2 in both the Zhujiang Hospital cohort and the TARGET cohort. Results: We used the Zhujiang Hospital cohort as a discovery cohort to determine the cutoff value of CRLF2 expression. CRLF2-high patients accounted for approximately 6%. In addition, the percentage of bone marrow blast cells and initial white blood cell count in CRLF2-high patients were higher than those in CRLF2-low patients, and MRD turned negative slower. The results were validated in the TARGET cohort and indicated that CRLF2 overexpression could be subdivided by CRLF2 expression levels into 2 categories: CRLF2-high with a poor survival and CRLF2-medium with a good OS and EFS. Such heterogeneity was attributed to the different molecular mechanisms leading to CLRF2 upregulation, where the CRLF2 overexpression level was high in Ph-like B-ALL and medium in high hyperdiploid B-ALL. Conclusion: This study highlights the importance of the molecular mechanisms of the upregulation of CRLF2 expression in predicting the prognosis of pediatric B-ALL patients.

Control of cholesterol biosynthesis, uptake and storage in hepatocytes by Cideb.

Cideb, a member of CIDE family proteins, has emerged as an important regulator in the development of obesity and diabetes by controlling fatty acid synthesis and VLDL secretion in hepatocytes. Here, we investigated the role of Cideb in cholesterol biosynthesis, uptake and storage in the liver by using Cideb-null mice as a model system. Cideb-null mice and wild-type mice were treated with normal diet (ND) or high cholesterol diet (HCD) for one month. The metabolic parameters of cholesterol metabolism and expression profiles of genes in cholesterol biosynthesis and storage were measured. Cideb-null mice had lower levels of plasma cholesterol and LDL when fed with both ND and HCD and increased rate of cholesterol absorption. Furthermore, the liver of Cideb-null mice has lower rates of cholesterol biosynthesis and reduced expression levels of sterol response element-binding protein (SREBP) cleavage-activation protein (SCAP), and lower levels of nuclear form of SREBP2 and its downstream target genes in cholesterol biosynthesis pathway under a normal diet treatment. On the contrary, hepatic cholesterol biosynthesis rate between wild-type and Cideb-null mice was similar after high cholesterol diet treatment. Interestingly, hepatic cholesterol storage in the liver of Cideb-null mice was significantly increased due to its increased LDL receptor (LDLR) and acyl-CoA cholesterol acyltransferase (ACAT) expression. Finally, we observed drastically reduced cholesterol levels in the heart of Cideb-null mice fed with a high cholesterol diet. Overall, our data suggest that Cideb is a novel regulator in controlling cholesterol homeostasis in the liver. Therefore, Cideb could serve as an important therapeutical target for the treatment of atherosclerosis and cardiovascular diseases.

Establishing knowledge on the sequence arrangement pattern of nucleated protein folding.

The heat-tolerance mechanisms of (hyper)thermophilic proteins provide a unique opportunity to investigate the unsolved protein folding problem. In an attempt to determine whether the interval between residues in sequence might play a role in determining thermostability, we constructed a sequence interval-dependent value function to calculate the residue pair frequency. Additionally, we identified a new sequence arrangement pattern, where like-charged residues tend to be adjacently assembled, while unlike-charged residues are distributed over longer intervals, using statistical analysis of a large sequence database. This finding indicated that increasing the intervals between unlike-charged residues can increase protein thermostability, with the arrangement patterns of these charged residues serving as thermodynamically favorable nucleation points for protein folding. Additionally, we identified that the residue pairs K-E, R-E, L-V and V-V involving long sequence intervals play important roles involving increased protein thermostability. This work demonstrated a novel approach for considering sequence intervals as keys to understanding protein folding. Our findings of novel relationships between residue arrangement and protein thermostability can be used in industry and academia to aid the design of thermostable proteins.

HIV coreceptor tropism determination and mutational pattern identification.

In the early stages of infection, Human Immunodeficiency Virus Type 1 (HIV-1) generally selects CCR5 as the primary coreceptor for entering the host cell. As infection progresses, the virus evolves and may exhibit a coreceptor-switch to CXCR4. Accurate determination coreceptor usage and identification key mutational patterns associated tropism switch are essential for selection of appropriate therapies and understanding mechanism of coreceptor change. We developed a classifier composed of two coreceptor-specific weight matrices (CMs) based on a full-scale dataset. For this classifier, we found an AUC of 0.97, an accuracy of 95.21% and an MCC of 0.885 (sensitivity 92.92%; specificity 95.54%) in a ten-fold cross-validation, outperforming all other methods on an independent dataset (13% higher MCC value than geno2pheno and 15% higher MCC value than PSSM). A web server (http://spg.med.tsinghua.edu.cn/CM.html) based on our classifier was provided. Patterns of genetic mutations that occur along with coreceptor transitions were further identified based on the score of each sequence. Six pairs of one-AA mutational patterns and three pairs of two-AA mutational patterns were identified to associate with increasing propensity for X4 tropism. These mutational patterns offered new insights into the mechanism of coreceptor switch and aided in monitoring coreceptor switch.

Cross-Species Analysis of Gene Expression and Function in Prefrontal Cortex, Hippocampus and Striatum.

BACKGROUND: Mouse has been extensively used as a tool for investigating the onset and development of human neurological disorders. As a first step to construct a transgenic mouse model of human brain lesions, it is of fundamental importance to clarify the similarity and divergence of genetic background between non-diseased human and mouse brain tissues. METHODS: We systematically compared, based on large scale integrated microarray data, the transcriptomes of three anatomically distinct brain regions; prefrontal cortex (PFC), hippocampus (HIP) and striatum (STR), across human and mouse. The widely used DAVID web server was used to decipher the biological functions of the highly expressed genes that were identified using a previously reported approach. Venn analysis was used to depict the overlapping ratios of the notably enriched biological process (BP) terms (one-tailed Fisher's exact test and Benjamini correction; adjusted p < 0.01) between two brain tissues. GOSemSim, an R package, was selected to perform GO semantic similarity analysis. Next, we adjusted signal intensities of orthologous genes by the total signals in all samples within species, and used one minus Pearson's correlation coefficient to assess the expression distance. Hierarchical clustering and principal component analysis (PCA) were selected for expression pattern analysis. Lineage specific expressed orthologous genes were identified by comparison of the most extreme sub-datasets across species and further verified using reverse transcription PCR (RT-PCR) and quantitative real-time PCR (qRT-PCR). RESULTS: We found that the number of the significantly enriched BP terms of the highly expressed genes in human brain regions is larger than that in mouse corresponding brain regions. The mainly involved BP terms in human brain tissues associated with protein-membrane targeting and selenium metabolism are species-specific. The overlapping ratios of all the significantly enriched BP terms between any two brain tissues across species are lower than that within species, but the pairwise semantic similarities are very high between any two brain tissues from either human or mouse. Hierarchical clustering analysis shows the biological functions of the highly expressed genes in brain tissues are more consistent within species than interspecies; whereas it shows the expression patterns of orthologous genes are evidently conserved between human and mouse equivalent brain tissues. In addition, we identified four orthologous genes (COX5B, WIF1, SLC4A10 and PLA2G7) that are species-specific, which have been widely studied and confirmed to be closely linked with neuro- physiological and pathological functions. CONCLUSION: Our study highlights the similarities and divergences in gene function and expression between human and mouse corresponding brain regions, including PFC, HIP and STR.

Gene expression profile in the fat tissue of Fsp27 deficient mice.

Fsp27 is a lipid droplet-associated protein almost exclusively expressed in adipocytes where it facilitates unilocular lipid droplet formation. In mice, Fsp27 deficiency is associated with increased basal lipolysis, browning of white fat and a healthy metabolic profile, whereas energetically challenged Fsp27 deficient mice (ob/ob/Fsp27 (-/-)) show dramatically reduced fat mass, hepatic steatosis and insulin resistance which represents a typical lipodystrophy phenotype. Here, we investigate the effect of Fsp27 depletion on the gene expression of gonadal white adipose tissue (GWAT) under normal or energetically challenged condition (Fsp27 (-/-) vs Wild type; ob/ob/Fsp27 (-/-) vs ob/ob). We systematically analyzed the change in signaling pathway in Fsp27 deficient mice. The raw data have been deposited into Gene Expression Omnibus (GEO): GSE59807 and GSE22693.

Insulin resistance and white adipose tissue inflammation are uncoupled in energetically challenged Fsp27-deficient mice.

Fsp27 is a lipid droplet-associated protein almost exclusively expressed in adipocytes where it facilitates unilocular lipid droplet formation. In mice, Fsp27 deficiency is associated with increased basal lipolysis, 'browning' of white fat and a healthy metabolic profile, whereas a patient with congenital CIDEC deficiency manifested an adverse lipodystrophic phenotype. Here we reconcile these data by showing that exposing Fsp27-null mice to a substantial energetic stress by crossing them with ob/ob mice or BATless mice, or feeding them a high-fat diet, results in hepatic steatosis and insulin resistance. We also observe a striking reduction in adipose inflammation and increase in adiponectin levels in all three models. This appears to reflect reduced activation of the inflammasome and less adipocyte death. These findings highlight the importance of Fsp27 in facilitating optimal energy storage in adipocytes and represent a rare example where adipose inflammation and hepatic insulin resistance are disassociated.

High-accuracy identification of incident HIV-1 infections using a sequence clustering based diversity measure.

Accurate estimates of HIV-1 incidence are essential for monitoring epidemic trends and evaluating intervention efforts. However, the long asymptomatic stage of HIV-1 infection makes it difficult to effectively distinguish incident infections from chronic ones. Current incidence assays based on serology or viral sequence diversity are both still lacking in accuracy. In the present work, a sequence clustering based diversity (SCBD) assay was devised by utilizing the fact that viral sequences derived from each transmitted/founder (T/F) strain tend to cluster together at early stage, and that only the intra-cluster diversity is correlated with the time since HIV-1 infection. The dot-matrix pairwise alignment was used to eliminate the disproportional impact of insertion/deletions (indels) and recombination events, and so was the proportion of clusterable sequences (Pc) as an index to identify late chronic infections with declined viral genetic diversity. Tested on a dataset containing 398 incident and 163 chronic infection cases collected from the Los Alamos HIV database (last modified 2/8/2012), our SCBD method achieved 99.5% sensitivity and 98.8% specificity, with an overall accuracy of 99.3%. Further analysis and evaluation also suggested its performance was not affected by host factors such as the viral subtypes and transmission routes. The SCBD method demonstrated the potential of sequencing based techniques to become useful for identifying incident infections. Its use may be most advantageous for settings with low to moderate incidence relative to available resources. The online service is available at http://www.bioinfo.tsinghua.edu.cn:8080/SCBD/index.jsp.

Accurate prediction of protein structural class.

Because of the increasing gap between the data from sequencing and structural genomics, the accurate prediction of the structural class of a protein domain solely from the primary sequence has remained a challenging problem in structural biology. Traditional sequence-based predictors generally select several sequence features and then feed them directly into a classification program to identify the structural class. The current best sequence-based predictor achieved an overall accuracy of 74.1% when tested on a widely used, non-homologous benchmark dataset 25PDB. In the present work, we built a multiple linear regression (MLR) model to convert the 440-dimensional (440D) sequence feature vector extracted from the Position Specific Scoring Matrix (PSSM) of a protein domain to a 4-dimensinal (4D) structural feature vector, which could then be used to predict the four major structural classes. We performed 10-fold cross-validation and jackknife tests of the method on a large non-homologous dataset containing 8,244 domains distributed among the four major classes. The performance of our approach outperformed all of the existing sequence-based methods and had an overall accuracy of 83.1%, which is even higher than the results of those predicted secondary structure-based methods.

A new pattern-based method for identifying recent HIV-1 infections from the viral env sequence.

The long asymptomatic stage of HIV infection poses a great challenge in identifying recent HIV infections. This is a bottleneck for monitoring HIV epidemic trends and evaluating the effectiveness of national AIDS control programs. Several serological methods were used to address this issue with some success. Because of high false-positive rates in patients with advanced infection or in ART treatment, UNAIDS still hesitates to recommend their use in routine surveillance. We developed a new pattern-based method for measuring intra-patient viral genetic diversity for determination of recent infections and estimation of population incidence. This method is verified by using several datasets (424 subtype B and 77 CRF07_BC samples) with clearly identified HIV-1 infection times. Pattern-based diversities of recent infections are significantly lower than that of chronic ones. With larger window periods varying from 200 to 350 days, a higher accuracy (90%-95%) not affected by advanced disease nor ART treatment could be obtained. The pattern-based genetic method is supplementary to the existing serology-based assays, both of which could be suitable for use in low and high epidemic regions, respectively.

Salt bridges in the hyperthermophilic protein Ssh10b are resilient to temperature increases.

A double mutant cycle (DMC) approach was employed to estimate the effect of temperature on the contribution of two highly conserved salt bridges to protein stability in the hyperthermophilic protein Ssh10b. The coupling free energy were 2.4 +/- 0.4 kJ/mol at 298 K and 2.2 +/- 0.4 kJ/mol at 353 K for Glu-54/Arg-57, and 6.0 +/- 0.2 kJ/mol at 298 K and 5.9 +/- 0.6 kJ/mol at 353 K for Glu-36/Lys-68. The stability free energy of Ssh10b decrease greatly with increasing temperature, while the direct contribution of these two salt bridges to protein stability remain almost constant, providing evidence supporting the theoretical prediction that salt bridges are extremely resilient to temperature increases and thus are specially suited to improving protein stability at high temperatures. The reason for the difference in coupling free energy between salt bridges Glu-54/Arg-57 and Glu-36/Lys-68 is discussed. Comparing our results with published DMC data for the contribution of salt bridges to stability in other proteins, we found that the energy contribution of a salt bridge formed by two charged residues far apart in the primary sequence is higher than that of those formed between two very close ones. Implications of this finding are useful for engineering proteins with enhanced thermostability.