Comparative analysis of tandem repeats from hundreds of species reveals unique insights into centromere evolution

Background

Centromeres are essential for chromosome segregation, yet their DNA sequences evolve rapidly. In most animals and plants that have been studied, centromeres contain megabase-scale arrays of tandem repeats. Despite their importance, very little is known about the degree to which centromere tandem repeats share common properties between different species across different phyla. We used bioinformatic methods to identify high-copy tandem repeats from 282 species using publicly available genomic sequence and our own data.

Results

Our methods are compatible with all current sequencing technologies. Long Pacific Biosciences sequence reads allowed us to find tandem repeat monomers up to 1,419 bp. We assumed that the most abundant tandem repeat is the centromere DNA, which was true for most species whose centromeres have been previously characterized, suggesting this is a general property of genomes. High-copy centromere tandem repeats were found in almost all animal and plant genomes, but repeat monomers were highly variable in sequence composition and length. Furthermore, phylogenetic analysis of sequence homology showed little evidence of sequence conservation beyond approximately 50 million years of divergence. We find that despite an overall lack of sequence conservation, centromere tandem repeats from diverse species showed similar modes of evolution.

Conclusions

While centromere position in most eukaryotes is epigenetically determined, our results indicate that tandem repeats are highly prevalent at centromeres of both animal and plant genomes. This suggests a functional role for such repeats, perhaps in promoting concerted evolution of centromere DNA across chromosomes.     

Validating clustering for gene expression data

MOTIVATION: Many clustering algorithms have been proposed for the analysis of gene expression data, but little guidance is available to help choose among them. We provide a systematic framework for assessing the results of clustering algorithms. Clustering algorithms attempt to partition the genes into groups exhibiting similar patterns of variation in expression level. Our methodology is to apply a clustering algorithm to the data from all but one experimental condition. The remaining condition is used to assess the predictive power of the resulting clusters-meaningful clusters should exhibit less variation in the remaining condition than clusters formed by chance. RESULTS: We successfully applied our methodology to compare six clustering algorithms on four gene expression data sets. We found our quantitative measures of cluster quality to be positively correlated with external standards of cluster quality.     

The Characterization of Twenty Sequenced Human Genomes

We present the analysis of twenty human genomes to evaluate the prospects for identifying rare functional variants that contribute to a phenotype of interest. We sequenced at high coverage ten “case” genomes from individuals with severe hemophilia A and ten “control” genomes. We summarize the number of genetic variants emerging from a study of this magnitude, and provide a proof of concept for the identification of rare and highly-penetrant functional variants by confirming that the cause of hemophilia A is easily recognizable in this data set. We also show that the number of novel single nucleotide variants (SNVs) discovered per genome seems to stabilize at about 144,000 new variants per genome, after the first 15 individuals have been sequenced. Finally, we find that, on average, each genome carries 165 homozygous protein-truncating or stop loss variants in genes representing a diverse set of pathways.     

How accurately is ncRNA aligned within whole-genome multiple alignments?

Background  

Multiple alignment of homologous DNA sequences is of great interest to biologists since it provides a window into evolutionary processes. At present, the accuracy of whole-genome multiple alignments, particularly in noncoding regions, has not been thoroughly evaluated.     

Anticancer immunotherapy by CTLA-4 blockade: obligatory contribution of IL-2 receptors and negative prognostic impact of soluble CD25

The cytotoxic T lymphocyte antigen-4 (CTLA-4)-blocking antibody ipilimumab induces immune-mediated long-term control of metastatic melanoma in a fraction of patients. Although ipilimumab undoubtedly exerts its therapeutic effects via immunostimulation, thus far clinically useful, immunologically relevant biomarkers that predict treatment efficiency have been elusive. Here, we show that neutralization of IL-2 or blocking the α and β subunits of the IL-2 receptor (CD25 and CD122, respectively) abolished the antitumor effects and the accompanying improvement of the ratio of intratumoral T effector versus regulatory cells (Tregs), which were otherwise induced by CTLA-4 blockade in preclinical mouse models. CTLA-4 blockade led to the reduction of a suppressive CD4⁺ T cell subset expressing Lag3, ICOS, IL-10 and Egr2 with a concomitant rise in IL-2-producing effector cells that lost FoxP3 expression and accumulated in regressing tumors. While recombinant IL-2 improved the therapeutic efficacy of CTLA-4 blockade, the decoy IL-2 receptor α (IL-2Rα, sCD25) inhibited the anticancer effects of CTLA-4 blockade. In 262 metastatic melanoma patients receiving ipilimumab, baseline serum concentrations of sCD25 represented an independent indicator of overall survival, with high levels predicting resistance to therapy. Altogether, these results unravel a role for IL-2 and IL-2 receptors in the anticancer activity of CTLA-4 blockade. Importantly, our study provides the first immunologically relevant biomarker, namely elevated serum sCD25, that predicts resistance to CTLA-4 blockade in patients with melanoma.     

Sequence-based heuristics for faster annotation of non-coding RNA families

MOTIVATION: Non-coding RNAs (ncRNAs) are functional RNA molecules that do not code for proteins. Covariance Models (CMs) are a useful statistical tool to find new members of an ncRNA gene family in a large genome database, using both sequence and, importantly, RNA secondary structure information. Unfortunately, CM searches are extremely slow. Previously, we created rigorous filters, which provably sacrifice none of a CM's accuracy, while making searches significantly faster for virtually all ncRNA families. However, these rigorous filters make searches slower than heuristics could be. RESULTS: In this paper we introduce profile HMM-based heuristic filters. We show that their accuracy is usually superior to heuristics based on BLAST. Moreover, we compared our heuristics with those used in tRNAscan-SE, whose heuristics incorporate a significant amount of work specific to tRNAs, where our heuristics are generic to any ncRNA. Performance was roughly comparable, so we expect that our heuristics provide a high-quality solution that--unlike family-specific solutions--can scale to hundreds of ncRNA families. AVAILABILITY: The source code is available under GNU Public License at the supplementary web site.     

Activation of hypothalamic RIP‐Cre neurons promotes beiging of WAT via sympathetic nervous system

Activation of brown adipose tissue (BAT) and beige fat by cold increases energy expenditure. Although their activation is known to be differentially regulated in part by hypothalamus, the underlying neural pathways and populations remain poorly characterized. Here, we show that activation of rat‐insulin‐promoter‐Cre (RIP‐Cre) neurons in ventromedial hypothalamus (VMH) preferentially promotes recruitment of beige fat via a selective control of sympathetic nervous system (SNS) outflow to subcutaneous white adipose tissue (sWAT), but has no effect on BAT. Genetic ablation of APPL2 in RIP‐Cre neurons diminishes beiging in sWAT without affecting BAT, leading to cold intolerance and obesity in mice. Such defects are reversed by activation of RIP‐Cre neurons, inactivation of VMH AMPK, or treatment with a β3‐adrenergic receptor agonist. Hypothalamic APPL2 enhances neuronal activation in VMH RIP‐Cre neurons and raphe pallidus, thereby eliciting SNS outflow to sWAT and subsequent beiging. These data suggest that beige fat can be selectively activated by VMH RIP‐Cre neurons, in which the APPL2–AMPK signaling axis is crucial for this defending mechanism to cold and obesity.     

<Emphasis Type="Italic">mab-31</Emphasis> and the TGF-β pathway act in the ray lineage to pattern <Emphasis Type="Italic">C. elegans</Emphasis>male sensory rays

Background  

C. elegans TGF-β-like Sma/Mab signaling pathway regulates both body size and sensory ray patterning. Most of the components in this pathway were initially identified by genetic screens based on the small body phenotype, and many of these mutants display sensory ray patterning defect. At the cellular level, little is known about how and where these components work although ray structural cell has been implicated as one of the targets. Based on the specific ray patterning abnormality, we aim to identify by RNAi approach additional components that function specifically in the ray lineage to elucidate the regulatory role of TGF-β signaling in ray differentiation.