Integration of a Systems Biological Network Analysis and QTL Results for Biomass Heterosis in Arabidopsis thaliana

To contribute to a further insight into heterosis we applied an integrative analysis to a systems biological network approach and a quantitative genetics analysis towards biomass heterosis in early Arabidopsis thaliana development. The study was performed on the parental accessions C24 and Col-0 and the reciprocal crosses. In an over-representation analysis it was tested if the overlap between the resulting gene lists of the two approaches is significantly larger than expected by chance. Top ranked genes in the results list of the systems biological analysis were significantly over-represented in the heterotic QTL candidate regions for either hybrid as well as regarding mid-parent and best-parent heterosis. This suggests that not only a few but rather several genes that influence biomass heterosis are located within each heterotic QTL region. Furthermore, the overlapping resulting genes of the two integrated approaches were particularly enriched in biomass related pathways. A chromosome-wise over-representation analysis gave rise to the hypothesis that chromosomes number 2 and 4 probably carry a majority of the genes involved in biomass heterosis in the early development of Arabidopsis thaliana.     

Heterosis in crosses between geographically separated populations of Drosophila melanogaster

Summary An experiment was performed to test the hypothesis that the genetic distance between populations estimated from enzyme loci could be used to predict the amount of heterosis that would occur in crosses between these populations. A partial diallel cross using 11 populations of Drosophila melanogaster from the AustralianPacific region and from England was carried out. Heterosis for larval viability, fecundity, cold shock mortality, and an index of these three traits was recorded. When two populations originating from the same location were crossed, no heterosis occurred, but otherwise heterosis was significant for all traits. For larval viability, a similar low level of heterosis occurred in all crosses. For cold shock mortality, the level of heterosis varied widely and fecundity showed a pattern intermediate between these two. The geographic distance between the sites from which populations originated was not correlated with the amount of heterosis in their crosses. There was a tendency for populations from ecologically different environments to show heterosis in crosses. Genetic distance based on ten enzyme loci was correlated with heterosis for cold shock mortality and the combined trait index. These results can be explained by the hypothesis that genes affecting larval viability are subject to strong, uniform selection in all populations, which limits the extent to which gene frequencies can drift apart. However, genes affecting cold shock mortality and the enzyme loci are subject to different selection pressures in different environments. This divergent selection combined with genetic drift causes divergence in gene frequency and heterosis.     

Cell number counts - The fw2.2 and CNR genes and implications for controlling plant fruit and organ size

Two key determinants of plant and organ size are cell number and cell size, and altering either one may affect the plant organ size, but cell number control often plays a predominant role in natural populations. Domesticated crops usually have larger fruit and harvested organ sizes than wild progenitors. Crop yields have increased significantly by breeding, often via heterosis, which is associated with increased plant and organ size primarily achieved by cell number increases. A small class of genes is now known that control plant and organ sizes though cell number or cell size. The fw2.2 gene was found to control a major QTL for tomato fruit size by negatively affecting cell numbers. Orthologs to these fw2.2 genes underlie QTLs for fruit sizes in other species, and their expression can be negatively correlated with increased cell number. In maize decreased or increased expression of the fw2.2 ortholog ZmCNR1, increases or decreases cell number, respectively, thereby affecting maize organ size throughout the plant and thus also whole plant size. Therefore, these genes should now be considered as more general regulators of plant cell number and organ size. The exact molecular function of these transmembrane domain proteins remains unknown, as does any clear relationship to the cell cycle. Because these genes control organ sizes in diverse plants and important crop species, and because they can affect whole plant size, interest arose into how effects of such genes could parallel agronomic crop improvements, in particular that by heterosis, as it also affects cell number. In joining these subjects here in discussion we speculate on how single gene cell number regulation and heterosis may cooperate in crop improvement.     

Targeted association analysis identified japonica rice varieties achieving Na+/K+ homeostasis without the allelic make-up of the salt tolerant indica variety Nona Bokra

During the last decade, a large number of QTLs and candidate genes for rice tolerance to salinity have been reported. Using 124 SNP and 52 SSR markers, we targeted 14 QTLs and 65 candidate genes for association mapping within the European Rice Core collection (ERCC) comprising 180 japonica accessions. Significant differences in phenotypic response to salinity were observed. Nineteen distinct loci significantly associated with one or more phenotypic response traits were detected. Linkage disequilibrium between these loci was extremely low, indicating a random distribution of favourable alleles in the ERCC. Analysis of the function of these loci indicated that all major tolerance mechanisms were present in the ERCC although the useful level of expression of the different mechanisms was scattered among different accessions. Under moderate salinity stress some accessions achieved the same level of control of Na(+) concentration and Na(+)/K(+) equilibrium as the indica reference variety for salinity tolerance Nona Bokra, although without sharing the same alleles at several loci associated with Na(+) concentration. This suggests (a) differences between indica and japonica subspecies in the effect of QTLs and genes involved in salinity tolerance and (b) further potential for the improvement of tolerance to salinity above the tolerance level of Nona Bokra, provided the underlying mechanisms are complementary at the whole plant level. No accession carried all favourable alleles, or showed the best phenotypic responses for all traits measured. At least nine accessions were needed to assemble the favourable alleles and all the best phenotypic responses. An effective strategy for the accumulation of the favourable alleles would be marker-assisted population improvement.     

Prediction of heterosis using genome-wide SNP-marker data: application to egg production traits in white Leghorn crosses

Prediction of heterosis has a long history with mixed success, partly due to low numbers of genetic markers and/or small data sets. We investigated the prediction of heterosis for egg number, egg weight and survival days in domestic white Leghorns, using ∼400 000 individuals from 47 crosses and allele frequencies on ∼53 000 genome-wide single nucleotide polymorphisms (SNPs). When heterosis is due to dominance, and dominance effects are independent of allele frequencies, heterosis is proportional to the squared difference in allele frequency (SDAF) between parental pure lines (not necessarily homozygous). Under these assumptions, a linear model including regression on SDAF partitions crossbred phenotypes into pure-line values and heterosis, even without pure-line phenotypes. We therefore used models where phenotypes of crossbreds were regressed on the SDAF between parental lines. Accuracy of prediction was determined using leave-one-out cross-validation. SDAF predicted heterosis for egg number and weight with an accuracy of ∼0.5, but did not predict heterosis for survival days. Heterosis predictions allowed preselection of pure lines before field-testing, saving ∼50% of field-testing cost with only 4% loss in heterosis. Accuracies from cross-validation were lower than from the model-fit, suggesting that accuracies previously reported in literature are overestimated. Cross-validation also indicated that dominance cannot fully explain heterosis. Nevertheless, the dominance model had considerable accuracy, clearly greater than that of a general/specific combining ability model. This work also showed that heterosis can be modelled even when pure-line phenotypes are unavailable. We concluded that SDAF is a useful predictor of heterosis in commercial layer breeding.     

Enriched partial correlations in genome-wide gene expression profiles of hybrids (A. thaliana): a systems biological approach towards the molecular basis of heterosis

Heterosis is a well-known phenomenon but the underlying molecular mechanisms are not yet established. To contribute to the understanding of heterosis at the molecular level, we analyzed genome-wide gene expression profile data of Arabidopsis thaliana in a systems biological approach. We used partial correlations to estimate the global interaction structure of regulatory networks. Our hypothesis states that heterosis comes with an increased number of partial correlations which we interpret as increased numbers of regulatory interactions leading to enlarged adaptability of the hybrids. This hypothesis is true for mid-parent heterosis for our dataset of gene expression in two homozygous parental lines and their reciprocal crosses. For the case of best-parent heterosis just one hybrid is significant regarding our hypothesis based on a resampling analysis. Summarizing, both metabolome and gene expression level of our illustrative dataset support our proposal of a systems biological approach towards a molecular basis of heterosis.     

Whole-exome sequencing reveals insights into genetic susceptibility to Congenital Zika Syndrome

Congenital Zika Syndrome (CZS) is a critical illness with a wide range of severity caused by Zika virus (ZIKV) infection during pregnancy. Life-threatening neurodevelopmental dysfunctions are among the most common phenotypes observed in affected newborns. Risk factors that contribute to susceptibility and response to ZIKV infection may be related to the virus itself, the environment, and maternal genetic background. Nevertheless, the newborn’s genetic contribution to the critical illness is still not elucidated. Here, we aimed to identify possible genetic variants as well as relevant biological pathways that might be associated with CZS phenotypes. For this purpose, we performed a whole-exome sequencing in 40 children born to women with confirmed exposure to ZIKV during pregnancy. We investigated the occurrence of rare harmful single-nucleotide variants (SNVs) possibly associated with inborn errors in genes ontologically related to CZS phenotypes. Moreover, an exome-wide association analysis was also performed using a case-control design (29 CZS cases and 11 controls), for both common and rare variants. Five out of the 29 CZS patients harbored known pathogenic variants likely to contribute to mild to severe manifestations observed. Approximately, 30% of affected individuals carried at least one pathogenic or likely pathogenic SNV in genes candidates to play a role in CZS. Our common variant association analysis detected a suggestive protective effect of the rs2076469 in DISP3 gene (p-value: 1.39 x 10⁻⁵). The IL12RB2 gene (p-value: 2.18x10^-11) also showed an unusual distribution of nonsynonymous rare SNVs in control samples. Finally, genes harboring harmful variants are involved in processes related to CZS phenotypes such as neurological development and immunity. Therefore, both rare and common variations may be likely to contribute as the underlying genetic cause of CZS susceptibility. The variations and pathways identified in this study may also have implications for the development of therapeutic strategies in the future.     

A systematic dissection in oilseed rape provides insights into the genetic architecture and molecular mechanism of yield heterosis

Heterosis refers to the better performance of cross progeny compared with inbred parents, and its utilization contributes greatly to agricultural production. Several hypotheses have been proposed to explain heterosis mainly including dominance, over-dominance (or pseudo-overdominance) and epistasis. However, systematic dissection and verification of these hypotheses are rarely documented. Here, comparison of heterosis level across different traits showed that the strong heterosis of composite traits (such as yield) could be attributed to the multiplicative effects of moderate heterosis of component traits, whether at the genome or locus level. Yield heterosis was regulated by a complex trait-QTL network that was characterized by obvious centre-periphery structure, hub QTL, complex up/downstream and positive/negative feedback relationships. More importantly, we showed that better-parent heterosis on yield could be produced in a cross of two near-isogenic lines by the pyramiding and complementation of two major heterotic QTL showing partial-dominance on yield components. The causal gene (BnaA9.CYP78A9) of QC14 was identified, and its heterotic effect results from the heterozygous status of a CACTA-like transposable element in its upstream regulatory region, which led to partial dominance at expression and auxin levels, thus resulting in non-additive expression of downstream responsive genes involved in cell cycle and proliferation, eventually leading to the heterosis of cell number. Taken together, the results at the phenotypic, genetic and molecular levels were highly consistent, which demonstrated that the pyramiding effect of heterotic QTL and the multiplicative effect of individual component traits could well explain substantial parts of yield heterosis in oilseed rape. These results provide in-depth insights into the genetic architecture and molecular mechanism of yield heterosis.     

On the genetic control of heterosis for fruit shape in melon (Cucumis melo L.)

The objective of the present work is to study the genetic basis of heterosis for fruit shape (FS) in melon observed in a cross between the Spanish cultivar "Piel de Sapo" (PS) and the Korean accession PI 161375 (Songwang Charmi [SC]) using a set of near-isogenic lines (NILs) with contrasting phenotypes for FS, each carrying a single chromosomal introgression from SC within the genetic background of PS. We investigated the FS of homozygous NILs, hybrids NIL x PS, and all 2-way crosses between NILs to test the main heterosis hypotheses (dominance, overdominance, and epistatic interactions). Gene action of alleles of quantitative trait loci inducing fruit enlargement was dominance, whereas those inducing rounder fruit were additive or recessive. Only minor epistatic interactions were found. Therefore, the most plausible explanation for FS heterosis in this cross is in agreement with the dominance complementation hypothesis. Over 70% of the hybrid heterosis could be achieved by combining just 2 loci, indicating that the genetic control of FS heterosis in this cross is relatively simple. FS is proposed as a reproductive trait in melon because of the high correlation to the number of seeds produced along the fruit longitudinal axis.     

Partitioning of genetic effects on lifetime performance of mice

Summary A total of 2,457 lifetime performance records of 29 genetic groups of mice was analyzed using multiple regression of records on the proportion of gene contribution from 6 lines (designated as Lines M_P, m_Q, W_P, w_Q, C_P and c_Q). Genetic effects were partitioned into line additive, line maternal, direct heterosis, maternal heterosis and paternal heterosis effects. The line additive and line maternal effects were expressed as deviations from Line c_Q. Seventeen of 25 line additive effects differed significantly (P<0.05) from Line c_Q whereas only 4 of 25 line maternal effects deviated significantly from Line c_Q. Deviations in line additive effects from c_Q were negative in all lines examined whereas deviations in line maternal effects from c_Q were all positive, indicating a negative relationship between line additive and line maternal effects. Direct heterosis effects were all positive and significant (P < 0.01) except in the M_PxW_P cross which was produced by mating Lines M_P and W_P of the same base population (P). Maternal heterosis effects were significant in 10 of 20 cases whereas paternal heterosis effects were significant in 13 of 20 cases. Although direct heterosis is a major component of total heterosis effects (sum of direct, maternal and paternal heterosis), the results suggest that parental heterosis may need to be considered in producing multiple way crosses. The fitting of line additive, line maternal, direct heterosis, maternal heterosis and paternal heterosis effects in the multiple regression model effectively accounted for all genetic effects in lifetime performance.Animal Research Centre Contribution No. 1326     

Use of a Controlled-Nutrient Experiment to Test Heterosis Hypotheses

A controlled-nutrient (CN) experiment was conducted to test three heterosis hypotheses with reference to tomato yield, and its components, for a set of two inbred lines and their hybrid that had previously exhibited considerable heterosis under field conditions. The CN treatments consisted of periodic applications of differential doses of nutrient solution to plants reared individually in containers filled with vermiculite. Ripe fruit were harvested, counted and weighed over a period of 340 days. The data permitted the partitioning of yield into a closed system of five component variables. Heterosis was not exhibited by yield, nor yield components, at any of the four nutrient levels. Hence the total heterosis phenomenon was classified as nutrient-dependent: heterosis occurring under field conditions, but not under the nutritional restrictions of the CN experiment. Three heterosis hypotheses were examined for their ability to explain all of the nutrient-dependent aspects of the heterosis phenomenon. Hypothesis 1: Heterosis is a consequence of a more efficient hybrid metabolic system in that it can produce more product with equal input. Hypothesis 2: Heterosis is a consequence of the somatic multiplication of additive component traits. Hypothesis 3: Heterosis is a consequence of a faster hybrid growth rate. Although none of the hypotheses are rejected by the field data, the first two are rejected by the CN experimental results. The third hypothesis fits all aspects of the nutrient-dependent heterosis phenomenon remarkably well. It is speculated that the indeterminate pattern of plant development responsible for yield and its components is due to two major gene systems: genes that determine morphogenetic, and genes that determine growth rate manifestations of growth. Under this hypothesis, the CN technique permits separation of the responses due to these two gene systems.     

Identification and characterization of a repertoire of genes differentially expressed in developing top ear shoots between a superior hybrid and its parental inbreds in Zea mays L.

Heterosis has been widely used in crop breeding and production; however, little is known about the genes controlling trait heterosis. The shortage of genes known to function in heterosis significantly limits our understanding of the molecular basis underlying heterosis. Here, we report 748 genes differentially expressed (DG) in the developing top ear shoots between a maize heterotic F₁ hybrid (Mo17 × B73) and its parental inbreds identified using maize microarrays containing 28,608 unigene features. Of the 748 DG, over 600 were new for the inbred and hybrid combination. The DG were enriched for 35 of the total 213 maize gene ontology (GO) terms, including those describing photosynthesis, respiration, DNA replication, metabolism, and hormone biosynthesis. From the DG, we identified six genes involved in glycolysis, three genes in the citrate cycle, and four genes in the C4-dicarboxylic acid cycle. We mapped 533 of the 748 DG to the maize B73 genome, 298 (55.9 %) of which mapped to the QTL intervals of 11 maize ear traits. Moreover, we compared the repertoire of the DG with that of 14-day seedlings of the same inbred and hybrid combination. Only approximately 5 % of the DG was shared between the two organs and developmental stages. Furthermore, we mapped 417 (55.7 %) of the 748 maize DG to the QTL intervals of 26 rice yield-related traits. Therefore, this study provides a repertoire of genes useful for identification of genes involved in maize ear trait heterosis and information for a better understanding of the molecular basis underlying heterosis in maize.     

DNA methylation in cotton hybrids and their parents

The possible role of methylation in the performance of heterosis has been analyzed in many crops. To further study this possibility, we investigated both the differences in cytosine methylation patterns between cotton heterotic hybrid/nonheterotic hybrids and their parental lines and the change in methylation level from seedling stage to flowering stage by using the methylation-sensitive amplified polymorphism (MSAP) method. The results showed that the number of demethylation loci in highly heterotic hybrids was greater that in lowly heterotic hybrids, and the level of DNA cytosine methylation in cotton at the seedling stage is higher than that at the flowering stage. The altered methylation patterns at low-copy genomic regions can be confirmed by DNA gel blot analysis. A total of 39 fragments that showed different methylation patterns were cloned and sequenced. The methylation status of these genes was modified differentially in hybrid and parents, suggesting that these genes might play a role in the performance of heterosis.     

Overdominant epistatic loci are the primary genetic basis of inbreeding depression and heterosis in rice. I. Biomass and grain yield

To understand the genetic basis of inbreeding depression and heterosis in rice, main-effect and epistatic QTL associated with inbreeding depression and heterosis for grain yield and biomass in five related rice mapping populations were investigated using a complete RFLP linkage map of 182 markers, replicated phenotyping experiments, and the mixed model approach. The mapping populations included 254 F(10) recombinant inbred lines derived from a cross between Lemont (japonica) and Teqing (indica) and two BC and two testcross hybrid populations derived from crosses between the RILs and their parents plus two testers (Zhong 413 and IR64). For both BY and GY, there was significant inbreeding depression detected in the RI population and a high level of heterosis in each of the BC and testcross hybrid populations. The mean performance of the BC or testcross hybrids was largely determined by their heterosis measurements. The hybrid breakdown (part of inbreeding depression) values of individual RILs were negatively associated with the heterosis measurements of their BC or testcross hybrids, indicating the partial genetic overlap of genes causing hybrid breakdown and heterosis in rice. A large number of epistatic QTL pairs and a few main-effect QTL were identified, which were responsible for >65% of the phenotypic variation of BY and GY in each of the populations with the former explaining a much greater portion of the variation. Two conclusions concerning the loci associated with inbreeding depression and heterosis in rice were reached from our results. First, most QTL associated with inbreeding depression and heterosis in rice appeared to be involved in epistasis. Second, most ( approximately 90%) QTL contributing to heterosis appeared to be overdominant. These observations tend to implicate epistasis and overdominance, rather than dominance, as the major genetic basis of heterosis in rice. The implications of our results in rice evolution and improvement are discussed.