Gene Genealogies Strongly Distorted by Weakly Interfering Mutations in Constant Environments

Neutral nucleotide diversity does not scale with population size as expected, and this “paradox of variation” is especially severe for animal mitochondria. Adaptive selective sweeps are often proposed as a major cause, but a plausible alternative is selection against large numbers of weakly deleterious mutations subject to Hill–Robertson interference. The mitochondrial genealogies of several species of whale lice (Amphipoda: Cyamus) are consistently too short relative to neutral-theory expectations, and they are also distorted in shape (branch-length proportions) and topology (relative sister-clade sizes). This pattern is not easily explained by adaptive sweeps or demographic history, but it can be reproduced in models of interference among forward and back mutations at large numbers of sites on a nonrecombining chromosome. A coalescent simulation algorithm was used to study this model over a wide range of parameter values. The genealogical distortions are all maximized when the selection coefficients are of critical intermediate sizes, such that Muller''s ratchet begins to turn. In this regime, linked neutral nucleotide diversity becomes nearly insensitive to N. Mutations of this size dominate the dynamics even if there are also large numbers of more strongly and more weakly selected sites in the genome. A genealogical perspective on Hill–Robertson interference leads directly to a generalized background-selection model in which the effective population size is progressively reduced going back in time from the present.OBSERVED levels of apparently neutral nucleotide diversity (π_n) are typically lower than expected under the assumptions of standard equilibrium theories, and they vary much less among species than do estimates of long-term effective population sizes (Nei and Grauer 1984; Bazin et al. 2006; Nabholz et al. 2008). Many explanations have been proposed for the apparent shortfalls and the lack of proportionality with population size, including (1) complex demographic histories (e.g., recurring population bottlenecks), (2) adaptive selective sweeps (Maynard Smith and Haigh 1974; Gillespie 1999), and (3) selection against deleterious mutations (Charlesworth et al. 1993, 1995; McVean and Charlesworth 2000; Comeron et al. 2008). Of these three possibilities, bottlenecks and sweeps are by far the most frequently mentioned, even though deleterious mutations occur at high rates in all species, regardless of ecological circumstances (Eyre-Walker and Keightley 2007). Here we show that weakly deleterious mutations can distort genealogies in three different ways and dramatically reduce nucleotide diversities in large populations of nonrecombining chromosomes. The mitochondrial genealogies of several species of whale lice (Kaliszewska et al. 2005) are distorted in exactly these ways, and several lines of evidence suggest that bottlenecks and adaptive sweeps are not likely to be the primary causes.Mitochondria have been proposed to be especially sensitive to selective sweeps. Animal mitochondrial genomes contain more than three dozen essential protein and structural RNA genes, so they are large targets for both mutation and selection (Ballard and Whitlock 2004). They do not undergo sexual recombination, so every advantageous mutation that fixes will reduce variation throughout the genome. Mitochondrial nucleotide diversity therefore could depend strongly on rates of environmental change, which could be similar for species with very different population sizes. Indeed, if rates of mitochondrial adaptation were mutation limited, then larger populations might actually experience higher rates of adaptive substitution and as a result show lower average levels of neutral diversity than smaller populations (Gillespie 2000, 2001). This idea was recently invoked to explain the remarkable similarity of average levels of mitochondrial nucleotide diversity among the major animal classes which appear to have very different average population sizes and substantially different average levels of nuclear nucleotide and amino acid diversity (Bazin et al. 2006).Unconditionally deleterious mutations can also depress linked neutral diversity by reducing the effective population size either through (1) background selection against relatively strongly selected mutations (Charlesworth et al. 1993, 1995) or (2) Hill–Robertson interference (Hill and Robertson 1966) among large numbers of relatively weakly selected mutations (reviewed by Comeron et al. 2008). The second of these processes, called “weak-selection Hill–Robertson interference” (wsHRi) by McVean and Charlesworth (2000) and “interference selection” (IS) by Comeron and Kreitman (2002), can shorten genealogies, give them strongly nonneutral branch-length proportions, and skew their topologies (Higgs and Woodcock 1995; Maia et al. 2004).To date, weak interference has mainly been studied by forward simulation, with the aim of assessing its possible effects on patterns of optimal synonymous codon use within eukaryotic nuclear genes and genomes, in the presence of recombination (Comeron and Guthrie 2005; Loewe and Charlesworth 2007; Comeron et al. 2008). In an attempt to understand the striking genealogical distortions seen in whale-louse mitochondria (Kaliszewska et al. 2005), we have developed a structured-coalescent algorithm that accurately models selection of arbitrary strength on a nonrecombining chromosome of finite length. All of the distortions seen in the whale-louse mitochondria are replicated under parameters that might plausibly apply to whale lice and many other animal species, and these distortions scale only weakly with population size.Whale lice are permanent, obligate ectoparasites of cetaceans. They feed on the dead outer surface of their host''s skin, and they appear to be harmless. They are amphipod Crustacea comprising a monophyletic family, Cyamidae, with ∼50 described species in several genera. Three of these species (Cyamus ovalis, C. gracilis, and C. erraticus) occur on right whales (Eubalaena spp.) but not regularly on any other hosts. Most adult right whales carry large populations of all three species.Right whales in the North Pacific, the North Atlantic, and the southern hemisphere have been separated for ∼5 million years, and so have their cyamids (Rosenbaum et al. 2000; Gaines et al. 2005; Kaliszewska et al. 2005). For this reason the right whales in different ocean systems are now considered distinct species (Eubalaena japonica, E. glacialis, and E. australis), and we refer to their cyamids as North Pacific C. ovalis, North Atlantic C. ovalis, southern C. ovalis, and so on, in anticipation that a future revision of the genus Cyamus will recognize them as “triplet” sibling species (3 × 3 = 9 species in all). We studied their mitochondrial population genetics with the initial aim of quantifying patterns of genetic differentiation among the cyamid populations on individual whales within local populations (Kaliszewska et al. 2005).We had reasoned (incorrectly) that the pattern of differentiation among whales might say something about their social interactions, since cyamids can transfer only between whales that are in direct physical contact with each other. We found very low levels of differentiation among whales and to our surprise literally no differentiation among the major southern hemisphere breeding aggregations that calve off the coasts of South America, South Africa, Australia, and New Zealand (Kaliszewska et al. 2005). This absence of population structure seems remarkable by terrestrial standards but is easily explained by modest rates of cyamid exchange among whales within local populations and between the major breeding aggregations, given the enormous sizes of cyamid populations. Right whales are highly gregarious (spending hours per day in social interactions), mobile (traveling thousands of kilometers per year on annual foraging migrations), and mortal (carrying their cyamid populations to the sea floor when they die). Thus cyamids have many opportunities to transfer between whales, and they might be expected to have evolved an inclination to do so when the opportunity presents itself (Hamilton and May 1977).The well-defined ecology of right-whale cyamids allows their population sizes to be estimated directly. The number of adult cyamids per whale (∼500–10,000, varying by species) times the number of whales per ocean (∼50,000–200,000, prior to human exploitation) equals the number of cyamids per species (Kaliszewska et al. 2005). Thus for all three nominal species of right-whale cyamids, long-term census population sizes are expected to have been in the range 2.5 × 10⁷–2 × 10⁹. Given conservative estimates of the per-generation mitochondrial mutation rate, even the lower end of this range predicts levels of synonymous nucleotide diversity at least an order of magnitude larger than those actually seen in the cyamids, which are consistently modest and similar to those seen in typical terrestrial arthropods (Kaliszewska et al. 2005). The three North Atlantic and southern hemisphere sibling-species pairs are strongly reciprocally monophyletic, as illustrated for C. ovalis in Figure 1. This is not expected at mutation–drift equilibrium, given their very large population sizes.Open in a separate windowFigure 1.—Mitochondrial gene genealogies for North Atlantic and southern hemisphere Cyamus ovalis, estimated by UPGMA from partial COI sequences. Left: The intraspecific genealogies coalesce globally at ∼0.5 and 1 MY and share an ancestor at ∼5 MY (Kaliszewska et al. 2005). Center and right: Each intraspecific genealogy is aligned with a generalized skyline plot (Strimmer and Pybus 2001) showing estimates of θ at different times in the past under a piecewise constant model of population size change fit by GENIE 3.0 (Pybus and Rambaut 2002). Three different point estimates of present-day θ are also indicated on the plots (synonymous-site nucleotide diversity π, Watterson''s θ estimated from synonymous sites, and θ estimated jointly with the apparent exponential growth rate by the MCMC coalescent algorithm in LAMARC). Values of Tajima''s D are lower (−1.5 to −1.6) when estimated from the sequences than when estimated from the branch lengths of the trees (−2.1 to −2.4), as expected because multiple substitutions occur at some sites. Similar values of D_T (−2.3 for both species) were obtained from the highest-likelihood trees found by BEAST with a fully parameterized GTR substitution model and a coalescent prior. Those trees have standardized imbalance statistics (−I_S) of −3.7 (southern hemisphere C. ovalis) and −2.2 (North Atlantic C. ovalis). The trees shown here have more extreme values of −I_S (−5.5 and −4.0, respectively), probably as a consequence of artificially pectinate branching orders induced by UPGMA among sets of identical sequences. Slightly different sets of sequences were used to make the two-species genealogy on the left and the single-species genealogies in the center. The long interspecific branches in the two-species tree are based on a multispecies maximum-likelihood analysis involving smaller numbers of much longer (4.1 kb) sequences (Kaliszewskaet al. 2005, Figure 3).These dramatic deficits of variation are not easily explained by population bottlenecks or adaptive selective sweeps. The bottleneck hypothesis is especially problematic because it requires the long-term near extinction of right whales in all three ocean systems. (Short bottlenecks such as those caused by human exploitation of right whales are not expected to have a noticeable effect on cyamid genetic diversity because cyamid populations remain large, and rates of genetic drift low, even when there are few whales.) That all three right whales have survived for millions of years suggests that they have maintained reasonably large population sizes, and the mitochondrial nucleotide diversity of southern right whales is consistent with this assumption (Kaliszewska et al. 2005), as is the nucleotide diversity of a cyamid nuclear gene (described below). Right whales eat copepods and krill, which are relatively close to the base of marine food webs, and right-whale populations are thought to be food limited. Thus a long-term, severe depression of their numbers would also seem to imply a collapse of marine ecosystems worldwide, for which there is no evidence.Several features of cyamid mitochondrial nucleotide diversity are also inconsistent with the bottleneck model and with adaptive sweeps as well. The most obvious of these features is the uniformity of cyamid mitochondrial diversity among species (π = 0.007–0.015 for COI sequences in the seven species surveyed by Kaliszewska et al. 2005). Gene genealogies estimated from these sequences also seem remarkably uniform in total depth, with last common ancestors differing in age by only a factor of 3 and in six of the seven species by less than a factor of 2 (see Kaliszewska et al. 2005, Figure 4). Adaptive sweeps might be expected to occur at roughly random intervals and not to be well coordinated in time among seven species in three different ocean systems. Interspecific coordination of such sweeps (over the whole globe) would seem to be required if they were to be a plausible primary cause of the genealogical shortening.Open in a separate windowFigure 4.—Apparent average effective sizes of ancestral populations, for models with different values of s. Values of N_e are given on the vertical axis (logarithmically scaled). They are estimated from the variance of expected contributions to the present, for the adults of any given generation. The parameters are those of Figure 2 (N = 65,536 = “64k,” μ = 1.5 × 10⁻⁶, L_s = 2048, U = 0.0031). In theory the curve for s = 0 should be perfectly horizontal. The discrepancy appears to be caused by subtle flaws in the shape of the very broad “idealized” simulated distribution that was used in this calculation. The curves for strong-selection cases are perfectly horizontal at times beyond a few thousand generations because the mutation-number distributions are compact, with little stochastic variation. The two lowest curves are those for the selection coefficients (s = 2⁻¹⁰, U/s = 3.2, and s = 2⁻¹¹, U/s = 6.3) that produce the most extreme values of the polymorphism and tree-shape statistics, given the other parameters (Figure 2).In addition to showing too little nucleotide variation, the cyamid mitochondrial genomes show strong and consistent excesses of rare nucleotide states, reflecting the “comb-like” or “star-like” shapes of the genealogies, in which deeper branches tend to be much too short relative to terminal branches (as if the trees had been “squished” from behind). This kind of distortion causes negative values of Tajima''s (1989) D and related statistics. It can be caused by population expansion from a bottleneck or by lineage expansion under positive selection (Kaplan et al. 1989; Slatkin and Hudson 1991; Rogers and Harpending 1992; Bamshad and Wooding 2003). However, the form of branch-length distortion seen in the cyamid genealogies suggests a slow, steady, roughly exponential form of population or lineage growth, not the relatively sudden increases suggested by the bottleneck and selective-sweep hypotheses. Generalized skyline plots (Strimmer and Pybus 2001) describing the histories of population size implied by the shapes of the northern and southern C. ovalis gene genealogies are shown in Figure 1. They are remarkably similar, as are the growth rates and estimates of present-day θ (= 2N_fμ_n) obtained by fitting exponential growth models using the coalescent algorithms in LAMARC (Kuhner et al. 1998, 2004) or BEAST (Drummond and Rambaut 2007).Cyamid populations cannot have grown in numbers as seemingly implied by these analyses. The number of cyamids on each whale appears to be set mainly by microhabitat limitations (e.g., by the area of rough callosity tissue on the head, where C. ovalis and C. gracilis live), and these features of their environment have hardly changed for millions of years, as demonstrated by the strong similarities of northern and southern right whales and their cyamids. Likewise, the numbers of right whales cannot have increased gradually from vanishingly small numbers over several hundred thousand years, for the reasons discussed above.The genealogical signals of “growth” therefore seem likely to be caused by selection. Environmental change is the most obvious potential cause of selection, but the apparent rate of growth seen here is strangely slow—in fact, slower than glacial. The orbitally forced Plio-Pleistocene glacial climate cycles have a major period of ∼100,000 years (Lambert et al. 2008 and references therein), but all seven of the right-whale cyamids for which we have mitochondrial population samples appear to have been “expanding,” more or less continuously, through at least several such cycles. The seemingly fairly consistent rate of branch-length foreshortening seen in the genealogies therefore suggests the action of a process that is relatively homogeneous in time, in addition to being very slow overall.The cyamid mitochondrial genealogies also appear to be topologically skewed, with sister clades too unequal in size, on average. In random bifurcating trees, the distribution of sister-clade sizes is uniform (Yule 1924; Heard 1992; Rogers 1994). Deviations from this null expectation can be quantified by statistics such as Colless''s (1982) index of tree imbalance (Shao and Sokal 1990; Rogers 1996). Our estimates of the cyamid genealogies tend to be excessively imbalanced (Figure 1). Strong topological imbalance is not caused by classic adaptive sweeps or by population growth following a bottleneck, but previous theoretical work has indicated that it can be caused by selection (Higgs and Woodcock 1995; Maia et al. 2004).To summarize, the cyamid mitochondrial genealogies are consistently much too short, too squished, and too skewed, relative to neutral-theory expectations. Owing to several special features of cyamid and right-whale biology, selection seems to be the only plausible explanation for this set of distortions, but conventional adaptive sweeps do not seem likely to be the primary cause. We therefore asked whether weakly deleterious mutations might be sufficient to generate the observed combination of patterns, in the absence of environmental change. Previous work (mentioned above) showed that interference among weak mutations at many sites can strongly affect linked neutral variation, but this work did not fully explore the parameter space relevant to our system or connect all the patterns in a genealogical setting.To address this question we first carried out forward simulations of populations of nonrecombining chromosomes with large numbers of nucleotide positions subject to forward and back mutations with unconditional fitness effects of size s. Large numbers of linked neutral sites were used to estimate genealogies and to calculate population statistics of interest. We found that for a range of intermediate values of s, considerable fitness variation was maintained and all three of the genealogical distortions (and the signal of apparent exponential growth) seen in the cyamid genealogies reached impressively large maxima. However, the computational burden of full forward simulation prevented us from considering realistic parameter values (i.e., large N and small μ), and it was not obvious that extrapolations based on compound parameters (e.g., Nμ and Ns) would work as hoped in all respects (see Comeron et al. 2008). We developed an equivalent coalescent algorithm that accurately reproduces all results of the forward simulations and allows for realistic parameter values. Under parameters relevant to cyamid mitochondria, the distortions of genealogical depth, proportions, and topology can be even more extreme than those seen in the cyamids, and the mean pairwise coalescence times (and resulting neutral nucleotide diversities) associated with maximally distorting intermediate values of s (U_s/s ∼ 5, where U_s is the total genomic mutation rate at sites with selection coefficents of size s) depend only weakly on N.All parameters of this model (including those of the environment) remain constant over time, yet in some respects it displays apparently nonequilibirum behavior. Under weak to intermediate selection (U_s/s > 10), the effective population size appears to become progressively smaller as time recedes into the past, giving rise to the illusion of growth. And in the maximally distorting range of intermediate selection coefficients, the distributions of deleterious mutation numbers and the shapes of genealogies show conspicuous dynamical instability of a form that could be taken to suggest “adaptive evolution” in response to episodes of environmental change. Adaptive mutations contribute importantly to this process, but they are reversions at some of the many sites previously mutated to mildly deleterious states. Subtle patterns of environmental change that converted previously optimal nucleotide states to slightly suboptimal states could give rise to a category of “virtual reversions” that would augment (or even outnumber) simple reversions, and the effects of such a process might well be consistent with the distortions seen in the cyamid mitochondrial genealogies. However, models with no environmental change of any kind appear to explain the observations surprisingly well.     

P2X4 Activation Modulates Volume-sensitive Outwardly Rectifying Chloride Channels in Rat Hepatoma Cells

Volume-sensitive outwardly rectifying (VSOR) Cl⁻ channels are critical for the regulatory volume decrease (RVD) response triggered upon cell swelling. Recent evidence indicates that H₂O₂ plays an essential role in the activation of these channels and that H₂O₂ per se activates the channels under isotonic isovolumic conditions. However, a significant difference in the time course for current onset between H₂O₂-induced and hypotonicity-mediated VSOR Cl⁻ activation is observed. In several cell types, cell swelling induced by hypotonic challenges triggers the release of ATP to the extracellular medium, which in turn, activates purinergic receptors and modulates cell volume regulation. In this study, we have addressed the effect of purinergic receptor activation on H₂O₂-induced and hypotonicity-mediated VSOR Cl⁻ current activation. Here we show that rat hepatoma cells (HTC) exposed to a 33% hypotonic solution responded by rapidly activating VSOR Cl⁻ current and releasing ATP to the extracellular medium. In contrast, cells exposed to 200 μm H₂O₂ VSOR Cl⁻ current onset was significantly slower, and ATP release was not detected. In cells exposed to either 11% hypotonicity or 200 μm H₂O₂, exogenous addition of ATP in the presence of extracellular Ca²⁺ resulted in a decrease in the half-time for VSOR Cl⁻ current onset. Conversely, in cells that overexpress a dominant-negative mutant of the ionotropic receptor P2X4 challenged with a 33% hypotonic solution, the half-time for VSOR Cl⁻ current onset was significantly slowed down. Our results indicate that, at high hypotonic imbalances, swelling-induced ATP release activates the purinergic receptor P2X4, which in turn modulates the time course of VSOR Cl⁻ current onset in a extracellular Ca²⁺-dependent manner.     

Phylogenetic relationships of Deschampsia antarctica (Poaceae): Insights from nuclear ribosomal ITS

Deschampsia antarctica E. Desv. is the only monocot in the Antarctic floristic zone. We evaluated the phylogenetic relationships of Deschampsia antarctica to other grasses using parsimony as the optimality criterion. Five different sets of gap, transversion and transitions costs were explored to analyze the effect of parameter choice on the phylogenetic results. Both internal transcribed spacers (ITS1 and ITS2) and the 5.8S subunit of nuclear ribosomal DNA were included in the analysis. A total of 43 species were analyzed including seven species of Deschampsia. Deschampsia antarctica forms a well supported group with five species of Deschampsia. Deschampsia does not appear monophyletic as D. flexuosa (L.) Trin. is not included in this clade. The clade to which D. antarctica belongs is sister to some Aveneae in all analyses. This study is the first contribution that evaluates the phylogenetic position of D. antarctica in relation to other species of Deschampsia.     

A preformed complex of postsynaptic proteins is involved in excitatory synapse development

Nonsynaptic clusters of postsynaptic proteins have been documented; however, their role remains elusive. We monitored the trafficking of several candidate proteins implicated in synaptogenesis, when nonsynaptic clusters of scaffold proteins are most abundant. We find a protein complex consisting of two populations that differ in their content, mobility, and involvement in synapse formation. One subpopulation is mobile and relies on actin transport for delivery to nascent and existing synapses. These mobile clusters contain the scaffolding proteins PSD-95, GKAP, and Shank. A proportion of mobile clusters that exhibits slow movement and travels short distances contains neuroligin-1. The second group consists of stationary nonsynaptic scaffold complexes that mainly contain neuroligin-1, can recruit synaptophysin-containing axonal transport vesicles, and are readily transformed to functional presynaptic contacts that recycle the vital dye FM 4-64. These results postulate a mechanism whereby preformed scaffold protein complexes serve as predetermined postsynaptic hotspots for establishment of new functional excitatory synapses.     

Reconciling Carbon-cycle Concepts, Terminology, and Methods

Recent projections of climatic change have focused a great deal of scientific and public attention on patterns of carbon (C) cycling as well as its controls, particularly the factors that determine whether an ecosystem is a net source or sink of atmospheric carbon dioxide (CO₂). Net ecosystem production (NEP), a central concept in C-cycling research, has been used by scientists to represent two different concepts. We propose that NEP be restricted to just one of its two original definitions—the imbalance between gross primary production (GPP) and ecosystem respiration (ER). We further propose that a new term—net ecosystem carbon balance (NECB)—be applied to the net rate of C accumulation in (or loss from [negative sign]) ecosystems. Net ecosystem carbon balance differs from NEP when C fluxes other than C fixation and respiration occur, or when inorganic C enters or leaves in dissolved form. These fluxes include the leaching loss or lateral transfer of C from the ecosystem; the emission of volatile organic C, methane, and carbon monoxide; and the release of soot and CO₂ from fire. Carbon fluxes in addition to NEP are particularly important determinants of NECB over long time scales. However, even over short time scales, they are important in ecosystems such as streams, estuaries, wetlands, and cities. Recent technological advances have led to a diversity of approaches to the measurement of C fluxes at different temporal and spatial scales. These approaches frequently capture different components of NEP or NECB and can therefore be compared across scales only by carefully specifying the fluxes included in the measurements. By explicitly identifying the fluxes that comprise NECB and other components of the C cycle, such as net ecosystem exchange (NEE) and net biome production (NBP), we can provide a less ambiguous framework for understanding and communicating recent changes in the global C cycle.     

Natural hybridization in Saxifraga callosa Sm

Saxifraga callosa Sm. is an evergreen perennial species distributed from Eastern Spain, through the Western Alps and the Apennines, to southern Italy. The existence of high morphological variation within different subspecies indicates that phenotypic characters are useful but not sufficient taxonomic tools. Indeed, available morphological data already suggested that S. callosa subentity lantoscana may be an outcross between S. callosa and S. cochlearis. In this work, by analyzing ITS (Internal Transcribed Sequences), AFLP (Amplified Fragment Length Polymorphisms), and cpDNA (chloroplast DNA) markers, a comprehensive study of the genomic relationships among S. callosa and related species has been carried out. The sequence of the ITS region of S. callosa subentity lantoscana gave no conclusive results on the taxonomy status of S. callosa subentity lantoscana. On the other hand, the use of the "NewHybrids" software to analyze an AFLP data-set (208 polymorphic amplified fragments) supported a significant posterior probability that S. callosa subentity lantoscana individuals are natural hybrids between S. callosa and S. cochlearis. The level of introgression of genes from alien genomes was confirmed by a simpler and quick methodology that analyze length variation in cpDNA sequences.     

Genetic diversity and origin of weedy rice (Oryza sativa f. spontanea) populations found in North-eastern China revealed by simple sequence repeat (SSR) markers

BACKGROUND AND AIMS: Weedy rice (Oryza sativa f. spontanea) is one of the most notorious weeds occurring in rice-planting areas worldwide. The objectives of this study are to determine the genetic diversity and differentiation of weedy rice populations from Liaoning Province in North-eastern China and to explore the possible origin of these weedy populations by comparing their genetic relationships with rice varieties (O. sativa) and wild rice (O. rufipogon) from different sources. METHODS: Simple sequence repeat (SSR) markers were used to estimate the genetic diversity of 30 weedy rice populations from Liaoning, each containing about 30 individuals, selected rice varieties and wild O. rufipogon. Genetic differentiation and the relationships of weedy rice populations were analysed using cluster analysis (UPGMA) and principle component analysis (PCA). KEY RESULTS: The overall genetic diversity of weedy rice populations from Liaoning was relatively high (H(e) = 0.313, I = 0.572), with about 35 % of the genetic variation found among regions. The Liaoning weedy rice populations were closely related to rice varieties from Liaoning and japonica varieties from other regions but distantly related to indica rice varieties and wild O. rufipogon. CONCLUSIONS: Weedy rice populations from Liaoning are considerably variable genetically and most probably originated from Liaoning rice varieties by mutation and intervarietal hybrids. Recent changes in farming practices and cultivation methods along with less weed management may have promoted the re-emergence and divergence of weedy rice in North-eastern China.     

Quantitative trait loci for grain moisture at harvest and field grain drying rate in maize (Zea mays, L.)

Hybrids with low grain moisture (GM) at harvest are specially required in mid- to short-season environments. One of the most important factors determining this trait is field grain drying rate (FDR). To produce hybrids with low GM at harvest, inbred lines can be obtained through selection for either GM or FDR. Thus, a single-cross population (181 F _2:3-generation plants) of two divergent inbred lines was evaluated to locate QTL affecting GM at harvest and FDR as a starting point for marker assisted selection (MAS). Moisture measurements were made with a hand-held moisture meter. Detection of QTL was facilitated with interval mapping in one and two dimensions including an interaction term, and a genetic linkage map of 122 SSR loci covering 1,557.8 cM. The markers were arranged in ten linkage groups. QTL mapping was made for the mean trait performance of the F _2:3 population across years. Ten QTL and an interaction were associated with GM. These QTL accounted for 54.8 and 65.2% of the phenotypic and genotypic variation, respectively. Eight QTL and two interactions were associated with FDR accounting for 35.7 and 45.2% of the phenotypic and genotypic variation, respectively. Two regions were in common between traits. The interaction between QTL for GM at harvest had practical implications for MAS. We conclude that MAS per se will not be an efficient method for reducing GM at harvest and/or increasing FDR. A selection index including both molecular marker information and phenotypic values, each appropriately weighted, would be the best selection strategy.