Displaying publications 1 - 20 of 208 in total

Abstract:
Sort:
  1. Chatrchyan S, Khachatryan V, Sirunyan AM, Tumasyan A, Adam W, Bergauer T, et al.
    Phys Rev Lett, 2014 Apr 25;112(16):161802.
    PMID: 24815637
    Results are presented of a search for a "natural" supersymmetry scenario with gauge mediated symmetry breaking. It is assumed that only the supersymmetric partners of the top quark (the top squark) and the Higgs boson (Higgsino) are accessible. Events are examined in which there are two photons forming a Higgs boson candidate, and at least two b-quark jets. In 19.7  fb-1 of proton-proton collision data at s=8  TeV, recorded in the CMS experiment, no evidence of a signal is found and lower limits at the 95% confidence level are set, excluding the top squark mass below 360 to 410 GeV, depending on the Higgsino mass.
  2. Khachatryan V, Sirunyan AM, Tumasyan A, Adam W, Asilar E, Bergauer T, et al.
    Eur Phys J C Part Fields, 2017;77(1):15.
    PMID: 28260978 DOI: 10.1140/epjc/s10052-016-4504-z
    A measurement of the top quark pair production ([Formula: see text]) cross section in proton-proton collisions at the centre-of-mass energy of 8[Formula: see text] is presented using data collected with the CMS detector at the LHC, corresponding to an integrated luminosity of 19.6[Formula: see text]. This analysis is performed in the [Formula: see text] decay channels with one isolated, high transverse momentum electron or muon and at least four jets, at least one of which is required to be identified as originating from hadronization of a b quark. The calibration of the jet energy scale and the efficiency of b jet identification are determined from data. The measured [Formula: see text] cross section is [Formula: see text]. This measurement is compared with an analysis of 7[Formula: see text] data, corresponding to an integrated luminosity of 5.0[Formula: see text], to determine the ratio of 8[Formula: see text] to 7[Formula: see text] cross sections, which is found to be [Formula: see text]. The measurements are in agreement with QCD predictions up to next-to-next-to-leading order.
  3. Davis HR, Chan KO, Das I, Brennan IG, Karin BR, Jackman TR, et al.
    Mol Phylogenet Evol, 2020 06;147:106785.
    PMID: 32135306 DOI: 10.1016/j.ympev.2020.106785
    The gekkonid genus Cyrtodactylus is a highly diverse group of lizards (280 + species), which covers an expansive geographic range. Although this genus has been the focus of many taxonomic and molecular systematic studies, species on the Southeast Asian island of Borneo have remained understudied, leading to an unclear evolutionary history with cascading effects on taxonomy and biogeographic inferences. We assembled the most comprehensive multilocus Bornean dataset (one mitochondrial and three nuclear loci) that included 129 novel sequences and representatives from each known Cyrtodactylus species on the island to validate taxonomic status, assess species diversity, and elucidate biogeographic patterns. Our results uncovered a high proportion of cryptic diversity and revealed numerous taxonomic complications, especially within the C. consobrinus, C. malayanus, and C. pubisulcus groups. Comparisons of pairwise genetic distances and a preliminary species delimitation analysis using the Automatic Barcode Gap Discovery (ABGD) method demonstrated that some wide-ranging species on Borneo likely comprise multiple distinct and deeply divergent lineages, each with more restricted distributional ranges. We also tested the prevailing biogeographic hypothesis of a single invasion from Borneo into the Philippines. Our analyses revealed that Philippine taxa were not monophyletic, but were likely derived from multiple separate invasions into the geopolitical areas comprising the Philippines. Although our investigation of Bornean Cyrtodactylus is the most comprehensive to-date, it highlights the need for expanded taxonomic sampling and suggests that our knowledge of the evolutionary history, systematics, and biogeography of Bornean Cyrtodactylus is far from complete.
  4. Chan KO, Grismer LL, Brown RM
    Zootaxa, 2014;3900(4):569-80.
    PMID: 25543757 DOI: 10.11646/zootaxa.3900.4.7
    A new species of Narrow-mouthed frog of the genus Kaloula is described from northern Peninsular Malaysia. Kaloula latidisca sp. nov. is genetically and morphologically most similar to K. baleata and K. indochinensis but differs from those and other congeners by the unique combination of the following characters: (1) adult males SVL 49.2-56.2 mm (x̅=53.5 ± 3.0; N=4); (2) finger tips expanded into large, transversely expanded discs (disc width 2.8-3.1 mm, x̅=3.0 ± 0.1); (3) inner metatarsal tubercle large, oval, distinctly raised, slightly shorter than first toe; (4) three subarticular tubercles on fourth toe; (5) toe webbing formula: I 1-2 II 1-3 III 2-3.5 IV 4-2 V; and (6) yellow to orange irregularly shaped patch on the axillary, inguinal and posterior region of thigh.
  5. Chan KO, Hutter CR, Wood PL, Grismer LL, Das I, Brown RM
    Mol Ecol, 2020 10;29(20):3970-3987.
    PMID: 32808335 DOI: 10.1111/mec.15603
    Most new cryptic species are described using conventional tree- and distance-based species delimitation methods (SDMs), which rely on phylogenetic arrangements and measures of genetic divergence. However, although numerous factors such as population structure and gene flow are known to confound phylogenetic inference and species delimitation, the influence of these processes is not frequently evaluated. Using large numbers of exons, introns, and ultraconserved elements obtained using the FrogCap sequence-capture protocol, we compared conventional SDMs with more robust genomic analyses that assess population structure and gene flow to characterize species boundaries in a Southeast Asian frog complex (Pulchrana picturata). Our results showed that gene flow and introgression can produce phylogenetic patterns and levels of divergence that resemble distinct species (up to 10% divergence in mitochondrial DNA). Hybrid populations were inferred as independent (singleton) clades that were highly divergent from adjacent populations (7%-10%) and unusually similar (<3%) to allopatric populations. Such anomalous patterns are not uncommon in Southeast Asian amphibians, which brings into question whether the high levels of cryptic diversity observed in other amphibian groups reflect distinct cryptic species-or, instead, highly admixed and structured metapopulation lineages. Our results also provide an alternative explanation to the conundrum of divergent (sometimes nonsister) sympatric lineages-a pattern that has been celebrated as indicative of true cryptic speciation. Based on these findings, we recommend that species delimitation of continuously distributed "cryptic" groups should not rely solely on conventional SDMs, but should necessarily examine population structure and gene flow to avoid taxonomic inflation.
  6. Chan KO, Grismer LL, Brown RM
    Mol Phylogenet Evol, 2018 10;127:1010-1019.
    PMID: 30030179 DOI: 10.1016/j.ympev.2018.07.005
    The family Rhacophoridae is one of the most diverse amphibian families in Asia, for which taxonomic understanding is rapidly-expanding, with new species being described steadily, and at increasingly finer genetic resolution. Distance-based methods frequently have been used to justify or at least to bolster the recognition of new species, particularly in complexes of "cryptic" species where obvious morphological differentiation does not accompany speciation. However, there is no universally-accepted threshold to distinguish intra- from interspecific genetic divergence. Moreover, indiscriminant use of divergence thresholds to delimit species can result in over- or underestimation of species diversity. To explore the range of variation in application of divergence scales, and to provide a family-wide assessment of species-level diversity in Old-World treefrogs (family Rhacophoridae), we assembled the most comprehensive multi-locus phylogeny to date, including all 18 genera and approximately 247 described species (∼60% coverage). We then used the Automatic Barcode Gap Discovery (ABGD) method to obtain different species-delimitation schemes over a range of prior intraspecific divergence limits to assess the consistency of divergence thresholds used to demarcate current species boundaries. The species-rich phylogeny was able to identify a number of taxonomic errors, namely the incorrect generic placement of Chiromantis inexpectatus, which we now move to the genus Feihyla, and the specific identity of Rhacophorus bipunctatus from Peninsular Malaysia, which we tentatively reassign to R. rhodopus. The ABGD analysis demonstrated overlap between intra- and interspecific divergence limits: genetic thresholds used in some studies to synonymize taxa have frequently been used in other studies to justify the recognition of new species. This analysis also highlighted numerous groups that could potentially be split or lumped, which we earmark for future examination. Our large-scale and en bloc approach to species-level phylogenetic systematics contributes to the resolution of taxonomic uncertainties, reveals possible new species, and identifies numerous groups that require critical examination. Overall, we demonstrate that the taxonomy and evolutionary history of Old-World tree frogs are far from resolved, stable or adequately characterized at the level of genus, species, and/or population.
  7. Grismer LL, Wood PLJ, Grismer JL, Quah ESH, Thy N, Phimmachak S, et al.
    Zootaxa, 2019 Jul 16;4638(2):zootaxa.4638.2.1.
    PMID: 31712473 DOI: 10.11646/zootaxa.4638.2.1
    An integrative taxonomic analysis of the Ptychozoon lionotum group across its range in Indochina and Sundaland recovers P. lionotum sensu lato Annandale, 1905 as paraphyletic with respect to P. popaense Grismer, Wood, Thura, Grismer, Brown, Stuart, 2018a and composed of four allopatric, genetically divergent, ND2 mitochondrial lineages. Multivariate and univariate analyses of continuous and discrete morphological and color pattern characters statistically and discretely diagnose each lineage from one another and together, with maximum likelihood and Bayesian inference analyses, provide the foundation for the recognition of each lineage as a new species-hypotheses corroborated with a Generalized Mixed Yule Coalescent species delimitation analysis. Ptychozoon cicakterbang sp. nov. ranges throughout Peninsular Malaysia to Pulau Natuna Besar, Indonesia; P. kabkaebin sp. nov. is endemic to northern and central Laos; and P. tokehos sp. nov. ranges from southern Thailand south of the Isthmus of Kra northward to Chiang Mai, fringing the Chao Phraya Basin and ranging southward through Cambodia to southern Vietnam. Ptychozoon lionotum sensu stricto ranges from northwestern Laos through southern Myanmar to eastern India. The phylogeographic structure within each species varies considerably with P. lionotum s.s. showing no genetic divergence across its 1,100 km range compared to P. cicakterbang sp. nov. showing upwards of 8.2% sequence divergence between syntopic individuals. Significant phylogeographic structure exists within P. tokehos sp. nov. and increased sampling throughout Thailand may require additional taxonomic changes within this species.
  8. Chan KO, Alexander AM, Grismer LL, Su YC, Grismer JL, Quah ESH, et al.
    Mol Ecol, 2017 Oct;26(20):5435-5450.
    PMID: 28802073 DOI: 10.1111/mec.14296
    Accurately delimiting species boundaries is a nontrivial undertaking that can have significant effects on downstream inferences. We compared the efficacy of commonly used species delimitation methods (SDMs) and a population genomics approach based on genomewide single-nucleotide polymorphisms (SNPs) to assess lineage separation in the Malaysian Torrent Frog Complex currently recognized as a single species (Amolops larutensis). First, we used morphological, mitochondrial DNA and genomewide SNPs to identify putative species boundaries by implementing noncoalescent and coalescent-based SDMs (mPTP, iBPP, BFD*). We then tested the validity of putative boundaries by estimating spatiotemporal gene flow (fastsimcoal2, ABBA-BABA) to assess the extent of genetic isolation among putative species. Our results show that the A. larutensis complex runs the gamut of the speciation continuum from highly divergent, genetically isolated lineages (mean Fst  = 0.9) to differentiating populations involving recent gene flow (mean Fst  = 0.05; Nm  > 5). As expected, SDMs were effective at delimiting divergent lineages in the absence of gene flow but overestimated species in the presence of marked population structure and gene flow. However, using a population genomics approach and the concept of species as separately evolving metapopulation lineages as the only necessary property of a species, we were able to objectively elucidate cryptic species boundaries in the presence of past and present gene flow. This study does not discount the utility of SDMs but highlights the danger of violating model assumptions and the importance of carefully considering methods that appropriately fit the diversification history of a particular system.
  9. Khachatryan V, Sirunyan AM, Tumasyan A, Adam W, Bergauer T, Dragicevic M, et al.
    Phys Rev Lett, 2015 Feb 13;114(6):061801.
    PMID: 25723204
    A search for new long-lived particles decaying to leptons is presented using proton-proton collisions produced by the LHC at √[s]=8  TeV. Data used for the analysis were collected by the CMS detector and correspond to an integrated luminosity of 19.7  fb(-1). Events are selected with an electron and muon with opposite charges that both have transverse impact parameter values between 0.02 and 2 cm. The search has been designed to be sensitive to a wide range of models with nonprompt e-μ final states. Limits are set on the "displaced supersymmetry" model, with pair production of top squarks decaying into an e-μ final state via R-parity-violating interactions. The results are the most restrictive to date on this model, with the most stringent limit being obtained for a top squark lifetime corresponding to cτ=2  cm, excluding masses below 790 GeV at 95% confidence level.
  10. Khachatryan V, Sirunyan AM, Tumasyan A, Adam W, Bergauer T, Dragicevic M, et al.
    Phys Rev Lett, 2015 Feb 6;114(5):051801.
    PMID: 25699433
    A study of vector boson scattering in pp collisions at a center-of-mass energy of 8 TeV is presented. The data sample corresponds to an integrated luminosity of 19.4  fb(-1) collected with the CMS detector. Candidate events are selected with exactly two leptons of the same charge, two jets with large rapidity separation and high dijet mass, and moderate missing transverse energy. The signal region is expected to be dominated by electroweak same-sign W-boson pair production. The observation agrees with the standard model prediction. The observed significance is 2.0 standard deviations, where a significance of 3.1 standard deviations is expected based on the standard model. Cross section measurements for W(±)W(±) and WZ processes in the fiducial region are reported. Bounds on the structure of quartic vector-boson interactions are given in the framework of dimension-eight effective field theory operators, as well as limits on the production of doubly charged Higgs bosons.
  11. Khachatryan V, Sirunyan AM, Tumasyan A, Adam W, Asilar E, Bergauer T, et al.
    Phys Rev Lett, 2016 Feb 19;116(7):071801.
    PMID: 26943527 DOI: 10.1103/PhysRevLett.116.071801
    A search for narrow resonances in proton-proton collisions at sqrt[s]=13  TeV is presented. The invariant mass distribution of the two leading jets is measured with the CMS detector using a data set corresponding to an integrated luminosity of 2.4  fb^{-1}. The highest observed dijet mass is 6.1 TeV. The distribution is smooth and no evidence for resonant particles is observed. Upper limits at 95% confidence level are set on the production cross section for narrow resonances with masses above 1.5 TeV. When interpreted in the context of specific models, the limits exclude string resonances with masses below 7.0 TeV, scalar diquarks below 6.0 TeV, axigluons and colorons below 5.1 TeV, excited quarks below 5.0 TeV, color-octet scalars below 3.1 TeV, and W^{'} bosons below 2.6 TeV. These results significantly extend previously published limits.
  12. Khachatryan V, Sirunyan AM, Tumasyan A, Adam W, Asilar E, Bergauer T, et al.
    Phys Rev Lett, 2016 Jan 22;116(3):032301.
    PMID: 26849587 DOI: 10.1103/PhysRevLett.116.032301
    The production cross sections of the B^{+}, B^{0}, and B_{s}^{0} mesons, and of their charge conjugates, are measured via exclusive hadronic decays in p+Pb collisions at the center-of-mass energy sqrt[s_{NN}]=5.02  TeV with the CMS detector at the CERN LHC. The data set used for this analysis corresponds to an integrated luminosity of 34.6  nb^{-1}. The production cross sections are measured in the transverse momentum range between 10 and 60  GeV/c. No significant modification is observed compared to proton-proton perturbative QCD calculations scaled by the number of incoherent nucleon-nucleon collisions. These results provide a baseline for the study of in-medium b quark energy loss in Pb+Pb collisions.
  13. Khachatryan V, Sirunyan AM, Tumasyan A, Adam W, Bergauer T, Dragicevic M, et al.
    Phys Rev Lett, 2015 Mar 13;114(10):101801.
    PMID: 25815923
    Results are presented from a search for new decaying massive particles whose presence is inferred from an imbalance in transverse momentum and which are produced in association with a single top quark that decays into a bottom quark and two light quarks. The measurement is performed using 19.7  fb^{-1} of data from proton-proton collisions at a center-of-mass energy of 8 TeV, collected with the CMS detector at the CERN LHC. No deviations from the standard model predictions are observed and lower limits are set on the masses of new invisible bosons. In particular, scalar and vector particles, with masses below 330 and 650 GeV, respectively, are excluded at 95% confidence level, thereby substantially extending a previous limit published by the CDF Collaboration.
  14. Sirunyan AM, Tumasyan A, Adam W, Asilar E, Bergauer T, Brandstetter J, et al.
    Phys Rev Lett, 2018 Apr 06;120(14):142301.
    PMID: 29694144 DOI: 10.1103/PhysRevLett.120.142301
    The relative yields of ϒ mesons produced in pp and Pb-Pb collisions at sqrt[s_{NN}]=5.02  TeV and reconstructed via the dimuon decay channel are measured using data collected by the CMS experiment. Double ratios are formed by comparing the yields of the excited states, ϒ(2S) and ϒ(3S), to the ground state, ϒ(1S), in both Pb-Pb and pp collisions at the same center-of-mass energy. The double ratios, [ϒ(nS)/ϒ(1S)]_{Pb-Pb}/[ϒ(nS)/ϒ(1S)]_{pp}, are measured to be 0.308±0.055(stat)±0.019(syst) for the ϒ(2S) and less than 0.26 at 95% confidence level for the ϒ(3S). No significant ϒ(3S) signal is found in the Pb-Pb data. The double ratios are studied as a function of collision centrality, as well as ϒ transverse momentum and rapidity. No significant dependencies are observed.
  15. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al.
    Phys Rev Lett, 2020 Sep 04;125(10):102001.
    PMID: 32955327 DOI: 10.1103/PhysRevLett.125.102001
    The first study of charm quark diffusion with respect to the jet axis in heavy ion collisions is presented. The measurement is performed using jets with p_{T}^{jet}>60  GeV/c and D^{0} mesons with p_{T}^{D}>4  GeV/c in lead-lead (Pb-Pb) and proton-proton (pp) collisions at a nucleon-nucleon center-of-mass energy of sqrt[s_{NN}]=5.02  TeV, recorded by the CMS detector at the LHC. The radial distribution of D^{0} mesons with respect to the jet axis is sensitive to the production mechanisms of the meson, as well as to the energy loss and diffusion processes undergone by its parent parton inside the strongly interacting medium produced in Pb-Pb collisions. When compared to Monte Carlo event generators, the radial distribution in pp collisions is found to be well described by pythia, while the slope of the distribution predicted by sherpa is steeper than that of the data. In Pb-Pb collisions, compared to the pp results, the D^{0} meson distribution for 4
  16. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Bergauer T, Dragicevic M, et al.
    Phys Rev Lett, 2020 Oct 09;125(15):152001.
    PMID: 33095627 DOI: 10.1103/PhysRevLett.125.152001
    Using a data sample of proton-proton collisions at sqrt[s]=13  TeV, corresponding to an integrated luminosity of 140  fb^{-1} collected by the CMS experiment in 2016-2018, the B_{s}^{0}→X(3872)ϕ decay is observed. Decays into J/ψπ^{+}π^{-} and K^{+}K^{-} are used to reconstruct, respectively, the X(3872) and ϕ. The ratio of the product of branching fractions B[B_{s}^{0}→X(3872)ϕ]B[X(3872)→J/ψπ^{+}π^{-}] to the product B[B_{s}^{0}→ψ(2S)ϕ]B[ψ(2S)→J/ψπ^{+}π^{-}] is measured to be [2.21±0.29(stat)±0.17(syst)]%. The ratio B[B_{s}^{0}→X(3872)ϕ]/B[B^{0}→X(3872)K^{0}] is found to be consistent with one, while the ratio B[B_{s}^{0}→X(3872)ϕ]/B[B^{+}→X(3872)K^{+}] is two times smaller. This suggests a difference in the production dynamics of the X(3872) in B^{0} and B_{s}^{0} meson decays compared to B^{+}. The reported observation may shed new light on the nature of the X(3872) particle.
  17. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al.
    Phys Rev Lett, 2018 May 18;120(20):202005.
    PMID: 29864318 DOI: 10.1103/PhysRevLett.120.202005
    A search for resonancelike structures in the B_{s}^{0}π^{±} invariant mass spectrum is performed using proton-proton collision data collected by the CMS experiment at the LHC at sqrt[s]=8  TeV, corresponding to an integrated luminosity of 19.7  fb^{-1}. The B_{s}^{0} mesons are reconstructed in the decay chain B_{s}^{0}→J/ψϕ, with J/ψ→μ^{+}μ^{-} and ϕ→K^{+}K^{-}. The B_{s}^{0}π^{±} invariant mass distribution shows no statistically significant peaks for different selection requirements on the reconstructed B_{s}^{0} and π^{±} candidates. Upper limits are set on the relative production rates of the X(5568) and B_{s}^{0} states times the branching fraction of the decay X(5568)^{±}→B_{s}^{0}π^{±}. In addition, upper limits are obtained as a function of the mass and the natural width of possible exotic states decaying into B_{s}^{0}π^{±}.
  18. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al.
    Phys Rev Lett, 2018 May 18;120(20):201801.
    PMID: 29864370 DOI: 10.1103/PhysRevLett.120.201801
    A search for narrow resonances decaying to bottom quark-antiquark pairs is presented, using a data sample of proton-proton collisions at sqrt[s]=8  TeV corresponding to an integrated luminosity of 19.7  fb^{-1}. The search is extended to masses lower than those reached in typical searches for resonances decaying into jet pairs at the LHC, by taking advantage of triggers that identify jets originating from bottom quarks. No significant excess of events is observed above the background predictions. Limits are set on the product of cross section and branching fraction to bottom quarks for spin 0, 1, and 2 resonances in the mass range of 325-1200 GeV. These results improve on the limits for resonances decaying into jet pairs in the 325-500 GeV mass range.
  19. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al.
    Phys Rev Lett, 2018 May 18;120(20):202301.
    PMID: 29864330 DOI: 10.1103/PhysRevLett.120.202301
    The prompt D^{0} meson azimuthal anisotropy coefficients, v_{2} and v_{3}, are measured at midrapidity (|y|<1.0) in Pb-Pb collisions at a center-of-mass energy sqrt[s_{NN}]=5.02  TeV per nucleon pair with data collected by the CMS experiment. The measurement is performed in the transverse momentum (p_{T}) range of 1 to 40  GeV/c, for central and midcentral collisions. The v_{2} coefficient is found to be positive throughout the p_{T} range studied. The first measurement of the prompt D^{0} meson v_{3} coefficient is performed, and values up to 0.07 are observed for p_{T} around 4  GeV/c. Compared to measurements of charged particles, a similar p_{T} dependence, but smaller magnitude for p_{T}<6  GeV/c, is found for prompt D^{0} meson v_{2} and v_{3} coefficients. The results are consistent with the presence of collective motion of charm quarks at low p_{T} and a path length dependence of charm quark energy loss at high p_{T}, thereby providing new constraints on the theoretical description of the interactions between charm quarks and the quark-gluon plasma.
  20. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al.
    Phys Rev Lett, 2018 Dec 14;121(24):241802.
    PMID: 30608761 DOI: 10.1103/PhysRevLett.121.241802
    Three of the most significant measured deviations from standard model predictions, the enhanced decay rate for B→D^{(*)}τν, hints of lepton universality violation in B→K^{(*)}ℓℓ decays, and the anomalous magnetic moment of the muon, can be explained by the existence of leptoquarks (LQs) with large couplings to third-generation quarks and masses at the TeV scale. The existence of these states can be probed at the LHC in high energy proton-proton collisions. A novel search is presented for pair production of LQs coupled to a top quark and a muon using data at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 35.9  fb^{-1}, recorded by the CMS experiment. No deviation from the standard model prediction has been observed and scalar LQs decaying exclusively into tμ are excluded up to masses of 1420 GeV. The results of this search are combined with those from previous searches for LQ decays into tτ and bν, which excluded scalar LQs below masses of 900 and 1080 GeV. Vector LQs are excluded up to masses of 1190 GeV for all possible combinations of branching fractions to tμ, tτ and bν. With this analysis, all relevant couplings of LQs with an electric charge of -1/3 to third-generation quarks are probed for the first time.
Filters
Contact Us

Please provide feedback to Administrator (afdal@afpm.org.my)

External Links