Sexual and reproductive behaviour of Drosophila melanogaster from a microclimatically interslope differentiated population of "Evolution Canyon" (Mount Carmel, Israel).

The strong microscale interslope environmental differences in "Evolution Canyon" provide an excellent natural model for sympatric speciation. Our previous studies revealed significant slope-specific differences for a fitness complex of Drosophila. This complex involved either adaptation traits (tolerance to high temperature, different viability and longevity pattern) or behavioural differentiation, manifested in habitat choice and non-random mating. This remarkable differentiation has evolved despite a very small interslope distance (a few hundred metres only). Our hypothesis is that strong interslope microclimatic contrast caused differential selection for fitness-related traits accompanied by behavioural differentiation and reinforced by some sexual isolation, which started incipient speciation. Here we describe the results of a systematic analysis of sexual behaviour in a non-choice situation and several reproductive parameters of D. melanogaster populations from the opposite slopes of "Evolution Canyon". The evidence indicates that: (i) mate choice derives from differences in mating propensity and discrimination; (ii) females from the milder north-facing slope discriminate strongly against males of the opposite slope; (iii) both sexes of the south-facing slope display distinct reproductive and behavioural patterns with females showing increased fecundity, shorter time before remating and relatively higher receptivity, and males showing higher mating propensity. These patterns represent adaptive life strategies contributing to higher fitness.     

Microsatellite analysis of Drosophila melanogaster populations along a microclimatic contrast at lower Nahel Oren canyon, Mount Carmel, Israel

Drosophila melanogaster populations collected at the south-facing slope (SFS) and north-facing slope (NFS) of lower Nahel Oren canyon, Mount Carmel, Israel display significant differences in survival and longevity at temperature, drought, and starvation stresses. Furthermore, significant assortative mating was previously observed between populations of the two slopes. We used a set of 48 microsatellite markers to analyze patterns of genetic differentiation between D. melanogaster populations from both slopes and D. simulans. Consistent with previous reports, we found D. simulans to be well differentiated from D. melanogaster. Genetic differentiation between SFS and NFS D. melanogaster populations was low (F(ST) = 0.0012). Also a tree of individuals based on the proportion of shared alleles and a model-based clustering method provided no evidence for population substructuring.     

Multidimensional analysis of Drosophila wing variation in Evolution Canyon

Environmental stress has been suggested to be a major evolutionary force, both through inducing strong selection and because of its direct impact on developmental buffering processes that alter the evolvability of organisms. In particular, temperature has attracted much attention because of its importance as an ecological feature and the relative ease with which it can be experimentally manipulated in the lab. Evolution Canyon, Lower Nahal Oren, Israel, is a well studied natural site where ecological parameters are suspected to drive evolutionary differentiation. In this study, using Drosophila melanogaster isofemale lines derived from wild flies collected on both slopes of the canyon, we investigated the effect of developmental temperature upon the different components of phenotypic variation of a complex trait: the wing. Combining geometric and traditional morphometrics, we find only limited evidence for a differentiation among slopes. Investigating simultaneously phenotypic plasticity, genetic variation among isofemale lines, variation among individuals and fluctuating asymmetry, we could not identify a consistent effect of the stressful conditions encountered on the south facing slope. The prevailing structuring effect is that of the experimentally manipulated temperature which clearly influences wing mean size and shape. Variability, in contrast, is not consistently affected by temperature. Finally, we investigated the specific relationship between individual variation and fluctuating asymmetry. Using metric multi-dimensional scaling we show that the related patterns of wing shape variation are not identical, supporting the view that the underlying developmental processes are to a certain extent different.     

Fluctuating Helical Asymmetry and Morphology of Snails (Gastropoda) in Divergent Microhabitats at 'Evolution Canyons I and II,' Israel

Background

Developmental instability of shelled gastropods is measured as deviations from a perfect equiangular (logarithmic) spiral. We studied six species of gastropods at ‘Evolution Canyons I and II’ in Carmel and the Galilee Mountains, Israel, respectively. The xeric, south-facing, ‘African’ slopes and the mesic, north-facing, ‘European’ slopes have dramatically different microclimates and plant communities. Moreover, ‘Evolution Canyon II’ receives more rainfall than ‘Evolution Canyon I.’

Methodology/Principal Findings

We examined fluctuating asymmetry, rate of whorl expansion, shell height, and number of rotations of the body suture in six species of terrestrial snails from the two ‘Evolution Canyons.’ The xeric ‘African’ slope should be more stressful to land snails than the ‘European’ slope, and ‘Evolution Canyon I’ should be more stressful than ‘Evolution Canyon II.’ Only Eopolita protensa jebusitica showed marginally significant differences in fluctuating helical asymmetry between the two slopes. Contrary to expectations, asymmetry was marginally greater on the ‘European’ slope. Shells of Levantina spiriplana caesareana at ‘Evolution Canyon I,’ were smaller and more asymmetric than those at ‘Evolution Canyon II.’ Moreover, shell height and number of rotations of the suture were greater on the north-facing slopes of both canyons.

Conclusions/Significance

Our data is consistent with a trade-off between drought resistance and thermoregulation in snails; Levantina was significantly smaller on the ‘African’ slope, for increasing surface area and thermoregulation, while Eopolita was larger on the ‘African’ slope, for reducing water evaporation. In addition, ‘Evolution Canyon I’ was more stressful than Evolution Canyon II’ for Levantina.     

Adaptive differentiation of thermotolerance in Drosophila along a microclimatic gradient

We examined whether a remarkable occurrence - the physiological evolution of two Drosophila melanogaster populations, despite a spatial separation of only 100-400 m, was idiosyncratic and temporary, or persisted over multiple years. We ascertained the high-temperature tolerance of Drosophila descended from populations on the north-facing slope (NFS) and south-facing slope (SFS) of 'Evolution Canyon' (Lower Nahal Oren, Mt Carmel, Israel), which were collected in 1997, 1999, and 2000. Results for these Drosophila uniformly resembled other studies in many respects: an inverse relationship between survival and heat-shock temperature, male-female differences in thermotolerance, and inducible thermotolerance. Importantly, for all years of collection, SFS flies consistently exceeded NFS flies in basal and inducible thermotolerance after diverse heat shocks, with and without thermal pretreatment, and whether isofemale lines, synthetic populations, or inbred lines were compared. Inbred lines, however, had lower thermotolerance than outbred lines. Several nonexclusive processes may explain the evolution of such physiological differentiation.     

Sequence polymorphism of candidate behavioural genes in Drosophila melanogaster flies from 'Evolution canyon'

This study focuses on the molecular features of three candidate behavioural genes in Drosophila melanogaster from the opposite slopes of Nahal Oren Canyon, Mount Carmel, Israel. These slopes display dramatic physical and biotic contrasts. Spatial variation of microclimatic conditions leads to adaptive differentiation and partial sexual isolation of populations, as suggested by our previous studies. The chosen candidate genes presumably contributing to genetic variation in sexual behaviour of Drosophila in the Canyon were desaturase, period, and no-on-transient A. These genes are known to include polymorphic repeated sequences, insertions/deletions, or nucleotide substitutions. The idea was that their polymorphism might be one of the determinants of behavioural peculiarities of flies derived from the opposite slopes. Indeed, interslope differences in the sequence encoding the (Thr-Gly)n repeat (exon 5) of the period gene were established, suggesting evolutionary functional importance. In particular, we unraveled variation in the length and composition of this region in different NFS (north-facing slope) and SFS (south-facing slope) lines. The 'European' allele (n = 20) was a 2.6-fold more abundant on the NFS compared to the SFS. This predominance probably gives some advantages for flies inhabiting wet and less warm conditions of the NFS. We suggest that repeat length/composition may influence the functional features of flies, i.e. habitat choice, nonrandom mating, and temperature adaptation. A series of female single-mate-choice tests show that females derived from NFS distinguish between males with specific per alleles (n = 17 vs. n = 20), as well as between males originated from the opposing slopes. Females from SFS were less discriminating and did not manifest significant deviation from random mating.     

Structure of variation in enzyme activity in natural Drosophila melanogaster populations

Enzyme activity variation was assessed in several isofemale lines originating from two Hungarian Drosophila melanogaster populations. Samples from each population were taken from from two villages; 8-9 isofemale lines were established from each village. The activities of ADH, alphaGPDH, IDH and 6PGDH were determined in the adults (in the F1 generation) and in the larvae (in the F3 generation) as well. Enzyme activities were measured on starch gel after electrophoresis. The activity of the enzyme was detected in a single individual and it was also possible to determine its genotype. The results showed that most of the variation occurred within sites for all four enzymes. This within site variation was more or less equally partitioned into within and between isofemale line (family) components. A smaller portion of variation was attributable to the differences between the populations. Nevertheless, adult alphaGPDH, and larval IDH and 6PGDH activities exhibited significant differences between the two populations. Variation in larval activities of all enzymes was higher than that of the adults, but 6PGDH had considerably higher variation in the adults. The greater variation in larval activities probably reflected the greater environmental variation in the microhabitat of the larvae compared to that of the adults. Larval activities of the investigated enzymes showed much stronger correlation than adult activities. The correlation pattern in the adults differed greatly between the two populations.     

Enzyme Activities and Alcohol Tolerance in Isofemale Lines of Drosophila melanogaster Originating from Different Habitats

Enzyme activity variation was studied in a Drosophila melanogaster population from two villages (Tiszafüred and Tiszaszolos) in Hungary. Two habitats (distillery and farmyard) were sampled in both villages and 8-9 isofemale lines were established from each sample with a total of 35 lines. The activities of ADH, alphaGPDH, IDH and 6PGDH were determined on starch gel after electrophoresis in 10 F1 females of each of the 35 isofemale lines. Three sublines were established from three selected isofemale lines of all four samples (altogether 36 sublines). Alcohol tolerance of the adult flies was assayed in these sublines. The activity of ADH was similar in the two habitats; so was the sensitivity to ethanol. Accordingly, no differences in adaptation to environmental ethanol were detected between the two habitats. The deviations between the two habitats in average activities and in the total variation of enzyme activities were not consistent in the two villages. These results suggest that founder effects and genetic drift are more pronounced in distilleries than selection. The association among enzyme activities varied greatly both between the two villages and between the two habitats. The two parameters of alcohol tolerance were not significantly different between the two habitats in any of the two villages.     

The Midas cichlid species complex: incipient sympatric speciation in Nicaraguan cichlid fishes?

Abstract Sympatric speciation is a contentious concept, although theoretical models as well as empirical evidence support its relevance in evolutionary biology. The Midas cichlid species complex (Amphilophus citrinellus, labiatus, zaliosus) from several crater lakes in Nicaragua fits several of the key characteristics of a sympatric speciation model. In particular, in A. citrinellus (i) strong assortative mating on the basis of colour polymorphism and (ii) ecological differentiation based on morphological polymorphisms involving the feeding apparatus and body shape might both be mechanisms of incipient speciation. Seven microsatellite markers and mtDNA control region sequences [836 base pairs (bp)] were used to study the population genetic structure of 519 specimens of Midas cichlid populations from the two Great Lakes Managua and Nicaragua, and three crater lakes in Nicaragua, Central America. The three named species of the species complex occupy different ecological niches, are morphologically distinct and can be distinguished genetically. We uncovered allopatric genetic differentiation of populations of A. citrinellus from different lakes and distant locations within Lake Managua and, more interestingly, incipient genetic differentiation of several sympatric populations based on colouration (in A. citrinellus and A. labiatus) but not on the morphology of the pharyngeal jaws (in A. citrinellus). Sexual selection and assortative mating might be the driven forces of diversification within named species. The Midas cichlid species complex in Nicaragua is an excellent model system for the study of the incipient stages of adaptation, speciation and the formation of species flocks.     

The endosymbiont Wolbachia in Eurasian populations of Drosophila melanogaster

The endosymbiotic [alpha]-proteobacteria Wolbachia is widely spread among arthropods and Filariidae nematodes. This bacterium is transmitted vertically via a transovarian route. Wolbachia is a cause of several reproductive abnormalities in the host species. We analyzed the isofemale lines created using flies collected from Drosophila melanogaster natural populations for infection with the endosymbiont Wolbachia. Wolbachia were genotyped according to five variable markers: the presence of insertion sequence IS5 in two loci, the copy number of two minisatellite repeats, and an inversion. Overall, 665 isofemale lines isolated from the populations of D. melanogaster from Ukraine, Belarus, Moldova, Caucasus, Central Asia, Ural, Udmurtia, Altai, West and East Siberia, and Far East in 1974 through 2005 were used in the work. The samples from Ukrainian, Altaian, and Middle Asian populations were largest. The infection rate of D. melanogaster populations from Middle Asia, Altaian, and Eastern Europe (Ukraine, Moldavia, and Belarus) with Wolbachia amounted to 64, 56, and 39%, respectively. The D. melanogaster population from the Caucasus displayed heterogeneity in the genotypes of this cytoplasmic infection. The Wolbachia genotype wMel, detected in all the populations studied, was the most abundant. The genotype wMelCS2 was always present in the populations from Middle Asia and Altai and was among the rare variants in the D. melanogaster populations from the Eastern Europe. Single instances of the Wolbachia genotype wMelCS occurred in a few flies from the Central Asian and Altai populations, but was not found this genotype in the other regions.     

Positive selection drives the evolution of the Acp29AB accessory gland protein in Drosophila.

Nucleotide sequence variation at the Acp29AB gene region has been surveyed in Drosophila melanogaster from Spain (12 lines), Ivory Coast (14 lines), and Malawi (13 lines) and in one line of D. simulans. The approximately 1.7-kb region studied encompasses the Acp29AB gene that codes for a male accessory gland protein and its flanking regions. Seventy-seven nucleotide and 8 length polymorphisms were detected. Nonsynonymous polymorphism was an order of magnitude lower than synonymous polymorphism, but still high relative to other non-sex-related genes. In D. melanogaster variation at this region revealed no major genetic differentiation between East and West African populations, while differentiation was highly significant between the European and the two African populations. Comparison of polymorphism and divergence at synonymous and nonsynonymous sites showed an excess of fixed nonsynonymous changes, which indicates that the evolution of the Acp29AB protein has been driven by directional selection at least after the split of the D. melanogaster and D. simulans lineages. The pattern of variation in extant populations of D. melanogaster favors a scenario where the fixation of advantageous replacement substitutions occurred in the early stages of speciation and balancing selection is maintaining variation in this species.     

Mitochondrial DNA variation and genetic structure in populations of Drosophila melanogaster

The understanding of the genetic structure of a species can be improved by considering together data from different types of genetic markers. In the past, a number of worldwide populations of Drosophila melanogaster have been extensively studied for several such markers, including allozymes, chromosomal inversions, and quantitative characters. Here we present results from a study of restriction- fragment-length polymorphisms of mitochondrial DNA (mtDNA) in 92 isofemale lines from many of the same geographic populations of D. melanogaster. Eleven restriction enzymes were used, of which four revealed restriction-site polymorphism. A total of 24 different haplotypes were observed, of which 18 were unique to single populations. In many populations, the unique haplotypes have reached high frequency without being observed in neighboring populations. A Wagner parsimony tree reveals that mutationally close variants show geographical clumping, suggesting local differentiation of mtDNA in populations. The Old-World and the New-World populations are differentiated, with the predominant Old-World haplotype being virtually absent from the New World. These results contrast with those for the nuclear genes, in which many loci show parallel clines in different continents, and suggest a common origin of D. melanogaster populations in North America.      

Evolvability of Hsp70 expression under artificial election for inducible thermotolerance in independent populations of Drosophila melanogaster

To test whether expression of the inducible heat-shock protein Hsp70 increases under selection for inducible thermotolerance in Drosophila melanogaster, we performed artificial selection on replicate sets of Drosophila lines founded from two independent populations. Selection entailed pretreatment at 36 degrees C to induce thermotolerance and Hsp70 expression, followed by a more severe heat shock, whose temperature varied between sexes and among generations to achieve 50% mortality. Inducible thermotolerance increased slowly and continuously in selected lines and was 37%-50% greater than in controls after 10-11 generations. Lines founded from the two populations differed in their coevolution of Hsp70 expression. In lines founded from Evolution Canyon, Israel, Hsp70 level initially increased and thereafter was unchanged; replicate lines exhibited two temporal patterns of response to selection. In lines founded from Australia, Hsp70 levels increased throughout selection. In both cases, however, the increase in Hsp70 level averaged only 15%, suggesting that pleiotropy in Hsp70 function constrains evolutionary increase in its expression.     

Forty years of solitude: life-history divergence and behavioural isolation between laboratory lines of Drosophila melanogaster

The study of the early stages of speciation can benefit from examination of differences between populations of known history that have been separated for a short time, such as a few thousands of generations. We asked whether two lines of Drosophila melanogaster that were isolated more than 40 years ago have evolved differences in life-history characters, or have begun to evolve behavioural or postzygotic isolation. One line, which is resistant to DDT, showed lower egg production and a shorter lifespan than a susceptible line. These differences are not a pleiotropic effect of resistance because they are not attributable to the chromosome that contains the resistance factors. The two lines have begun to become behaviourally isolated. Again, the isolation is not attributable to genes on the chromosome that contains resistance factors. The lines show only prezygotic isolation; there is no evidence of reduced fitness of F1 or F2 hybrids. These lines and others like them, should be excellent subjects for analyses of genetic changes that could lead to speciation.     

Adaptive microclimatic evolution of the dehydrin 6 gene in wild barley at ��Evolution Canyon��, Israel

Dehydrins are one of the major stress-induced gene families, and the expression of dehydrin 6 (Dhn6) is strictly related to drought in barley. In order to investigate how the evolution of the Dhn6 gene is associated with adaptation to environmental changes, we examined 48 genotypes of wild barley, Hordeum spontaneum, from "Evolution Canyon" at Mount Carmel, Israel. The Dhn6 sequences of the 48 genotypes were identified, and a recent insertion of 342?bp at 5'UTR was found in the sequences of 11 genotypes. Both nucleotide and haplotype diversity of single nucleotide polymorphism in Dhn6 coding regions were higher on the AS ("African" slope or dry slope) than on the ES ("European" slope or humid slope), and the applied Tajima D and Fu-Li test rejected neutrality of SNP diversity. Expression analysis indicated that the 342?bp insertion at 5'UTR was associated with the earlier up-regulation of Dhn6 after dehydration. The genetic divergence of amino acids sequences indicated significant positive selection of Dhn6 among the wild barley populations. The diversity of Dhn6 in microclimatic divergence slopes suggested that Dhn6 has been subjected to natural selection and adaptively associated with drought resistance of wild barley at "Evolution Canyon".     

Intrinsic reproductive isolation between Trinidadian populations of the guppy, Poecilia reticulata

Although Trinidadian populations of the guppy, Poecilia reticulata, show considerable adaptive genetic differentiation, they have been assumed to show little or no reproductive isolation. We tested this assumption by crossing Caroni (Tacarigua River) and Oropuche (Oropuche R.) drainage populations from Trinidad's Northern Range, and by examining multiple aspects of reproductive compatibility in the F1, F2 and BC1 generations. In open-aquarium experiments, F1 males performed fewer numbers of mating behaviours relative to parental population controls. This is the first documentation of hybrid behavioural sterility within a species, and it suggests that such sterility may feasibly be involved in causing speciation. The crosses also uncovered hybrid breakdown for embryo viability, brood size and sperm counts. In contrast, no reductions in female fertility were detected, indicating that guppies obey Haldane's rule for sterility. Intrinsic isolation currently presents a much stronger obstacle to gene flow than behavioural isolation, and our results indicate that Trinidadian populations constitute a useful model for investigating incipient speciation.     

Divergence of larval resource acquisition for water conservation and starvation resistance in Drosophila melanogaster

Laboratory selection experiments have evidenced storage of energy metabolites in adult flies of desiccation and starvation resistant strains of D. melanogaster but resource acquisition during larval stages has received lesser attention. For wild populations of D. melanogaster, it is not clear whether larvae acquire similar or different energy metabolites for desiccation and starvation resistance. We tested the hypothesis whether larval acquisition of energy metabolites is consistent with divergence of desiccation and starvation resistance in darker and lighter isofemale lines of D. melanogaster. Our results are interesting in several respects. First, we found contrasting patterns of larval resource acquisition, i.e., accumulation of higher carbohydrates during 3rd instar larval stage of darker flies versus higher levels of triglycerides in 1st and 2nd larval instars of lighter flies. Second, 3rd instar larvae of darker flies showed ~40?h longer duration of development at 21°C; and greater accumulation of carbohydrates (trehalose and glycogen) in fed larvae as compared with larvae non-fed after 150?h of egg laying. Third, darker isofemale lines have shown significant increase in total water content (18%); hemolymph (86%) and dehydration tolerance (11%) as compared to lighter isofemale lines. Loss of hemolymph water under desiccation stress until death was significantly higher in darker as compared to lighter isofemale lines but tissue water loss was similar. Fourth, for larvae of darker flies, about 65% energy content is contributed by carbohydrates for conferring greater desiccation resistance while the larvae of lighter flies acquire 2/3 energy from lipids for sustaining starvation resistance; and such energy differences persist in the newly eclosed flies. Thus, larval stages of wild-caught darker and lighter flies have evolved independent physiological processes for the accumulation of energy metabolites to cope with desiccation or starvation stress.     

Microsatellite variation among diverging populations of Drosophila mojavensis

Divergence and speciation may occur by various means, depending on the particular history, selective environments, and genetic composition of populations. In Drosophila mojavensis, a good model of incipient speciation, understanding the population genetic structure within this group facilitates our ability to understand the context in which reproductive isolation among populations is developing. Here we report the genetic structure and relationships of D. mojavensis populations at nuclear loci. We surveyed 29 populations throughout the distribution of D. mojavensis for four microsatellite loci to differentiation among populations of this species. These loci reveal four distinct geographical regions of D. mojavensis populations in the south-western United States and north-western Mexico--(i) Baja California peninsula (Baja), (ii) Sonora, Mexico-southern Arizona, United States (Sonora), (iii) Mojave Desert and Grand Canyon (Mojave), and (iv) Santa Catalina Island (Catalina). While all regions show strong isolation, Mojave and Catalina are highly diverged from other regions. Within any region, populations are largely homogenous over broad geographical distances. Based on the population structure, we find clear geographical barriers to gene flow appear to have a strong effect in isolating populations across regions for this species.     

Sexual isolation and extreme morphological divergence in the Cumana guppy: a possible case of incipient speciation

Theory predicts that sexual selection can promote the evolution of reproductive isolation and speciation. Those cases in which sexual selection has led to speciation should be characterized by significant differentiation in male display traits and correlated female preferences in the absence of post-zygotic isolation, accompanied by little genetic or other morphological differentiation. Previous evidence indicates that a cluster of populations of the guppy (Poecilia reticulata Peters) from Cumana, Venezuela, the 'Cumana guppy', differs significantly in female preferences from a nearby guppy population (A. Lindholm & F. Breden, Am. Nat., 160: 2002, S214). Here, we further document sexual isolation between these populations. In addition, these populations exhibit significant divergence in male display traits correlated to differences in between-population mating success, little mitochondrial genetic differentiation, and we find no evidence for genetic incompatibility between a Cumana population and several geographically isolated populations. These results suggest that divergent sexual selection has contributed to differentiation of the Cumana guppy, and this may be the first example of incipient speciation in the guppy.     

Vip3Aa tolerance response of Helicoverpa armigera populations from a Cry1Ac cotton planting region

Transgenic cotton, Gossypium hirsutum L., that expresses the Bacillus thuringiensis (Bt) Cry1Ac toxin, holds great promise in controlling target insect pests. Evolution of resistance by target pests is the primary threat to the continued efficacy of Bt cotton. To thwart pest resistance evolution, a transgenic cotton culitvar that produces two different Bt toxins, cry1Ac and vip3A genes, was proposed as a successor of cry1Ac cotton. This article reports on levels of Vip3Aa tolerance in Helicoverpa armigera (Hübner) (Lepidoptera: Noctuidae) populations from the Cry1Ac cotton planting region in China based on bioassays of the F1 generation of isofemale lines. In total, 80 isofemale families of H. armigera from Xiajin county of Shandong Province (an intensive Bt cotton planting area) and 93 families from Anci county of Hebei Province (a multiple-crop system including corn [Zea mays L.] , soybean [Glycine max (L.) Merr.], peanut (Arachis hypogaea L.), and Bt cotton) were screened with a discriminating concentration of both Cry1Ac- and Vip3A-containing diets in 2009. From data on the relative average development rates and percentage of larval weight inhibition of F1 full-sib families tested simultaneously on Cry1Ac and Vip3Aa, results indicate that responses to Cry1Ac and Vip3Aa were not genetically correlated in field population ofH. armigera. This indicates that the threat of cross-resistance between Cry1Ac and Vip3A is low in field populations of H. armigera. Thus, the introduction of Vip3Aa/Cry1Ac-producing lines could delay resistance evolution in H. armigera in Bt cotton planting area of China.