Genotype-by-management interactions for grain yield and grain protein concentration of wheat

The magnitude of genotype-by-management (G x M) interactions for grain yield and grain protein concentration was examined in a multi-environment trial (MET) involving a diverse set of 272 advanced breeding lines from the Queensland wheat breeding program. The MET was structured as a series of management-regimes imposed at 3 sites for 2 years. The management-regimes were generated at each site-year as separate trials in which planting time, N fertiliser application rate, cropping history, and irrigation were manipulated. irrigation was used to simulate different rainfall regimes. From the combined analysis of variance, the G x M interaction variance components were found to be the largest source of G x E interaction variation for both grain yield (0.117 +/- 0.005 t(2) ha(-2); 49% of total G x E 0.238 +/- 0.028 t(2) ha(-2)) and grain protein concentration (0.445 +/- 0.020%(2); 82% of total G x E 0.546 +/- 0.057%(2)), and in both cases this source of variation was larger than the genotypic variance component (grain yield 0.068 +/- 0.014 t(2) ha(-2) and grain protein 0.203 +/- 0.026%(2)). The genotypic correlation between the traits varied considerably with management-regime, ranging from -0.98 to -0.31, with an estimate of 0.0 for one trial. Pattern analysis identified advanced breeding lines with improved grain yield and grain protein concentration relative to the cultivars Hartog, Sunco and Meteor. It is likely that a large component of the previously documented G x E interactions for grain yield of wheat in the northern grains region are in part a result of G x M interactions. The implications of the strong influence of G x M interactions for the conduct of wheat breeding METs in the northern region are discussed. (C) 2001 Elsevier Science B.V. All rights reserved.

Power of the joint segregation analysis method for testing mixed major-gene and polygene inheritance models of quantitative traits

Understanding the genetic architecture of quantitative traits can greatly assist the design of strategies for their manipulation in plant-breeding programs. For a number of traits, genetic variation can be the result of segregation of a few major genes and many polygenes (minor genes). The joint segregation analysis (JSA) is a maximum-likelihood approach for fitting segregation models through the simultaneous use of phenotypic information from multiple generations. Our objective in this paper was to use computer simulation to quantify the power of the JSA method for testing the mixed-inheritance model for quantitative traits when it was applied to the six basic generations: both parents (P-1 and P-2), F-1, F-2, and both backcross generations (B-1 and B-2) derived from crossing the F-1 to each parent. A total of 1968 genetic model-experiment scenarios were considered in the simulation study to quantify the power of the method. Factors that interacted to influence the power of the JSA method to correctly detect genetic models were: (1) whether there were one or two major genes in combination with polygenes, (2) the heritability of the major genes and polygenes, (3) the level of dispersion of the major genes and polygenes between the two parents, and (4) the number of individuals examined in each generation (population size). The greatest levels of power were observed for the genetic models defined with simple inheritance; e.g., the power was greater than 90% for the one major gene model, regardless of the population size and major-gene heritability. Lower levels of power were observed for the genetic models with complex inheritance (major genes and polygenes), low heritability, small population sizes and a large dispersion of favourable genes among the two parents; e.g., the power was less than 5% for the two major-gene model with a heritability value of 0.3 and population sizes of 100 individuals. The JSA methodology was then applied to a previously studied sorghum data-set to investigate the genetic control of the putative drought resistance-trait osmotic adjustment in three crosses. The previous study concluded that there were two major genes segregating for osmotic adjustment in the three crosses. Application of the JSA method resulted in a change in the proposed genetic model. The presence of the two major genes was confirmed with the addition of an unspecified number of polygenes.

The effect of channel geometry and wall boundary conditions on the formation of extrusion surface instabilities for LLDPE

It is believed that surface instabilities can occur during the extrusion of linear low density polyethylene due to high extensional stresses at the exit of the die. Local crack development can occur at a critical stress level when melt rupture is reached. This high extensional stress results from the rearrangement of the flow at the boundary transition between the wall exit and the free surface. The stress is highest at the extrudate surface and decreases into the bulk of the material. The location of the region where the critical level is reached can determine the amplitude of the extrudate surface distortion, This paper studies the effect of wall slip on the numerically simulated extensional stress level at the die exit and correlates this to the experimentally determined amplitude of the surface instability. The effect of die exit radius and die wall roughness on extrusion surface instabilities is also correlated to the exit stress level in the same way. Whereas full slip may completely suppress the surface instability, a reduction in the exit stress level and instability amplitude is also shown for a rounded die exit and a slight increase in instability is shown to result from a rough die wall. A surface instability map demonstrates how the shear rate for onset of extrusion surface instabilities can be predicted on the basis of melt strength measurements and simulated stress peaks at the exit of the die. (C) 2001 Elsevier Science B.V. All rights reserved.

Stay-green: A consequence of the balance between supply and demand for nitrogen during grain filling?

Retention of green leaf area in grain sorghum under post-anthesis drought, known as stay-green, is associated with greater biomass production, lodging resistance and yield. The stay-green phenomenon can be examined at a cell, leaf, or whole plant level. At a cell level, the retention of chloroplast proteins such as LHCP2, OEC33 and Rubisco until late in senescence has been reported in sorghum containing the KS19 source of stay-green, indicating that photosynthesis may be maintained for longer during senescence in these genotypes. At a leaf level, longevity of photosynthetic apparatus is intimately related to nitrogen (N) status. At a whole plant level, stay-green can be viewed as a consequence of the balance between N demand by the grain and N supply during grain filling. To examine some of these concepts, nine hybrids varying in the B35 and KS19 sources of stay-green were grown under a postanthesis water deficit. Genotypic variation in delayed onset and reduced rate of leaf senescence were explained by differences in specific leaf nitrogen (SLN) and N uptake during grain filling. Matching N supply from age-related senescence and N uptake during grain tilling with grain N demand found that the shortfall in N supply for grain filling was greater in the senescent than stay-green hybrids, resulting in more accelerated leaf senescence in the former. We hypothesise that increased N uptake by stay-green hybrids is a result of greater biomass accumulation during grain filling in response to increased sink demand (higher grain numbers) which, in turn, is the result of increased radiation use efficiency and transpiration efficiency due to higher SLN. Delayed leaf senescence resulting from higher SLN should, in turn, allow snore carbon and nitrogen to be allocated to the roots of stay-green hybrids during grain filling, thereby maintaining a greater capacity to extract N from the soil compared with senescent hybrids.

Shoot control of hypernodulation and aberrant root formation in the har1-1 mutant of Lotus japonicus

The har1-1 mutant of Lotus japonicus B-129-S9 Gifu is characterized by two phenotypes: greater than normal nodulation (hypernodulation) and significantly inhibited root growth in the presence of its microsymbiont Mesorhizobium loti strain NZP2235. We demonstrate that the two traits co-segregate, suggesting a single genetic alteration involving developmental pleiotropy. A cross between the mutant and genotype Funakura (with wild-type root and nodule morphology) demonstrated Mendelian recessive segregation of both phenotypes (root and nodule) in 216 F2 individuals. Using DNA-amplification fingerprinting polymorphisms in Gifu har1-1 and Funakura, the mutant locus was positioned between two markers at about 7 and 13 cM distance. Reciprocal hypocotyl grafting of shoots and roots showed that the hypernodulation and reduced root phenotypes are both predominantly controlled by the shoot.

Genetic linkage analysis of the lesser grain borer Rhyzopertha dominica identifies two loci that confer high-level resistance to the fumigant phosphine

High levels of inheritable resistance to phosphine in Rhyzopertha dominica have recently, been detected in Australia and hi art effort to isolate the genes responsible For resistance we have used random amplified DNA fingerprinting (RAF) to produce a genetic linkage map of R. dominica. The map consists of 94 dominant DNA markers with art average distance between markers of 4.6 cM and defines nine linkage groups with a total recombination distance of 390.1 cM. We have identified two loci that are responsible for high-level resistance. One provides similar to50x resistance to phosphine while the other provides 12.5x resistance and in combination, the two genes act synergistically to provide a resistance level 250 x greater than that of fully susceptible beetles. The haploid genome size has been determined to be 4.76 x 10(8) bp, resulting in an average physical distance of 1.2 Mbp per map unit. No recombination has been observed between either of the two resistance loci and their adjacent DNA markers in a population of 44 fully resistant F-5 individuals, which indicates that the genes are likely to reside within 0.91 cM (1.1 Mbp) of the DNA markers.

Boundary value problems for systems of difference equations associated with systems of second-order ordinary differential equations

We study difference equations which arise as discrete approximations to two-point boundary value problems for systems of second-order ordinary differential equations. We formulate conditions which guarantee a priori bounds on first differences of solutions to the discretized problem. We establish existence results for solutions to the discretized boundary value problems subject to nonlinear boundary conditions. We apply our results to show that solutions to the discrete problem converge to solutions of the continuous problem in an aggregate sense. (C) 2002 Elsevier Science Ltd. All rights reserved.

Existence of multiple solutions for second-order discrete boundary value problems

We give conditions on f involving pairs of discrete lower and discrete upper solutions which lead to the existence of at least three solutions of the discrete two-point boundary value problem yk+1 - 2yk + yk-1 + f (k, yk, vk) = 0, for k = 1,..., n - 1, y0 = 0 = yn,, where f is continuous and vk = yk - yk-1, for k = 1,..., n. In the special case f (k, t, p) = f (t) greater than or equal to 0, we give growth conditions on f and apply our general result to show the existence of three positive solutions. We give an example showing this latter result is sharp. Our results extend those of Avery and Peterson and are in the spirit of our results for the continuous analogue. (C) 2002 Elsevier Science Ltd. All rights reserved.

Difference equations associated with fully nonlinear boundary value problems for second order ordinary differential equations

We study the continuous problem y"=f(x,y,y'), xc[0,1], 0=G((y(0),y(1)),(y'(0), y'(1))), and its discrete approximation (y(k+1)-2y(k)+y(k-1))/h(2) =f(t(k), y(k), v(k)), k = 1,..., n-1, 0 = G((y(0), y(n)), (v(1), v(n))), where f and G = (g(0), g(1)) are continuous and fully nonlinear, h = 1/n, v(k) = (y(k) - y(k-1))/h, for k =1,..., n, and t(k) = kh, for k = 0,...,n. We assume there exist strict lower and strict upper solutions and impose additional conditions on f and G which are known to yield a priori bounds on, and to guarantee the existence of solutions of the continuous problem. We show that the discrete approximation also has solutions which approximate solutions of the continuous problem and converge to the solution of the continuous problem when it is unique, as the grid size goes to 0. Homotopy methods can be used to compute the solution of the discrete approximation. Our results were motivated by those of Gaines.

Shear stress and sediment transport calculations for sheet flow under waves

A simple method is provided for calculating transport rates of not too fine (d(50) greater than or equal to 0.20 mm) sand under sheet flow conditions. The method consists of a Meyer-Peter-type transport formula operating on a time-varying Shields parameter, which accounts for both acceleration-asymmetry and boundary layer streaming. While velocity moment formulae, e.g.., = Constant x calibrated against U-tube measurements, fail spectacularly under some real waves (Ribberink, J.S., Dohmen-Janssen, C.M., Hanes, D.M., McLean, S.R., Vincent, C., 2000. Near-bed sand transport mechanisms under waves. Proc. 27th Int. Conf. Coastal Engineering, Sydney, ASCE, New York, pp. 3263-3276, Fig. 12), the new method predicts the real wave observations equally well. The reason that the velocity moment formulae fail under these waves is partly the presence of boundary layer streaming and partly the saw-tooth asymmetry, i.e., the front of the waves being steeper than the back. Waves with saw-tooth asymmetry may generate a net landward sediment transport even if = 0, because of the more abrupt acceleration under the steep front. More abrupt accelerations are associated with thinner boundary layers and greater pressure gradients for a given velocity magnitude. The two real wave effects are incorporated in a model of the form Q(s)(t) = Q(s)[theta(t)] rather than Q(S)(t) = Q(S)[u(infinity)(t)], i.e., by expressing the transport rate in terms of an instantaneous Shields parameter rather than in terms of the free stream velocity, and accounting for both streaming and accelerations in the 0(t) calculations. The instantaneous friction velocities u(*)(t) and subsequently theta(t) are calculated as follows. Firstly, a linear filter incorporating the grain roughness friction factor f(2.5) and a phase angle phi(tau) is applied to u(infinity)(t). This delivers u(*)(t) which is used to calculate an instantaneous grain roughness Shields parameter theta(2.5)(t). Secondly, a constant bed shear stress is added which corresponds to the streaming related bed shear stress -rho ($) over bar((u) over tilde(w) over tilde)(infinity) . The method can be applied to any u(infinity)(t) time series, but further experimental validation is recommended before application to conditions that differ strongly from the ones considered below. The method is not recommended for rippled beds or for sheet flow with typical prototype wave periods and d(50) < 0.20 turn. In such scenarios, time lags related to vertical sediment movement become important, and these are not considered by the present model. (C) 2002 Elsevier Science B.V. All rights reserved.

Vertical fluxes of sediment in oscillatory sheet flow

Time series of vertical sediment fluxes are derived from concentration time series in sheet flow under waves. While the concentrations C(z,t) vary very little with time for \z\ < 10d(50), the measured vertical sediment fluxes Q(zs)(z,t) vary strongly with time in this vertical band and their time variation follows, to some extent, the variation of the grain roughness Shields parameter 02,5(t). Thus, sediment distribution models based on the pickup function boundary condition are in some qualitative agreement with the measurements. However, the pickup function models are only able to model the upward bursts of sediment during the accelerating phases of the flow. They are, so far, unable to model the following strong downward sediment fluxes, which are observed during the periods of flow deceleration. Classical pickup functions, which essentially depend on the Shields parameter, are also incapable of modelling the secondary entrainment fluxes, which sometimes occur at free stream velocity reversal. The measured vertical fluxes indicate that the effective sediment settling velocity in the high [(0.3 < C(z,t) < 0.4] concentration area is typically only a few percent of the clear water settling velocity, while the measurements of Richardson and Jeronimo [Chem. Eng. Sci. 34 (1979) 1419], from a different physical setting, lead to estimates of the order 20%. The data does not support gradient diffusion as a model for sediment entrainment from the bed. That is, detailed modelling of the observed near-bed fluxes would require diffusivities that go negative during periods of flow deceleration. An observed general trend for concentration variability to increase with elevation close to the bed is also irreconcilable with diffusion models driven by a bottom boundary condition. (C) 2002 Published by Elsevier Science B.V.

A real symmetric Lanczos subspace implementation of quantum scattering using boundary inhomogeneities

We present an efficient and robust method for calculating state-to-state reaction probabilities utilising the Lanczos algorithm for a real symmetric Hamiltonian. The method recasts the time-independent Artificial Boundary Inhomogeneity technique recently introduced by Jang and Light (J. Chem. Phys. 102 (1995) 3262) into a tridiagonal (Lanczos) representation. The calculation proceeds at the cost of a single Lanczos propagation for each boundary inhomogeneity function and yields all state-to-state probabilities (elastic, inelastic and reactive) over an arbitrary energy range. The method is applied to the collinear H + H-2 reaction and the results demonstrate it is accurate and efficient in comparison with previous calculations. (C) 2002 Elsevier Science B.V. All rights reserved.

Boundary-layer transition measurements in hypervelocity flows in a shock tunnel

Experiments to investigate the transition process in hypervelocity boundary layers were performed in the T4 free-piston shock tunnel. An array of thin-film heat-transfer gauges was used to detect the location and extent of the transitional region on a 1500 mm long x 120 turn wide flat plate, which formed one of the walls of a duct. The experiments were performed in a Mach 6 flow of air with 6- and 12-MJ/kg nozzle-supply enthalpies at unit Reynolds numbers ranging from 1.6 x 10(6) to 4.9 x 10(6) m(-1). The results show that the characteristics typical of transition taking place through the initiation, growth, and merger of turbulent spots are evident in the heat-transfer signals. A 2-mm-high excrescence located 440 turn from the leading edge was found to be capable of generating a turbulent wedge within an otherwise laminar boundary layer at a unit Reynolds number of 2.6 x 10(6) m(-1) at the 6-MJ/kg condition. A tripping strip, located 100 mm from the leading edge and consisting of a line 37 teeth of 2 rum height equally spaced and spanning the test surface, was also found to be capable of advancing the transition location at the same condition and at the higher enthalpy condition.

Columnar to equiaxed transition of eutectic in hypoeutectic aluminium-silicon alloys

Directional solidification of unmodified and strontium modified binary, high-purity aluminium-7 wt% silicon and commercial A356 alloys has been carried out to investigate the mechanism of eutectic solidification. The microstructure of the eutectic growth inter-face was investigated with optical microscopy and Electron Backscattering Diffraction (EBSD). In the commercial alloys, the eutectic solidification inter-face extends in the growth direction and creates a eutectic mushy zone. A planar eutectic growth front is observed in the high-purity alloys. The eutectic aluminium has mainly the same crystallographic orientation as the dendrites in the unmodified alloys and the strontium modified high-purity alloy. A more complex eutectic grain structure is found in the strontium modified commercial alloy. A mechanism involving constitutional undercooling and a columnar to equiaxed transition explains the differences between pure and commercial alloys. It is probably caused by the segregation of iron and magnesium and the activation of nucleants in the commercial alloy. (C) 2002 Acta Materialia Inc. Published by Elsevier Science Ltd. All rights reserved.