Analysis of MGE Dm412 localization pattern in Drosophila melanogaster lines undergoing long-term selection for a quantitative character]

The subpopulation composed of the mixture of Drosophila isogenic lines with interrupted wing radial vein (mutation radius incompletus, ri) was subjected to long-term selection in different directions for increase or decrease in expression of the ri gene. As a result, the lines with contrasting different values of mean character phenotype were developed. The isogenic lines of mean character phenotype were developed. The isogenic lines and F2 from their crosses with selected lines were analysed by the pattern of copia-like MGE DM412 localization. The isogenic lines were shown to have individual pattern, the selected lines differing strongly from them. Selection led to the loss of Dm412 localization sites during negative selection, while positive selection results both in loss and acquisition of sites. Correlation between the phenotype of the quantitative character and the pattern of MGE Dm412 was revealed.     

Change in distribution of mobile genetic elements in genome of <Emphasis Type="Italic">Drosophila melanogaster</Emphasis>: Cause or consequence of selection for quantitative traits?

The distribution pattern of the hobo transposon and Dm412 retrotransposon hybridization sites on the salivary gland polytene chromosomes from larvae of the Drosophila melanogaster isogenic strain 51 used to analyze the effect of the transposition of transposable elements (TEs) on selection for quantitative traits was studied. It was shown that no more than half of the Dm412 hybridization sites were retained 15 years after isogenization; the frequency of the Dm412 transposition varied from 2.0 × 10⁻⁴ to 8.8 × 10⁻⁵ sites per genome for generation depending on whether the appearance of the same hybridization sites in a part of individuals was considered as independent events or as the manifestation of the appearing sample heterogeneity. The distribution patterns of hobo hybridization sites in two isofemale strains derived from the isogenic strain 51 differed much more markedly; the number of the hobo sites in one of the derivative strains was threefold smaller than in the other one and only some of the sites were common. Within each derivative strain, the TE distribution was uniform, which suggests that inbreeding had no effect on Dm412 activity in this strain. The rates of change in the distribution patterns of various TE in the strain 51 corresponded to their spontaneous transposition rates. Since the isogenic strain accumulates polymorphism in the TE distribution without selection, the TEs are more likely to be the markers of selection events rather than their inducers. Thus, when studying the effects of various environmental factors on TE transposition even in isogenic strains, it is necessary to perform additional close inbreeding to reduce the potential polymorphism.     

Differential gene expression of wheat progeny with contrasting levels of transpiration efficiency

High water use efficiency or transpiration efficiency (TE) in wheat is a desirable physiological trait for increasing grain yield under water-limited environments. The identification of genes associated with this trait would facilitate the selection for genotypes with higher TE using molecular markers. We performed an expression profiling (microarray) analysis of approximately 16,000 unique wheat ESTs to identify genes that were differentially expressed between wheat progeny lines with contrasting TE levels from a cross between Quarrion (high TE) and Genaro 81 (low TE). We also conducted a second microarray analysis to identify genes responsive to drought stress in wheat leaves. Ninety-three genes that were differentially expressed between high and low TE progeny lines were identified. One fifth of these genes were markedly responsive to drought stress. Several potential growth-related regulatory genes, which were down-regulated by drought, were expressed at a higher level in the high TE lines than the low TE lines and are potentially associated with a biomass production component of the Quarrion-derived high TE trait. Eighteen of the TE differentially expressed genes were further analysed using quantitative RT-PCR on a separate set of plant samples from those used for microarray analysis. The expression levels of 11 of the 18 genes were positively correlated with the high TE trait, measured as carbon isotope discrimination (Δ¹³C). These data indicate that some of these TE differentially expressed genes are candidates for investigating processes that underlie the high TE trait or for use as expression quantitative trait loci (eQTLs) for TE. Electronic Supplementary Material Supplementary material is available for this article at     

Response of the MGE 412 pattern to truncation selection of a quantitative trait in an isogenic line of drosophila after severe heat shock (SHS)

Positive and negative selection on the total length of two fragments of an interrupted longitudinal wing vein in an isogenic line of Drosophila melanogaster was accompanied by changes in the genomic localization pattern of MGE 412. Strong truncation selection was conducted in the population of effective size Ne = 160 for 50 generations. Twenty-six out of 35 polymorphic HHS-induced segments of MGE localization behaved as independent copies and markers, whereas 9 segments proved to be selective. The second group included "hot" segments of HHS transposition induction (43B, 97E, etc.). Thus, final consensus patterns of induced MGE transpositions have a random and an adaptive component in generation 50 of positive and negative selection. Selective patterns probably include modifier MGEs, which generate induced genetic regulatory variation of polygenes controlling the selected quantitative trait in the isogenic line after HHS.     

Computer Simulation of Concerted Evolution of Patterns of Polygenes, Transposable Elements, and Origin Identity Labels: Negative Selection

Computer simulation of the population dynamics of the genomic patterns of polygenes, transposable elements (TEs), and origin identity labels (OILs) in the course of negative selection for an additive quantitative trait has been performed. It was demonstrated that active polygene alleles disappear very rapidly, whereas the patterns of TEs and OILs continue their evolution determined by strict selective inbreeding and gene drift. Dendrograms of the patterns of polygenes, TEs, and OILs were constructed for all generations. It was demonstrated that the final consensus pattern of OILs consists of the fragments of the original patterns, which contain neither active polygene alleles nor modifier or marker TEs. Neutral TE copies were present in the final pattern, as should be expected in the case of gene drift. Inbreeding coefficient increased steadily but by generation 100 reached values higher than 0.9. All other parameters and initial conditions being the same, the responses to negative and positive selections were asymmetric.     

Negative selection and computer models of the joint evolution of the patterns of polygenes,transposable elements,and origin identity labels

Computer simulation of the population dynamics of the genomic patterns of polygenes, transposable elements (TEs), and origin identity labels (OILs) in the course of negative selection for an additive quantitative trait has been performed. It was demonstrated that active polygene alleles disappear very rapidly, whereas the patterns of TEs and OILs continue their evolution determined by strict selective inbreeding and gene drift. Dendrograms of the patterns of polygenes, TEs, and OILs were constructed for all generations. It was demonstrated that the final consensus pattern of OILs consists of the fragments of the original patterns, which contain neither active polygene alleles nor modifier or marker TEs. Neutral TE copies were present in the final pattern, as should be expected in the case of gene drift. Inbreeding coefficient increased steadily but by generation 100 reached values higher than 0.9. All other parameters and initial conditions being the same, the responses to negative and positive selections were asymmetric.     

DENSITY-DEPENDENT SELECTION ON CONTINUOUS CHARACTERS: A QUANTITATIVE GENETIC MODEL

A quantitative genetic model of density-dependent selection is presented and analysed with parameter values obtained from laboratory selection experiments conducted by Mueller and his coworkers. The ecological concept of r- and K-selection is formulated in terms of selection gradients on underlying phenotypic characters that influence the density-dependent measure of fitness. Hence the selection gradients on traits are decomposed into two components, one that changes in the direction to increase r, and one that changes in the direction to increase K. The relative importance of the two components is determined by temporal fluctuations in population density. The evolutionary rate of r and K (per-generation changes in r and K due to the genetic responses of the underlying traits) is also formulated. Numerical simulation has shown that with moderate genetic variances of the underlying characters, r and K can evolve rapidly and the evolutionary rate is influenced by synergistic interaction between characters that contribute to r and K. But strong r-selection can occur only with severe and continuous disturbances of populations so that the population density is kept low enough to prevent K-selection.     

Maintenance of transposable element copy number in natural populations of Drosophila melanogaster and D. simulans

To investigate the main forces controlling the containment of transposable elements (TE) in natural populations, we analyzed the copia, mdg1, and 412 elements in various populations of Drosophila melanogaster and D. simulans. A lower proportion of insertion sites on the X chromosome in comparison with the autosomes suggests that selection against the detrimental effects of TE insertions is the major force containing TE copies in populations of Drosophila. This selection effect hypothesis is strengthened by the absence of the negative correlation between recombination rate and TE copy number along the chromosomes, which was expected under the alternative ectopic exchange model (selection against the deleterious rearrangements promoted by recombination between TE insertions). A cline in 412 copy number in relation to latitude was observed among the natural populations of D. simulans, with very high numbers existing in some local populations (around 60 copies in a sample from Canberra, Australia). An apparent absence of selection effects in this Canberra sample and a value of transposition rate equal to 1–2 × 10^-3 whatever the population and its copy number agree with the idea of recent but temporarily drastic TE movements in local populations. The high values of transposition rate in D. simulans clearly disfavor the hypothesis that the low amount of transposable elements in this species could result from a low transposition rate. This revised version was published online in August 2006 with corrections to the Cover Date.     

Change in the system of wing venation of Drosophila melanogaster under heat shock and selection

A change in the system of wing venation of Drosophila melanogaster appeared in response to heat shock and positive selection pressure directed to restoring the normal formation of wing radial vein, L2, that had been violated earlier by the recessive mutation of radius incompletes. Positive selection was effective, L2 having been formed correctly and completely to 35 generation. Besides, (+)-selection was accompanied by appearance of a small fragment of an additional vein at the wing tip. Selection directed to increase of size of this fragment resulted in the sufficient changes in the system of wing venation as a whole. It is suggested that, during evolution, transformation of wing venation of Drosophila was effected by the change of the way of prepattern realization, whereas the cells of wing plate continued to allow the formation of veins practically over a whole area.     

Truncation family selection and nonsystematic inbreeding leads to a rapid fixation of a pattern of mobile genetic elements in a computer model

A computer simulation model of the population dynamics of a polygenic system and a pattern of mobile genetic elements (MGEs) under directional truncation selection for a quantitative trait was developed. Modifier MGEs were shown to be rapidly and adaptively fixed (or lost) together with the modified polygenes. Marker MGEs and independent MGE copies were fixed and lost just as rapidly but in a random manner. Using specific marking of initial haploid genomes and direct computing of the mean proportion of identical encounters at each locus in each generation, it was shown that the mean nonselective inbreeding coefficient F(n) dramatically increases in the course of selection, reaching values 0.7-0.9 in 15-20 generations. As a result, adaptive homozygotization of polygenes and modifier MGEs and random homozygotization of marker MGEs, independent MGE copies, and all other genes of the genome occurs. These results confirm the hypothesis on the "champion" polygene pattern advanced earlier to explain the data of selection experiments.     

The role of the von Willebrand factor in renal diseases and haemodialysis patients

During a period of twenty years, the von Willebrand factor (VWf) biological activity was evaluated in 805 patients with vein thrombosis, diabetes mellitus, chronic renal failure and ischemic heart disease. The examined patients were 168 with vein thrombosis, 129 with diabetes mellitus, 412 with chronic renal failure (CRF), and 96 with ischemic heart disease. The biological activity was also determined in 104 haemodialysis patients using four different haemodialytic membranes: 30 on cuprophan membrane, 30 on polymethylmetacrylate membrane (PMMA), 24 on hemophane and 20 patients on polysulphone (PS) membrane. In 42 patients with arterio-venous fistula prone to thrombosis, the biological activity of the von Willebrand Factor was 178% in comparison to 106% in the control group. The biological activity of VWF was increased in patients with vein thrombosis (p < 0.02), in patients with diabetes mellitus (p < 0.01), CRF (p < 0.05), and in patients with ischemic heart disease (p < 0.01). The highest biological activity was found in patients on PMMA (p < 0.001), then cuprophan (p < 0.05) and hemophane membrane (p < 0.01), while the lowest increase of its concentration was noticed in patients on PS without statistical significance. In arteriovenous fistula prone to thrombosis patients biological activity of the von Willebrand Factor was significantly increased (p < 0.01). Our investigations show the importance of VWF as a marker of endothelial disfunction, a possible predictor of A-V fistula thrombosis, and a possible marker of haemodialysis membranes biocompatibility.     

Spatial differences in patterns of modification: selection on hairy in Drosophila melanogaster wings

Artificial selection was carried out for over 45 generations to enhance and suppress expression of the mutation hairy on the Drosophila melanogaster wing. Whole chromosome mapping of X‐linked and autosomal modifiers of sense organ number displayed regional differences in magnitude and direction of their effects. Regional specificity of modifier effects was also seen in some interchromosomal interactions. Scanning electron microscopy allowed precise measurement of sense organ size and position along the L3 longitudinal wing vein. Sense organ size varied in a predictable fashion along the proximal–distal axis, and the dorsal pattern differed from the ventral pattern. The high and low selection lines differed most in the proximal portion of the L3 vein. Extra sense organs in the High line were often associated with vein fragments at locations predicted from ancestral vein patterns. Thus, regional specificity of polygenic or quantitative trait locus modifier effects was identified in several different parts of the wing. This revised version was published online in July 2006 with corrections to the Cover Date.     

A QUANTITATIVE-GENETIC MODEL FOR SELECTION ON DEVELOPMENTAL NOISE

We propose a simple model for analyzing the effects of microenvironmental variation in quantitative genetics. Our model assumes that the sensitivity of the phenotype to fluctuations in microenvironment has a genetic basis and allows for genetic correlation between trait value and microenvironmental sensitivity. We analyze the effects of short-term stabilizing and directional selection on the genotypic and microenvironmental components of phenotypic variance. Our model predicts that stabilizing selection on a quantitative trait increases developmental canalization. We show that stabilizing selection can result in an increase in the heritability. Our findings may provide an explanation for the results of selection experiments in which artificial stabilizing selection did not change the heritability coefficient or increased it.     

MHC polymorphism and disease resistance to vibrio anguillarum in 8 families of half-smooth tongue sole (Cynoglossus semilaevis)

Background

Most quantitative traits are controlled by multiple quantitative trait loci (QTL). The contribution of each locus may be negligible but the collective contribution of all loci is usually significant. Genome selection that uses markers of the entire genome to predict the genomic values of individual plants or animals can be more efficient than selection on phenotypic values and pedigree information alone for genetic improvement. When a quantitative trait is contributed by epistatic effects, using all markers (main effects) and marker pairs (epistatic effects) to predict the genomic values of plants can achieve the maximum efficiency for genetic improvement.

Results

In this study, we created 126 recombinant inbred lines of soybean and genotyped 80 makers across the genome. We applied the genome selection technique to predict the genomic value of somatic embryo number (a quantitative trait) for each line. Cross validation analysis showed that the squared correlation coefficient between the observed and predicted embryo numbers was 0.33 when only main (additive) effects were used for prediction. When the interaction (epistatic) effects were also included in the model, the squared correlation coefficient reached 0.78.

Conclusions

This study provided an excellent example for the application of genome selection to plant breeding.     

Selective sweep at a quantitative trait locus in the presence of background genetic variation

We model selection at a locus affecting a quantitative trait (QTL) in the presence of genetic variance due to other loci. The dynamics at the QTL are related to the initial genotypic value and to the background genetic variance of the trait, assuming that background genetic values are normally distributed, under three different forms of selection on the trait. Approximate dynamics are derived under the assumption of small mutation effect. For similar strengths of selection on the trait (i.e, gradient of directional selection beta) the way background variation affects the dynamics at the QTL critically depends on the shape of the fitness function. It generally causes the strength of selection on the QTL to decrease with time. The resulting neutral heterozygosity pattern resembles that of a selective sweep with a constant selection coefficient corresponding to the early conditions. The signature of selection may also be blurred by mutation and recombination in the later part of the sweep. We also study the race between the QTL and its genetic background toward a new optimum and find the conditions for a complete sweep. Overall, our results suggest that phenotypic traits exhibiting clear-cut molecular signatures of selection may represent a biased subset of all adaptive traits.     

Deleterious Mutations, Apparent Stabilizing Selection and the Maintenance of Quantitative Variation

Apparent stabilizing selection on a quantitative trait that is not causally connected to fitness can result from the pleiotropic effects of unconditionally deleterious mutations, because as N. Barton noted, "...individuals with extreme values of the trait will tend to carry more deleterious alleles...." We use a simple model to investigate the dependence of this apparent selection on the genomic deleterious mutation rate, U; the equilibrium distribution of K, the number of deleterious mutations per genome; and the parameters describing directional selection against deleterious mutations. Unlike previous analyses, we allow for epistatic selection against deleterious alleles. For various selection functions and realistic parameter values, the distribution of K, the distribution of breeding values for a pleiotropically affected trait, and the apparent stabilizing selection function are all nearly Gaussian. The additive genetic variance for the quantitative trait is kQa2, where k is the average number of deleterious mutations per genome, Q is the proportion of deleterious mutations that affect the trait, and a2 is the variance of pleiotropic effects for individual mutations that do affect the trait. In contrast, when the trait is measured in units of its additive standard deviation, the apparent fitness function is essentially independent of Q and a2; and beta, the intensity of selection, measured as the ratio of additive genetic variance to the "variance" of the fitness curve, is very close to s = U/k, the selection coefficient against individual deleterious mutations at equilibrium. Therefore, this model predicts appreciable apparent stabilizing selection if s exceeds about 0.03, which is consistent with various data. However, the model also predicts that beta must equal Vm/VG, the ratio of new additive variance for the trait introduced each generation by mutation to the standing additive variance. Most, although not all, estimates of this ratio imply apparent stabilizing selection weaker than generally observed. A qualitative argument suggests that even when direct selection is responsible for most of the selection observed on a character, it may be essentially irrelevant to the maintenance of variation for the character by mutation-selection balance. Simple experiments can indicate the fraction of observed stabilizing selection attributable to the pleiotropic effects of deleterious mutations.     

DENSITY-DEPENDENT EVOLUTION OF LIFE-HISTORY TRAITS IN DROSOPHILA MELANOGASTER

Populations of Drosophila melanogaster were maintained for 36 generations in r- and K-selected environments in order to test the life-history predictions of theories on density-dependent selection. In the r-selection environment, populations were reduced to low densities by density-independent adult mortality, whereas populations in the K-selection environment were maintained at their carrying capacity. Some of the experimental results support the predictions or r- and K-selection theory; relative to the r-selected populations, the K-selected populations evolved an increased larval-to-adult viability, larger body size, and longer development time at high larval densities. Mueller and Ayala (1981) found that K-selected populations also have a higher rate of population growth at high densities. Other predictions of the thoery are contradicted by the lack of differences between the r and K populations in adult longevity and fecundity and a slower rate of development for r-selected individuals at low densities. The differences between selected populations in larval survivorship, larval-to-adult development time, and adult body size are strongly dependent on larval density, and there is a significant interaction between populations and larval density for each trait. This manifests an inadequacy of the theory on r- and K-selection, which does not take into account such interactions between genotypes and environments. We describe mechanisms that may explain the evolution of preadult life-history traits in our experiment and discuss the need for changes in theories of density-dependent selection.     

Genetic basis of kernel nutritional traits during maize domestication and improvement

The nutritional traits of maize kernels are important for human and animal nutrition, and these traits have undergone selection to meet the diverse nutritional needs of humans. However, our knowledge of the genetic basis of selecting for kernel nutritional traits is limited. Here, we identified both single and epistatic quantitative trait loci (QTLs) that contributed to the differences of oil and carotenoid traits between maize and teosinte. Over half of teosinte alleles of single QTLs increased the values of the detected oil and carotenoid traits. Based on the pleiotropism or linkage information of the identified single QTLs, we constructed a trait–locus network to help clarify the genetic basis of correlations among oil and carotenoid traits. Furthermore, the selection features and evolutionary trajectories of the genes or loci underlying variations in oil and carotenoid traits revealed that these nutritional traits produced diverse selection events during maize domestication and improvement. To illustrate more, a mutator distance–relative transposable element (TE) in intron 1 of DXS2, which encoded a rate‐limiting enzyme in the methylerythritol phosphate pathway, was identified to increase carotenoid biosynthesis by enhancing DXS2 expression. This TE occurs in the grass teosinte, and has been found to have undergone selection during maize domestication and improvement, and is almost fixed in yellow maize. Our findings not only provide important insights into evolutionary changes in nutritional traits, but also highlight the feasibility of reintroducing back into commercial agricultural germplasm those nutritionally important genes hidden in wild relatives.     

Calcium activation of heart mitochondrial oxidative phosphorylation: rapid kinetics of mVO2, NADH, AND light scattering

Parallel activation of heart mitochondria NADH and ATP production by Ca(2+) has been shown to involve the Ca(2+)-sensitive dehydrogenases and the F(0)F(1)-ATPase. In the current study we hypothesize that the response time of Ca(2+)-activated ATP production is rapid enough to support step changes in myocardial workload ( approximately 100 ms). To test this hypothesis, the rapid kinetics of Ca(2+) activation of mV(O(2)), [NADH], and light scattering were evaluated in isolated porcine heart mitochondria at 37 degrees C using a variety of optical techniques. The addition of Ca(2+) was associated with an initial response time (IRT) of mV(O(2)) that was dose-dependent with a minimum IRT of 0.27 +/- 0.02 s (n = 41) at 535 nm Ca(2+). The IRTs for NADH fluorescence and light scattering in response to Ca(2+) additions were similar to mV(O(2)). The Ca(2+) IRT for mV(O(2)) was significantly shorter than 1.6 mm ADP (2.36 +/- 0.47 s; p < or = 0.001, n = 13), 2.2 mm P(i) (2.32 +/- 0.29, p < or = 0.001, n = 13), or 10 mm creatine (15.6.+/-1.18 s, p < or = 0.001, n = 18) under similar experimental conditions. Calcium effects were inhibited with 8 microm ruthenium red (2.4 +/- 0.31 s; p < or = 0.001, n = 16) and reversed with EGTA (1.6 +/- 0.44; p < or = 0.01, n = 6). Estimates of Ca(2+) uptake into mitochondria using optical Ca(2+) indicators trapped in the matrix revealed a sufficiently rapid uptake to cause the metabolic effects observed. These data are consistent with the notion that extramitochondrial Ca(2+) can modify ATP production, via an increase in matrix Ca(2+) content, rapidly enough to support cardiac work transitions in vivo.     

The Effect of Linkage Disequilibrium and Family Relationships on the Reliability of Genomic Prediction

Although the concept of genomic selection relies on linkage disequilibrium (LD) between quantitative trait loci and markers, reliability of genomic predictions is strongly influenced by family relationships. In this study, we investigated the effects of LD and family relationships on reliability of genomic predictions and the potential of deterministic formulas to predict reliability using population parameters in populations with complex family structures. Five groups of selection candidates were simulated by taking different information sources from the reference population into account: (1) allele frequencies, (2) LD pattern, (3) haplotypes, (4) haploid chromosomes, and (5) individuals from the reference population, thereby having real family relationships with reference individuals. Reliabilities were predicted using genomic relationships among 529 reference individuals and their relationships with selection candidates and with a deterministic formula where the number of effective chromosome segments (M_e) was estimated based on genomic and additive relationship matrices for each scenario. At a heritability of 0.6, reliabilities based on genomic relationships were 0.002 ± 0.0001 (allele frequencies), 0.022 ± 0.001 (LD pattern), 0.018 ± 0.001 (haplotypes), 0.100 ± 0.008 (haploid chromosomes), and 0.318 ± 0.077 (family relationships). At a heritability of 0.1, relative differences among groups were similar. For all scenarios, reliabilities were similar to predictions with a deterministic formula using estimated M_e. So, reliabilities can be predicted accurately using empirically estimated M_e and level of relationship with reference individuals has a much higher effect on the reliability than linkage disequilibrium per se. Furthermore, accumulated length of shared haplotypes is more important in determining the reliability of genomic prediction than the individual shared haplotype length.