Stratospheric temperature trends: impact of ozone variability and the QBO

In most climate simulations used by the Intergovernmental Panel on Climate Change 2007 fourth assessment report, stratospheric processes are only poorly represented. For example, climatological or simple specifications of time-varying ozone concentrations are imposed and the quasi-biennial oscillation (QBO) of equatorial stratospheric zonal wind is absent. Here we investigate the impact of an improved stratospheric representation using two sets of perturbed simulations with the Hadley Centre coupled ocean atmosphere model HadGEM1 with natural and anthropogenic forcings for the 1979–2003 period. In the first set of simulations, the usual zonal mean ozone climatology with superimposed trends is replaced with a time series of observed zonal mean ozone distributions that includes interannual variability associated with the solar cycle, QBO and volcanic eruptions. In addition to this, the second set of perturbed simulations includes a scheme in which the stratospheric zonal wind in the tropics is relaxed to appropriate zonal mean values obtained from the ERA-40 re-analysis, thus forcing a QBO. Both of these changes are applied strictly to the stratosphere only. The improved ozone field results in an improved simulation of the stepwise temperature transitions observed in the lower stratosphere in the aftermath of the two major recent volcanic eruptions. The contribution of the solar cycle signal in the ozone field to this improved representation of the stepwise cooling is discussed. The improved ozone field and also the QBO result in an improved simulation of observed trends, both globally and at tropical latitudes. The Eulerian upwelling in the lower stratosphere in the equatorial region is enhanced by the improved ozone field and is affected by the QBO relaxation, yet neither induces a significant change in the upwelling trend.

Electronic marking of statistics assessments for bioscience students

We describe the main features of a program written to perform electronic marking of quantitative or simple text questions. One of the main benefits is that it can check answers for being consistent with earlier errors, so can cope with a range of numerical questions. We summarise our experience of using it in a statistics course taught to 200 bioscience students.

Selection of candidate coding DNA barcoding regions for use on land plants

An in silico screen of 41 of the 81 coding regions of the Nicotiana plastid genome generated a shortlist of 12 candidates as DNA barcoding loci for land plants. These loci were evaluated for amplification and sequence variation against a reference set of 98 land plant taxa. The deployment of multiple primers and a modified multiplexed tandem polymerase chain reaction yielded 85–94% amplification across taxa, and mean sequence differences between sister taxa of 6.1 from 156 bases of accD to 22 from 493 bases of matK. We conclude that loci should be combined for effective diagnosis, and recommend further investigation of the following six loci: matK, rpoB, rpoC1, ndhJ, ycf5 and accD.

Calculating the exclusion probability and paternity index for X-chromosomal loci in the presence of substructure

There has been recent interest in the use of X-chromosomal loci for forensic and relatedness testing casework, with many authors developing new X-linked short tandem repeat (STR) loci suitable for forensic use. Here we present formulae for two key quantities in paternity testing, the average probability of exclusion and the paternity index, which are suitable for Xchromosomal loci in the presence of population substructure.

Gene conversion and evolution of Xq28 duplicons involved in recurring inversions causing severe hemophilia A

Inversions breaking the 1041 bp int1h-1 or the 9.5-kb int22h-1 sequence of the F8 gene cause hemophilia A in 1/30,000 males. These inversions are due to homologous recombination between the above sequences and their inverted copies on the same DNA molecule, respectively, int1h-2 and int22h-2 or int22h-3. We find that (1) int1h and int22h duplicated more than 25 million years ago; (2) the identity of the copies (>99%) of these sequences in humans and other primates is due to gene conversion; (3) gene conversion is most frequent in the internal regions of int22h; (4) breakpoints of int22h-related inversions also tend to involve the internal regions of int22h; (5) sequence variations in a sample of human X chromosomes defined eight haplotypes of int22h-1 and 27 of int22h-2 plus int22h-3; (6) the latter two sequences, which lie, respectively, 500 and 600 kb telomeric to int22h-1 are five-fold more identical when in cis than when in trans, thus suggesting that gene conversion may be predominantly intrachromosomal; (7) int1h, int22h, and flanking sequences evolved at a rate of about 0.1% substitutions per million years during the divergence between humans and other primates, except for int1h during the human-chimpanzee divergence, when its rate of evolution was significantly lower. This is reminiscent of the slower evolution of palindrome arms in the male specific regions of the Y chromosome and we propose, as an explanation, that intrachromosomal gene conversion and cosegregation of the duplicated regions favors retention of the ancestral sequence and thus reduces the evolution rate.

Population genetic models can be used to study the evolution of the interacting behaviours of parents and their progeny

An example of the evolution of the interacting behaviours of parents and progeny is studied using iterative equations linking the frequencies of the gametes produced by the progeny to the frequencies of the gametes in the parental generation. This population genetics approach shows that a model in which both behaviours are determined by a single locus can lead to a stable equilibrium in which the two behaviours continue to segregate. A model in which the behaviours are determined by genes at two separate loci leads eventually to fixation of the alleles at both loci but this can take many generations of selection. Models of the type described in this paper will be needed to understand the evolution of complex behaviour when genomic or experimental information is available about the genetic determinants of behaviour and the selective values of different genomes. (c) 2007 Elsevier Inc. All rights reserved.

Effect of viral load on the outcome of herpes zoster

Varicella-zoster virus (VZV) is a member of the Herpesviridae family, primary infection with which causes varicella, more commonly known as chicken pox. Characteristic of members of the alphaherpesvirus subfamily, VZV is neurotropic and establishes latency in sensory neurons. Reactivation of VZV causes herpes zoster, also known as shingles. The most frequent complication following zoster is chronic and often debilitating pain called postherpetic neuralgia (PHN), which can last for months after the disappearance of a rash. During episodes of acute zoster, VZV viremia occurs in some, but not all, patients; however, the effect of the viral load on the disease outcome is not known. Here we describe the development of a highly specific, sensitive, and reproducible real-time PCR assay to investigate the factors that may contribute to the presence and levels of baseline viremia in patients with zoster and to determine the relationship between viremia and the development and persistence of PHN. VZV DNA was detected in the peripheral blood mononuclear cells (PBMCs) of 78% of patients with acute zoster and in 9% of healthy asymptomatic blood donors. The presence of VZV in the PBMCs of patients with acute zoster was independently associated with age and being on antivirals but not with gender, immune status, extent of rash, the age of the rash at the time of blood sampling, having a history of prodromal pain, or the extent of acute pain. Prodromal pain was significantly associated with higher baseline viral loads. Viral load levels were not associated with the development or persistence of PHN at 6, 12, or 26 weeks.

The expected performance of single nucleotide polymorphism loci in paternity testing

We discuss the utility of single nucleotide polymorphism loci for full trio and mother-unavailable paternity testing cases, in the presence of population substructure and relatedness of putative and actual fathers. We focus primarily on the expected number of loci required to gain specified probabilities of mismatches, and report the expected proportion of paternity indices greater than three threshold values for these loci. (c) 2004 Elsevier Ireland Ltd. All rights reserved.

Detecting non-multiplicative genotype relative risks from transmissions of parental alleles to affected children

The differential transmission of alleles from parents to affected children indicates that the locus under investigation is either directly involved in the occurrence of the disease or that there are allelic associations with other loci that are directly involved. Conditional logistic regression applied to a diallelic locus leads to a test with two degrees of freedom. The power of a single degree of freedom test to detect non-multiplicative allelic effects is discussed here.

API-CALC 1.0: a computer program for calculating the average probability of identity allowing for substructure, inbreeding and the presence of close relatives

Individual identification via DNA profiling is important in molecular ecology, particularly in the case of noninvasive sampling. A key quantity in determining the number of loci required is the probability of identity (PIave), the probability of observing two copies of any profile in the population. Previously this has been calculated assuming no inbreeding or population structure. Here we introduce formulae that account for these factors, whilst also accounting for relatedness structure in the population. These formulae are implemented in API-CALC 1.0, which calculates PIave for either a specified value, or a range of values, for F-IS and F-ST.

DNA variation in a 13-Mb region including the F9 gene: inferring the genealogical history and causal role of a hemophilia B mutation (IVS 5+13 A -> G)

About 5.5% of all UK hemophilia B patients have the base substitution IVS 5+13 A-->G as the only change in their factor (F)IX gene (F9). This generates a novel donor splice site which fits the consensus better than the normal intron 5 donor splice. Use of the novel splice site should result in a missense mutation followed by the abnormal addition of four amino acids to the patients' FIX. In order to explain the prevalence of this mutation, its genealogical history is examined. Analysis of restriction fragment length polymorphism in the 21 reference UK individuals (from different families) with the above mutation showed identical haplotypes in 19 while two differed from the rest and from each other. In order to investigate the history of the mutation and to verify that it had occurred independently more than once, the sequence variation in 1.5-kb segments scattered over a 13-Mb region including F9 was examined in 18 patients and 15 controls. This variation was then analyzed with a recently developed Bayesian approach that reconstructs the genealogy of the gene investigated while providing evidence of independent mutations that contribute disconnected branches to the genealogical tree. The method also provides minimum estimates of the age of the mutation inherited by the members of coherent trees. This revealed that 17 or 18 mutant genes descend from a founder who probably lived 450 years ago, while one patient carries an independent mutation. The independent recurrence of the IVS5+13 A-->G mutation strongly supports the conclusion that it is the cause of these patients' mild hemophilia.

Linkage disequilibrium assessment via log-linear modeling of SNP haplotype frequencies

Analyses of high-density single-nucleotide polymorphism (SNP) data, such as genetic mapping and linkage disequilibrium (LD) studies, require phase-known haplotypes to allow for the correlation between tightly linked loci. However, current SNP genotyping technology cannot determine phase, which must be inferred statistically. In this paper, we present a new Bayesian Markov chain Monte Carlo (MCMC) algorithm for population haplotype frequency estimation, particulary in the context of LD assessment. The novel feature of the method is the incorporation of a log-linear prior model for population haplotype frequencies. We present simulations to suggest that 1) the log-linear prior model is more appropriate than the standard coalescent process in the presence of recombination (>0.02cM between adjacent loci), and 2) there is substantial inflation in measures of LD obtained by a "two-stage" approach to the analysis by treating the "best" haplotype configuration as correct, without regard to uncertainty in the recombination process. Genet Epidemiol 25:106-114, 2003. (C) 2003 Wiley-Liss, Inc.