Statistical methods for examining genetic influences of resistance to anti-epileptic drugs

It is generally accepted that genetics may be an important factor in explaining the variation between patients’ responses to certain drugs. However, identification and confirmation of the responsible genetic variants is proving to be a challenge in many cases. A number of difficulties that maybe encountered in pursuit of these variants, such as non-replication of a true effect, population structure and selection bias, can be mitigated or at least reduced by appropriate statistical methodology. Another major statistical challenge facing pharmacogenetics studies is trying to detect possibly small polygenic effects using large volumes of genetic data, while controlling the number of false positive signals. Here we review statistical design and analysis options available for investigations of genetic resistance to anti-epileptic drugs.

Pesticides, preference heterogeneity and environmental taxes

In this paper we present results from two choice experiments (CE), designed to take account of the different negative externalities associated with pesticide use in agricultural production. For cereal production, the most probable impact of pesticide use is a reduction in environmental quality. For fruit and vegetable production, the negative externality is on consumer health. Using latent class models we find evidence of the presence of preference heterogeneity in addition to reasonably high willingness to pay (WTP) estimates for a reduction in the use of pesticides for both environmental quality and consumer health. To place our WTP estimates in a policy context we convert them into an equivalent pesticide tax by type of externality. Our tax estimates suggest that pesticide taxes based on the primary externality resulting from a particular mode of agricultural production are a credible policy option that warrants further consideration.

Population genetic diversity of the endemic Sardinian newt Euproctus platycephalus: implications for conservation

The Sardinian mountain newt Euproctus platycephalus, endemic to the island of Sardinia, (Italy), is considered a rare and threatened species and is classed as critically endangered by IUCN. It inhabits streams, small lakes and pools on the main mountain systems of the island. Threats from climatic and anthropogenic factors have raised concerns for the long-term survival of newt populations on the island. MtDNA sequencing was used to investigate the genetic population structure and phylogeography of this endemic species. Patterns of genetic variation were assessed by sequencing the complete Dloop region and part of the 12SrRNA, from 74 individuals representing four different populations. Analyses of molecular variance suggest that populations are significantly differentiated, and the distribution of haplotypes across the island shows strong geographical structuring. However, phylogenetic analyses also suggest that the Sardinian population consists of two distinct mtDNA groups, which may reflect ancient isolation and expansion events. Population structure, evolutionary history of the species and implications for the conservation of newt populations are discussed.

Genomic and genetic analyses of diversity and plant interactions of Pseudomonas fluorescens

Background: Pseudomonas fluorescens are common soil bacteria that can improve plant health through nutrient cycling, pathogen antagonism and induction of plant defenses. The genome sequences of strains SBW25 and Pf0-1 were determined and compared to each other and with P. fluorescens Pf-5. A functional genomic in vivo expression technology (IVET) screen provided insight into genes used by P. fluorescens in its natural environment and an improved understanding of the ecological significance of diversity within this species. Results: Comparisons of three P. fluorescens genomes (SBW25, Pf0-1, Pf-5) revealed considerable divergence: 61% of genes are shared, the majority located near the replication origin. Phylogenetic and average amino acid identity analyses showed a low overall relationship. A functional screen of SBW25 defined 125 plant-induced genes including a range of functions specific to the plant environment. Orthologues of 83 of these exist in Pf0-1 and Pf-5, with 73 shared by both strains. The P. fluorescens genomes carry numerous complex repetitive DNA sequences, some resembling Miniature Inverted-repeat Transposable Elements (MITEs). In SBW25, repeat density and distribution revealed 'repeat deserts' lacking repeats, covering approximately 40% of the genome. Conclusions: P. fluorescens genomes are highly diverse. Strain-specific regions around the replication terminus suggest genome compartmentalization. The genomic heterogeneity among the three strains is reminiscent of a species complex rather than a single species. That 42% of plant-inducible genes were not shared by all strains reinforces this conclusion and shows that ecological success requires specialized and core functions. The diversity also indicates the significant size of genetic information within the Pseudomonas pan genome.

Evaluation of cacao (Theobroma cacao L.) clones against seven Colombian isolates of Moniliophthora roreri (Cif.) Evans et al., representing four major genetic groupings of the pathogen in cacao

Artificial pod inoculation was used to compare the relative aggressiveness of seven Colombian isolates of Moniliophthora roreri (the causal agent of moniliasis or frosty pod disease), representing four major genetic groupings of the pathogen in cacao (cocoa), when applied to five diverse cacao genotypes (ICS-1, ICS-95, TSH-565, SCC-61 and CAP-34) at La Suiza Experimental Farm, Santander Department, Colombia. The following variables were evaluated 9 weeks after inoculation of 2- to 3-month-old pods with spore suspensions (1.2 x 10(5) spores mL(-1)): (i) disease incidence (DI); (ii) external severity (ES); and (iii) internal severity (IS). IS was found to be of greatest value in classifying the reaction of the host genotype against M. roreri. Genetic variation reported between isolates and cacao genotypes was not matched by similar diversity in their aggressiveness. All isolates were generally highly aggressive against most cacao genotypes, with only two isolates showing reduced IS and ES reactions. There was considerable variation between clones in the IS and ES scores, but one cultivated clone (ICS-95) displayed a significant level of resistance against all seven isolates. This clone may be useful in cacao breeding initiatives for resistance to moniliasis of cacao.

Influence of environmental stress on distributions of times to first division in Escherichia coli populations, as determined by digital-image analysis of individual cells

The distributions of times to first cell division were determined for populations of Escherichia coli stationary-phase cells inoculated onto agar media. This was accomplished by using automated analysis of digital images of individual cells growing on agar and calculation of the "box area ratio." Using approximately 300 cells per experiment, the mean time to first division and standard deviation for cells grown in liquid medium at 37 degrees C and inoculated on agar and incubated at 20 degrees C were determined as 3.0 h and 0.7 h, respectively. Distributions were observed to tail toward the higher values, but no definitive model distribution was identified. Both preinoculation stress by heating cultures at 50 degrees C and postinoculation stress by growth in the presence of higher concentrations of NaCl increased mean times to first division. Both stresses also resulted in an increase in the spread of the distributions that was proportional to the mean division time, the coefficient of variation being constant at approximately 0.2 in all cases. The "relative division time," which is the time to first division for individual cells expressed in terms of the cell size doubling time, was used as measure of the "work to be done" to prepare for cell division. Relative division times were greater for heat-stressed cells than for those growing under osmotic stress.

Genetic analysis of heat tolerance at anthesis in rice

Genetic analysis of heat tolerance will help breeders produce rice (Oryza sativa L.) varieties adapted to future climates. An F6 population of 181 recombinant inbred lines of Bala (tolerant) × Azucena (susceptible) was screened for heat tolerance at anthesis by measuring spikelet fertility at 30°C (control) and 38°C (high temperature) in experiments conducted in the Philippines and the United Kingdom. The parents varied significantly for absolute spikelet fertility under control (79–87%) and at high temperature (2.9–47.1%), and for relative spikelet fertility (high temperature/control) at high temperature (3.7–54.9%). There was no correlation between spikelet fertility in control and high-temperature conditions and no common quantitative trait loci (QTLs) were identified. Two QTLs for spikelet fertility under control conditions were identified on chromosomes 2 and 4. Eight QTLs for spikelet fertility under high-temperature conditions were identified on chromosomes 1, 2, 3, 8, 10, and 11. The most significant heat-responsive QTL, contributed by Bala and explaining up to 18% of the phenotypic variation, was identified on chromosome 1 (38.35 mega base pairs on the rice physical genome map). This QTL was also found to influence plant height, explaining 36.6% of the phenotypic variation. A comparison with other studies of abiotic (drought, cold, salinity) stresses showed QTLs at similar positions on chromosomes 1, 3, 8, and 10, suggesting common underlying stress-responsive regions of the genome.

Detection of somaclonal variation during cocoa somatic embryogenesis characterised using cleaved amplified polymorphic sequence and the new freeware Artbio

The scarcity and stochastic nature of genetic mutations presents a significant challenge for scientists seeking to characterise de novo mutation frequency at specific loci. Such mutations can be particularly numerous during regeneration of plants from in vitro culture and can undermine the value of germplasm conservation efforts. We used cleaved amplified polymorphic sequence (CAPS) analysis to characterise new mutations amongst a clonal population of cocoa plants regenerated via a somatic embryogenesis protocol used previously for cocoa cryopreservation. Efficacy of the CAPS system for mutation detection was greatly improved after an ‘a priori’ in silico screen of reference target sequences for actual and potential restriction enzyme recognition sites using a new freely available software called Artbio. Artbio surveys known sequences for existing restriction enzyme recognition sites but also identifies all single nucleotide polymorphism (SNP) deviations from such motifs. Using this software, we performed an in silico screen of seven loci for restriction sites and their potential mutant SNP variants that were possible from 21 restriction enzymes. The four most informative locus-enzyme combinations were then used to survey the regenerant populations for de novo mutants. We characterised the pattern of point mutations and, using the outputs of Artbio, calculated the ratio of base substitution in 114 somatic embryo-derived cocoa regenerants originating from two explant genotypes. We found 49 polymorphisms, comprising 26.3% of the samples screened, with an inferred rate of 2.8 × 10−3 substitutions/screened base. This elevated rate is of a similar order of magnitude to previous reports of de novo microsatellite length mutations arising in the crop and suggests caution should be exercised when applying somatic embryogenesis for the conservation of plant germplasm.

The compositional and configurational heterogeneity of matrix habitats shape woodland carabid communities in wooded-agricultural landscapes

Landscape heterogeneity (the composition and configuration of matrix habitats) plays a major role in shaping species communities in wooded-agricultural landscapes. However, few studies consider the influence of different types of semi-natural and linear habitats in the matrix, despite their known ecological value for biodiversity. Objective To investigate the importance of the composition and configuration of matrix habitats for woodland carabid communities and identify whether specific landscape features can help to maintain long-term populations in wooded-agricultural environments. Methods Carabids were sampled from woodlands in 36 tetrads of 4 km2 across southern Britain. Landscape heterogeneity including an innovative representation of linear habitats was quantified for each tetrad. Carabid community response was analysed using ordination methods combined with variation partitioning and additional response trait analyses. Results Woodland carabid community response was trait-specific and better explained by simultaneously considering the composition and configuration of matrix habitats. Semi-natural and linear features provided significant refuge habitat and functional connectivity. Mature hedgerows were essential for slow-dispersing carabids in fragmented landscapes. Species commonly associated with heathland were correlated with inland water and woodland patches despite widespread heathland conversion to agricultural land, suggesting that species may persist for some decades when elements representative of the original habitat are retained following landscape modification. Conclusions Semi-natural and linear habitats have high biodiversity value. Landowners should identify features that can provide additional resources or functional connectivity for species relative to other habitat types in the landscape matrix. Agri-environment options should consider landscape heterogeneity to identify the most efficacious changes for biodiversity.

A genome-wide test of the differential susceptibility hypothesis reveals a genetic predictor of differential response to psychological treatments for child anxiety disorders

Background: The differential susceptibly hypothesis suggests that certain genetic variants moderate the effects of both negative and positive environments on mental health and may therefore be important predictors of response to psychological treatments. Nevertheless, the identification of such variants has so far been limited to preselected candidate genes. In this study we extended the differential susceptibility hypothesis from a candidate gene to a genome-wide approach to test whether a polygenic score of environmental sensitivity predicted response to Cognitive Behavioural Therapy (CBT) in children with anxiety disorders. Methods: We identified variants associated with environmental sensitivity using a novel method in which within-pair variability in emotional problems in 1026 monozygotic (MZ) twin pairs was examined as a function of the pairs’ genotype. We created a polygenic score of environmental sensitivity based on the whole-genome findings and tested the score as a moderator of parenting on emotional problems in 1,406 children and response to individual, group and brief parent-led CBT in 973 children with anxiety disorders. Results: The polygenic score significantly moderated the effects of parenting on emotional problems and the effects of treatment. Individuals with a high score responded significantly better to individual CBT than group CBT or brief parent-led CBT (remission rates: 70.9%, 55.5% and 41.6% respectively). Conclusions: Pending successful replication, our results should be considered exploratory. Nevertheless, if replicated, they suggest that individuals with the greatest environmental sensitivity may be more likely to develop emotional problems in adverse environments, but also benefit more from the most intensive types of treatment.

The heterogeneity of wooded-agricultural landscape mosaics influences woodland bird community assemblages

Context Landscape heterogeneity (the composition and configuration of different landcover types) plays a key role in shaping woodland bird assemblages in wooded-agricultural mosaics. Understanding how species respond to landscape factors could contribute to preventing further decline of woodland bird populations. Objective To investigate how woodland birds with different species traits respond to landscape heterogeneity, and to identify whether specific landcover types are important for maintaining diverse populations in wooded-agricultural environments. Methods Birds were sampled from woodlands in 58 2 x 2 km tetrads across southern Britain. Landscape heterogeneity was quantified for each tetrad. Bird assemblage response was determined using redundancy analysis combined with variation partitioning and response trait analyses. Results For woodland bird assemblages, the independent explanatory importance of landscape composition and landscape configuration variables were closely interrelated. When considered simultaneously during variation partitioning, the community response was better represented by compositional variables. Different species responded to different landscape features and this could be explained by traits relating to woodland association, foraging strata and nest location. Ubiquitous, generalist species, many of which were hole-nesters or ground foragers, correlated positively with urban landcover while specialists of broadleaved woodland avoided landscapes containing urban areas. Species typical of coniferous woodland correlated with large conifer plantations. Conclusions At the 2 x 2 km scale, there was evidence that the availability of resources provided by proximate landcover types was highly important for shaping woodland bird assemblages. Further research to disentangle the effects of composition and configuration at different spatial scales is advocated.

Genetic and phenotypic characterization of complex hereditary spastic paraplegia

The hereditary spastic paraplegias are a heterogeneous group of degenerative disorders that are clinically classified as either pure with predominant lower limb spasticity, or complex where spastic paraplegia is complicated with additional neurological features, and are inherited in autosomal dominant, autosomal recessive or X-linked patterns. Genetic defects have been identified in over 40 different genes, with more than 70 loci in total. Complex recessive spastic paraplegias have in the past been frequently associated with mutations in SPG11 (spatacsin), ZFYVE26/SPG15, SPG7 (paraplegin) and a handful of other rare genes, but many cases remain genetically undefined. The overlap with other neurodegenerative disorders has been implied in a small number of reports, but not in larger disease series. This deficiency has been largely due to the lack of suitable high throughput techniques to investigate the genetic basis of disease, but the recent availability of next generation sequencing can facilitate the identification of disease- causing mutations even in extremely heterogeneous disorders. We investigated a series of 97 index cases with complex spastic paraplegia referred to a tertiary referral neurology centre in London for diagnosis or management. The mean age of onset was 16 years (range 3 to 39). The SPG11 gene was first analysed, revealing homozygous or compound heterozygous mutations in 30/97 (30.9%) of probands, the largest SPG11 series reported to date, and by far the most common cause of complex spastic paraplegia in the UK, with severe and progressive clinical features and other neurological manifestations, linked with magnetic resonance imaging defects. Given the high frequency of SPG11 mutations, we studied the autophagic response to starvation in eight affected SPG11 cases and control fibroblast cell lines, but in our restricted study we did not observe correlations between disease status and autophagic or lysosomal markers. In the remaining cases, next generation sequencing was carried out revealing variants in a number of other known complex spastic paraplegia genes, including five in SPG7 (5/97), four in FA2H (also known as SPG35) (4/97) and two in ZFYVE26/SPG15. Variants were identified in genes usually associated with pure spastic paraplegia and also in the Parkinson’s disease-associated gene ATP13A2, neuronal ceroid lipofuscinosis gene TPP1 and the hereditary motor and sensory neuropathy DNMT1 gene, highlighting the genetic heterogeneity of spastic paraplegia. No plausible genetic cause was identified in 51% of probands, likely indicating the existence of as yet unidentified genes.

Genetic dissection of vitamin E biosynthesis in tomato

Vegetables are critical for human health as they are a source of multiple vitamins including vitamin E (VTE). In plants, the synthesis of VTE compounds, tocopherol and tocotrienol, derives from precursors of the shikimate and methylerythritol phosphate pathways. Quantitative trait loci (QTL) for alpha-tocopherol content in ripe fruit have previously been determined in an Solanum pennellii tomato introgression line population. In this work, variations of tocopherol isoforms (alpha, beta, gamma, and delta) in ripe fruits of these lines were studied. In parallel all tomato genes structurally associated with VTE biosynthesis were identified and mapped. Previously identified VTE QTL on chromosomes 6 and 9 were confirmed whilst novel ones were identified on chromosomes 7 and 8. Integrated analysis at the metabolic, genetic and genomic levels allowed us to propose 16 candidate loci putatively affecting tocopherol content in tomato. A comparative analysis revealed polymorphisms at nucleotide and amino acid levels between Solanum lycopersicum and S. pennellii candidate alleles. Moreover, evolutionary analyses showed the presence of codons evolving under both neutral and positive selection, which may explain the phenotypic differences between species. These data represent an important step in understanding the genetic determinants of VTE natural variation in tomato fruit and as such in the ability to improve the content of this important nutriceutical.

Analysis of genotypic variation in genes associated with virulence in Aggregatibacter actinomycetemcomitans clinical isolates

Background and Objective: Although certain serotypes of Aggregatibacter actinomycetemcomitans are associated more with aggressive periodontitis than are other serotypes, the correlation between distinct lineages and virulence traits in this species is poorly understood. This study aimed to evaluate the polymorphism of genes encoding putative virulence factors of clinical isolates, and to correlate these findings with A. actinomycetemcomitans serotypes, genotypes and periodontal status of the hosts. Material and Methods: Twenty-six clinical isolates from diverse geographic populations with different periodontal conditions were evaluated. Genotyping was performed using pulse-field gel electrophoresis. Polymorphisms in the genes encoding leukotoxin, Aae, ApaH and determinants for serotype-specific O polysaccharide were investigated. Results: The isolates were classified into serotypes a-f, and exhibited three apaH genotypes, five aae alleles and 25 macrorestriction profiles. Two serotype b isolates (7.7%), obtained from Brazilian patients with aggressive periodontitis, were associated with the highly leukotoxic genotype; these isolates showed identical fingerprint patterns and aae and apaH genotypes. Serotype c, obtained from various periodontal conditions, was the most prevalent among Brazilian isolates, and isolates were distributed in two aae alleles, but formed a genetically distinct group based on apaH analysis. Cluster analysis showed a close relationship between fingerprinting genotypes and serotypes/apaH genotypes, but not with aae genotypes. Conclusion: Apart from the deletion in the ltx promoter region, no disease-associated markers were identified. Non-JP2-like strains recovered from individuals with periodontal disease exhibited considerable genetic variation regarding aae/apaH genotypes, serotypes and XhoI DNA fingerprints.

Genetic Control of Residual Variance for Teat Number in Pigs

The genetic improvement in litter size in pigs has been substantial during the last 10-15 years. The number of teats on the sow must increase as well to meet the needs of the piglets, because each piglet needs access to its own teat. We applied a genetic heterogeneity model on teat numberin sows, and estimated medium-high heritability for teat number (0.5), but low heritability for residual variance (0.05), indicating that selection for reduced variance might have very limited effect. A numerically positive correlation (0.8) between additive genetic breeding values for mean and for variance was found, but because of the low heritability for residual variance, the variance will increase very slowly with the mean.