Here we briefly review possible indirect effects of dark matter (DM) of the Universe. It includes effects in cosmic rays (CR): first of all, the positron excess at $\sim$ 500 GeV and possible electron-positron excess at 1-1.5 TeV. We tell that the main and least model-dependent constraint on such possible interpretation of CR effects goes from gamma-ray background. Even ordinary $e^+e^-$ mode of DM decay or annihilation produces prompt photons (FSR) so much that leads to contradiction with data on cosmic gamma-rays. We present our attempts to possibly avoid gamma-ray constraint. They concern peculiarities of both space distribution of DM and their physics. The latter involves complications of decay/annihilation modes of DM, modifications of Lagrangian of DM-ordinary matter interaction, and inclusion of mode with identical fermions in final state. In this way, no possibilities to suppress were found except, possibly, mode with identical fermions. While the case of spatial distribution variation allows achieving consistency between different data. Also we consider stable form of dark matter which can interact with baryons. We show which constraint such DM candidate can get from damping effect in plasma during large scale structure formation in comparison with other existing constraints.

The Primordial Black Holes (PBHs) are gradually involved into consideration as the phenomenon having reliable basis. We discuss here the possibility of their agglomeration into clusters that may have several prominent observable features. The clusters can form due to closed domain walls appearance in the natural and the hybrid inflation with subsequent evolution and gravitational collapse. Early dustlike stages of dominance of heavy metastable dissipative particles, at which star-like objects are formed, can also naturally lead to formation of black hole clusters, remaining in the Universe after decay of particles, from which they have originated. The dynamical evolution of such clusters discussed here is of the crucial importance. Such a model inherits all the advantages of the single PBHs like possible explanation of existence of supermassive black holes (origin of the early quasars), binary BH merges registered by LIGO/Virgo through gravitational waves, contribution to reionization of the Universe, but also has additional benefits. The cluster could alleviate or completely avoid existing constraints on the single PBH abundance making PBHs a real dark matter candidate. The most of existing constraints on (single) PBH density should be re-considered as applied to the clusters. Also unidentified cosmic gamma-ray point-like sources could be (partially) accounted for by them. One can conclude, that it seems really to be much more viable model with respect to the single PBHs.

Cold dark matter (DM) scenario may be cured of several problems by involving self-interaction of dark matter. Viability of the models of long-range interacting DM crucially depends on the effectiveness of recombination of the DM particles, making thereby their interaction short-range. Usually in numeric calculations, recombination is described by cross section obtained on a feasible quantum level. However in a wide range of parameter values, a classical treatment, where the particles are bound due to dipole radiation, is applicable. The cross sections, obtained in both approaches, are very different and lead to diverse consequences. Classical cross section has a steeper dependence on relative velocity, what leads to the fact that, after decoupling of DM particles from thermal background of "dark photons" (carriers of DM long-range interaction), recombination process does not "freeze out", diminishing gradually density of unbound DM particles. Our simplified estimates show, that at the taken parameter values (the mass of DM particle is $100$ GeV, interaction constant is $100^{-1}$, and quite natural assumptions on initial conditions, from which the result is very weakly dependent) the difference in residual density reaches about $5$ orders of magnitude on pre-galactic stage. This estimate takes into account thermal effects induced by dipole radiation and recombination, which resulted in the increase of both temperature and density of DM particles by a half order of magnitude.

The nonbaryonic dark matter of the Universe is assumed to consist of new stable forms of matter. Their stability reflects symmetry of micro world and mechanisms of its symmetry breaking. In the early Universe heavy metastable particles can dominate, leaving primordial black holes (PBHs) after their decay, as well as the structure of particle symmetry breaking gives rise to cosmological phase transitions, from which massive black holes and/or their clusters can originate. PBHs can be formed in such transitions within a narrow interval of masses about $10^{17}$g and, avoiding severe observational constraints on PBHs, can be a candidate for the dominant form of dark matter. PBHs in this range of mass can give solution of the problem of reionization in the Universe at the redshift $z\sim 5... 10$. Clusters of massive PBHs can serve as a nonlinear seeds for galaxy formation, while PBHs evaporating in such clusters can provide an interesting interpretation for the observations of point-like gamma-ray sources. Analysis of possible PBH signatures represents a universal probe for super-high energy physics in the early Universe in studies of indirect effects of the dark matter.