Estimating Maximal Symmetries of Regression Functions via Subgroup Lattices

We present a method for estimating the maximal symmetry of a regression function. Knowledge of such a symmetry can be used to significantly improve modelling by removing the modes of variation resulting from the symmetries. Symmetry estimation is carried out using hypothesis testing for invariance strategically over the subgroup lattice of a search group G acting on the feature space. We show that the estimation of the unique maximal invariant subgroup of G can be achieved by testing on only a finite portion of the subgroup lattice when G_max is a compact subgroup of G, even for infinite search groups and lattices (such as for the 3D rotation group SO(3)). We then show that the estimation is consistent when G is finite. We demonstrate the performance of this estimator in low dimensional simulations, on a synthetic image classification on MNIST data, and apply the methods to an application using satellite measurements of the earth's magnetic field.Comment: 35 Pages, 11 figure

The genomes of two key bumblebee species with primitive eusocial organization

Background: The shift from solitary to social behavior is one of the major evolutionary transitions. Primitively eusocial bumblebees are uniquely placed to illuminate the evolution of highly eusocial insect societies. Bumblebees are also invaluable natural and agricultural pollinators, and there is widespread concern over recent population declines in some species. High-quality genomic data will inform key aspects of bumblebee biology, including susceptibility to implicated population viability threats. Results: We report the high quality draft genome sequences of Bombus terrestris and Bombus impatiens, two ecologically dominant bumblebees and widely utilized study species. Comparing these new genomes to those of the highly eusocial honeybee Apis mellifera and other Hymenoptera, we identify deeply conserved similarities, as well as novelties key to the biology of these organisms. Some honeybee genome features thought to underpin advanced eusociality are also present in bumblebees, indicating an earlier evolution in the bee lineage. Xenobiotic detoxification and immune genes are similarly depauperate in bumblebees and honeybees, and multiple categories of genes linked to social organization, including development and behavior, show high conservation. Key differences identified include a bias in bumblebee chemoreception towards gustation from olfaction, and striking differences in microRNAs, potentially responsible for gene regulation underlying social and other traits. Conclusions: These two bumblebee genomes provide a foundation for post-genomic research on these key pollinators and insect societies. Overall, gene repertoires suggest that the route to advanced eusociality in bees was mediated by many small changes in many genes and processes, and not by notable expansion or depauperation

Autocrine Activation of the MET Receptor Tyrosine Kinase in Acute Myeloid Leukemia

Although the treatment of acute myeloid leukemia (AML) has improved significantly, more than half of all patients develop disease that is refractory to intensive chemotherapy. Functional genomics approaches offer a means to discover specific molecules mediating aberrant growth and survival of cancer cells. Thus, using a loss-of-function RNA interference genomic screen, we identified aberrant expression of the hepatocyte growth factor (HGF) as a critical factor in AML pathogenesis. We found HGF expression leading to autocrine activation of its receptor tyrosine kinase, MET, in nearly half of the AML cell lines and clinical samples studied. Genetic depletion of HGF or MET potently inhibited the growth and survival of HGF-expressing AML cells. However, leukemic cells treated with the specific MET kinase inhibitor crizotinib developed resistance due to compensatory upregulation of HGF expression, leading to restoration of MET signaling. In cases of AML where MET is coactivated with other tyrosine kinases, such as fibroblast growth factor receptor 1 (FGFR1), concomitant inhibition of FGFR1 and MET blocked compensatory HGF upregulation, resulting in sustained logarithmic cell kill both in vitro and in xenograft models in vivo. Our results demonstrate widespread dependence of AML cells on autocrine activation of MET, as well as the importance of compensatory upregulation of HGF expression in maintaining leukemogenic signaling by this receptor. We anticipate that these findings will lead to the design of additional strategies to block adaptive cellular responses that drive compensatory ligand expression as an essential component of the targeted inhibition of oncogenic receptors in human cancers

A 100%-complete sequence reveals unusually simple genomic features in the hot-spring red alga Cyanidioschyzon merolae

<p>Abstract</p> <p>Background</p> <p>All previously reported eukaryotic nuclear genome sequences have been incomplete, especially in highly repeated units and chromosomal ends. Because repetitive DNA is important for many aspects of biology, complete chromosomal structures are fundamental for understanding eukaryotic cells. Our earlier, nearly complete genome sequence of the hot-spring red alga <it>Cyanidioschyzon merolae </it>revealed several unique features, including just three ribosomal DNA copies, very few introns, and a small total number of genes. However, because the exact structures of certain functionally important repeated elements remained ambiguous, that sequence was not complete. Obviously, those ambiguities needed to be resolved before the unique features of the <it>C. merolae </it>genome could be summarized, and the ambiguities could only be resolved by completing the sequence. Therefore, we aimed to complete all previous gaps and sequence all remaining chromosomal ends, and now report the first nuclear-genome sequence for any eukaryote that is 100% complete.</p> <p>Results</p> <p>Our present complete sequence consists of 16546747 nucleotides covering 100% of the 20 linear chromosomes from telomere to telomere, representing the simple and unique chromosomal structures of the eukaryotic cell. We have unambiguously established that the <it>C. merolae </it>genome contains the smallest known histone-gene cluster, a unique telomeric repeat for all chromosomal ends, and an extremely low number of transposons.</p> <p>Conclusion</p> <p>By virtue of these attributes and others that we had discovered previously, <it>C. merolae </it>appears to have the simplest nuclear genome of the non-symbiotic eukaryotes. These unusually simple genomic features in the 100% complete genome sequence of <it>C. merolae </it>are extremely useful for further studies of eukaryotic cells.</p

Effects of the high-density lipoprotein mimetic agent CER-001 on coronary atherosclerosis in patients with acute coronary syndromes: a randomized trial†

Aim High-density lipoproteins (HDLs) have several potentially protective vascular effects. Most clinical studies of therapies targeting HDL have failed to show benefits vs. placebo. Objective To investigate the effects of an HDL-mimetic agent on atherosclerosis by intravascular ultrasonography (IVUS) and quantitative coronary angiography (QCA). Design and setting A prospective, double-blinded, randomized trial was conducted at 51 centres in the USA, the Netherlands, Canada, and France. Intravascular ultrasonography and QCA were performed to assess coronary atherosclerosis at baseline and 3 (2-5) weeks after the last study infusion. Patients Five hundred and seven patients were randomized; 417 and 461 had paired IVUS and QCA measurements, respectively. Intervention Patients were randomized to receive 6 weekly infusions of placebo, 3 mg/kg, 6 mg/kg, or 12 mg/kg CER-001. Main outcome measures The primary efficacy parameter was the nominal change in the total atheroma volume. Nominal changes in per cent atheroma volume on IVUS and coronary scores on QCA were also pre-specified endpoints. Results The nominal change in the total atheroma volume (adjusted means) was −2.71, −3.13, −1.50, and −3.05 mm3 with placebo, CER-001 3 mg/kg, 6 mg/kg, and 12 mg/kg, respectively (primary analysis of 12 mg/kg vs. placebo: P = 0.81). There was also no difference among groups for the nominal change in per cent atheroma volume (0.02, −0.02, 0.01, and 0.19%; nominal P = 0.53 for 12 mg/kg vs. placebo). Change in the coronary artery score was −0.022, −0.036, −0.022, and −0.015 mm (nominal P = 0.25, 0.99, 0.55), and change in the cumulative coronary stenosis score was −0.51, 2.65, 0.71, and −0.77% (compared with placebo, nominal P = 0.85 for 12 mg/kg and nominal P = 0.01 for 3 mg/kg). The number of patients with major cardiovascular events was 10 (8.3%), 16 (13.3%), 17 (13.7%), and 12 (9.8%) in the four groups. Conclusion CER-001 infusions did not reduce coronary atherosclerosis on IVUS and QCA when compared with placebo. Whether CER-001 administered in other regimens or to other populations could favourably affect atherosclerosis must await further study. Name of the trial registry: Clinicaltrials.gov; Registry's URL: http://clinicaltrials.gov/ct2/show/NCT01201837?term=cer-001&rank=2; Trial registration number: NCT0120183

The first myriapod genome sequence reveals conservative arthropod gene content and genome organisation in the centipede Strigamia maritima.

Myriapods (e.g., centipedes and millipedes) display a simple homonomous body plan relative to other arthropods. All members of the class are terrestrial, but they attained terrestriality independently of insects. Myriapoda is the only arthropod class not represented by a sequenced genome. We present an analysis of the genome of the centipede Strigamia maritima. It retains a compact genome that has undergone less gene loss and shuffling than previously sequenced arthropods, and many orthologues of genes conserved from the bilaterian ancestor that have been lost in insects. Our analysis locates many genes in conserved macro-synteny contexts, and many small-scale examples of gene clustering. We describe several examples where S. maritima shows different solutions from insects to similar problems. The insect olfactory receptor gene family is absent from S. maritima, and olfaction in air is likely effected by expansion of other receptor gene families. For some genes S. maritima has evolved paralogues to generate coding sequence diversity, where insects use alternate splicing. This is most striking for the Dscam gene, which in Drosophila generates more than 100,000 alternate splice forms, but in S. maritima is encoded by over 100 paralogues. We see an intriguing linkage between the absence of any known photosensory proteins in a blind organism and the additional absence of canonical circadian clock genes. The phylogenetic position of myriapods allows us to identify where in arthropod phylogeny several particular molecular mechanisms and traits emerged. For example, we conclude that juvenile hormone signalling evolved with the emergence of the exoskeleton in the arthropods and that RR-1 containing cuticle proteins evolved in the lineage leading to Mandibulata. We also identify when various gene expansions and losses occurred. The genome of S. maritima offers us a unique glimpse into the ancestral arthropod genome, while also displaying many adaptations to its specific life history.This work was supported by the following grants: NHGRIU54HG003273 to R.A.G; EU Marie Curie ITN #215781 “Evonet” to M.A.; a Wellcome Trust Value in People (VIP) award to C.B. and Wellcome Trust graduate studentship WT089615MA to J.E.G; Marine rhythms of Life” of the University of Vienna, an FWF (http://www.fwf.ac.at/) START award (#AY0041321) and HFSP (http://www.hfsp.org/) research grant (#RGY0082/2010) to KT-­‐R; MFPL Vienna International PostDoctoral Program for Molecular Life Sciences (funded by Austrian Ministry of Science and Research and City of Vienna, Cultural Department -­‐Science and Research to T.K; Direct Grant (4053034) of the Chinese University of Hong Kong to J.H.L.H.; NHGRI HG004164 to G.M.; Danish Research Agency (FNU), Carlsberg Foundation, and Lundbeck Foundation to C.J.P.G.; U.S. National Institutes of Health R01AI55624 to J.H.W.; Royal Society University Research fellowship to F.M.J.; P.D.E. was supported by the BBSRC via the Babraham Institute;This is the final version of the article. It first appeared from PLOS via http://dx.doi.org/10.1371/journal.pbio.100200

Pan-cancer analysis of whole genomes

Cancer is driven by genetic change, and the advent of massively parallel sequencing has enabled systematic documentation of this variation at the whole-genome scale(1-3). Here we report the integrative analysis of 2,658 whole-cancer genomes and their matching normal tissues across 38 tumour types from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). We describe the generation of the PCAWG resource, facilitated by international data sharing using compute clouds. On average, cancer genomes contained 4-5 driver mutations when combining coding and non-coding genomic elements; however, in around 5% of cases no drivers were identified, suggesting that cancer driver discovery is not yet complete. Chromothripsis, in which many clustered structural variants arise in a single catastrophic event, is frequently an early event in tumour evolution; in acral melanoma, for example, these events precede most somatic point mutations and affect several cancer-associated genes simultaneously. Cancers with abnormal telomere maintenance often originate from tissues with low replicative activity and show several mechanisms of preventing telomere attrition to critical levels. Common and rare germline variants affect patterns of somatic mutation, including point mutations, structural variants and somatic retrotransposition. A collection of papers from the PCAWG Consortium describes non-coding mutations that drive cancer beyond those in the TERT promoter(4); identifies new signatures of mutational processes that cause base substitutions, small insertions and deletions and structural variation(5,6); analyses timings and patterns of tumour evolution(7); describes the diverse transcriptional consequences of somatic mutation on splicing, expression levels, fusion genes and promoter activity(8,9); and evaluates a range of more-specialized features of cancer genomes(8,10-18).Peer reviewe

New genetic loci link adipose and insulin biology to body fat distribution.

Body fat distribution is a heritable trait and a well-established predictor of adverse metabolic outcomes, independent of overall adiposity. To increase our understanding of the genetic basis of body fat distribution and its molecular links to cardiometabolic traits, here we conduct genome-wide association meta-analyses of traits related to waist and hip circumferences in up to 224,459 individuals. We identify 49 loci (33 new) associated with waist-to-hip ratio adjusted for body mass index (BMI), and an additional 19 loci newly associated with related waist and hip circumference measures (P < 5 × 10(-8)). In total, 20 of the 49 waist-to-hip ratio adjusted for BMI loci show significant sexual dimorphism, 19 of which display a stronger effect in women. The identified loci were enriched for genes expressed in adipose tissue and for putative regulatory elements in adipocytes. Pathway analyses implicated adipogenesis, angiogenesis, transcriptional regulation and insulin resistance as processes affecting fat distribution, providing insight into potential pathophysiological mechanisms

Finishing the euchromatic sequence of the human genome

The sequence of the human genome encodes the genetic instructions for human physiology, as well as rich information about human evolution. In 2001, the International Human Genome Sequencing Consortium reported a draft sequence of the euchromatic portion of the human genome. Since then, the international collaboration has worked to convert this draft into a genome sequence with high accuracy and nearly complete coverage. Here, we report the result of this finishing process. The current genome sequence (Build 35) contains 2.85 billion nucleotides interrupted by only 341 gaps. It covers ∼99% of the euchromatic genome and is accurate to an error rate of ∼1 event per 100,000 bases. Many of the remaining euchromatic gaps are associated with segmental duplications and will require focused work with new methods. The near-complete sequence, the first for a vertebrate, greatly improves the precision of biological analyses of the human genome including studies of gene number, birth and death. Notably, the human enome seems to encode only 20,000-25,000 protein-coding genes. The genome sequence reported here should serve as a firm foundation for biomedical research in the decades ahead