Rapid genome wide mapping of phosphine resistance loci by a simple regional averaging analysis in the red flour beetle, Tribolium castaneum

Background Next-generation sequencing technology is an important tool for the rapid, genome-wide identification of genetic variations. However, it is difficult to resolve the ‘signal’ of variations of interest and the ‘noise’ of stochastic sequencing and bioinformatic errors in the large datasets that are generated. We report a simple approach to identify regional linkage to a trait that requires only two pools of DNA to be sequenced from progeny of a defined genetic cross (i.e. bulk segregant analysis) at low coverage (<10×) and without parentage assignment of individual SNPs. The analysis relies on regional averaging of pooled SNP frequencies to rapidly scan polymorphisms across the genome for differential regional homozygosity, which is then displayed graphically. Results Progeny from defined genetic crosses of Tribolium castaneum (F4 and F19) segregating for the phosphine resistance trait were exposed to phosphine to select for the resistance trait while the remainders were left unexposed. Next generation sequencing was then carried out on the genomic DNA from each pool of selected and unselected insects from each generation. The reads were mapped against the annotated T. castaneum genome from NCBI (v3.0) and analysed for SNP variations. Since it is difficult to accurately call individual SNP frequencies when the depth of sequence coverage is low, variant frequencies were averaged across larger regions. Results from regional SNP frequency averaging identified two loci, tc_rph1 on chromosome 8 and tc_rph2 on chromosome 9, which together are responsible for high level resistance. Identification of the two loci was possible with only 5-7× average coverage of the genome per dataset. These loci were subsequently confirmed by direct SNP marker analysis and fine-scale mapping. Individually, homozygosity of tc_rph1 or tc_rph2 results in only weak resistance to phosphine (estimated at up to 1.5-2.5× and 3-5× respectively), whereas in combination they interact synergistically to provide a high-level resistance >200×. The tc_rph2 resistance allele resulted in a significant fitness cost relative to the wild type allele in unselected beetles over eighteen generations. Conclusion We have validated the technique of linkage mapping by low-coverage sequencing of progeny from a simple genetic cross. The approach relied on regional averaging of SNP frequencies and was used to successfully identify candidate gene loci for phosphine resistance in T. castaneum. This is a relatively simple and rapid approach to identifying genomic regions associated with traits in defined genetic crosses that does not require any specialised statistical analysis.

Resistance to Fusarium oxysporum f. sp fragariae and Predicted Breeding Values in Strawberry

Fusarium wilt of strawberry, incited by Fusarium oxysporum f. sp. fragariae (Fof), is a major disease of the cultivated strawberry (Fragaria xananassa) worldwide. An increase in disease outbreaks of the pathogen in Western Australia and Queensland plus the search for alternative disease management strategies place emphasis on the development of resistant cultivars. In response, a partial incomplete diallel cross involving four parents was performed for use in glasshouse resistance screenings. The resulting progeny were evaluated for their susceptibility to Fof. Best-performing progeny and suitability of progenies as parents were determined using data from disease severity ratings and analyzed using a linear mixed model incorporating a pedigree to produce best linear unbiased predictions of breeding values. Variation in disease response, ranging from highly susceptible to resistant, indicates a quantitative effect. The estimate of the narrow-sense heritability was 0.49 +/- 0.04 (SE), suggesting the population should be responsive to phenotypic recurrent selection. Several progeny genotypes have predicted breeding values higher than any of the parents. Knowledge of Fof resistance derived from this study can help select best parents for future crosses for the development of new strawberry cultivars with Fof resistance.

Induced Variations in Brassinosteroid Genes Define Barley Height and Sturdiness, and Expand the Green Revolution Genetic Toolkit

Reduced plant height and culm robustness are quantitative characteristics important for assuring cereal crop yield and quality under adverse weather conditions. A very limited number of short-culm mutant alleles were introduced into commercial crop cultivars during the Green Revolution. We identified phenotypic traits, including sturdy culm, specific for deficiencies in brassinosteroid biosynthesis and signaling in semidwarf mutants of barley (Hordeum vulgare). This set of characteristic traits was explored to perform a phenotypic screen of near-isogenic short-culm mutant lines from the brachytic, breviaristatum, dense spike, erectoides, semibrachytic, semidwarf, and slender dwarf mutant groups. In silico mapping of brassinosteroid-related genes in the barley genome in combination with sequencing of barley mutant lines assigned more than 20 historic mutants to three brassinosteroid-biosynthesis genes (BRASSINOSTEROID-6-OXIDASE, CONSTITUTIVE PHOTOMORPHOGENIC DWARF, and DIMINUTO) and one brassinosteroid-signaling gene (BRASSINOSTEROID-INSENSITIVE1 [HvBRI1]). Analyses of F2 and M2 populations, allelic crosses, and modeling of nonsynonymous amino acid exchanges in protein crystal structures gave a further understanding of the control of barley plant architecture and sturdiness by brassinosteroid-related genes. Alternatives to the widely used but highly temperature-sensitive uzu1.a allele of HvBRI1 represent potential genetic building blocks for breeding strategies with sturdy and climate-tolerant barley cultivars.

Beasts on Fields. Human-Wildlife Conflicts in Nature-Culture Borderlands

Human-wildlife conflicts are today an integral part of the rural development discourse. In this research, the main focus is on the spatial explanation which is not a very common approach in the reviewed literature. My research hypothesis is based on the assumption that human-wildlife conflicts occur when a wild animal crosses a perceived borderline between the nature and culture and enters into the realms of the other. The borderline between nature and culture marks a perceived division of spatial content in our senses of place. The animal subject that crosses this border becomes a subject out of place meaning that the animal is then spatially located in a space where it should not be or where it does not belong according to tradition, custom, rules, law, public opinion, prevailing discourse or some other criteria set by human beings. An appearance of a wild animal in a domesticated space brings an uncontrolled subject into that space where humans have previously commanded total control of all other natural elements. A wild animal out of place may also threaten the biosecurity of the place in question. I carried out a case study in the Liwale district in south-eastern Tanzania to test my hypothesis during June and July 2002. I also collected documents and carried out interviews in Dar es Salaam in 2003. I studied the human-wildlife conflicts in six rural villages, where a total of 183 persons participated in the village meetings. My research methods included semi-structured interviews, participatory mapping, questionnaire survey and Q- methodology. The rural communities in the Liwale district have a long-history of co-existing with wildlife and they still have traditional knowledge of wildlife management and hunting. Wildlife conservation through the establishment of game reserves during the colonial era has escalated human-wildlife conflicts in the Liwale district. This study shows that the villagers perceive some wild animals differently in their images of the African countryside than the district and regional level civil servants do. From the small scale subsistence farmers point of views, wild animals continue to challenge the separation of the wild (the forests) and the domestics spaces (the cultivated fields) by moving across the perceived borders in search of food and shelter. As a result, the farmers may loose their crops, livestock or even their own lives in the confrontations of wild animals. Human-wildlife conflicts in the Liwale district are manifold and cannot be explained simply on the basis of attitudes or perceived images of landscapes. However, the spatial explanation of these conflicts provides us some more understanding of why human-wildlife conflicts are so widely found across the world.

Genetic characterization of field-evolved resistance to phosphine in the rusty grain beetle, Cryptolestes ferrugineus (Laemophloeidae: Coleoptera)

Inheritance of resistance to phosphine fumigant was investigated in three field-collected strains of rusty grain beetle, Cryptolestes ferrugineus, Susceptible (S-strain), Weakly Resistant (Weak-R) and Strongly Resistant (Strong-R). The strains were purified for susceptibility, weak resistance and strong resistance to phosphine, respectively, to ensure homozygosity of resistance genotype. Crosses were established between S-strain × Weak-R, S-strain × Strong-R and Weak-R × Strong-R, and the dose mortality responses to phosphine of these strains and their F1, F2 and F1-backcross progeny were obtained. The fumigations were undertaken at 25 °C and 55% RH for 72 h. Weak-R and Strong-R showed resistance factors of 6.3 × and 505 × compared with S-strain at the LC50. Both weak and strong resistances were expressed as incompletely recessive with degrees of dominance of − 0.48 and − 0.43 at the LC50, respectively. Responses of F2 and F1-backcross progeny indicated the existence of one major gene in Weak-R, and at least two major genes in Strong-R, one of which was allelic with the major factor in Weak-R. Phenotypic variance analyses also estimated that the number of independently segregating genes conferring weak resistance was 1 (nE = 0.89) whereas there were two genes controlling strong resistance (nE = 1.2). The second gene, unique to Strong-R, interacted synergistically with the first gene to confer a very high level of resistance (~ 80 ×). Neither of the two major resistance genes was sex linked. Despite the similarity of the genetics of resistance to that previously observed in other pest species, a significant proportion (~ 15 to 30%) of F1 individuals survived at phosphine concentrations higher than predicted. Thus it is likely that additional dominant heritable factors, present in some individuals in the population, also influenced the resistance phenotype. Our results will help in understanding the process of selection for phosphine resistance in the field which will inform resistance management strategies. In addition, this information will provide a basis for the identification of the resistance genes.

Genetic characterization of field-evolved resistance to phosphine in the rusty grain beetle, Cryptolestes ferrugineus (Laemophloeidae: Coleoptera)

Inheritance of resistance to phosphine fumigant was investigated in three field-collected strains of rusty grain beetle, Cryptolestes ferrugineus, Susceptible (S-strain), Weakly Resistant (Weak-R) and Strongly Resistant (Strong-R). The strains were purified for susceptibility, weak resistance and strong resistance to phosphine, respectively, to ensure homozygosity of resistance genotype. Crosses were established between S-strain × Weak-R, S-strain × Strong-R and Weak-R × Strong-R, and the dose mortality responses to phosphine of these strains and their F1, F2 and F1-backcross progeny were obtained. The fumigations were undertaken at 25 °C and 55% RH for 72 h. Weak-R and Strong-R showed resistance factors of 6.3 × and 505 × compared with S-strain at the LC50. Both weak and strong resistances were expressed as incompletely recessive with degrees of dominance of − 0.48 and − 0.43 at the LC50, respectively. Responses of F2 and F1-backcross progeny indicated the existence of one major gene in Weak-R, and at least two major genes in Strong-R, one of which was allelic with the major factor in Weak-R. Phenotypic variance analyses also estimated that the number of independently segregating genes conferring weak resistance was 1 (nE = 0.89) whereas there were two genes controlling strong resistance (nE = 1.2). The second gene, unique to Strong-R, interacted synergistically with the first gene to confer a very high level of resistance (~ 80 ×). Neither of the two major resistance genes was sex linked. Despite the similarity of the genetics of resistance to that previously observed in other pest species, a significant proportion (~ 15 to 30%) of F1 individuals survived at phosphine concentrations higher than predicted. Thus it is likely that additional dominant heritable factors, present in some individuals in the population, also influenced the resistance phenotype. Our results will help in understanding the process of selection for phosphine resistance in the field which will inform resistance management strategies. In addition, this information will provide a basis for the identification of the resistance genes.

Genetic conservation of Phosphine Resistance in the Rice Weevil Sitophilus oryzae (L.)

High levels of resistance to phosphine in the rice weevil Sitophilus oryzae have been detected in Asian countries including China and Vietnam, however there is limited knowledge of the genetic mechanism of resistance in these strains. We find that the genetic basis of strong phosphine resistance is conserved between strains of S. oryzae from China, Vietnam and Australia. Each of four strongly resistant strains has an identical amino acid variant in the encoded dihydrolipoamide dehydrogenase (DLD) enzyme that was previously identified as a resistance factor in Rhyzopertha dominica and Tribolium castaneum. The unique amino acid substitution, Asparagine > Threonine (N505T) of all strongly resistant S. oryzae corresponds to the position of an Asparagine > Histidine variant (N506H) that was previously reported in strongly resistant R. dominica. Progeny (F16 and F18) from two independent crosses showed absolute linkage of N505T to the strong resistance phenotype, indicating that if N505T was not itself the resistance variant that it resided within 1 or 2 genes of the resistance factor. Non-complementation between the strains confirmed the shared genetic basis of strong resistance, which was supported by the very similar level of resistance between the strains, with LC50 values ranging from 0.20 to 0.36 mgL-1 for a 48 hour exposure at 25°C. Thus, the mechanism of high level resistance to phosphine is strongly conserved between R. dominica, T. castaneum and S. oryzae. A fitness cost associated with strongly resistant allele was observed in segregating populations in the absence of selection.

Behavioural, population, and genetic processes affecting metapopulation dynamics of the Glanville fritillary butterfly

In my thesis I have been studying the effects of population fragmentation and extinction-recolonization dynamics on genetic and evolutionary processes in the Glanville fritillary butterfly (Melitaea cinxia). By conducting crosses within and among newly-colonized populations and using several fitness measures, I found a strong decrease in fitness following colonization by a few related individuals, and a strong negative relationship between parental relatedness and offspring fitness. Thereafter, I was interested in determining the number and relatedness of individuals colonizing new populations, which I did using a set of microsatellites I had previously developed for this species. Additionally, I am interested in the evolution of key life-history traits. By following the lifetime reproductive success of males emerging at different times in a semi-natural setup, I demonstrated that protandry is adaptive in males, and I was able to rule out, for M. cinxia, alternative incidental hypotheses evoked to explain the evolution of protandry in insects. Finally, in work I did together with Prof. Hanna Kokko, I am proposing bet-hedging as a new mechanism that could explain the evolution of polyandry in M. cinxia.

Phase diagram of vortex matter in layered superconductors with tilted columnar pinning centers

We study the vortex matter phase diagram of a layered superconductor in the presence of columnar pinning defects, tilted with respect to the normal to the layers. We use numerical minimization of the free energy written as a functional of the time-averaged vortex density of the Ramakrishnan-Yussouff form, supplemented by the appropriate pinning potential. We study the case where the pin density is smaller than the areal vortex density. At lower pin concentrations, we find, for temperatures of the order of the melting temperature of the unpinned lattice, a Bose glass type phase which at lower temperatures converts, via a first-order transition, to a Bragg glass, while, at higher temperatures, it crosses over to an interstitial liquid. At somewhat higher concentrations, no transition to a Bragg glass is found even at the lowest temperatures studied. While qualitatively the behavior we find is similar to that obtained using the same procedures for columnar pins normal to the layers, there are important and observable quantitative differences, which we discuss.

Steady-state dendritic cells expressing cognate antigen terminate memory CD8+ T-cell responses

Antigen stimulation of naive T cells in conjunction with strong costimulatory signals elicits the generation of effector and memory populations. Such terminal differentiation transforms naive T cells capable of differentiating along several terminal pathways in response to pertinent environmental cues into cells that have lost developmental plasticity and exhibit heightened responsiveness. Because these cells exhibit little or no need for the strong costimulatory signals required for full activation of naive T cells, it is generally considered memory and effector T cells are released from the capacity to be inactivated. Here, we show that steadystate dendritic cells constitutively presenting an endogenously expressed antigen inactivate fully differentiated memory and effector CD8+ T cells in vivo through deletion and inactivation. These findings indicate that fully differentiated effector and memory T cells exhibit a previously unappreciated level of plasticity and provide insight into how memory and effector T-cell populations may be regulated. © 2008 by The American Society of Hematology.

Reproduction and population structure in European aspen

The European aspen (Populus tremula) is a keystone species for biodiversity in boreal forests. However, the future of aspen may be threatened, because large aspens have mostly been removed from managed forests, whereas regeneration and the long-term persistence of mature trees are subjects of concern in protected areas. Aspen is a pioneer tree, and it can reproduce both sexually by seed and asexually by root suckers. Through asexual reproduction aspen forms clones, groups of genetically identical trees (ramets). In my thesis, I have studied the structure of aspen populations in terms of number, size, clonal and demographic properties. Additionally, I have investigated the emergence and survival of seedlings as well as the seed quantity and quality in crosses between the European and hybrid aspen. To study the regeneration and population structure, mature aspens were recorded in old-growth and managed forests in eastern Finland based on a large-scale inventory (11 400 ha). In addition, small aspen trees were surveyed on sample plots. Clonal structure was investigated both by morphological characters and by DNA-based markers (microsatellites). Seedling emergence and survival was studied with two sowing experiments. With crosses between European and hybrid aspens we wanted to study whether elevated temperatures due to climate change would benefit the different crosses of European and hybrid aspen unequally and thus affect the gene flow between the two species. The average volumes of mature aspen were 5.3 m3/ha in continuous old-growth, and 0.8 m3/ha in managed forests. Results indicate also that large aspen trees in managed forests are a legacy of the past less intensively managed forest landscapes. Long-term persistence of aspen in protected areas can only be secured by restoration measures creating sufficiently large gaps for regeneration. More emphasis should be given to sparing aspens in thinnings and to retaining of mature aspens in regeneration cutting in managed forests. Aspen was found to be spatially aggregated in the landscape. This could be explained by site type, disturbance history and / or limitations in seed dispersal. Clonal structure does not explain the spatial aggregation, since average size of the clones was only 2.3 ramets, and most clones (70 %) consisted of just one ramet. The small size of the clones suggests that most of them are relatively young. Therefore, sexual reproduction may be more common than has previously been thought. Seedling emergence was most successful in mineral soil especially, when the site had been burned. Only few seedlings occurred on humus. Survival of the seedlings was low, and strongly dependent on moisture, but also on seedbed conditions. The seeds were found to maintain their germinability longer than has earlier been thought to be possible. Interspecific crosses produced more seeds with higher quality than intraspecific crosses. When temperature was elevated, germination of hybrid aspen seeds increased more than seeds from P. tremula x P. tremula crosses. These results suggest that hybrid aspen may have a significant genetic impact on the European aspen, and this effect may become strengthened by climate warming.

First passage probability during random vibration

The probability that a random process crosses an arbitrary level for the first time is expressed as a Gram—Charlier series, the leading term of which is the Poisson approximation. The coefficients of this series are related to the moments of the number of level crossings. The results are applicable to both stationary and non-stationary processes. Some numerical results are presented for the response process of a linear single-degree-of-freedom oscillator under Gaussian white noise excitation.

Surface freshwater from Bay of Bengal runoff and Indonesian Throughflow in the Tropical Indian Ocean

According to recent estimates, the annual total continental runoff into the Bay of Bengal (BoB) is about 2950 km 3, which is more than half that into the entire tropical Indian Ocean (IO). Here we use climatological observations to trace the seasonal pathways of near surface freshwater from BoB runoff and Indonesian Throughflow (ITF) by removing the net contribution from precipitation minus evaporation. North of 20 degrees S, the amount of freshwater from BoB runoff and ITF changes with season in a manner consistent with surface currents from drifters. BoB runoff reaches remote regions of the Arabian Sea; it also crosses the equator in the east to join the ITF. This freshwater subsequently flows west across the southern tropical IO in the South Equatorial Current.

Estimation of cutting volume with three inventory methods for harvest planning in Canadian boreal forests.

Two methods of pre-harvest inventory were designed and tested on three cutting sites containing a total of 197 500 m3 of wood. These sites were located on flat-ground boreal forests located in northwestern Quebec. Both methods studied involved scaling of trees harvested to clear the road path one year (or more) prior to harvest of adjacent cut-blocks. The first method (ROAD) considers the total road right-of-way volume divided by the total road area cleared. The resulting volume per hectare is then multiplied by the total cut-block area scheduled for harvest during the following year to obtain the total estimated cutting volume. The second method (STRATIFIED) also involves scaling of trees cleared from the road. However, in STRATIFIED, log scaling data are stratified by forest stand location. A volume per hectare is calculated for each stretch of road that crosses a single forest stand. This volume per hectare is then multiplied by the remaining area of the same forest stand scheduled for harvest one year later. The sum of all resulting estimated volumes per stand gives the total estimated cutting-volume for all cut-blocks adjacent to the studied road. A third method (MNR) was also used to estimate cut-volumes of the sites studied. This method represents the actual existing technique for estimating cutting volume in the province of Quebec. It involves summing the cut volume for all forest stands. The cut volume is estimated by multiplying the area of each stand by its estimated volume per hectare obtained from standard stock tables provided by the governement. The resulting total estimated volume per cut-block for all three methods was then compared with the actual measured cut-block volume (MEASURED). This analysis revealed a significant difference between MEASURED and MNR methods with the MNR volume estimate being 30 % higher than MEASURED. However, no significant difference from MEASURED was observed for volume estimates for the ROAD and STRATIFIED methods which respectively had estimated cutting volumes 19 % and 5 % lower than MEASURED. Thus the ROAD and STRATIFIED methods are good ways to estimate cut-block volumes after road right-of-way harvest for conditions similar to those examined in this study.

Optical conductivity of perovskite manganites

We calculate the optical conductivity sigma(omega) for doped rare-earth manganites based on the recently proposed microscopic ``two fluid'' l-b model. We study the temperature dependence of sigma(omega) for La(0.825)Sr(0.175)MnO(3), which has a metallic ground state. At low temperatures, the calculated sigma(omega) shows a ``two-peak'' structure consisting of a far-infrared coherent Drude peak and a broad mid-infrared ``polaron'' peak, as observed in experiments. Upon heating, the Drude peak rapidly loses spectral weight, and sigma(omega) crosses over to having just a single broad mid-infrared peak. The temperature dependence of the mid-infrared peak and the spectral weight transfer between the two peaks are also in agreement with experimental findings. We also study the doping dependence of sigma(omega) for the same compound. The integrated spectral weight under the Drude peak increases rapidly as the doping level is increased from an underdoped, insulating state (x = 0.1) to a highly doped, metallic state (x = 0.3), again in agreement with trends seen experimentally.