Genetic analysis of complex demographic scenarios: Spatially expanding populations of the cane toad, Bufo marinus

Inferring the spatial expansion dynamics of invading species from molecular data is notoriously difficult due to the complexity of the processes involved. For these demographic scenarios, genetic data obtained from highly variable markers may be profitably combined with specific sampling schemes and information from other sources using a Bayesian approach. The geographic range of the introduced toad Bufo marinus is still expanding in eastern and northern Australia, in each case from isolates established around 1960. A large amount of demographic and historical information is available on both expansion areas. In each area, samples were collected along a transect representing populations of different ages and genotyped at 10 microsatellite loci. Five demographic models of expansion, differing in the dispersal pattern for migrants and founders and in the number of founders, were considered. Because the demographic history is complex, we used an approximate Bayesian method, based on a rejection-regression algorithm. to formally test the relative likelihoods of the five models of expansion and to infer demographic parameters. A stepwise migration-foundation model with founder events was statistically better supported than other four models in both expansion areas. Posterior distributions supported different dynamics of expansion in the studied areas. Populations in the eastern expansion area have a lower stable effective population size and have been founded by a smaller number of individuals than those in the northern expansion area. Once demographically stabilized, populations exchange a substantial number of effective migrants per generation in both expansion areas, and such exchanges are larger in northern than in eastern Australia. The effective number of migrants appears to be considerably lower than that of founders in both expansion areas. We found our inferences to be relatively robust to various assumptions on marker. demographic, and historical features. The method presented here is the only robust, model-based method available so far, which allows inferring complex population dynamics over a short time scale. It also provides the basis for investigating the interplay between population dynamics, drift, and selection in invasive species.

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.

Genetic analysis of potassium use efficiency in Brassica oleracea

Potassium (K) fertilizers are used in intensive and extensive agricultural systems to maximize production. However, there are both financial and environmental costs to K-fertilization. It is therefore important to optimize the efficiency with which K-fertilizers are used. Cultivating crops that acquire and/or utilize K more effectively can reduce the use of K-fertilizers. The aim of the present study was to determine the genetic factors affecting K utilization efficiency (KUtE), defined as the reciprocal of shoot K concentration (1/K(shoot)), and K acquisition efficiency (KUpE), defined as shoot K content, in Brassica oleracea. Genetic variation in K(shoot) was estimated using a structured diversity foundation set (DFS) of 376 accessions and in 74 commercial genotypes grown in glasshouse and field experiments that included phosphorus (P) supply as a treatment factor. Chromosomal quantitative trait loci (QTL) associated with K(shoot) and KUpE were identified using a genetic mapping population grown in the glasshouse and field. Putative QTL were tested using recurrent backcross substitution lines in the glasshouse. More than two-fold variation in K(shoot) was observed among DFS accessions grown in the glasshouse, a significant proportion of which could be attributed to genetic factors. Several QTL associated with K(shoot) were identified, which, despite a significant correlation in K(shoot) among genotypes grown in the glasshouse and field, differed between these two environments. A QTL associated with K(shoot) in glasshouse-grown plants (chromosome C7 at 62 center dot 2 cM) was confirmed using substitution lines. This QTL corresponds to a segment of arabidopsis chromosome 4 containing genes encoding the K(+) transporters AtKUP9, AtAKT2, AtKAT2 and AtTPK3. There is sufficient genetic variation in B. oleracea to breed for both KUtE and KUpE. However, as QTL associated with these traits differ between glasshouse and field environments, marker-assisted breeding programmes must consider carefully the conditions under which the crop will be grown.

Genetic analysis reveals a novel locus in Vicia faba decoupling pigmentation in the flower from that in the extra-floral nectaries

The genetics of the stipule spot pigmentation (SSP) in faba bean (Vicia faba L.) was studied using four inbred lines, of which Disco/2 was zero-tannin (zt2) with colourless stipule spots, ILB938/2 was normal-tannin (ZT2) with colourless stipule spots, and both Aurora/2 and Mélodie/2 were ZT2 with coloured stipule spots. Crosses Mélodie/2 × ILB 938/2, Mélodie/2 × Disco/2, ILB 938/2 × Aurora/2 and ILB 938/2 × Disco/2 (A, B, C and D, respectively) were prepared, along with reciprocals and backcrosses, and advanced through single-seed descent. All F1 hybrid plants had pigmented stipule spots, and in the F2 generation, the segregation ratio fit 3 coloured:1 colourless in crosses A, B and C and 9:7 in cross D. In the F3 generation, the ratio fit 5:3 in crosses A and C and 25:39 in cross D, and in the F4 generation, 9:7 in cross A. SSP was linked to the zero-tannin characteristics (white flower) only in cross B. The results show that coloured stipule spot is dominant to colourless and that colouration is determined by two unlinked complementary recessive genes. We propose the symbols ssp2 for the gene associated with zt2 in Disco/2 and ssp1 for the gene not associated with tannin content in ILB938/2. The novel ssp1 locus was mapped at F5 in cross ‘A’ using Medicago truncatula-derived single-nucleotide polymorphism and was on chromosome 1 of faba bean, in a well-conserved region of M. truncatula chromosome 5 containing some candidate Myb and basic helix–loop–helix transcription factor genes.

Recent developments in genetic data analysis: what can they tell us about human demographic history?

Over the last decade, a number of new methods of population genetic analysis based on likelihood have been introduced. This review describes and explains the general statistical techniques that have recently been used, and discusses the underlying population genetic models. Experimental papers that use these methods to infer human demographic and phylogeographic history are reviewed. It appears that the use of likelihood has hitherto had little impact in the field of human population genetics, which is still primarily driven by more traditional approaches. However, with the current uncertainty about the effects of natural selection, population structure and ascertainment of single-nucleotide polymorphism markers, it is suggested that likelihood-based methods may have a greater impact in the future.

High Resolution Melt (HRM) analysis is an efficient tool to genotype EMS mutants in complex crop genomes

Background: Targeted Induced Loci Lesions IN Genomes (TILLING) is increasingly being used to generate and identify mutations in target genes of crop genomes. TILLING populations of several thousand lines have been generated in a number of crop species including Brassica rapa. Genetic analysis of mutants identified by TILLING requires an efficient, high-throughput and cost effective genotyping method to track the mutations through numerous generations. High resolution melt (HRM) analysis has been used in a number of systems to identify single nucleotide polymorphisms (SNPs) and insertion/deletions (IN/DELs) enabling the genotyping of different types of samples. HRM is ideally suited to high-throughput genotyping of multiple TILLING mutants in complex crop genomes. To date it has been used to identify mutants and genotype single mutations. The aim of this study was to determine if HRM can facilitate downstream analysis of multiple mutant lines identified by TILLING in order to characterise allelic series of EMS induced mutations in target genes across a number of generations in complex crop genomes. Results: We demonstrate that HRM can be used to genotype allelic series of mutations in two genes, BraA.CAX1a and BraA.MET1.a in Brassica rapa. We analysed 12 mutations in BraA.CAX1.a and five in BraA.MET1.a over two generations including a back-cross to the wild-type. Using a commercially available HRM kit and the Lightscanner™ system we were able to detect mutations in heterozygous and homozygous states for both genes. Conclusions: Using HRM genotyping on TILLING derived mutants, it is possible to generate an allelic series of mutations within multiple target genes rapidly. Lines suitable for phenotypic analysis can be isolated approximately 8-9 months (3 generations) from receiving M3 seed of Brassica rapa from the RevGenUK TILLING service.

Analysis of DNA polymorphism in ancient barley herbarium material: validation of the KASP SNP genotyping platform

The availability of crop specimens archived in herbaria and old seed collections represent valuable resources for the analysis of plant genetic diversity and crop domestication. The ability to extract ancient DNA (aDNA) from such samples has recently allowed molecular genetic investigations to be undertaken in ancient materials. While analyses of aDNA initially focused on the use of markers which occur in multiple copies such as the internal transcribed spacer region (ITS) within ribosomal DNA and those requiring amplification of short DNA regions of variable length such as simple sequence repeats (SSRs), emphasis is now moving towards the genotyping of single nucleotide polymorphisms (SNPs), traditionally undertaken in aDNA by Sanger sequencing. Here, using a panel of barley aDNA samples previously surveyed by Sanger sequencing for putative causative SNPs within the flowering-time gene PPD-H1, we assess the utility of the Kompetitive Allele Specific PCR (KASP) genotyping platform for aDNA analysis. We find KASP to out-perform Sanger sequencing in the genotyping of aDNA samples (78% versus 61% success, respectively), as well as being robust to contamination. The small template size (≥46 bp) and one-step, closed-tube amplification/genotyping process make this platform ideally suited to the genotypic analysis of aDNA, a process which is often hampered by template DNA degradation and sample cross-contamination. Such attributes, as well as its flexibility of use and relatively low cost, make KASP particularly relevant to the genetic analysis of aDNA samples. Furthermore, KASP provides a common platform for the genotyping and analysis of corresponding SNPs in ancient, landrace and modern plant materials. The extended haplotype analysis of PPD-H1 undertaken here (allelic variation at which is thought to be important for the spread of domestication and local adaptation) provides further resolution to the previously identified geographic cline of flowering-time allele distribution, illustrating how KASP can be used to aid genetic analyses of aDNA from plant species. We further demonstrate the utility of KASP by genotyping ten additional genetic markers diagnostic for morphological traits in barley, shedding light on the phenotypic traits, alleles and allele combinations present in these unviable ancient specimens, as well as their geographic distributions.

Geographical variation in the clouded leopard, Neofelis nebulosa, reveals two species

The clouded leopard, Neofelis nebulosa, is an endangered semiarboreal felid with a wide distribution in tropical forests of southern and southeast Asia, including the islands of Sumatra and Borneo in the Indonesian archipelago [1]. In common with many larger animal species, it displays morphological variation within its wide geographical range and is currently regarded as comprising of up to four subspecies [2-4]. It is widely recognized that taxonomic designation has a major impact on conservation planning and action [5-8]. Given that the last taxonomic revision was made over 50 years ago [2], a more detailed examination of geographical variation is needed. We describe here the results of a morphometric analysis of the pelages of 57 clouded leopards sampled throughout the species' range. We conclude that there are two distinct morphological groups, which differ primarily in the size of their cloud markings. These results are supported by a recent genetic analysis [9]. On that basis, we give diagnoses for the distinction of two species, one in mainland Asia (N. nebulosa) and the other in Indonesia (N. diardi). The implications for conservation that arise from this new taxonomic arrangement are discussed.

The development of ISSR-derived SCAR markers around the Seasonal Flowering Locus (SFL) in fragaria vesca

Fragaria vesca is a short-lived perennial with a seasonal-flowering habit. Seasonality of flowering is widespread in the Rosaceae and is also found in the majority of temperate polycarpic perennials. Genetic analysis has shown that seasonal flowering is controlled by a single gene in F. vesca, the SEASONAL FLOWERING LOCUS (SFL). Here, we report progress towards the marker-assisted selection and positional cloning of SFL, in which three ISSR markers linked to SFL were converted to locus-specific sequence-characterized amplified region (SCAR1–SCAR3) markers to allow large-scale screening of mapping progenies. We believe this is the first study describing the development of SCAR markers from ISSR profiles. The work also provides useful insight into the nature of polymorphisms generated by the ISSR marker system. Our results indicate that the ISSR polymorphisms originally detected were probably caused by point mutations in the positions targeted by primer anchors (causing differential PCR failure), by indels within the amplicon (leading to variation in amplicon size) and by internal sequence differences (leading to variation in DNA folding and so in band mobility). The cause of the original ISSR polymorphism was important in the selection of appropriate strategies for SCAR-marker development. The SCAR markers produced were mapped using a F. vesca f. vesca × F. vesca f. semperflorens testcross population. Marker SCAR2 was inseparable from the SFL, whereas SCAR1 mapped 3.0 cM to the north of the gene and SCAR3 1.7 cM to its south.

Mining gene networks with application to GAW15 problem 1

The genetic analysis workshop 15 (GAW15) problem 1 contained baseline expression levels of 8793 genes in immortalised B cells from 194 individuals in 14 Centre d’Etude du Polymorphisme Humane (CEPH) Utah pedigrees. Previous analysis of the data showed linkage and association and evidence of substantial individual variations. In particular, correlation was examined on expression levels of 31 genes and 25 target genes corresponding to two master regulatory regions. In this analysis, we apply Bayesian network analysis to gain further insight into these findings. We identify strong dependences and therefore provide additional insight into the underlying relationships between the genes involved. More generally, the approach is expected to be applicable for integrated analysis of genes on biological pathways.

Adaptive divergence and speciation among sexual and pseudoviviparous populations of Festuca

Pseudovivipary is an environmentally induced flowering abnormality in which vegetative shoots replace seminiferous (sexual) inflorescences. Pseudovivipary is usually retained in transplantation experiments, indicating that the trait is not solely induced by the growing environment. Pseudovivipary is the defining characteristic of Festuca vivipara, and arguably the only feature separating this species from its closest seminiferous relative, Festuca ovina. We performed phylogenetic and population genetic analysis on sympatric F. ovina and F. vivipara samples to establish whether pseudovivipary is an adaptive trait that accurately defines the separation of genetically distinct Festuca species. Chloroplast and nuclear marker-based analyses revealed that variation at a geographical level can exceed that between F. vivipara and F. ovina. We deduced that F. vivipara is a recent species that frequently arises independently within F. ovina populations and has not accumulated significant genetic differentiation from its progenitor. We inferred local gene flow between the species. We identified one amplified fragment length polymorphism marker that may be linked to a pseudovivipary-related region of the genome, and several other markers provide evidence of regional local adaptation in Festuca populations. We conclude that F. vivipara can only be appropriately recognized as a morphologically and ecologically distinct species; it lacks genetic differentiation from its relatives. This is the first report of a ‘failure in normal flowering development’ that repeatedly appears to be adaptive, such that the trait responsible for species recognition constantly reappears on a local basis.

Integration of retrotransposons-based markers in a linkage map of barley

A deeper understanding of random markers is important if they are to be employed for a range of objectives. The sequence specific amplified polymorphism (S-SAP) technique is a powerful genetic analysis tool which exploits the high copy number of retrotransposon long terminal repeats (LTRs) in the plant genome. The distribution and inheritance of S-SAP bands in the barley genome was studied using the Steptoe × Morex (S × M) double haploid (DH) population. Six S-SAP primer combinations generated 98 polymorphic bands, and map positions were assigned to all but one band. Eight putative co-dominant loci were detected, representing 16 of the mapped markers. Thus at least 81 of the mapped S-SAP loci were dominant. The markers were distributed along all of the seven chromosomes and a tendency to cluster was observed. The distribution of S-SAP markers over the barley genome concurred with the knowledge of the high copy number of retrotransposons in plants. This experiment has demonstrated the potential for the S-SAP technique to be applied in a range of analyses such as genetic fingerprinting, marker assisted breeding, biodiversity assessment and phylogenetic analyses.

The biosurfactant viscosin produced by Pseudomonas fluorescens SBW25 aids spreading motility and plant growth promotion

Food security depends on enhancing production and reducing loss to pests and pathogens. A promising alternative to agrochemicals is the use of plant growth-promoting rhizobacteria (PGPR), which are commonly associated with many, if not all, plant species. However, exploiting the benefits of PGPRs requires knowledge of bacterial function and an in-depth understanding of plant-bacteria associations. Motility is important for colonization efficiency and microbial fitness in the plant environment, but the mechanisms employed by bacteria on and around plants are not well understood. We describe and investigate an atypical mode of motility in Pseudomonas fluorescens SBW25 that was revealed only after flagellum production was eliminated by deletion of the master regulator fleQ. Our results suggest that this ‘spidery spreading’ is a type of surface motility. Transposon mutagenesis of SBW25ΔfleQ (SBW25Q) produced mutants, defective in viscosin production, and surface spreading was also abolished. Genetic analysis indicated growth-dependency, production of viscosin, and several potential regulatory and secretory systems involved in the spidery spreading phenotype. Moreover, viscosin both increases efficiency of surface spreading over the plant root and protects germinating seedlings in soil infected with the plant pathogen Pythium. Thus, viscosin could be a useful target for biotechnological development of plant growth promotion agents.

Alterations in predicted regulatory and coding regions of the sterol 14α-demethylase gene (CYP51) confer decreased azole sensitivity in the oilseed rape pathogen Pyrenopeziza brassicae

The incidence and severity of light leaf spot epidemics caused by the ascomycete fungus Pyrenopeziza brassicae on UK oilseed rape crops is increasing. The disease is currently controlled by a combination of host resistance, cultural practices and fungicide applications. We report decreases in sensitivities of modern UK P. brassicae isolates to the azole (imidazole and triazole) class of fungicides. By cloning and sequencing the P. brassicae CYP51 (PbCYP51) gene, encoding the azole target sterol 14α-demethylase, we identified two non-synonymous mutations encoding substitutions G460S and S508T associated with reduced azole sensitivity. We confirmed the impact of the encoded PbCYP51 changes on azole sensitivity and protein activity by heterologous expression in a Saccharomyces cerevisiae mutant YUG37::erg11 carrying a controllable promoter of native CYP51 expression. In addition, we identified insertions in the predicted regulatory regions of PbCYP51 in isolates with reduced azole sensitivity. The presence of these insertions was associated with enhanced transcription of PbCYP51 in response to sub-inhibitory concentrations of the azole fungicide tebuconazole. Genetic analysis of in vitro crosses of sensitive and resistant isolates confirmed the impact of PbCYP51 alterations in coding and regulatory sequences on a reduced sensitivity phenotype, as well as identifying a second major gene at another locus contributing to resistance in some isolates. The least sensitive field isolates carry combinations of upstream insertions and non-synonymous mutations, suggesting PbCYP51 evolution is on-going and the progressive decline in azole sensitivity of UK P. brassicae populations will continue. The implications for the future control of light leaf spot are discussed.

Fighting in fig wasps: do males avoid killing brothers or do they never meet them?

1. In many fig wasp species, armoured wingless males regularly engage in lethal fights for access to females inside figs, which act as discrete mating patches. 2. Kin selection generally opposes killing brothers, because their reproductive success provides indirect genetic benefits (inclusive fitness). However, siblicide may be avoided if (i) brothers do not occur in the same figs, or (ii) males avoid fighting brothers in the same fig. Alternatively, (iii) siblicide may occur because intense mate competition between brothers at the local scale overcomes kin selection effects, or (iv) males do not recognise kin. 3. A fig may also contain wasps from other closely related species and it is not known if males also fight with these individuals. 4. Nine microsatellite loci were used in the first genetic analysis of fighting in fig wasps. We assigned species and sibling identities to males and tested alternative fighting scenarios for three Sycoscapter wasp species in figs of Ficus rubiginosa. 5. Approximately 60% of figs contained males frommore than one Sycoscapter species and approximately 80% of fights were between conspecifics, but a surprising 20% were between heterospecific males. 6.Within species, fewfigs contained brothers, suggesting that females typically lay one son per fig. Overall, most males do not compete with brothers and all fights observed were between unrelated males. Key words:Competition, fighting, genetics, kin selection, microsatellites, relatedness.