Identifying QTLs involved in hybrid performance and heterotic group complementarity: new GWAS models applied to factorial and admixed diallel maize hybrid panels.

KEY MESSAGE: An original GWAS model integrating the ancestry of alleles was proposed and allowed the detection of background specific additive and dominance QTLs involved in heterotic group complementarity and hybrid performance. Maize genetic diversity is structured into genetic groups selected and improved relative to each other. This process increases group complementarity and differentiation over time and ensures that the hybrids produced from inter-group crosses exhibit high performances and heterosis. To identify loci involved in hybrid performance and heterotic group complementarity, we introduced an original association study model that disentangles allelic effects from the heterotic group origin of the alleles and compared it with a conventional additive/dominance model. This new model was applied on a factorial between Dent and Flint lines and a diallel between Dent-Flint admixed lines with two different layers of analysis: within each environment and in a multiple-environment context. We identified several strong additive QTLs for all traits, including some well-known additive QTLs for flowering time (in the region of Vgt1/2 on chromosome 8). Yield trait displayed significant non-additive effects in the diallel panel. Most of the detected Yield QTLs exhibited overdominance or, more likely, pseudo-overdominance effects. Apparent overdominance at these QTLs contributed to a part of the genetic group complementarity. The comparison between environments revealed a higher stability of additive QTL effects than non-additive ones. Several QTLs showed variations of effects according to the local heterotic group origin. We also revealed large chromosomic regions that display genetic group origin effects. Altogether, our results illustrate how admixed panels combined with dedicated GWAS modeling allow the identification of new QTLs that could not be revealed by a classical hybrid panel analyzed with traditional modeling.

Improving genomic predictions with inbreeding and nonadditive effects in two admixed maize hybrid populations in single and multienvironment contexts.

Genetic admixture, resulting from the recombination between structural groups, is frequently encountered in breeding populations. In hybrid breeding, crossing admixed lines can generate substantial nonadditive genetic variance and contrasted levels of inbreeding which can impact trait variation. This study aimed at testing recent methodological developments for the modeling of inbreeding and nonadditive effects in order to increase prediction accuracy in admixed populations. Using two maize (Zea mays L.) populations of hybrids admixed between dent and flint heterotic groups, we compared a suite of five genomic prediction models incorporating (or not) parameters accounting for inbreeding and nonadditive effects with the natural and orthogonal interaction approach in single and multienvironment contexts. In both populations, variance decompositions showed the strong impact of inbreeding on plant yield, height, and flowering time which was supported by the superiority of prediction models incorporating this effect (+0.038 in predictive ability for mean yield). In most cases dominance variance was reduced when inbreeding was accounted for. The model including additivity, dominance, epistasis, and inbreeding effects appeared to be the most robust for prediction across traits and populations (+0.054 in predictive ability for mean yield). In a multienvironment context, we found that the inclusion of nonadditive and inbreeding effects was advantageous when predicting hybrids not yet observed in any environment. Overall, comparing variance decompositions was helpful to guide model selection for genomic prediction. Finally, we recommend the use of models including inbreeding and nonadditive parameters following the natural and orthogonal interaction approach to increase prediction accuracy in admixed populations.

The Pitfalls of Heterosis Coefficients.

Heterosis (hybrid vigour) is a universal phenomenon of crucial agro-economic and evolutionary importance. We show that the most common heterosis coefficients do not properly measure deviation from additivity because they include both a component accounting for "real" heterosis and a term that is not related to heterosis, since it is derived solely from parental values. Therefore, these coefficients are inadequate whenever the aim of the study is to compare heterosis levels between different traits, environments, genetic backgrounds, or developmental stages, as these factors may affect not only the level of non-additivity, but also parental values. The only relevant coefficient for such comparisons is the so-called "potence ratio". Because most heterosis studies consider several traits/stages/environmental conditions, our observations support the use of the potence ratio, at least in non-agronomic contexts, because it is the only non-ambiguous heterosis coefficient.

Heterosis Is a Systemic Property Emerging From Non-linear Genotype-Phenotype Relationships: Evidence From in Vitro Genetics and Computer Simulations.

Heterosis, the superiority of hybrids over their parents for quantitative traits, represents a crucial issue in plant and animal breeding as well as evolutionary biology. Heterosis has given rise to countless genetic, genomic and molecular studies, but has rarely been investigated from the point of view of systems biology. We hypothesized that heterosis is an emergent property of living systems resulting from frequent concave relationships between genotypic variables and phenotypes, or between different phenotypic levels. We chose the enzyme-flux relationship as a model of the concave genotype-phenotype (GP) relationship, and showed that heterosis can be easily created in the laboratory. First, we reconstituted in vitro the upper part of glycolysis. We simulated genetic variability of enzyme activity by varying enzyme concentrations in test tubes. Mixing the content of "parental" tubes resulted in "hybrids," whose fluxes were compared to the parental fluxes. Frequent heterotic fluxes were observed, under conditions that were determined analytically and confirmed by computer simulation. Second, to test this model in a more realistic situation, we modeled the glycolysis/fermentation network in yeast by considering one input flux, glucose, and two output fluxes, glycerol and acetaldehyde. We simulated genetic variability by randomly drawing parental enzyme concentrations under various conditions, and computed the parental and hybrid fluxes using a system of differential equations. Again we found that a majority of hybrids exhibited positive heterosis for metabolic fluxes. Cases of negative heterosis were due to local convexity between certain enzyme concentrations and fluxes. In both approaches, heterosis was maximized when the parents were phenotypically close and when the distributions of parental enzyme concentrations were contrasted and constrained. These conclusions are not restricted to metabolic systems: they only depend on the concavity of the GP relationship, which is commonly observed at various levels of the phenotypic hierarchy, and could account for the pervasiveness of heterosis.

Influence of water temperature and food on the last stages of cultured pearl mineralization from the black-lip pearl oyster Pinctada margaritifera.

Environmental parameters, such as food level and water temperature, have been shown to be major factors influencing pearl oyster shell growth and molecular mechanisms involved in this biomineralization process. The present study investigates the effect of food level (i.e., microalgal concentration) and water temperature, in laboratory controlled conditions, on the last stages of pearl mineralization in order to assess their impact on pearl quality. To this end, grafted pearl oysters were fed at different levels of food and subjected to different water temperatures one month prior to harvest to evaluate the effect of these factors on 1) pearl and shell deposition rate, 2) expression of genes involved in biomineralization in pearl sacs, 3) nacre ultrastructure (tablet thickness and number of tablets deposited per day) and 4) pearl quality traits. Our results revealed that high water temperature stimulates both shell and pearl deposition rates. However, low water temperature led to thinner nacre tablets, a lower number of tablets deposited per day and impacted pearl quality with better luster and fewer defects. Conversely, the two tested food level had no significant effects on shell and pearl growth, pearl nacre ultrastructure or pearl quality. However, one gene, Aspein, was significantly downregulated in high food levels. These results will be helpful for the pearl industry. A wise strategy to increase pearl quality would be to rear pearl oysters at a high water temperature to increase pearl growth and consequently pearl size; and to harvest pearls after a period of low water temperature to enhance luster and to reduce the number of defects.

Uranium perturbs signaling and iron uptake response in Arabidopsis thaliana roots.

Uranium is a natural element which is mainly redistributed in the environment due to human activity, including accidents and spillages. Plants may be useful in cleaning up after incidents, although little is yet known about the relationship between metal speciation and plant response. Here, J-Chess modeling was used to predict U speciation and exposure conditions affecting U bioavailability for plants. The model was confirmed by exposing Arabidopsis thaliana plants to U under hydroponic conditions. The early root response was characterized using complete Arabidopsis transcriptome microarrays (CATMA). Expression of 111 genes was modified at the three timepoints studied. The associated biological processes were further examined by real-time quantitative RT-PCR. Annotation revealed that oxidative stress, cell wall and hormone biosynthesis, and signaling pathways (including phosphate signaling) were affected by U exposure. The main actors in iron uptake and signaling (IRT1, FRO2, AHA2, AHA7 and FIT1) were strongly down-regulated upon exposure to uranyl. A network calculated using IRT1, FRO2 and FIT1 as bait revealed a set of genes whose expression levels change under U stress. Hypotheses are presented to explain how U perturbs the iron uptake and signaling response. These results give preliminary insights into the pathways affected by uranyl uptake, which will be of interest for engineering plants to help clean areas contaminated with U.

Functional divergence in shrimp anti-lipopolysaccharide factors (ALFs): from recognition of cell wall components to antimicrobial activity.

Antilipopolysaccharide factors (ALFs) have been described as highly cationic polypeptides with a broad spectrum of potent antimicrobial activities. In addition, ALFs have been shown to recognize LPS, a major component of the Gram-negative bacteria cell wall, through conserved amino acid residues exposed in the four-stranded ß-sheet of their three dimensional structure. In penaeid shrimp, ALFs form a diverse family of antimicrobial peptides composed by three main variants, classified as ALF Groups A to C. Here, we identified a novel group of ALFs in shrimp (Group D ALFs), which corresponds to anionic polypeptides in which many residues of the LPS binding site are lacking. Both Group B (cationic) and Group D (anionic) shrimp ALFs were produced in a heterologous expression system. Group D ALFs were found to have impaired LPS-binding activities and only limited antimicrobial activity compared to Group B ALFs. Interestingly, all four ALF groups were shown to be simultaneously expressed in an individual shrimp and to follow different patterns of gene expression in response to a microbial infection. Group B was by far the more expressed of the ALF genes. From our results, nucleotide sequence variations in shrimp ALFs result in functional divergence, with significant differences in LPS-binding and antimicrobial activities. To our knowledge, this is the first functional characterization of the sequence diversity found in the ALF family.

Big defensins, a diverse family of antimicrobial peptides that follows different patterns of expression in hemocytes of the oyster Crassostrea gigas.

BACKGROUND: Big defensin is an antimicrobial peptide composed of a highly hydrophobic N-terminal region and a cationic C-terminal region containing six cysteine residues involved in three internal disulfide bridges. While big defensin sequences have been reported in various mollusk species, few studies have been devoted to their sequence diversity, gene organization and their expression in response to microbial infections. FINDINGS: Using the high-throughput Digital Gene Expression approach, we have identified in Crassostrea gigas oysters several sequences coding for big defensins induced in response to a Vibrio infection. We showed that the oyster big defensin family is composed of three members (named Cg-BigDef1, Cg-BigDef2 and Cg-BigDef3) that are encoded by distinct genomic sequences. All Cg-BigDefs contain a hydrophobic N-terminal domain and a cationic C-terminal domain that resembles vertebrate ß-defensins. Both domains are encoded by separate exons. We found that big defensins form a group predominantly present in mollusks and closer to vertebrate defensins than to invertebrate and fungi CSαß-containing defensins. Moreover, we showed that Cg-BigDefs are expressed in oyster hemocytes only and follow different patterns of gene expression. While Cg-BigDef3 is non-regulated, both Cg-BigDef1 and Cg-BigDef2 transcripts are strongly induced in response to bacterial challenge. Induction was dependent on pathogen associated molecular patterns but not damage-dependent. The inducibility of Cg-BigDef1 was confirmed by HPLC and mass spectrometry, since ions with a molecular mass compatible with mature Cg-BigDef1 (10.7 kDa) were present in immune-challenged oysters only. From our biochemical data, native Cg-BigDef1 would result from the elimination of a prepropeptide sequence and the cyclization of the resulting N-terminal glutamine residue into a pyroglutamic acid. CONCLUSIONS: We provide here the first report showing that big defensins form a family of antimicrobial peptides diverse not only in terms of sequences but also in terms of genomic organization and regulation of gene expression.

Systemic properties of metabolic networks lead to an epistasis-based model for heterosis.

The genetic and molecular approaches to heterosis usually do not rely on any model of the genotype-phenotype relationship. From the generalization of Kacser and Burns' biochemical model for dominance and epistasis to networks with several variable enzymes, we hypothesized that metabolic heterosis could be observed because the response of the flux towards enzyme activities and/or concentrations follows a multi-dimensional hyperbolic-like relationship. To corroborate this, we used the values of systemic parameters accounting for the kinetic behaviour of four enzymes of the upstream part of glycolysis, and simulated genetic variability by varying in silico enzyme concentrations. Then we "crossed" virtual parents to get 1,000 hybrids, and showed that best-parent heterosis was frequently observed. The decomposition of the flux value into genetic effects, with the help of a novel multilocus epistasis index, revealed that antagonistic additive-by-additive epistasis effects play the major role in this framework of the genotype-phenotype relationship. This result is consistent with various observations in quantitative and evolutionary genetics, and provides a model unifying the genetic effects underlying heterosis.

A relationship between antimicrobial peptide gene expression and capacity of a selected shrimp line to survive a Vibrio infection.

Understanding of antimicrobial defence mechanisms of penaeid shrimp should help in the design of efficient strategies for the management and disease control in aquaculture. In this study, we have specifically analysed the expression in circulating hemocytes of antimicrobial peptides (AMPs) encoding genes, such as PEN2 and PEN3, ALF, crustin, lysozyme and a putative cysteine-rich peptide. We evidenced a relationship between the level of expression of some AMPs and the successful response of the shrimp, Litopenaeus stylirostris, to circumvent a pathogenic Vibrio penaeicida infection. Additionally, significant differences in some AMP transcript amounts are evidenced between control, non-selected shrimp line and the third generation breeding of shrimp selected for their survival to natural V. penaeicida infections. On the basis of these results, it will now be of great interest to determine if these AMPs are directly involved in the resistance of shrimp to infection or if they only reflect other acquired defence mechanisms which can confer a resistance.

Evidence of a bactericidal permeability increasing protein in an invertebrate, the Crassostrea gigas Cg-BPI.

A cDNA sequence with homologies to members of the LPS-binding protein and bactericidal/permeability-increasing protein (BPI) family was identified in the oyster Crassostrea gigas. The recombinant protein was found to bind LPS, to display bactericidal activity against Escherichia coli, and to increase the permeability of the bacterial cytoplasmic membrane. This indicated that it is a BPI rather than an LPS-binding protein. By in situ hybridization, the expression of the C. gigas BPI (Cg-bpi) was found to be induced in hemocytes after oyster bacterial challenge and to be constitutive in various epithelia of unchallenged oysters. Thus, Cg-bpi transcripts were detected in the epithelial cells of tissues/organs in contact with the external environment (mantle, gills, digestive tract, digestive gland diverticula, and gonad follicles). Therefore, Cg-BPI, whose expression profile and biological properties are reminiscent of mammalian BPIs, may provide a first line of defense against potential bacterial invasion. To our knowledge, this is the first characterization of a BPI in an invertebrate.

Localization of deformed wing virus infection in queen and drone Apis mellifera L.

The distribution of deformed wing virus infection within the honey bee reproductive castes (queens, drones) was investigated by in situ hybridization and immunohistology from paraffin embedded sections. Digoxygenin or CY5.5 fluorochrome end-labelled nucleotide probes hybridizing to the 3' portion of the DWV genome were used to identify DWV RNA, while a monospecific antibody to the DWV-VP1 structural protein was used to identify viral proteins and particles. The histological data were confirmed by quantitative RT-PCR of dissected organs. Results showed that DWV infection is not restricted to the digestive tract of the bee but spread in the whole body, including queen ovaries, queen fat body and drone seminal vesicles.

The early responses of Arabidopsis thaliana cells to cadmium exposure explored by protein and metabolite profiling analyses.

To get more insight into plant cell response to cadmium (Cd) stress, both proteomic and metabolomic "differential display" analyses were performed on Arabidopsis thaliana cells exposed to different concentrations of the toxic chemical. After a 24 h treatment, soluble proteins extracted from untreated and treated cells were separated by 2-D-PAGE and image analyses were performed to quantify and compare protein levels. Proteins up- and down-regulated in response to Cd were identified by MS and mapped into specific metabolic pathways and cellular processes, highlighting probable activation of the carbon, nitrogen, and sulfur metabolic pathways. For some of these proteins, Northern blot and RT-PCR analyses were performed to test transcript accumulation in response to Cd. In parallel, metabolite profiling analyses by LC coupled to ESI MS were initiated to better characterize the metabolic adaptation to the chemical stress. This study revealed that the main variation at the metabolite level came from the presence of six different families of phytochelatins, in A. thaliana cells treated with Cd, whose accumulation increases with Cd concentrations. Taken together these data provide an overview of the molecular and cellular changes elicited by Cd exposure.

Simplified modelling of metabolic pathways for flux prediction and optimization: lessons from an in vitro reconstruction of the upper part of glycolysis.

Explicit modelling of metabolic networks relies on well-known mathematical tools and specialized computer programs. However, identifying and estimating the values of the very numerous enzyme parameters inherent to the models remain a tedious and difficult task, and the rate equations of the reactions are usually not known in sufficient detail. A way to circumvent this problem is to use 'non-mechanistic' models, which may account for the behaviour of the systems with a limited number of parameters. Working on the first part of glycolysis reconstituted in vitro, we showed how to derive, from titration experiments, values of effective enzyme activity parameters that do not include explicitly any of the classical kinetic constants. With a maximum of only two parameters per enzyme, this approach produced very good estimates for the flux values, and enabled us to determine the optimization conditions of the system, i.e. to calculate the set of enzyme concentrations that maximizes the flux. This fast and easy method should be valuable in the context of integrative biology or for metabolic engineering, where the challenge is to deal with the dramatic increase in the number of parameters when the systems become complex.

Characterization of a defensin from the oyster Crassostrea gigas. Recombinant production, folding, solution structure, antimicrobial activities, and gene expression.

In invertebrates, defensins were found in arthropods and in the mussels. Here, we report for the first time the identification and characterization of a defensin (Cg-Def) from an oyster. Cg-def mRNA was isolated from Crassostrea gigas mantle using an expressed sequence tag approach. To gain insight into potential roles of Cg-Def in oyster immunity, we produced the recombinant peptide in Escherichia coli, characterized its antimicrobial activities, determined its solution structure by NMR spectroscopy, and quantified its gene expression in vivo following bacterial challenge of oysters. Recombinant Cg-Def was active in vitro against Gram-positive bacteria but showed no or limited activities against Gram-negative bacteria and fungi. The activity of Cg-Def was retained in vitro at a salt concentration similar to that of seawater. The Cg-Def structure shares the so-called cystine-stabilized alpha-beta motif (CS-alphabeta) with arthropod defensins but is characterized by the presence of an additional disulfide bond, as previously observed in the mussel defensin (MGD-1). Nevertheless, despite a similar global fold, the Cg-Def and MGD-1 structures mainly differ by the size of their loops and by the presence of two aspartic residues in Cg-Def. Distribution of Cg-def mRNA in various oyster tissues revealed that Cg-def is mainly expressed in mantle edge where it was detected by mass spectrometry analyses. Furthermore, we observed that the Cg-def messenger concentration was unchanged after bacterial challenge. Our results suggest that Cg-def gene is continuously expressed in the mantle and would play a key role in oyster by providing a first line of defense against pathogen colonization.

Evidence in oyster of a plasma extracellular superoxide dismutase which binds LPS.

We have characterized in the oyster Crassostrea gigas an extracellular superoxide dismutase (Cg-EcSOD) which appears to bind lipopolysaccharides (LPS). The protein has been purified from the oyster plasma and identified as a Cu/ZnSOD according to its N-terminal sequencing and biological activity. Cg-EcSOD expression and synthesis are restricted to hemocytes as revealed by in situ hybridization and immunocytochemistry. Cg-EcSOD-expressing hemocytes were seen in blood circulation, in connective tissues, and closely associated to endothelium blood vessels. Cg-EcSOD presents in its amino acid sequence a LPS-binding motif found in the endotoxin receptor CD14 and we show that the protein displays an affinity to Escherichia coli bacteria and with LPS and Lipid A. Additionally, an RGD motif known to be implicated in the association to membrane integrin receptor is present in the amino acid sequence. The purified Cg-EcSOD was shown to bind to oyster hemocytes and to be immunocolocalized with a beta-integrin-like receptor.

An extension to the metabolic control theory taking into account correlations between enzyme concentrations.

The classical metabolic control theory [Kacser, H. & Burns, J.A. (1973) Symp. Soc. Exp. Biol.27, 65-104; Heinrich, R. & Rapoport, T. (1974) Eur. J. Biochem.42, 89-95.] does not take into account experimental evidence for correlations between enzyme concentrations in the cell. We investigated the implications of two causes of linear correlations: competition between enzymes, which is a mere physical adaptation of the cell to the limitation of resources and space, and regulatory correlations, which result from the existence of regulatory networks. These correlations generate redistribution of enzyme concentrations when the concentration of an enzyme varies; this may dramatically alter the flux and metabolite concentration curves. In particular, negative correlations cause the flux to have a maximum value for a defined distribution of enzyme concentrations. Redistribution coefficients of enzyme concentrations allowed us to calculate the 'combined response coefficient' that quantifies the response of flux or metabolite concentration to a perturbation of enzyme concentration.

Assessing factors for reliable quantitative proteomics based on two-dimensional gel electrophoresis.

We statistically analysed various factors to get accurate estimates of protein quantities from two-dimensional gels. Yeast proteins were labelled with (35)S or stained with Coomassie Brilliant Blue G-250, and spots were automatically quantified with software packages Kepler, ImageQuaNT, Melanie 3.0 and Progenesis. The different software packages proved to have very similar performances. With (35)S-labelled actin spot as a reference, we studied the staining efficiency of colloidal Coomassie blue as a function of amino acid composition of the protein, and derived an equation to estimate the number of molecules per cell from blue-stained proteins. Absolute quantification of most glycolytic enzymes was carried out in two yeast strains.

A new class (penaeidin class 4) of antimicrobial peptides from the Atlantic white shrimp (Litopenaeus setiferus) exhibits target specificity and an independent proline-rich-domain function.

A highly pure, chemically defined representative of a new class of antimicrobial peptide from the Atlantic white shrimp (Litopenaeus setiferus), penaeidin class 4 [Pen4-1 (penaeidin class 4 isoform 1)], was produced synthetically. Chemical synthesis was achieved by native ligation from two separate domains yielding a bioactive peptide that reflected the characteristics of native penaeidin. Synthetic Pen4-1 proved to be an effective antimicrobial peptide, particularly against the broad-spectrum pathogen Fusarium oxysporum, exhibiting a complex effect on reproductive growth at inhibitory concentrations resulting in the suppression of spore formation. Pen4-1 exhibits unique features [not previously observed for penaeidins from the Pacific white shrimp (L. vannamei)], including target-species specificity against Gram-positive bacteria, indicating a potential partitioning of antimicrobial function among this family of peptides. The proline-rich domain of penaeidin class 4 alone was an active antimicrobial peptide, having the same target range as the full-length Pen4-1. These findings indicate that the proline-rich domain of penaeidin is sufficient to confer target specificity and that divergence in this domain between classes can result in a gain in antimicrobial function as observed for the proline-rich domain of Pen4-1.