Spatial and temporal analysis of bottlenose dolphin strandings in South Carolina, 1992-2005

This CD contains summary data of bottlenose dolphins stranded in South Carolina using a Geographical Information System (GIS) and contains two published manuscripts in .pdf files. The intent of this CD is to provide data on bottlenose dolphin strandings in South Carolina to marine mammal researchers and managers. This CD is an accumulation of 14 years of stranding data collected through the collaborations of the National Ocean Service, Center for Coastal Environmental Health and Biomolecular Research (CCEHBR), the South Carolina Department of Natural Resources, and numerous volunteers and veterinarians that comprised the South Carolina Marine Mammal Stranding Network. Spatial and temporal information can be visually represented on maps using GIS. For this CD, maps were created to show relationships of stranding densities with land use, human population density, human interaction with dolphins, high geographical regions of live strandings, and seasonal changes. Point maps were also created to show individual strandings within South Carolina. In summary, spatial analysis revealed higher densities of bottlenose dolphin strandings in Charleston and Beaufort Counties, which consist of urban land with agricultural input. This trend was positively correlated with higher human population levels in these coastal counties as compared with other coastal counties. However, spatial analysis revealed that certain areas within a county may have low human population levels but high stranding density, suggesting that the level of effort to respond to strandings is not necessarily positively correlated with the density of strandings in South Carolina. Temporal analysis revealed a significantly higher density of bottlenose dolphin strandings in the northern portion of the State in the fall, mostly due to an increase of neonate strandings. On a finer geographic scale, seasonal stranding densities may fluctuate depending on the region of interest. Charleston Harbor had the highest density of live bottlenose dolphin strandings compared to the rest of the State. This was due in large part to the number of live dolphin entanglements in the crab pot fishery, the largest source of fishery-related mortality for bottlenose dolphins in South Carolina (Burdett and McFee 2004). Spatial density calculations also revealed that Charleston and Beaufort accounted for the majority of dolphins that were involved with human activities. 1

Assessment of cattle genetic introgression into domestic yak populations using mitochondrial and microsatellite DNA markers

Hybridization between yak Poephagus grunniens and taurine Bos taurus or indicine B. indicus cattle has been widely practiced throughout the yak geographical range, and gene flow is expected to have occurred between these species. To assess the impact of cattle admixture on domestic yak, we examined 1076 domestic yak from 29 populations collected in China, Bhutan, Nepal, India, Pakistan, Kyrgyzstan, Mongolia and Russia using mitochondrial DNA and 17 autosomal microsatellite loci. A cattle diagnostic marker-based analysis reveals cattle-specific mtDNA and/or autosomal microsatellite allele introgression in 127 yak individuals from 22 populations. The mean level of cattle admixture across the populations, calculated using allelic information at 17 autosomal microsatellite loci, remains relatively low (mY(cattle) = 2.66 +/- 0.53% and Q(cattle) = 0.69 +/- 2.58%), although it varies a lot across populations as well as among individuals within population. Although the level of cattle admixture shows a clear geographical structure, with higher levels of admixture in the Qinghai-Tibetan Plateau and Mongolian and Russian regions, and lower levels in the Himalayan and Pamir Plateau region, our results indicate that the level of cattle admixture is not significantly correlated with the altitude across geographical regions as well as within geographical region. Although yak-cattle hybridization is primarily driven to produce F-1 hybrids, our results show that the subsequent gene flow between yak and cattle took place and has affected contemporary genetic make-up of domestic yak. To protect yak genetic integrity, hybridization between yak and cattle should be tightly controlled.

Assessing seasonal precipitation trends in India using parametric and non-parametric statistical techniques

This paper provides an insight into the long-term trends of the four seasonal and annual precipitations in various climatological regions and sub-regions in India. The trends were useful to investigate whether Indian seasonal rainfall is changing in terms of magnitude or location-wise. Trends were assessed over the period of 1954-2003 using parametric ordinary least square fits and non-parametric Mann-Kendall technique. The trend significance was tested at the 95% confidence level. Apart from the trends for individual climatological regions in India and the average for the whole of India, trends were also specifically determined for the possible smaller geographical areas in order to understand how different the trends would be from the bigger spatial scales. The smaller geographical regions consist of the whole southwestern continental state of Kerala. It was shown that there are decreasing trends in the spring and monsoon rainfall and increasing trends in the autumn and winter rainfalls. These changes are not always homogeneous over various regions, even in the very short scales implying a careful regional analysis would be necessary for drawing conclusions regarding agro-ecological or other local projects requiring change in rainfall information. Furthermore, the differences between the trend magnitudes and directions from the two different methods are significantly small and fall well within the significance limit for all the cases investigated in Indian regions (except where noted). © 2010 Springer-Verlag.

Molecular phylogeography and population structure of a mid-elevation montane frog Leptobrachium ailaonicum in a fragmented habitat of southwest China

Leptobrachium ailaonicum is a vulnerable anuran restricted to a patchy distribution associated with small mountain streams surrounded by forested slopes at mid-elevations (approximately 2000-2600 m) in the subtropical Mount Wuliang and Mount Ailao ranges in southwest China (Yunnan Province) and northern Vietnam. Given high habitat specificity and lack of suitable habitat in lower elevations between these ranges, we hypothesized limited gene flow between populations throughout its range. We used two mitochondrial genes to construct a phylogeographic pattern within this species in order to test our hypothesis. We also examined whether this phylogeographic pattern is a response to past geological events and/or climatic oscillations. A total of 1989 base pairs were obtained from 81 individuals of nine populations yielding 51 unique haplotypes. Both Bayesian and maximum parsimony phylogenetic analyses revealed four deeply divergent and reciprocally monophyletic mtDNA lineages that approximately correspond to four geographical regions separated by deep river valleys. These results suggest a long history of allopatric separation by vicariance. The distinct geographic distributions of four major clades and the estimated divergence time suggest spatial and temporal separations that coincide with climatic and paleogeographic changes following the orogeny and uplift of Mount Ailao during the late Miocene to mid Pliocene in southwest China. At the southern distribution, the presence of two sympatric yet differentiated clades in two areas are interpreted as a result of secondary contact between previously allopatric populations during cooler Pleistocene glacial cycles. Analysis of molecular variance indicates that most of the observed genetic variation occurs among the four regions implying long-term interruption of maternal gene flow, suggesting that L ailaonicum may represent more than one distinct species and should at least be separated into four management units corresponding to these four geographic lineages for conservation. (C) 2009 Elsevier Inc. All rights reserved.

不同小麦品质类型与生态环境关系的研究

小麦（Triticum aestivum L.）是世界上种植面积最大，总产量最高，食物加工种类最丰富的粮食作物，占世界人口35 ％－40 ％的人们以此为主要食物。因此小麦产量的高低和品质的优劣直接影响人们对食物需求的安全和满意程度，也影响着人类的营养平衡以及面粉和食品加工业的发展。随着生活水平的提高，人们对于小麦的品质越来越重视。培育优质专用小麦新品种，制定优质专用小麦品种品质生态区划，从而在不同程度上实现小麦的区域化种植和产业化经营具有重要的意义。 影响小麦品质的因素主要是遗传因素和环境因素，其中环境因素又包括各种自然生态因素和人为因素。研究表明，小麦品质的环境间的差异大于品种间的差异，气候条件是影响小麦品质的最重要的因子，小麦品质的地域间的差异反映出了小麦的品质区域分布规律。为了满足市场对不同品质小麦的需求，对小麦进行区域化研究具有重要的理论和现实意义。本研究结合四川的地理、气候特点，研究不同品质类型与生态环境的关系，为在复杂的生态环境内进行品质区划提供依据。 本研究首先根据四川省小麦种植区域的生态特点，在四川省多个典型生态区：川南丘陵的荣县、川西南高原的西昌、川西平原的双流布点种植，采用的小麦试验材料为不同品质类型：中筋小麦川育12、川育14、川育16由本所提供；弱筋小麦川麦32和强筋小麦川麦36由四川省农科院作物所提供。通过研究品质性状与品种及各个生态因子包括地点、土壤土质差异等的关系，明确不同生态环境中适宜种植的小麦品种类型，强筋小麦、中筋小麦更适合于在荣县、双流地区种植，弱筋小麦更适合于在西昌地区种植，为品种品质区划奠定基础。 其次，选择了本课题组育成的稳定中间品系，对其品质性状SDS沉降值进行了多年测定。分析了品质性状SDS沉降值与多种气候因子的相关性，结果表明SDS沉降值与日均温、日照时数成正相关，与降水量成负相关，为品质育种提供了理论依据。 此外，以中筋小麦新品种小麦川育14为材料，应用三元二次正交旋转回归模型设计试验，研究主要栽培因子播期、密度和施肥量对产量的影响，并建立函数模型。经计算机模拟寻优，筛选出了高产高效栽培组合措施，并确定了置信域。结合四川省不同的地理情况，在平原和丘陵地区分别进行实验，并各自建立了高产高效栽培组合措施，为川育14品种的推广提供了理论指导。 Wheat is one of the most important crops in the world. About 35%-40% people all over the world, take the wheat as their most important food. So the quality, as well as the quantity of the wheat makes a direct effect on people’s demands of food and their satisfaction. It also effects on human’s healthy, and the development of the Food processing industry. With the development of the living standard, people pay more attention to the quality of wheat. So, we set a special ecology zoning for wheat. It is significant to carry out planting the wheat in special zoning in varying degrees. The main factors affecting wheat quality are heredity and environment including many ecological factors and the factors in cultivation. As to the quality，the difference between ecology and cultivation is more important than the difference between special wheat. In so many factors, climate is the most important one. From the difference in quality between different zones，we can conclude the rule of distribution abort quality of wheat. Finding out the intersection of numerous wheat not only can meet the demand of food production，but also has important signification in theory and realism。In our research, according to the complex geography in Sichuan province, we study the relationship between numerous kinds of quality characters in wheat and the ecology. So, we can set a foundation for more research. In this research, firstly, we plant wheat in some typical ecological regions of SICHUAN province: RONGXIAN（south of SC）、XICHANG（south of SC）, SHUANGLIU(west of SC). The materials of the experiments: ChuanYu12, ChuanYu14, ChuanYu16(from our institute), Chuanmai32, Chuanmai36 (from the Chinese academy of agriculture sciences of Sichuan. Through the research on the relationship between the quality of wheat and those ecology factors, we can make a definition that which area is perfect matched with which kind of wheat. And it can satisfy the demand of people. Secondly, select many sorts of wheat from our research group. All of them are selected and bred more than 3 years(2003-2005). And we make every-year determination as well. We’ve gotten SDS value from those 9, and various data on factors of climate. We also got to know the relation ship between those numbers. Thirdly, use Chuanyu14 as material, the mathematical model of the relation between the production of wheat and main agricultural measures such as date, density and fertilizer. The model was established by association of three elements two return, rotate and regression. We set a suitable model and get a suitable method which can make high harvest. Based on various kinds of geographical regions in Sichuan province, we set different models which can be used in plain and hill. So, we can plant Chuan Yu 14 in Sichuan province under the result in research.

Population genetic structure of crimson snapper Lutjanus erythropterus in East Asia, revealed by analysis of the mitochondrial control region

The population genetic structure of the crimson snapper Lutjanus erythropterus in East Asia was examined with a 427-bp hypervariable portion of the mtDNA control region. A total of 262 samples were collected and 75 haplotypes were obtained. Neutrality tests (Tajima's and Fu's) suggested that Lutjanus erythropterus in East Asia had experienced a bottleneck followed by population expansion since the late Pleistocene. Despite the low phylogeographic structures in mtDNA haplotypes, a hierarchical examination of populations in 11 localities from four geographical regions using analysis of molecular variance (AMOVA) indicated significant genetic differentiation among regions (Phi(CT) = 0.08564, p < 0.01). Limited gene flow between the eastern region (including a locality in the western Pacific Ocean and two localities in the East Sea) and three geographic regions of the South China Sea largely contributed to the genetic subdivision. However, comparisons among three geographic regions of the South China Sea showed little to no genetic difference. Populations of Lutjanus erythropterus in East Asia are inferred to be divided into two major groups: an eastern group, including populations of the western Pacific Ocean and the East Sea, and a South China Sea group, consisting of populations from northern Malaysia to South China. The results suggest that fishery management should reflect the genetic differentiation and diversity in East Asia. (c) 2006 International Council for the Exploration of the Sea. Published by Elsevier Ltd. All rights reserved.

Distribution patterns of infection with multiple types of human papillomaviruses and their association with risk factors

Background: Infection with multiple types of human papillomavirus (HPV) is one of the main risk factors associated with the development of cervical lesions. In this study, cervical samples collected from 1, 810 women with diverse sociocultural backgrounds, who attended to their cervical screening program in different geographical regions of Colombia, were examined for the presence of cervical lesions and HPV by Papanicolau testing and DNA PCR detection, respectively. Principal Findings: The negative binomial distribution model used in this study showed differences between the observed and expected values within some risk factor categories analyzed. Particularly in the case of single infection and coinfection with more than 4 HPV types, observed frequencies were smaller than expected, while the number of women infected with 2 to 4 viral types were higher than expected. Data analysis according to a negative binomial regression showed an increase in the risk of acquiring more HPV types in women who were of indigenous ethnicity (+37.8%), while this risk decreased in women who had given birth more than 4 times (-31.1%), or were of mestizo (-24.6%) or black (-40.9%) ethnicity. Conclusions: According to a theoretical probability distribution, the observed number of women having either a single infection or more than 4 viral types was smaller than expected, while for those infected with 2-4 HPV types it was larger than expected. Taking into account that this study showed a higher HPV coinfection rate in the indigenous ethnicity, the role of underlying factors should be assessed in detail in future studies.

Midlatitude daytime D region ionosphere variations measured from radio atmospherics

We measured the midlatitude daytime ionospheric D region electron density profile height variations in July and August 2005 near Duke University by using radio atmospherics (or sferics for short), which are the high-power, broadband very low frequency (VLF) signals launched by lightning discharges. As expected, the measured daytime D region electron density profile heights showed temporal variations quantitatively correlated with solar zenith angle changes. In the midlatitude geographical regions near Duke University, the observed quiet time heights decreased from ∼80 km near sunrise to ∼71 km near noon when the solar zenith angle was minimum. The measured height quantitative dependence on the solar zenith angle was slightly different from the low-latitude measurement given in a previous work. We also observed unexpected spatial variations not linked to the solar zenith angle on some days, with 15% of days exhibiting regional differences larger than 0.5 km. In these 2 months, 14 days had sudden height drops caused by solar flare X-rays, with a minimum height of 63.4 km observed. The induced height change during a solar flare event was approximately proportional to the logarithm of the X-ray flux. In the long waveband (wavelength, 1-8 Å), an increase in flux by a factor of 10 resulted in 6.3 km decrease of the height at the flux peak time, nearly a perfect agreement with the previous measurement. During the rising and decaying phases of the solar flare, the height changes correlated more consistently with the short, rather than the long, wavelength X-ray flux changes. © 2010 by the American Geophysical Union.

Are Prorocentrum hoffmannianum and Prorocentrum belizeanum (Dinophyceae, Prorocentrales), the same species? An integration of morphological and molecular data.

The taxonomic assignment of Prorocentrum species is based on morphological characteristics; however, morphological variability has been found for several taxa isolated from different geographical regions. In this study, we evaluated species boundaries of Prorocentrum hoffmannianum and Prorocentrum belizeanum based on morphological and molecular data. A detailed morphological analysis was done, concentrating on the periflagellar architecture. Molecular analyses were performed on partial Small Sub-Unit (SSU) rDNA, partial Large Sub-Unit (LSU) rDNA, complete Internal Transcribed Spacer Regions (ITS1-5.8S-ITS2), and partial cytochrome b (cob) sequences. We concatenated the SSU-ITS-LSU fragments and constructed a phylogenetic tree using Bayesian Inference (BI) and maximum likelihood (ML) methods. Morphological analyses indicated that the main characters, such as cell size and number of depressions per valve, normally used to distinguish P. hoffmannianum from P. belizeanum, overlapped. No clear differences were found in the periflagellar area architecture. Prorocentrum hoffmannianum and P. belizeanum were a highly supported monophyletic clade separated into three subclades, which broadly corresponded to the sample collection regions. Subtle morphological overlaps found in cell shape, size, and ornamentation lead us to conclude that P. hoffmanianum and P. belizeanum might be considered conspecific. The molecular data analyses did not separate P. hoffmannianum and P. belizeanum into two morphospecies, and thus, we considered them to be the P. hoffmannianum species complex because their clades are separated by their geographic origin. These geographic and genetically distinct clades could be referred to as ribotypes: (A) Belize, (B) Florida-Cuba, (C1) India, and (C2) Australia.

Incidence trends for childhood type 1 diabetes in Europe during 1989-2003 and predicted new cases 2005-20: a multicentre prospective registration study

Background: The incidence of type 1 diabetes in children younger than 15 years is increasing. Prediction of future incidence of this disease will enable adequate fund allocation for delivery of care to be planned. We aimed to establish 15-year incidence trends for childhood type 1 diabetes in European centres, and thereby predict the future burden of childhood diabetes in Europe.
Methods: 20 population-based EURODIAB registers in 17 countries registered 29 311 new cases of type 1 diabetes, diagnosed in children before their 15th birthday during a 15-year period, 1989–2003. Age-specific log linear rates of increase were estimated in five geographical regions, and used in conjunction with published incidence rates and population projections to predict numbers of new cases throughout Europe in 2005, 2010, 2015, and 2020.
Findings: Ascertainment was better than 90% in most registers. All but two registers showed significant yearly increases in incidence, ranging from 0·6% to 9·3%. The overall annual increase was 3·9% (95% CI 3·6–4·2), and the increases in the age groups 0–4 years, 5–9 years, and 10–14 years were 5·4% (4·8–6·1), 4·3% (3·8–4·8), and 2·9% (2·5–3·3), respectively. The number of new cases in Europe in 2005, is estimated as 15 000, divided between the 0–4 year, 5–9 year, and 10–14 year age-groups in the ratio 24%, 35%, and 41%, respectively. In 2020, the predicted number of new cases is 24 000, with a doubling in numbers in children younger than 5 years and a more even distribution across age-groups than at present (29%, 37%, and 34%, respectively). Prevalence under age 15 years is predicted to rise from 94 000 in 2005, to 160 000 in 2020.
Interpretation: If present trends continue, doubling of new cases of type 1 diabetes in European children younger than 5 years is predicted between 2005 and 2020, and prevalent cases younger than 15 years will rise by 70%. Adequate health-care resources to meet these children’s needs should be made available.

Genetic diversity of liver flukes (Fasciola hepatica) from Eastern Europe.

The genetic diversity of liver fluke populations in three different countries from Eastern Europe (Greece, Bulgaria, and Poland) in comparison with available data from other countries was determined. Specifically, SNPs from regions of two nuclear genes, 28S rDNA, ß-tubulin 3 and an informative region of the mitochondrial genome were examined. Two major lineages for the 28S rDNA gene based on the highly polymorphic 105th nucleotide position were found. These lineages were widely and almost equally spread not only through the countries studied but also in other investigated geographical areas. Two basic lineages and additional haplotypes were defined for the mtDNA gene region, consisting of the cytochrome c oxidase subunit III gene, transfer RNA histidine gene and cytochome b gene. The basic lineages were observed within Greek, Bulgarian, and Polish Fasciola hepatica populations but the distribution of additional haplotypes differed between the populations from the three countries. For the ß-tubulin 3 gene multiple polymorphic sites were revealed but no explicit clades. The SNPs were spread unequally in all studied geographical regions with an evident distinction between the Greek and Polish specimens. Additional genotypes for the 28S rDNA region as well as haplotypes of the mtDNA region that were typical for the Greek or Polish populations were observed. Significant polymorphisms for ß-tubulin 3 gene were displayed with decreasing percentage of presence within populations from Greece to Poland. There was an amino acid substitution in ß-tubulin 3 protein found only among Polish specimens. It is hypothesized that genotypic differences between Greek, Bulgarian, and Polish liver fluke populations are due to territorial division and genetic drift in past epochs.

Hibernation and non-shivering thermogenesis in the Hottentot golden mole (Amblysomus hottentottus longiceps)

Although heterothermy (hibernation and torpor) is a common feature among mammals, there is debate over whether it is a derived or ancestral trait relative to endothermic homeothermy. Determination of the physiological characteristics of primitive mammals is central to understanding the evolution of endothermy. Moreover, evaluation of physiological mechanisms responsible for endothermic heat production [e.g. non-shivering thermogenesis (NST)] is key to understanding how early mammals responded to historical climate changes and colonised different geographical regions. Here we investigated the capacity for NST and heterothermy in the Hottentot golden mole, a basal eutherian mammal. NST was measured as the metabolic response to injections of noradrenalin and heterothermy by recording body temperature in free-ranging animals. We found that hibernation and torpor occurred and that the seasonal phenotypic adjustment of NST capacity was similar to that found in other placental mammals. Using phylogenetically independent contrasts, we compared measured values of NST with those obtained from the literature. This showed that all variation in NST was accounted for by differences in phylogeny and not zoogeography. These findings lend support to the observation that NST and heterothermy occur in the Afrotheria, the basal placental mammalian clade. Furthermore, this work suggests that heterothermy, rather than homeothermy is a plesiomorphic trait in mammals and supports the notion that NST mechanisms are phylogenetically ancient.

Using simulation to assess cardiac first-responder schemes exhibiting stochastic and spatial complexities

A Monte-Carlo simulation-based model has been constructed to assess a public health scheme involving mobile-volunteer cardiac First-Responders. The scheme being assessed aims to improve survival of Sudden-Cardiac-Arrest (SCA) patients, through reducing the time until administration of life-saving defibrillation treatment, with volunteers being paged to respond to possible SCA incidents alongside the Emergency Medical Services. The need for a model, for example, to assess the impact of the scheme in different geographical regions, was apparent upon collection of observational trial data (given it exhibited stochastic and spatial complexities). The simulation-based model developed has been validated and then used to assess the scheme's benefits in an alternative rural region (not a part of the original trial). These illustrative results conclude that the scheme may not be the most efficient use of National Health Service resources in this geographical region, thus demonstrating the importance and usefulness of simulation modelling in aiding decision making.

A unique homologue of the eukaryotic protein-modifier ubiquitin present in the bacterium Bacteroides fragilis, a predominant resident of the human gastro-intestinal tract

In the complete genome sequences of Bacteroides fragilis NCTC9343 and 638R, we have discovered a gene, ubb, the product of which has 63% identity to human ubiquitin and cross-reacts with antibodies raised against bovine ubiquitin. The sequence of ubb is closest in identity (76%) to the ubiquitin gene from a Migratory Grasshopper entomopoxvirus, suggesting acquisition by inter-kingdom horizontal gene transfer. We have screened clinical isolates of B. fragilis from diverse geographical regions and found that ubb is present in some, but not all strains. The gene is transcribed and the mRNA translated in B. fragilis, but deletion of ubb did not have a detrimental effect on growth. BfUbb has a predicted signal sequence; both full length and processed forms were detected in whole cell extracts, while the processed form was found in concentrated culture supernatants. Purified recombinant BfUbb inhibited in vitro ubiquitination and was able to covalently bind the human E1 activating enzyme, suggesting it could act as a suicide substrate in vivo. B. fragilis is one of the predominant members of the normal human resident gastro-intestinal microbiota with estimates up to >1011 cells g-1 of faeces by culture. These data indicate that the gastro-intestinal tract of some individuals could contain a significant amount of aberrant ubiquitin with the potential to inappropriately activate the host immune system and/or interfere with eukaryotic ubiquitin activity. This discovery could have profound implications in relation to our understanding of human diseases such as inflammatory bowel and autoimmune diseases.

The effects of new topical treatments on management of glaucoma in Scotland:An examination of ophthalmological health care

Background: The management of glaucoma has been changed in the past decade by the introduction of new drugs. The impact of these changes on clinical care of patients was examined by examining operation and prescribing rates for glaucoma in four geographical areas of Scotland for the years 1994 to 1999. Methods: A retrospective analysis of national health statistics: primary care prescribing data, hospital derived operation rates, consultant numbers, optometrist numbers, and eye test data, expressed by estimated population at risk of glaucoma. The outcome measures were prescribing volume and cost for glaucoma medications, and operation rates, corrected for population estimated to be at risk of glaucoma (PEG), for trabeculectomy, for Scotland as a whole, and for four geographical "regions" (north east, south east, central, and south west Scotland). Results: Prescribed items per 1000 population estimated to have glaucoma (PEG) increased by 24.9% between 1994 and 1999. This was above the general increase in prescribing in Scotland (17.8%). This increase varied in the four health regions evaluated (14.3% to 31.9%). Prescribing of topical ß blockers increased little (6.4%), but there was a large increase in the use of new products (topical prostaglandins, carbonic anhydrase inhibitors, and a agonists), at the expense of miotics (47.7% fall), and older sympathomimetics. This change in prescribing pattern was accompanied by a 61.5% increase in cost (range 42.2% to 73.4% in the four regions). New drugs accounted for more than half of total glaucoma expenditure in 1999. Operation rates (corrected for PEG) fell by 45.9% (range 43.1 to 58.6%) between 1994 and 1999. Other indicators suggested increased activity in ophthalmic areas (for example, cataract operations, eye tests, numbers of optometrists and ophthalmic surgeons all increased). Within north east Scotland operation rates decreased and prescribing increased less than in other regions, both from lowest regional baseline in 1994. Conclusions: The introduction of new drug classes has had dramatic effects on the prescribing of glaucoma treatments. There has been a decline in older treatments and an increase in new agents, which has been associated with a large reduction in operation rates for glaucoma in Scotland over 6 years. Comparison of prescribing and operation data indicates regional differences in healthcare delivery for glaucoma.