Sparsest factor analysis for clustering variables: a matrix decomposition approach

We propose a new procedure for sparse factor analysis (FA) such that each variable loads only one common factor. Thus, the loading matrix has a single nonzero element in each row and zeros elsewhere. Such a loading matrix is the sparsest possible for certain number of variables and common factors. For this reason, the proposed method is named sparsest FA (SSFA). It may also be called FA-based variable clustering, since the variables loading the same common factor can be classified into a cluster. In SSFA, all model parts of FA (common factors, their correlations, loadings, unique factors, and unique variances) are treated as fixed unknown parameter matrices and their least squares function is minimized through specific data matrix decomposition. A useful feature of the algorithm is that the matrix of common factor scores is re-parameterized using QR decomposition in order to efficiently estimate factor correlations. A simulation study shows that the proposed procedure can exactly identify the true sparsest models. Real data examples demonstrate the usefulness of the variable clustering performed by SSFA

Insights from Amphioxus into the Evolution of Vertebrate Cartilage

Central to the story of vertebrate evolution is the origin of the vertebrate head, a problem difficult to approach using paleontology and comparative morphology due to a lack of unambiguous intermediate forms. Embryologically, much of the vertebrate head is derived from two ectodermal tissues, the neural crest and cranial placodes. Recent work in protochordates suggests the first chordates possessed migratory neural tube cells with some features of neural crest cells. However, it is unclear how and when these cells acquired the ability to form cellular cartilage, a cell type unique to vertebrates. It has been variously proposed that the neural crest acquired chondrogenic ability by recruiting proto-chondrogenic gene programs deployed in the neural tube, pharynx, and notochord. To test these hypotheses we examined the expression of 11 amphioxus orthologs of genes involved in neural crest chondrogenesis. Consistent with cellular cartilage as a vertebrate novelty, we find that no single amphioxus tissue co-expresses all or most of these genes. However, most are variously co-expressed in mesodermal derivatives. Our results suggest that neural crest-derived cartilage evolved by serial cooption of genes which functioned primitively in mesoderm

A flexible framework for sparse simultaneous component based data integration

<p>Abstract</p> <p>1 Background</p> <p>High throughput data are complex and methods that reveal structure underlying the data are most useful. Principal component analysis, frequently implemented as a singular value decomposition, is a popular technique in this respect. Nowadays often the challenge is to reveal structure in several sources of information (e.g., transcriptomics, proteomics) that are available for the same biological entities under study. Simultaneous component methods are most promising in this respect. However, the interpretation of the principal and simultaneous components is often daunting because contributions of each of the biomolecules (transcripts, proteins) have to be taken into account.</p> <p>2 Results</p> <p>We propose a sparse simultaneous component method that makes many of the parameters redundant by shrinking them to zero. It includes principal component analysis, sparse principal component analysis, and ordinary simultaneous component analysis as special cases. Several penalties can be tuned that account in different ways for the block structure present in the integrated data. This yields known sparse approaches as the lasso, the ridge penalty, the elastic net, the group lasso, sparse group lasso, and elitist lasso. In addition, the algorithmic results can be easily transposed to the context of regression. Metabolomics data obtained with two measurement platforms for the same set of <it>Escherichia coli </it>samples are used to illustrate the proposed methodology and the properties of different penalties with respect to sparseness across and within data blocks.</p> <p>3 Conclusion</p> <p>Sparse simultaneous component analysis is a useful method for data integration: First, simultaneous analyses of multiple blocks offer advantages over sequential and separate analyses and second, interpretation of the results is highly facilitated by their sparseness. The approach offered is flexible and allows to take the block structure in different ways into account. As such, structures can be found that are exclusively tied to one data platform (group lasso approach) as well as structures that involve all data platforms (Elitist lasso approach).</p> <p>4 Availability</p> <p>The additional file contains a MATLAB implementation of the sparse simultaneous component method.</p

Sequential and Coordinated Actions of c-Myc and N-Myc Control Appendicular Skeletal Development

BACKGROUND: During limb development, chondrocytes and osteoblasts emerge from condensations of limb bud mesenchyme. These cells then proliferate and differentiate in separate but adjacent compartments and function cooperatively to promote bone growth through the process of endochondral ossification. While many aspects of limb skeletal formation are understood, little is known about the mechanisms that link the development of undifferentiated limb bud mesenchyme with formation of the precartilaginous condensation and subsequent proliferative expansion of chondrocyte and osteoblast lineages. The aim of this study was to gain insight into these processes by examining the roles of c-Myc and N-Myc in morphogenesis of the limb skeleton. METHODOLOGY/PRINCIPAL FINDINGS: To investigate c-Myc function in skeletal development, we characterized mice in which floxed c-Myc alleles were deleted in undifferentiated limb bud mesenchyme with Prx1-Cre, in chondro-osteoprogenitors with Sox9-Cre and in osteoblasts with Osx1-Cre. We show that c-Myc promotes the proliferative expansion of both chondrocytes and osteoblasts and as a consequence controls the process of endochondral growth and ossification and determines bone size. The control of proliferation by c-Myc was related to its effects on global gene transcription, as phosphorylation of the C-Terminal Domain (pCTD) of RNA Polymerase II, a marker of general transcription initiation, was tightly coupled to cell proliferation of growth plate chondrocytes where c-Myc is expressed and severely downregulated in the absence of c-Myc. Finally, we show that combined deletion of N-Myc and c-Myc in early limb bud mesenchyme gives rise to a severely hypoplastic limb skeleton that exhibits features characteristic of individual c-Myc and N-Myc mutants. CONCLUSIONS/SIGNIFICANCE: Our results show that N-Myc and c-Myc act sequentially during limb development to coordinate the expansion of key progenitor populations responsible for forming the limb skeleton

Jet energy measurement with the ATLAS detector in proton-proton collisions at root s=7 TeV

The jet energy scale and its systematic uncertainty are determined for jets measured with the ATLAS detector at the LHC in proton-proton collision data at a centre-of-mass energy of √s = 7TeV corresponding to an integrated luminosity of 38 pb-1. Jets are reconstructed with the anti-kt algorithm with distance parameters R=0. 4 or R=0. 6. Jet energy and angle corrections are determined from Monte Carlo simulations to calibrate jets with transverse momenta pT≥20 GeV and pseudorapidities {pipe}η{pipe}<4. 5. The jet energy systematic uncertainty is estimated using the single isolated hadron response measured in situ and in test-beams, exploiting the transverse momentum balance between central and forward jets in events with dijet topologies and studying systematic variations in Monte Carlo simulations. The jet energy uncertainty is less than 2. 5 % in the central calorimeter region ({pipe}η{pipe}<0. 8) for jets with 60≤pT<800 GeV, and is maximally 14 % for pT<30 GeV in the most forward region 3. 2≤{pipe}η{pipe}<4. 5. The jet energy is validated for jet transverse momenta up to 1 TeV to the level of a few percent using several in situ techniques by comparing a well-known reference such as the recoiling photon pT, the sum of the transverse momenta of tracks associated to the jet, or a system of low-pT jets recoiling against a high-pT jet. More sophisticated jet calibration schemes are presented based on calorimeter cell energy density weighting or hadronic properties of jets, aiming for an improved jet energy resolution and a reduced flavour dependence of the jet response. The systematic uncertainty of the jet energy determined from a combination of in situ techniques is consistent with the one derived from single hadron response measurements over a wide kinematic range. The nominal corrections and uncertainties are derived for isolated jets in an inclusive sample of high-pT jets. Special cases such as event topologies with close-by jets, or selections of samples with an enhanced content of jets originating from light quarks, heavy quarks or gluons are also discussed and the corresponding uncertainties are determined. © 2013 CERN for the benefit of the ATLAS collaboration

Measurement of the inclusive and dijet cross-sections of b-jets in pp collisions at sqrt(s) = 7 TeV with the ATLAS detector

The inclusive and dijet production cross-sections have been measured for jets containing b-hadrons (b-jets) in proton-proton collisions at a centre-of-mass energy of sqrt(s) = 7 TeV, using the ATLAS detector at the LHC. The measurements use data corresponding to an integrated luminosity of 34 pb^-1. The b-jets are identified using either a lifetime-based method, where secondary decay vertices of b-hadrons in jets are reconstructed using information from the tracking detectors, or a muon-based method where the presence of a muon is used to identify semileptonic decays of b-hadrons inside jets. The inclusive b-jet cross-section is measured as a function of transverse momentum in the range 20 < pT < 400 GeV and rapidity in the range |y| < 2.1. The bbbar-dijet cross-section is measured as a function of the dijet invariant mass in the range 110 < m_jj < 760 GeV, the azimuthal angle difference between the two jets and the angular variable chi in two dijet mass regions. The results are compared with next-to-leading-order QCD predictions. Good agreement is observed between the measured cross-sections and the predictions obtained using POWHEG + Pythia. MC@NLO + Herwig shows good agreement with the measured bbbar-dijet cross-section. However, it does not reproduce the measured inclusive cross-section well, particularly for central b-jets with large transverse momenta.Comment: 10 pages plus author list (21 pages total), 8 figures, 1 table, final version published in European Physical Journal

Search for new phenomena in final states with an energetic jet and large missing transverse momentum in pp collisions at √ s = 8 TeV with the ATLAS detector

Results of a search for new phenomena in final states with an energetic jet and large missing transverse momentum are reported. The search uses 20.3 fb−1 of √ s = 8 TeV data collected in 2012 with the ATLAS detector at the LHC. Events are required to have at least one jet with pT > 120 GeV and no leptons. Nine signal regions are considered with increasing missing transverse momentum requirements between Emiss T > 150 GeV and Emiss T > 700 GeV. Good agreement is observed between the number of events in data and Standard Model expectations. The results are translated into exclusion limits on models with either large extra spatial dimensions, pair production of weakly interacting dark matter candidates, or production of very light gravitinos in a gauge-mediated supersymmetric model. In addition, limits on the production of an invisibly decaying Higgs-like boson leading to similar topologies in the final state are presente

Observation of associated near-side and away-side long-range correlations in √sNN=5.02 TeV proton-lead collisions with the ATLAS detector

Two-particle correlations in relative azimuthal angle (Δϕ) and pseudorapidity (Δη) are measured in √sNN=5.02  TeV p+Pb collisions using the ATLAS detector at the LHC. The measurements are performed using approximately 1  μb-1 of data as a function of transverse momentum (pT) and the transverse energy (ΣETPb) summed over 3.1<η<4.9 in the direction of the Pb beam. The correlation function, constructed from charged particles, exhibits a long-range (2<|Δη|<5) “near-side” (Δϕ∼0) correlation that grows rapidly with increasing ΣETPb. A long-range “away-side” (Δϕ∼π) correlation, obtained by subtracting the expected contributions from recoiling dijets and other sources estimated using events with small ΣETPb, is found to match the near-side correlation in magnitude, shape (in Δη and Δϕ) and ΣETPb dependence. The resultant Δϕ correlation is approximately symmetric about π/2, and is consistent with a dominant cos⁡2Δϕ modulation for all ΣETPb ranges and particle pT

A Late Role for bmp2b in the Morphogenesis of Semicircular Canal Ducts in the Zebrafish Inner Ear

BACKGROUND:The Bone Morphogenetic Protein (BMP) genes bmp2 and bmp4 are expressed in highly conserved patterns in the developing vertebrate inner ear. It has, however, proved difficult to elucidate the function of BMPs during ear development as mutations in these genes cause early embryonic lethality. Previous studies using conditional approaches in mouse and chicken have shown that Bmp4 has a role in semicircular canal and crista development, but there is currently no direct evidence for the role of Bmp2 in the developing inner ear. METHODOLOGY/PRINCIPAL FINDINGS:We have used an RNA rescue strategy to test the role of bmp2b in the zebrafish inner ear directly. Injection of bmp2b or smad5 mRNA into homozygous mutant swirl (bmp2b(-/-)) embryos rescues the early patterning defects in these mutants and the fish survive to adulthood. As injected RNA will only last, at most, for the first few days of embryogenesis, all later development occurs in the absence of bmp2b function. Although rescued swirl adult fish are viable, they have balance defects suggestive of vestibular dysfunction. Analysis of the inner ears of these fish reveals a total absence of semicircular canal ducts, structures involved in the detection of angular motion. All other regions of the ear, including the ampullae and cristae, are present and appear normal. Early stages of otic development in rescued swirl embryos are also normal. CONCLUSIONS/SIGNIFICANCE:Our findings demonstrate a critical late role for bmp2b in the morphogenesis of semicircular canals in the zebrafish inner ear. This is the first demonstration of a developmental role for any gene during post-embryonic stages of otic morphogenesis in the zebrafish. Despite differences in the early stages of semicircular canal formation between zebrafish and amniotes, the role of Bmp2 in semicircular canal duct outgrowth is likely to be conserved between different vertebrate species

Measurement of the cross-section of high transverse momentum vector bosons reconstructed as single jets and studies of jet substructure in pp collisions at √s = 7 TeV with the ATLAS detector

This paper presents a measurement of the cross-section for high transverse momentum W and Z bosons produced in pp collisions and decaying to all-hadronic final states. The data used in the analysis were recorded by the ATLAS detector at the CERN Large Hadron Collider at a centre-of-mass energy of √s = 7 TeV;{\rm Te}{\rm V}$and correspond to an integrated luminosity of$4.6\;{\rm f}{{{\rm b}}^{-1}}$. The measurement is performed by reconstructing the boosted W or Z bosons in single jets. The reconstructed jet mass is used to identify the W and Z bosons, and a jet substructure method based on energy cluster information in the jet centre-of-mass frame is used to suppress the large multi-jet background. The cross-section for events with a hadronically decaying W or Z boson, with transverse momentum${{p}_{{\rm T}}}\gt 320\;{\rm Ge}{\rm V}$and pseudorapidity$|\eta |\lt 1.9$, is measured to be${{\sigma }_{W+Z}}=8.5\pm 1.7$ pb and is compared to next-to-leading-order calculations. The selected events are further used to study jet grooming techniques