Population Size Expansion Research Articles

Empirical Analysis on the Mechanism of Industrial Park Driving Urban Expansion: A Case Study of Xining City

Taking Xining City as an example, this article analyzes the mechanism by which industrial park construction drives the expansion of urban population size and built-up area, based on a review of the process of industrial park development and urban population growth. It also discusses future urban governance models in light of urban development trends. The research finds: (1) In the process of urban development, industrial park construction is often the initial factor in the cumulative and cyclical development of a city; (2) As the level of development improves and the mode of economic growth changes, the government should timely adjust its strategies, shifting from the expansion of industrial park construction towards structural optimization and quality improvement. The most significant difference from previous research is that this paper emphasizes the importance of government planning. This study can not only demonstrate the general process of industrial parks promoting urban expansion, but more importantly, it explains the fundamental reasons for the transition of urban expansion to adjustment from a mechanism perspective, thereby eliminating the drawbacks of simply predicting urban scale evolution through data models.

Origin of the Laurentian Great Lakes fish fauna through upward adaptive radiation cascade prior to the Last Glacial Maximum

The evolutionary histories of adaptive radiations can be marked by dramatic demographic fluctuations. However, the demographic histories of ecologically-linked co-diversifying lineages remain understudied. The Laurentian Great Lakes provide a unique system of two such radiations that are dispersed across depth gradients with a predator-prey relationship. We show that the North American Coregonus species complex (“ciscoes”) radiated rapidly prior to the Last Glacial Maximum (80–90 kya), a globally warm period, followed by rapid expansion in population size. Similar patterns of demographic expansion were observed in the predator species, Lake Charr (Salvelinus namaycush), following a brief time lag, which we hypothesize was driven by predator-prey dynamics. Diversification of prey into deep water created ecological opportunities for the predators, facilitating their demographic expansion, which is consistent with an upward adaptive radiation cascade. This study provides a new timeline and environmental context for the origin of the Laurentian Great Lakes fish fauna, and firmly establishes this system as drivers of ecological diversification and rapid speciation through cyclical glaciation.

Digital Economy Mechanisms and Tests for Promoting Common Wealth

This paper focuses on the two development strategies of digital economy and common wealth, explores the effect, mechanism and spatial spillover effect of digital economy on common wealth, and further examines the heterogeneity of cities of different sizes and locations. The study finds that the digital economy can increase the level of common prosperity development by 37.1%, and effectively reduce the ratio of urban and rural per capita disposable income, and positively affect the level of common prosperity development in neighboring areas through spatial spillover effects. The mechanism test shows that public management, human capital, R&D and innovation play a significant positive mediating role in the development process of the digital economy affecting the common wealth. The study also confirms that the impact of the digital economy on shared prosperity decreases with the expansion of urban population size, and plays a greater advantage in the location of non-core cities.

Towards zero-carbon cities: how does environmental information disclosure impact the risk of carbon locking?

Utilizing panel data from 119 cities spanning from 2013 to 2019, the study rigorously investigates the influence of environmental information disclosure (EID) on mitigating the risk of urban carbon lock-in through the application of both two-way fixed effects and spatial econometric models. The findings robustly indicate that EID plays a significant role in facilitating carbon unlocking within cities, a conclusion that remains steadfast even after addressing potential endogeneity issues. Mechanistically, it is revealed that EID can effectively decrease the risk of carbon lock-in by enhancing urbanization levels and reducing energy intensity. Furthermore, the study unveils that the efficacy of EID in curbing urban carbon lock-in risk is also contingent upon factors such as carbon intensity, population size, and environmental quality, highlighting the multifaceted nature of EID's impact on urban carbon management. Specifically, the increase in carbon intensity, the expansion of population size, and the deterioration of environmental quality will reduce the carbon unlocking effect, which may not be conducive to the achievement of the goal of "carbon neutrality" in cities. The spatial effect test found that EID and carbon locking risk are both spatially significant positive autocorrelation, the deepening of local EID can significantly cut the risk of carbon locking in neighboring cities, and the urbanization pathway and energy saving pathway exist both locally and in neighboring places.

Evolution of urban land and population system coupling micro–dynamics and macro-stability: Trends and paths

Cities are constantly dynamics, but emerge the stability of urban systems as evidenced by robust urban scaling. Both properties serve as necessary perspectives for understanding the evolution of cities and the complexity of urban systems. However, the evolution trends and mechanisms behind the above relationships, particularly concerning geographic processes, remain to be explored in depth. In this study, the dynamic fluctuations within individual cities and the macro-stability of the evolution of urban systems over the past two decades are systematically demonstrated. Subsequently, associations between urban systems and individual cities are established on the basis of micro- and macro-dynamic mechanisms. The trajectory of individual cities can be divided into scaled–driven growth (SG) and urban excess growth (EG). The former is the result of scaling effects within the system, while the latter is affected by global external effects. Taking as examples urban built-up area and social productions, the evolution pathways that vary from different attributes are analyzed. The results indicate that urban built-up area is more affected by scaling effects whereas social productions is driven more by global external effects. Moreover, we examine the driving forces and mechanism of EG from the aspects of economic development level, technological and economic progress, industrial structure and regional policy through an empirical analysis. This study provides a nuanced understanding of the trends and pathways of urban land–population system from the perspective of complex urban systems and holism to deepen the understanding for general properties of cities and co-evolution for urban population size and land expansion.

Dipterocarpoidae genomics reveal their demography and adaptations to Asian rainforests

Dipterocarpoideae species form the emergent layer of Asian rainforests. They are the indicator species for Asian rainforest distribution, but they are severely threatened. Here, to understand their adaptation and population decline, we assemble high-quality genomes of seven Dipterocarpoideae species including two autotetraploid species. We estimate the divergence time between Dipterocarpoideae and Malvaceae and within Dipterocarpoideae to be 108.2 (97.8‒118.2) and 88.4 (77.7‒102.9) million years ago, and we identify a whole genome duplication event preceding dipterocarp lineage diversification. We find several genes that showed a signature of selection, likely associated with the adaptation to Asian rainforests. By resequencing of two endangered species, we detect an expansion of effective population size after the last glacial period and a recent sharp decline coinciding with the history of local human activities. Our findings contribute to understanding the diversification and adaptation of dipterocarps and highlight anthropogenic disturbances as a major factor in their endangered status.

Go west: Population genomics reveals unexpected population fluctuations and little gene flow in Western hemisphere populations of the predatory lady beetle, Hippodamia convergens.

Hippodamia convergens-the convergent lady beetle, has been used extensively in augmentative biological control of aphids, thrips, and whiteflies across its native range in North America, and was introduced into South America in the 1950s. Overwintering H. convergens populations from its native western range in the United States are commercially collected and released across its current range in the eastern USA, with little knowledge of the effectiveness of its augmentative biological control. Here we use a novel ddRADseq-based SNP/haplotype discovery approach to estimate its range-wide population diversity, differentiation, and recent evolutionary history. Our results indicate (1) significant population differentiation among eastern USA, western USA, and South American populations of H. convergens, with (2) little to no detectable recent admixture between them, despite repeated population augmentation, and (3) continued recent population size expansion across its range. These results contradict previous findings using microsatellite markers. In light of these new findings, the implications for the effectiveness of augmentative biological control using H. convergens are discussed. Additionally, because quantifying the non-target effects of augmentative biological control is a difficult problem in migratory beetles, our results could serve as a cornerstone in improving and predicting the efficacy of future releases of H. convergens across its range.

Genomic analysis of lineage-specific transmission of multidrug resistance tuberculosis in China

ABSTRACT Objectives We investigated the genetic diversities and lineage-specific transmission dynamics of multidrug-resistant tuberculosis (MDR-TB), with the goal of determining the potential factors driving the MDR epidemics in China. Methods We curated a large nationwide Mycobacterium tuberculosis (M. tuberculosis) whole genome sequence data set, including 1313 MDR strains. We reconstructed the phylogeny and mapped the transmission networks of MDR-TB across China using Bayesian inference. To identify drug-resistance variants linked to enhanced transmissibility, we employed ordinary least-squares (OLS) regression analysis. Result The majority of MDR-TB strains in China belong to lineage 2.2.1. Transmission chain analysis has indicated that the repeated and frequent transmission of L2.2.1 plays a central role in the establishment of MDR epidemic in China, but no occurrence of a large predominant MDR outbreak was detected. Using OLS regression, the most common single nucleotide polymorphisms (SNPs) associated with resistance to isoniazid (katG_p.Ser315Thr and katG_p.Ser315Asn) and rifampicin (rpoB_p.Ser450Leu, rpoB_p.His445Tyr, rpoB_p.His445Arg, rpoB_p.His445Asp, and rpoB_p.His445Asn) were more likely to be found in L2 clustered strains. Several putative compensatory mutations in rpoA, rpoC, and katG were significantly associated with clustering. The eastern, central, and southern regions of China had a high level of connectivity for the migration of L2 MDR strains throughout the country. The skyline plot showed distinct population size expansion dynamics for MDR-TB lineages in China. Conclusion MDR-TB epidemic in China is predominantly driven by the spread of highly transmissible Beijing strains. A range of drug-resistance mutations of L2 MDR-TB strains displayed minimal fitness costs and may facilitate their transmission.

PSO‐GA‐SVRmodel for S‐parameters of radio‐frequency power amplifier under different temperatures

AbstractSupport vector machine (SVR) has been introduced into the modeling of S‐parameters in radio‐frequency (RF) power amplifiers (PA). The modeling accuracy and speed of SVR are primarily affected by the penalty parametercand the kernel function coefficientγ. Using the traditional grid search technique to determine these two parameters is time‐consuming and labor‐intensive, and ensuring the model's accuracy is not easy. This article proposes an S‐parameters modeling method based on PSO‐GA‐SVR to improve the SVR's modeling accuracy and speed. The model mainly focuses on particle swarm optimization (PSO) and combines selection, crossover, and mutation operations in genetic algorithms (GA). The fitness values are arranged from small to large in each iteration process, and the first 1/3 are selected for crossover and mutation. Then, the resulting new particle swarms are introduced into the original particle swarm population for searching. On the one hand, PSO‐GA extends the population size and reduces the possibility of falling into local optimization. On the other hand, due to population size expansion, the number of iteration rounds is reduced, and the modeling speed is also increased. The experimental results show that compared to SVR, GA‐SVR, and PSO‐SVR, the proposed PSO‐GA‐SVR can improve the modeling accuracy by more than one magnitude or more while also increasing modeling speed by one magnitude or more. Furthermore, compared with the classical machine learning algorithms such as gradient boosting, random forest, and gcForset, the proposed PSO‐GA‐SVR improves the modeling accuracy by one order of magnitude and the modeling speed by two orders of magnitude more than gradient boosting, random forest, and improves the modeling speed by one order of magnitude more than gcForset.

Not Asian Anymore: Reconstruction of the History, Evolution, and Dispersal of the "Asian" Lineage of CPV-2c.

Variability has been one of the hallmarks of canine parvovirus type 2 (CPV-2) since its discovery, and several lineages and antigenic variants have emerged. Among these, a group of viruses commonly called Asian CPV-2c has recently been reported with increasing frequency in different regions. Currently, its global epidemiology and evolution are essentially unknown. The present work deals with this information gap by evaluating, via sequence, phylodynamic, and phylogeographic analyses, all the complete coding sequences of strains classified as Asian CPV-2c based on a combination of amino acid markers and phylogenetic analysis. After its estimated origin around 2008, this lineage circulated undetected in Asia until approximately 2012, when an expansion in viral population size and geographical distribution occurred, involving Africa, Europe, and North America. Asia was predicted to be the main nucleus of viral dispersal, leading to multiple introduction events in other continents/countries, where infection establishment, persistence, and rapid evolution occurred. Although the dog is the main host, other non-canine species were also involved, demonstrating the host plasticity of this lineage. Finally, although most of the strains showed an amino acid motif considered characteristic of this lineage, several exceptions were observed, potentially due to convergent evolution or reversion phenomena.

Cladistics of Echinococcus granulosus Sensu Stricto Genotypes Infecting the Slaughtered Pigs.

Echinococcus granulosus sensu stricto (s. s.) consists of the most widespread genotypes (G1, G3) implicated in human cystic echinococcosis worldwide. The present study aimed to evaluate the role of pigs in the transmission dynamics of E. granulosus s. s. genotypes, including the phylogenetics, evolutionary divergence and haplotype network analyses of north Indian pig isolates along with GenBank™ archived sequences. In totality, 223 slaughtered pigs were thoroughly screened for the presence of hydatid cysts. The amplification of the partialmitochondrial cytochrome C oxidase subunit 1 gene was performed for the molecular confirmation and phylogenetics of the retrieved metacestodes. The molecularly confirmed and sequenced present study isolates (n = 2) clustered with the E. granulosus genotype G3. The very low evolutionary divergence values (< 0.0027) of the present study isolates with E. ganulosus genotype G3 indicated the absence of diverged lineages. The significantly negative values of Tajima's D (-2.598) and Fu and Li's D (-4.054) of the overall dataset and for the Asian sequences signified an expansion of population size. The overall dataset exhibited low nucleotide (0.067 ± 0.055) and high haplotype (0.574 ± 0.015) diversities, also alluding to demographic expansion. The haplotype network showed that the pig isolates from South America and Europe constituted the predominant haplotype, Hap_2 along with Hap_3 and Hap_6, primarily associated to E. granulosus genotype G1; whereas, the Indian isolates formed different haplotypes (Hap_1 and Hap_5) belonging to genotype G3. The present study highlighted the important role of pigs in the transmission of E. granulosus s. s., which is of paramount significance given the public health and economic impact of cystic echinococcosis.

Phylogeography and demographic dynamics of Rhipicephalus microplus from North India.

Rhipicephalus microplus, a hematophagous vector prevalent in the tropics and subtropics, is responsible for huge economic losses throughout the globe. However, the taxonomy of the tick species, especially prevalent in north India and south China has been challenged in the recent past. The present study attempted to assess the cryptic status of R. microplus ticks of north India based on two mitochondrial markers; the 16S rRNA and cox1 gene sequences. The phylogenetic tree corresponding to both markers demonstrated the presence of three distinct genetic assemblages/ clades of R. microplus. The present study isolates (n=five and seven for the cox1 and 16S rRNA gene sequences, respectively) from north India along with other isolates from India assorted in the R. microplus clade C sensu. Based on the median joining network analysis corresponding to the 16S rRNA gene sequence, 18 haplotypes were recorded, exhibiting a stellate shape, which was indicative of rapid population expansion. For the cox1 gene, the haplotypes corresponding to clades A, B and C were distantly placed with two exceptions. While performing the population structure analysis, low nucleotide (0.04745±0.00416 and 0.01021±0.00146) and high haplotype diversities (0.913±0.032 and 0.794±0.058) were recorded for the different clades of R. microplus based on the cox1 and 16S rRNA mitochondrial markers, respectively. Eventually, high genetic differentiation and low gene flow were recorded among the different clades. A negative value for the neutrality indices (Tajima's D=-1.44125, Fu's Fs=-4.879, Fu and Li's D=-2.78031 and Fu and Li's F=-2.75229) corresponding to the 16S rRNA gene for the overall dataset evinced an expansion of population size. Based on the detailed studies, it was inferred that the R. microplus tick species circulating in north India belonged to clade C sensu, similar to that of the species prevalent in the other parts of the country as well as in the Indian subcontinent.

Exceptional population genomic homogeneity in the black brittle star Ophiocomina nigra (Ophiuroidea, Echinodermata) along the Atlantic-Mediterranean coast

The Atlantic-Mediterranean marine transition is characterised by strong oceanographic barriers and steep environmental gradients that generally result in connectivity breaks between populations from both basins and may lead to local adaptation. Here, we performed a population genomic study of the black brittle star, Ophiocomina nigra, covering most of its distribution range along the Atlantic-Mediterranean region. Interestingly, O. nigra is extremely variable in its coloration, with individuals ranging from black to yellow-orange, and different colour morphs inhabiting different depths and habitats. In this work, we used a fragment of the mitochondrial COI gene and 2,374 genome-wide ddRADseq-derived SNPs to explore: (a) whether the different colour morphs of O. nigra represent different evolutionary units; (b) the disruptive effects of major oceanographic fronts on its population structure; and (c) genomic signals of local adaptation to divergent environments. Our results revealed exceptional population homogeneity, barely affected by oceanographic fronts, with no signals of local adaptation nor genetic differentiation between colour morphs. This remarkable panmixia likely results from a long pelagic larval duration, a large effective population size and recent demographic expansions. Our study unveils an extraordinary phenotypic plasticity in O. nigra, opening further research questions on the ecological and molecular mechanisms underpinning coloration in Ophiuroidea.

Evolutionarily diverse origins of deformed wing viruses in western honey bees

Novel transmission routes can allow infectious diseases to spread, often with devastating consequences. Ectoparasitic varroa mites vector a diversity of RNA viruses, having switched hosts from the eastern to western honey bees (Apis cerana to Apis mellifera). They provide an opportunity to explore how novel transmission routes shape disease epidemiology. As the principal driver of the spread of deformed wing viruses (mainly DWV-A and DWV-B), varroa infestation has also driven global honey bee health declines. The more virulent DWV-B strain has been replacing the original DWV-A strain in many regions over the past two decades. Yet, how these viruses originated and spread remains poorly understood. Here, we use a phylogeographic analysis based on whole-genome data to reconstruct the origins and demography of DWV spread. We found that, rather than reemerging in western honey bees after varroa switched hosts, as suggested by previous work, DWV-A most likely originated in East Asia and spread in the mid-20th century. It also showed a massive population size expansion following the varroa host switch. By contrast, DWV-B was most likely acquired more recently from a source outside East Asia and appears absent from the original varroa host. These results highlight the dynamic nature of viral adaptation, whereby a vector's host switch can give rise to competing and increasingly virulent disease pandemics. The evolutionary novelty and rapid global spread of these host-virus interactions, together with observed spillover into other species, illustrate how increasing globalization poses urgent threats to biodiversity and food security.

New Records of Callinectes sapidus (Crustacea, Portunidae) from Malta and the San Leonardo River Estuary in Sicily (Central Mediterranean)

The current study provides new data on the known records of the Atlantic blue crab, Callinectes sapidus Rathbun, 1896, from the Maltese Islands through two ovigerous female specimens collected from Salini (Malta), a Natura 2000 site, in August 2020. Additional new records of the species were also presented from the San Leonardo River estuary, Syracuse (Sicily, Italy), in August 2022. The specimens collected from both Malta and Sicily were identified using both morphological and molecular analyses. The latter has shown that all the specimens analysed share their mitochondrial DNA barcode region with the most recorded haplotype of the same species from the Mediterranean Sea. Knowing the new distribution and range expansion records of alien species is important for the continuous monitoring of marine alien species, which is essential for the implementation of the best possible mitigation measures. This is especially relevant for C. sapidus, as it is considered one of the worst invaders of the Mediterranean Sea. Therefore, given that this alien species is edible, fishing for it may allow some control over its population size and further expansion. Knowing its distribution is crucial to inform fishers about its exact location, making such a mitigation measure more effective.

Phylogeographic Patterns of the Greater Cane Rat (Thryonomys Swinderianus) Populations from Eastern, Western and Southern Africa and Implications for Wildlife Conservation

Background and aim: African Greater Cane (ACG) rat populations are widely distributed across Africa but the ancestral origin of populations inhabiting different geographical regions is unknown. There is currently no information available on the genetic differentiation of these populations resulting from long-time geographical isolation, environmental and anthropogenic pressures they face in each region. This limits our ability to understand potential speciation processes in AGC and their implications for their conservation and of other small wild mammals exposed to similar ecological conditions in the savannah biomes in Africa. Methods: We analyzed mitochondrial genes targeting the D-loop region of AGC samples from two blocks of the Eastern Arc Mountains in Tanzania (eastern), three agro-ecological zones in Ghana (western) and four sites in South Africa (southern) to characterize the phylogeographical patterns of the species. Results: AGC populations from Tanzania revealed higher haplotype diversity than those from the two other régions. Analysis of variance (AMOVA) revealed higher genetic variations within than between populations in all geographic regions. Demographic history analysis revealed a negative and significant Tajima’s D for a single southern African population as well as close relatedness with the Eastern Africa populations, suggesting a common ancestral origin. Conclusion: This is the first study to compare maternal lineages of AGC populations from eastern, western and southern Africa and provides a basis for future genetic studies of the species, its evolutionary ecology and the conservation of these populations across their range. Conservation implications: Evidence of recent population size expansion underpinned by high genetic diversity observed in the Uluguru urban (Tanzania) Guinea savannah (Ghana) and Kwazulu Natal (South Africa) populations suggest that the AGC populations are not currently threatened, consistent with the current IUCN status. Also, the distinct haplotypes observed in each region suggest that the populations can be managed as meta-populations, thus providing opportunities for potential local game-farming programs.

Possible species discrimination of a blotched nerite Nerita albicilla with their distribution pattern and demographic history in the Indo-Pacific

The blotched nerite Nerita albicilla (Linnaeus 1758) is distributed in intertidal areas of the Indo-Pacific. In South Korea, it has been found only in the southernmost region of Jeju Island so far. Owing to its limited distribution, it can be a promising intertidal species helpful for monitoring global climate change effects in the Korean Peninsula. We performed population genetic analyses based on 393 COI haplotypes from 697 N. albicilla, including 167 from this study and 530 from public databases. The results showed that there are two distinct genetic lineages in N. albicilla: PAIO (Palearctic, Australasia, Indo-Malay, and Oceania) and Afrotropic lineages. DNA barcoding gap analyses indicated that the two lineages could be differentiated into two different species: N. albicilla (PAIO) and N. originalis sp. nov. (Afrotropic) (3.96%). Additionally, it was revealed that their divergence time was ca. 5.96 Ma and dramatic diversification of COI haplotypes occurred during the late Pliocene and Pleistocene. The results of MDA, BSP, and neutrality test implied recent population size expansion, which was estimated to be ca. 250 Ka. Finally, we discussed whether the observation of N. originalis sp. nov. in South Korea is due to the northward migration through ocean currents caused by global warming or due to artificial activity through marine transportation.

Genetic evidence for widespread population size expansion in North American boreal birds prior to the Last Glacial Maximum.

Pleistocene climate cycles are well documented to have shaped contemporary species distributions and genetic diversity. Northward range expansions in response to deglaciation following the Last Glacial Maximum (LGM; approximately 21 000 years ago) are surmised to have led to population size expansions in terrestrial taxa and changes in seasonal migratory behaviour. Recent findings, however, suggest that some northern temperate populations may have been more stable than expected through the LGM. We modelled the demographic history of 19 co-distributed boreal-breeding North American bird species from full mitochondrial gene sets and species-specific molecular rates. We used these demographic reconstructions to test how species with different migratory strategies were affected by glacial cycles. Our results suggest that effective population sizes increased in response to Pleistocene deglaciation earlier than the LGM, whereas genetic diversity was maintained throughout the LGM despite shifts in geographical range. We conclude that glacial cycles prior to the LGM have most strongly shaped contemporary genetic diversity in these species. We did not find a relationship between historic population dynamics and migratory strategy, contributing to growing evidence that major switches in migratory strategy during the LGM are unnecessary to explain contemporary migratory patterns.

Phylogeography and transmission of Mycobacterium tuberculosis spanning prisons and surrounding communities in Paraguay

Recent rises in incident tuberculosis (TB) cases in Paraguay and the increasing concentration of TB within prisons highlight the urgency of targeting strategies to interrupt transmission and prevent new infections. However, whether specific cities or carceral institutions play a disproportionate role in transmission remains unknown. We conducted prospective genomic surveillance, sequencing 471 Mycobacterium tuberculosis complex genomes, from inside and outside prisons in Paraguay’s two largest urban areas, Asunción and Ciudad del Este, from 2016 to 2021. We found genomic evidence of frequent recent transmission within prisons and transmission linkages spanning prisons and surrounding populations. We identified a signal of frequent M. tuberculosis spread between urban areas and marked recent population size expansion of the three largest genomic transmission clusters. Together, our findings highlight the urgency of strengthening TB control programs to reduce transmission risk within prisons in Paraguay, where incidence was 70 times that outside prisons in 2021.

Numbers and distribution of eagles in the Czech Republic: reason for optimism

At present, four eagle species are regular breeders in the Czech Republic: White-Tailed Sea-Eagle (Haliaeetus albicilla), Eastern Imperial Eagle (Aquila heliaca), Lesser Spotted Eagle (Aquila [Clanga] pomarina), and Golden Eagle (Aquila chrysaetos). There is also a historical record of confirmed breeding of Greater Spotted Eagle (Aquila [Clanga] clanga) in 1847 in the region of Eastern Bohemia. At present, this species is only a rare visitor on autumn and spring passage. White-Tailed Sea Eagle (Haliaeetus albicilla) has re-colonized the country since the mid-1980s, after approx. one hundred years absence as a breeder. The recolonization started at two separate areas: South Bohemia, with birds coming mainly from German and Poland source population, and South Moravia, with birds coming from Pannonia, i.e. mainly from Hungary. Later the species occupied nearly the whole area of suitable habitats. The population is still increasing and currently reaches ca. 160–180 breeding pairs. Eastern Imperial Eagle (Aquila heliaca) has bred in the Check Republic for the first time in 1998 (there is one unconfirmed note about its breeding also in 1930s). Since then, the population is slowly growing, similarly as the whole Pannonian population (with source populations mainly in Hungary and Eastern Slovakia). So far, the breeding is limited to Pannonian region of the country (i.e. south-eastern part) but the first pairs have spread north and west during the last years. 23 territorial pairs were monitored in 2023. The Lesser Spotted Eagle (LSE, Aquila [Clanga] pomarina) was never an abundant species in the Czech Republic. Scarce population existed especially in submountainous regions and in the mountains throughout the CR. The last breeding pairs were documented at the end of 1980’s and the beginning of 1990’s in Šumava Mts. After more than 20 years, breeding of LSE was surprisingly confirmed in western Bohemia. Since 2012, one pair breeds there regularly. In 2018, a small new population was discovered in the north of Morava region and currently three pairs breed there. Since 2020, another pair breeds in eastern Bohemia. Presently, two pairs breed in this region but we presume that another 1–2 breeding pairs could be hiding in this region. Therefore, in 2023 we estimate the breeding population of the LSE in the Czech Republic as 6–8 pairs. Golden Eagle (Aquila chrysaetos) started to breed in Northern Moravia (NE part of the country) in 2013 (similarly as White-Tailed Sea Eagle after more than one hundred years of absence), partly as a result of re-introduction project carried out during 2006–2017 (27 birds originating from Slovakia were released). Along with the growing population size and range expansion in the neighbouring Slovakia, there are now min 5–6 occupied territories in Czech Republic. Osprey (Pandion haliaetus) is a regular passage migrant (both spring and autumn), with some birds oversummering and very rarely also wintering. Historical data about its regular breeding in 19th or even early 20th century are not fully confirmed. Several breeding attempts were recorded also during the last 20 years. These included nest building but the clutch was never laid. Short-Toed Snake Eagle (Circaetus gallicus) is a very rare visitor in the Czech Republic, but the number of observations seems to increase rapidly during the last 10 years. With the increasing population in Hungary and possibly also Slovakia, its breeding cannot be excluded in the near future.