Laser interaction with an ultra-thin pre-structured target is investigated with the help of both two-dimensional and three-dimensional particle-in-cell simulations. With the existence of a periodic structure on the target surface, the laser seems to penetrate through the target at its fundamental frequency even if the plasma density of the target is much higher than the laser’s relativistically critical density. The particle-in-cell simulations show that the transmitted laser energy behind the pre-structured target is increased by about two orders of magnitude compared to that behind the flat target. Theoretical analyses show that the transmitted energy behind the pre-structured target is actually re-emitted by electron ‘islands’ formed by the surface plasma waves on the target surfaces. In other words, the radiation with the fundamental frequency is actually ‘surface emission’ on the target rear surface. Besides the intensity of the component with the fundamental frequency, the intensity of the high-order harmonics behind the pre-structured target is also much enhanced compared to that behind the flat target. The enhancement of the high-order harmonics is also related to the surface plasma waves generated on the target surfaces.

Triploid and pentaploid breeding is of great importance in agricultural production, but it is not always easy to obtain double ploidy parents. However, in fishes, chromosome ploidy is diversiform, which may provide natural parental resources for triploid and pentaploid breeding. Both tetraploid and hexaploid exist in Schizothorax fishes, which were thought to belong to different subfamilies with tetraploid Percocypris fishes in morphology, but they are sister genera in molecule. Fortunately, the pentaploid hybrid fishes have been successfully obtained by hybridization of Schizothorax wangchiachii (♀, 2n = 6X = 148) × Percocypris pingi (♂, 2n = 4X = 98). To understand the genetic and morphological difference among the hybrid fishes and their parents, four methods were used in this study: morphology, karyotype, red blood cell (RBC) DNA content determination and inter-simple sequence repeat (ISSR). In morphology, the hybrid fishes were steady, and between their parents with no obvious preference. The chromosome numbers of P. pingi have been reported as 2n = 4X = 98. In this study, the karyotype of S. wangchiachii was 2n = 6X = 148 = 36m + 34sm + 12st + 66t, while that the hybrid fishes was 2n = 5X = 123 = 39m + 28sm + 5st + 51t. Similarly, the RBC DNA content of the hybrid fishes was intermediate among their parents. In ISSR, the within-group genetic diversity of hybrid fishes was higher than that of their parents. Moreover, the genetic distance of hybrid fishes between P. pingi and S.wangchiachii was closely related to that of their parental ploidy, suggesting that parental genetic material stably coexisted in the hybrid fishes. This is the first report to show a stable pentaploid F1 hybrids produced by hybridization of a hexaploid and a tetraploid in aquaculture.

To compare the epidemiologic features (e.g. settings and transmission mode) and patient clinical characteristics associated with outbreaks of different norovirus (Nov) strains, we retrospectively analysed data of Nov outbreaks occurring in Guangzhou, China from 2012 to 2018. The results suggested that outbreaks of Nov GII.2, GII.17 and GII.4 Sydney exhibited different outbreak settings, transmission modes and symptoms. GII.2 outbreaks mainly occurred in kindergartens, elementary and high schools and were transmitted mainly through person-to-person contact. By contrast, GII.4 Sydney outbreaks frequently occurred in colleges and were primarily associated with foodborne transmission. Cases from GII.2 and GII.17 outbreaks reported vomiting more frequently than those from outbreaks associated with GII.4 Sydney.

Generally, earth rotating and non-spherical perturbation of the earth in re-entry motion model are simplified using the standard trajectory guidance method. The re-entry motion is also simplified to horizontal motion and vertical motion and controlled, respectively. The simplification of re-entry motion model will lead to loss of motion accuracy and location accuracy. The direct decomposition will lead to the reduction of control accuracy because the horizontal motion and the vertical motion are coupled in re-entry. To improve the standard trajectory guidance method, the standard trajectory guidance method based on decoupling control variables and waypoint is proposed in this paper. The proposed guidance method will not simplify earth rotating and non-spherical perturbation of the earth in motion equation or decompose the re-entry motion to horizontal motion and vertical motion. Trajectory waypoint is adopted to reduce the change frequency of tracking states, because tracking states change frequently if the entire standard trajectory is tracked in real time.

Thermoacoustic instability – self-sustained pressure oscillations triggered by temperature gradients – has become an increasingly studied topic in the context of energy conversion. Generally, the process relies on conductive heat transfer between a solid and the fluid in which the generated pressure oscillations are sustained. In the present study, the thermoacoustic theory is extended to include mass transfer; specifically, the working fluid is modified so as to incorporate a ‘reactive’ gas, able to exchange phase with a solid/liquid boundary through a sorption process (or through evaporation/condensation), such that most heat is transferred in the form of latent heat rather than through conduction. A set of differential equations is derived, accounting for phase-exchange heat and mass transfer, and de-coupled via a small-amplitude asymptotic expansion. These equations are solved and subsequently manipulated into the form of a wave equation, representing the small perturbation on the pressure field, and used to derive expressions for the time-averaged, second-order heat and mass fluxes. A stability analysis is performed on the wave equation, from which the marginal stability curve is calculated in terms of the temperature difference, $\unicode[STIX]{x0394}T_{onset}$ , required for initiation of self-sustained oscillations. Calculated stability curves are compared with published experimental results, showing good agreement. Effects of gas mixture composition are studied, indicating that a lower heat capacity of the inert component, combined with a low boiling temperature and high latent heat of the reactive component substantially lower $\unicode[STIX]{x0394}T_{onset}$ . Furthermore, an increase in the average mole fraction of the reactive gas, $C_{m}$ strongly affects onset conditions, leading to $\unicode[STIX]{x0394}T_{onset}\sim 5\,^{\circ }\text{C}$ at the highest value of $C_{m}$ achievable under atmospheric pressure. An analysis of the system limit cycle is performed for a wide range of parameters, indicating a systematic decrease in the temperature difference capable of sustaining the limit cycle, as well as a significant distortion of the acoustic wave form as the phase-exchange mechanism becomes dominant. These findings, combined, reveal the underlying mechanisms by which a phase-exchange engine may produce more acoustic power than its counterpart ‘classical’ thermoacoustic system, while its temperature difference is substantially lower.

The seasonality of individual influenza subtypes/lineages and the association of influenza epidemics with meteorological factors in the tropics/subtropics have not been well understood. The impact of the 2009 H1N1 pandemic on the prevalence of seasonal influenza virus remains to be explored. Using wavelet analysis, the periodicities of A/H3N2, seasonal A/H1N1, A/H1N1pdm09, Victoria and Yamagata were identified, respectively, in Panzhihua during 2006–2015. As a subtropical city in southwestern China, Panzhihua is the first industrial city in the upper reaches of the Yangtze River. The relationship between influenza epidemics and local climatic variables was examined based on regression models. The temporal distribution of influenza subtypes/lineages during the pre-pandemic (2006–2009), pandemic (2009) and post-pandemic (2010–2015) years was described and compared. A total of 6892 respiratory specimens were collected and 737 influenza viruses were isolated. A/H3N2 showed an annual cycle with a peak in summer–autumn, while A/H1N1pdm09, Victoria and Yamagata exhibited an annual cycle with a peak in winter–spring. Regression analyses demonstrated that relative humidity was positively associated with A/H3N2 activity while negatively associated with Victoria activity. Higher prevalence of A/H1N1pdm09 and Yamagata was driven by lower absolute humidity. The role of weather conditions in regulating influenza epidemics could be complicated since the diverse viral transmission modes and mechanism. Differences in seasonality and different associations with meteorological factors by influenza subtypes/lineages should be considered in epidemiological studies in the tropics/subtropics. The development of subtype- and lineage-specific prevention and control measures is of significant importance.

Hypoimmunity and numerous stresses are two major challenges in broiler industry. Nutrient intervention at the specific time of embryonic stage is a feasible way to improve animal performance. This study was conducted to investigate the possible effects of in ovo feeding (IOF) of vitamin C at embryonic age 15th day (E15) on growth performance, antioxidation and immune function of broilers. A total of 240 broiler fertile eggs were randomly divided into two groups (0 and 3 mg injected dose of vitamin C at E15), and new-hatched chicks from each treatment were randomly allocated into six replicates with 10 chicks per replicate after incubation. The results indicated that in ovo vitamin C injection improved the hatchability (P < 0.05) and increased immunoglobulin M (IgM) (at the broiler’s age 1st day, D1), IgG and IgM concentrations (D21), as well as lysozyme activity (D21, P < 0.05) and total antioxidant capacity (D42, P < 0.01) in plasma of broilers. On D21, the splenic expression level of DNA methyltransferase 1 (DNMT1) was up-regulated in vitamin C (VC) group, whereas interleukin (IL)-6, interferon-γ, ten-eleven translocation protein 1 and thymine-DNA glycosylase were down-regulated (P < 0.05). On D42, in ovo vitamin C injection up-regulated splenic expression levels of DNMT1, DNA methyltransferase 3B (DNMT3B) and growth arrest and DNA-damage-inducible protein beta (P < 0.05), whereas down-regulated splenic expression levels of IL-6, tumour necrosis factor-α and methyl-CpG-binding domain protein 4 (P < 0.05). Our findings suggested that IOF of 3 mg vitamin C at E15 could improve, to some extent, the antioxidant activity and immune function in plasma, corresponding with the lower expression of pro-inflammatory cytokines in spleen. However, IOF of vitamin C leading to the changes in the expression of DNA methyltransferases and demethylases may suggest an increased trend of DNA methylation level in spleen and whether DNA methylation variation is associated with the lower expression of pro-inflammatory cytokines in spleen warrants future study.

Heading date (HD) and flowering date (FD) are critical for yield potential and stability, so understanding their genetic foundation is of great significance in wheat breeding. Three related recombinant inbred line populations with a common female parent were developed to identify quantitative trait loci (QTL) for HD and FD in four environments. In total, 25 putative additive QTL and 20 pairwise epistatic effect QTL were detected in four environments. The additive QTL were distributed across 17 wheat chromosomes. Of these, QHd-1A, QHd-1D, QHd-2B, QHd-3B, QHd-4A, QHd-4B and QHd-6D were major and stable QTL for HD. QFd-1A, QFd-2B, QFd-4A and QFd-4B were major and stable QTL for FD. In addition, an epistatic interaction test showed that epistasis played important roles in controlling wheat HD and FD. Genetic relationships between HD/FD and five yield-related traits (YRTs) were characterized and ten QTL clusters (C1–C10) simultaneously controlling YRTs and HD/FD were identified. The present work laid a genetic foundation for improving yield potential in wheat molecular breeding programmes.

Muons produced by the Bethe–Heitler process from laser wakefield accelerated electrons interacting with high $Z$ materials have velocities close to the laser wakefield. It is possible to accelerate those muons with laser wakefield directly. Therefore for the first time we propose an all-optical ‘Generator and Booster’ scheme to accelerate the produced muons by another laser wakefield to supply a prompt, compact, low cost and controllable muon source in laser laboratories. The trapping and acceleration of muons are analyzed by one-dimensional analytic model and verified by two-dimensional particle-in-cell (PIC) simulation. It is shown that muons can be trapped in a broad energy range and accelerated to higher energy than that of electrons for longer dephasing length. We further extrapolate the dependence of the maximum acceleration energy of muons with the laser wakefield relativistic factor $\unicode[STIX]{x1D6FE}$ and the relevant initial energy $E_{0}$ . It is shown that a maximum energy up to 15.2 GeV is promising with $\unicode[STIX]{x1D6FE}=46$ and $E_{0}=1.45~\text{GeV}$ on the existing short pulse laser facilities.

The distribution of human papillomaviruses (HPVs) must be understood for the control and prevention of cervical cancer. Community-based Papanicolaou and HPV DNA tests were performed on 41 578 women. The prevalences of HPV genotypes 16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, 66 and 68 were assessed. In total, 10% women were infected/co-infected by these HPVs. The infection rate increased from 7.1% in women aged ⩽30 years to 10.4% in those aged 50–60 years, and then decreased slightly to 9.9% in those aged >60 years. The HPV 16 and 58 positivity rates were significantly higher among women with high-grade squamous intraepithelial lesions (HSILs) than among those with cervicitis/negativity for intraepithelial lesion or malignancy (NILM) or low-grade SILs (LSILs). The HPV 18, 52 and 68 infection rates were significantly lower in women with HSILs than in those with NILM or LSILs. The proportion of women infected by multiple HPV strains was higher among those with HSILs. The proportions of the five most common genotypes, HPV 16, 18, 33, 52 and 58, increased with the number of co-infecting strains. HPV 16 and 58 were the high-risk HPVs in the Shanghai community and should be the focus in HPV screening and vaccination.

Full-scale simulations of a (Magnetorheological) MR damper are carried out for revealing its hysteretic behaviors associated with implementation of semi-active control using the routine of computational fluid dynamics. By virtue of the structural symmetry of the MR damper, a two-dimensional configuration for finite element simulation is built up. Herschel-Bulkley model is employed to represent the property of the MR fluid, of which the control parameters and their relevances to the input current are addressed. Typical cases involving sinusoidal and irregular displacements, steady and transient currents loaded upon the MR damper are investigated. Numerical investigations reveal that the damper force has a positive correlation with input current, excitation amplitude and excitation frequency. The full-scale simulation is proved to exhibit a sound accuracy through the validation of experimental data. It provides a logical manner revealing the true performance of MR dampers under desirable operating modes in practice, and can be readily integrated with the gain design of the associated semi-actively controlled structure. This progress bypasses the technical challenge inherent in the traditional tests with low-frequency cyclic loadings due to the limitation of experimental setup. Besides, comparative study between two-dimensional and three-dimensional configuration simulations of the MR damper shows that former has a better applicability, which can be carried out on a low-cost platform.

Good canopy structure is essential for optimal maize (Zea mays L.) production. However, creating appropriate maize canopy structure can be difficult, because the characteristics of individual plants are altered by changes in plant age, density and interactions with neighbouring plants. The objective of the current study was to find a reliable method for building good maize canopy structure by analysing changes in canopy structure, light distribution and grain yield (GY). A modern maize cultivar (ZhengDan958) was planted at 12 densities ranging from 1.5 to 18 plants/m2 at two field locations in Xinjiang, China. At the silking stage (R1), plant and ear height increased with plant density as well as leaf area index (LAI), whereas leaf area per plant decreased logarithmically. The fraction of light intercepted by the plant (F) increased with increasing plant density, but the light extinction coefficient (K) decreased linearly from 0.61 to 0.39. Taking the optimum value of F (95%) as an example, and using measured values of K for each plant density at R1 and the equation from Beer's law, the corresponding (theoretical) LAI for each plant density was calculated and optimum plant density (9.72 plants/m2) obtained by calculating the difference between theoretical LAIs and actual observations. Further analysis showed that plant density ranging from 10.64 to 11.55 plants/m2 yielded a stable GY range. Therefore, taking into account the persistence time for maximum LAI, the plant density required to obtain an ideal GY maize canopy structure should be increased by 10–18% from 9.72 plants/m2.

The gastrointestinal tract (GIT) of animals is capable of sensing various kinds of nutrients via G-protein coupled receptor-mediated signaling transduction pathways, and the process is known as ‘gut nutrient chemosensing’. GPR40, GPR41, GPR43 and GPR119 are chemoreceptors for free fatty acids (FFAs) and lipid derivatives, but they are not well studied in small ruminants. The objective of this study is to determine the expression of GPR40, GPR41, GPR43 and GPR119 along the GIT of kid goats under supplemental feeding (S) v. grazing (G) during early development. In total, 44 kid goats (initial weight 1.35±0.12 kg) were slaughtered for sampling (rumen, abomasum, duodenum, jejunum, ileum, cecum, colon and rectum) between days 0 and 70. The expression of GPR41 and GPR43 were measured at both mRNA and protein levels, whereas GPR40 and GPR119 were assayed at protein level only. The effects of age and feeding system on their expression were variable depending upon GIT segments, chemoreceptors and expression level (mRNA or protein), and sometimes feeding system × age interactions (P<0.05) were observed. Supplemental feeding enhanced expression of GPR40, GPR41 and GPR43 in most segments of the GIT of goats, whereas G enhanced expression of GPR119. GPR41 and GPR43 were mainly expressed in rumen, abomasum and cecum, with different responses to age and feeding system. GPR41 and GPR43 expression in abomasum at mRNA level was greatly (P<0.01) affected by both age and feeding system; whereas their expression in rumen and abomasum at protein level were different, feeding system greatly (P<0.05) affected GPR41 expression, but had no effect (P>0.05) on GPR43 expression; and there were no feeding system×age interactions (P>0.05) on GPR41 and GPR43 protein expression. The expression of GPR41 and GPR43 in rumen and abomasum linearly (P<0.01) increased with increasing age (from days 0 to 70). Meanwhile, age was the main factor affecting GPR40 expression throughout the GIT. These outcomes indicate that age and feeding system are the two factors affecting chemoreceptors for FFAs and lipid derivatives expression in the GIT of kids goats, and S enhanced the expression of chemoreceptors for FFAs, whereas G gave rise to greater expression of chemoreceptors for lipid derivatives. Our results suggest that enhanced expression of chemoreceptors for FFAs might be one of the benefits of early supplemental feeding offered to young ruminants during early development.

Compensatory gain describes an accelerated growth seen in animals following a period of nutrient restriction. Methionine (Met) is the second limiting amino acid in typical swine diets and is essential for muscle growth. This study was conducted to determine (1) if a Met-deficient diet can cause growth retardation in growing pigs, (2) if returning to a normal feeding can yield compensatory gain in the pigs previously fed the Met-deficient diet, and (3) if this Met-deficiency followed by the normal feeding program affects carcass characteristics. Twenty individually-penned crossbred young barrows were randomly allotted to two dietary treatments (n = 10). One Met-deficient (D1) and one Met-adequate (D2) diets were formulated based on corn and soybean meal and fed to respective pigs for 31 days. After that, all pigs were fed the same commercial grower-finisher diet until market weight (around 125 kg), then slaughtered, and carcass characteristics measured. The D1 and D2 pigs began with similar body weights (23.5 vs. 23.6 kg; P = 0.935), but after 31-days on the dietary treatments, D1 pigs were lighter than D2 pigs (51.6 vs. 55.0 kg; P = 0.102). After feeding the normal diet for 55 days, D1 and D2 pigs had similar body weights (122.7 vs. 122.6 kg; P = 0.989). In terms of carcass characteristics, however, D1 pigs had thicker back-fat (at 10th rib; 2.95 vs. 2.51 cm; P = 0.015), heavier belly weight (11.0 vs. 9.6 kg; P = 0.005), lighter ham weights (untrimmed: 20.8 vs. 21.6 kg; P = 0.043; trimmed: 19.6 vs. 20.6 kg; P = 0.016), lighter picnic shoulder weight (8.72 vs. 9.80 kg; P = 0.041), lighter total lean cut weight (51.8 vs. 53.8 kg; P = 0.055), and lower lean cut percentage (56.4 vs. 59.0%; P = 0.012). These results indicate that the Met-deficient diet produced growth-retarded pigs, which showed compensatory gain after the normal feeding. At slaughter, the pigs previously fed the Met-deficient diet had more fat and less lean tissue than their non-deficient counterparts.

China has been the world's largest producer of eggs for the last 30 years. There have been considerable recent changes in the structure of the egg industry due to rapid economic growth, improved supply chains and favourable prices. Since 2012, revenue from poultry has increased by 8.7% on average. Annual growth rate of egg production has been approximately 0.6 million tonnes per year from 2000 to 2016 (FAO, 2017), and 2016 saw the total eggs produced reaching a peak of 31 million tonnes. Egg production in China continues to provide the population with a significant proportion of their intake of high quality dietary protein. Egg consumption is expected to continue to increase with the rise in the urban population. It is probable that the highest increase in demand for eggs will be for ‘out of home’ consumption and in processed food products. These levels of production are associated with the spread of intensive systems. The poultry sector is no longer dominated by hundreds of millions of smallholders keeping birds as a side-line activity and many small farmers have ceased production. Chinese consumers are becoming more focused on the quality and safety of eggs. Future developments in the egg production are expected to concentrate on quality, safety and traceability of eggs.

Cereal grains treated with organic acids were proved to increase ruminal resistant starch and can relieve the risk of ruminal acidosis. However, previous study mainly focussed on acid-treated barley, the effects of organic acid-treated corn is still unknown. The objectives of this study were to evaluate whether feeding ground corn steeped in citric acid (CA) would affect ruminal pH and fermentation patterns, milk production and innate immunity responses in dairy goats. Eight ruminally cannulated Saanen dairy goats were used in a crossover designed experiment. Each experimental period was 21 day long including 14 days for adaption to new diet and 7 days for sampling and data collection. The goats were fed high-grain diet contained 30% hay and 70% corn-based concentrate. The corn was steeped either in water (control) or in 0.5% (wt/vol) CA solution for 48 h. Goats fed CA diet showed improved ruminal pH status with greater mean and minimum ruminal pH, and shorter (P<0.05) duration of ruminal pH<5.6 and less area of ruminal pH<5.6, 5.8 and 6.0. Concentration of total volatile fatty acid and molar proportion of propionate were less but the molar proportion of acetate was greater (P<0.05) in goats fed the CA diet than the control diet. Concentration of ruminal lipopolysaccharide (LPS) was lower (P<0.05) and that of lactic acid also tended (P<0.10) to be lower in goats fed CA than the control. Although dry matter intake, actual milk yield, yield and content of milk protein and lactose were not affected, the milk fat content and 4% fat-corrected milk tended (P<0.10) to be greater in goats fed CA diet. For the inflammatory responses, peripheral LPS did not differ, whereas the concentration of LPS binding protein and serum amyloid A tended (P<0.10) to be less in goats fed CA diet. Similarly, goats fed CA diet had less (P<0.05) concentration of haptoglobin and tumour necrosis factor. These results indicated that feeding ground corn treated with CA effectively improved ruminal pH status, thus alleviated the risk of ruminal acidosis, reduced inflammatory response, and tend to improve milk yield and milk fat test.

Dengue is the fastest spreading mosquito-transmitted disease in the world. In China, Guangzhou City is believed to be the most important epicenter of dengue outbreaks although the transmission patterns are still poorly understood. We developed an autoregressive integrated moving average model incorporating external regressors to examine the association between the monthly number of locally acquired dengue infections and imported cases, mosquito densities, temperature and precipitation in Guangzhou. In multivariate analysis, imported cases and minimum temperature (both at lag 0) were both associated with the number of locally acquired infections (P < 0.05). This multivariate model performed best, featuring the lowest fitting root mean squared error (RMSE) (0.7520), AIC (393.7854) and test RMSE (0.6445), as well as the best effect in model validation for testing outbreak with a sensitivity of 1.0000, a specificity of 0.7368 and a consistency rate of 0.7917. Our findings suggest that imported cases and minimum temperature are two key determinants of dengue local transmission in Guangzhou. The modelling method can be used to predict dengue transmission in non-endemic countries and to inform dengue prevention and control strategies.

In vivo and in vitro trials were conducted to assess the effects of tributyrin (TB) supplementation on short-chain fatty acid (SFCA) concentrations, fibrolytic enzyme activity, nutrient digestibility and methanogenesis in adult sheep. Nine 12-month-old ruminally cannulated Small Tail ewes (initial body weight 55 ± 5.0 kg) without pregnancy were used for the in vitro trial. In vitro substrate made to offer TB at 0, 2, 4, 6 and 8 g/kg on a dry matter (DM) basis was incubated by ruminal microbes for 72 h at 39°C. Forty-five adult Small Tail ewes used for the in vivo trial were randomly assigned to five treatments with nine animals each for an 18-d period according to body weight (55 ± 5.0 kg). Total mixed ration fed to ewes was also used to offer TB at 0, 2, 4, 6 and 8 g/kg on a DM basis. The in vitro trial showed that TB supplementation linearly increased apparent digestibility of DM, crude protein, neutral detergent fibre and acid detergent fibre, and enhanced gas production and methane emissions. The in vivo trial showed that TB supplementation decreased DM intake, but enhanced ruminal fermentation efficiency. Both in vitro and in vivo trials showed that TB supplementation enhanced total SFCA concentrations and carboxymethyl cellulase activity. The results indicate that TB supplementation might exert advantage effects on rumen microbial metabolism, despite having an enhancing effect on methanogenesis.

Multiple human immunodeficiency virus (HIV)-1 genotypes in China were first discovered in Yunnan Province before disseminating throughout the country. As the HIV-1 epidemic continues to expand in Yunnan, genetic characteristics and transmitted drug resistance (TDR) should be further investigated among the recently infected population. Among 2828 HIV-positive samples newly reported in the first quarter of 2014, 347 were identified as recent infections with BED-captured enzyme immunoassay (CEIA). Of them, 291 were successfully genotyped and identified as circulating recombinant form (CRF)08_BC (47.4%), unique recombinant forms (URFs) (18.2%), CRF01_AE (15.8%), CRF07_BC (14.4%), subtype C (2.7%), CRF55_01B (0.7%), subtype B (0.3%) and CRF64_BC (0.3%). CRF08_BC and CRF01_AE were the predominant genotypes among heterosexual and homosexual infections, respectively. CRF08_BC, URFs, CRF01_AE and CRF07_BC expanded with higher prevalence in central and eastern Yunnan. The recent common ancestor of CRF01_AE, CRF07_BC and CRF08_BC dated back to 1983.1, 1992.1 and 1989.5, respectively. The effective population sizes (EPS) for CRF01_AE and CRF07_BC increased exponentially during 1991–1999 and 1994–1999, respectively. The EPS for CRF08_BC underwent two exponential growth phases in 1994–1998 and 2001–2002. Lastly, TDR-associated mutations were identified in 1.8% of individuals. These findings not only enhance our understanding of HIV-1 evolution in Yunnan but also have implications for vaccine design and patient management strategies.