Heterotic trait locus (HTL) mapping identifies intra-locus interactions that underlie reproductive hybrid vigor in Sorghum bicolor

Genomic prediction of hybrid performance in grain sorghum (Sorghum bicolor L.)

Genomic selection is expected to improve selection efficiency and genetic gain in breeding programs. The objective of this study was to assess the efficacy of predicting the performance of grain sorghum hybrids using genomic information of parental genotypes. One hundred and two public sorghum inbred parents were genotyped using genotyping-by-sequencing. Ninty-nine of the inbreds were crossed to three tester female parents generating a total of 204 hybrids for evaluation at two environments. The hybrids were sorted in to three sets of 77,59 and 68 and evaluated along with two commercial checks using a randomized complete block design in three replications. The sequence analysis generated 66,265 SNP markers that were used to predict the performance of 204 F1 hybrids resulted from crosses between the parents. Both additive (partial model) and additive and dominance (full model) were constructed and tested using various training population (TP) sizes and cross-validation procedures. Increasing TP size from 41 to 163 increased prediction accuracies for all traits. With the partial model, the five-fold cross validated prediction accuracies ranged from 0.03 for thousand kernel weight (TKW) to 0.58 for grain yield (GY) while it ranged from 0.06 for TKW to 0.67 for GY with the full model. The results suggest that genomic prediction could become an effective tool for predicting the performance of sorghum hybrids based on parental genotypes.

Association mapping and genomic selection for sorghum adaptation to tropical soils of Brazil in a sorghum multiparental random mating population

A multiparental random mating population used in sorghum breeding is amenable for the detection of QTLs related to tropical soil adaptation, fine mapping of underlying genes and genomic selection approaches. Tropical soils where low phosphorus (P) and aluminum (Al) toxicity limit sorghum [Sorghum bicolor (L.) Moench] production are widespread in the developing world. We report on BRP13R, a multiparental random mating population (MP-RMP), which is commonly used in sorghum recurrent selection targeting tropical soil adaptation. Recombination dissipated much of BRP13R's likely original population structure and average linkage disequilibrium (LD) persisted up to 2.5 Mb, establishing BRP13R as a middle ground between biparental populations and sorghum association panels. Genome-wide association mapping (GWAS) identified conserved QTL from previous studies, such as for root morphology and grain yield under low-P, and indicated the importance of dominance in the genetic architecture of grain yield. By overlapping consensus QTL regions, we mapped two candidate P efficiency genes to a ~ 5 Mb region on chromosomes 6 (ALMT) and 9 (PHO2). Remarkably, we find that only 200 progeny genotyped with ~ 45,000 markers in BRP13R can lead to GWAS-based positional cloning of naturally rare, subpopulation-specific alleles, such as for SbMATE-conditioned Al tolerance. Genomic selection was found to be useful in such MP-RMP, particularly if markers in LD with major genes are fitted as fixed effects into GBLUP models accommodating dominance. Shifts in allele frequencies in progeny contrasting for grain yield indicated that intermediate to minor-effect genes on P efficiency, such as SbPSTOL1 genes, can be employed in pre-breeding via allele mining in the base population. Therefore, MP-RMPs such as BRP13R emerge as multipurpose resources for efficient gene discovery and deployment for breeding sorghum cultivars adapted to tropical soils.

Dissecting Adaptive Traits with Nested Association Mapping: Genetic Architecture of Inflorescence Morphology in Sorghum

In the cereal crop sorghum (Sorghum bicolor) inflorescence morphology variation underlies yield variation and confers adaptation across precipitation gradients, but its genetic basis is poorly understood. We characterized the genetic architecture of sorghum inflorescence morphology using a global nested association mapping (NAM) population (2200 recombinant inbred lines) and 198,000 phenotypic observations from multi-environment trials for four inflorescence morphology traits (upper branch length, lower branch length, rachis length, and rachis diameter). Trait correlations suggest that lower and upper branch length are under somewhat independent control, while lower branch length and rachis diameter are highly pleiotropic. Joint linkage and genome-wide association mapping revealed an oligogenic architecture with 1–22 QTL per trait, each explaining 0.1–5.0% of variation across the entire NAM population. There is a significant enrichment (2.twofold) of QTL colocalizing with grass inflorescence gene homologs, notably with orthologs of maize Ramosa2 and rice Aberrant Panicle Organization1 and TAWAWA1. Still, many QTL do not colocalize with inflorescence gene homologs. In global georeferenced germplasm, allelic variation at the major inflorescence QTL is geographically patterned but only weakly associated with the gradient of annual precipitation. Comparison of NAM with diversity panel association suggests that naive association models may capture some true associations not identified by mixed linear models. Overall, the findings suggest that global inflorescence diversity in sorghum is largely controlled by oligogenic, epistatic, and pleiotropic variation in ancestral regulatory networks. The findings also provide a basis for genomics-enabled breeding of locally-adapted inflorescence morphology.

Identification of differentially expressed transcripts at critical developmental stages in sorghum [Sorghum bicolor (L.) Moench] in relation to grain yield heterosis

Evaluation of a set of 10 F₁ hybrids along with their female (27A and 7A) and male parents (C 43, RS 673, RS 627, CB 26, and CB 29) for grain yield and its component traits revealed that grain yield/plant followed by panicle weight, primary branches/panicle, and 100-seed weight exhibited high levels of heterosis. Eight hybrids exhibited 50% or more mid-parent heterosis for grain yield/plant, of which, one hybrid (27A × RS673) recorded heterobeltiosis above 50% (73.61%). Differential display analysis generated about 2995 reproducible transcripts, which were categorized as UPF₁-expressed in any one of the parents and F₁ (10.53-14.76%), BPnF₁-expressed in both parents but not in F₁ (4.56-11.44%), UPnF₁-expressed in either of the parents and not in F₁ (17.95-27.40%), F₁nBP-expressed only in F₁ but not in either of the parents (14.39-20.54%), and UET-expressed in both parents and F₁ (34.52-42.43%). A comparison between high and low heterotic hybrids revealed that the proportions of UPF₁ and F₁nBP transcript patterns were much higher in the former (21.31% and 45.24%) as compared to the latter (16.67% and 32.14%) at the booting and flowering stage, respectively, indicating the role of over-dominance and dominance in the manifestation of grain yield heterosis. Significant positive correlations were observed for differential transcript patterns with mid-parent and better-parent heterosis for the components of grain yield such as primary branches (0.63 and 0.61 at p < 0.01) and 100-seed weight (0.64 and 0.52 at p < 0.01). Cloning and sequence analysis of 16 transcripts that were differentially expressed in hybrids and their parental lines revealed that they code for genes involved in basic cellular processes, cellulose biosynthesis, and assimilate partitioning between various organs and allocation between various pathways, pyrimidine, and polyamine biosynthesis, enhancing ATP production and regulation of plant growth and development.

Genetic architecture of nonadditive inheritance in Arabidopsis thaliana hybrids

The ubiquity of nonparental hybrid phenotypes, such as hybrid vigor and hybrid inferiority, has interested biologists for over a century and is of considerable agricultural importance. Although examples of both phenomena have been subject to intense investigation, no general model for the molecular basis of nonadditive genetic variance has emerged, and prediction of hybrid phenotypes from parental information continues to be a challenge. Here we explore the genetics of hybrid phenotype in 435 Arabidopsis thaliana individuals derived from intercrosses of 30 parents in a half diallel mating scheme. We find that nonadditive genetic effects are a major component of genetic variation in this population and that the genetic basis of hybrid phenotype can be mapped using genome-wide association (GWA) techniques. Significant loci together can explain as much as 20% of phenotypic variation in the surveyed population and include examples that have both classical dominant and overdominant effects. One candidate region inherited dominantly in the half diallel contains the gene for the MADS-box transcription factor AGAMOUS-LIKE 50 (AGL50), which we show directly to alter flowering time in the predicted manner. Our study not only illustrates the promise of GWA approaches to dissect the genetic architecture underpinning hybrid performance but also demonstrates the contribution of classical dominance to genetic variance.

Genes with a Combination of Over-Dominant and Epistatic Effects Underlie Heterosis in Growth of Saccharomyces cerevisiae at High Temperature

Heterosis describes a phenotypic phenomenon of hybrid superiority over its homozygous parents. It is a genetically intriguing phenomenon with great importance for food production. Also called hybrid-vigor, heterosis is created by non-additive effects of genes in a heterozygous hybrid made by crossing two distinct homozygous parents. Few models have been proposed to explain how the combination of parental genes creates an exceptional hybrid performance. Over-dominant mode of inheritance is an attractive model since a single gene can potentially create the heterotic effect, but only a few such loci have been identified. To a collection of 120 hybrids, made by crossing 16 divergent Saccharomyces cerevisiae yeast strains, we applied a method for mapping heterozygous loci that non-additively contribute to heterotic growth at 37°. Among 803 candidate loci that were mapped, five were tested for their heterotic effect by analyzing backcrosses and F2 populations in a specific hybrid background. Consistently with the many mapped loci, specific analyses confirmed the minor heterotic effect of the tested candidate loci. Allele-replacement analyses of one gene, AEP3, further supported its heterotic effect. In addition to over-dominant effects, the contribution of epistasis to heterosis was evident from F2 population and allele-replacement analyses. Pairs of over-dominant genes contributed synergistically to heterosis. We show that minor over-dominant effects of multiple genes can combine to condition heterosis, similarly to loci affecting other quantitative traits. Furthermore, by finding of epistatic interactions between loci that each of them individually has an over-dominant effect on heterosis, we demonstrate how hybrid advantage could benefit from a synergistic combination of two interaction types (over-dominant and synergistic epistatic). Thus, by portraying the underlying genetic complexity, these findings advance our understanding of heterosis.

Genome scan for nonadditive heterotic trait loci reveals mainly underdominant effects in Saccharomyces cerevisiae

The overdominant model of heterosis explains the superior phenotype of hybrids by synergistic allelic interaction within heterozygous loci. To map such genetic variation in yeast, we used a population doubling time dataset of Saccharomyces cerevisiae 16 × 16 diallel and searched for major contributing heterotic trait loci (HTL). Heterosis was observed for the majority of hybrids, as they surpassed their best parent growth rate. However, most of the local heterozygous loci identified by genome scan were surprisingly underdominant, i.e., reduced growth. We speculated that in these loci adverse effects on growth resulted from incompatible allelic interactions. To test this assumption, we eliminated these allelic interactions by creating hybrids with local hemizygosity for the underdominant HTLs, as well as for control random loci. Growth of hybrids was indeed elevated for most hemizygous to HTL genes but not for control genes, hence validating the results of our genome scan. Assessing the consequences of local heterozygosity by reciprocal hemizygosity and allele replacement assays revealed the influence of genetic background on the underdominant effects of HTLs. Overall, this genome-wide study on a multi-parental hybrid population provides a strong argument against single gene overdominance as a major contributor to heterosis, and favors the dominance complementation model.

Dissecting repulsion linkage in the dwarfing gene Dw3 region for sorghum plant height provides insights into heterosis

Heterosis is a main contributor to yield increase in many crop species. Different mechanisms have been proposed for heterosis: dominance, overdominance, epistasis, epigenetics, and protein metabolite changes. However, only limited examples of molecular dissection and validation of these mechanisms are available. Here, we present an example of discovery and validation of heterosis generated by a combination of repulsion linkage and dominance. Using a recombinant inbred line population, a separate quantitative trait locus (QTL) for plant height (qHT7.1) was identified near the genomic region harboring the known auxin transporter Dw3 gene. With two loci having repulsion linkage between two inbreds, heterosis in the hybrid can appear as a single locus with an overdominance mode of inheritance (i.e., pseudo-overdominance). Individually, alleles conferring taller plant height exhibited complete dominance over alleles conferring shorter height. Detailed analyses of different height components demonstrated that qHT7.1 affects both the upper and lower parts of the plant, whereas Dw3 affects only the part below the flag leaf. Computer simulations show that repulsion linkage could influence QTL detection and estimation of effect in segregating populations. Guided by findings in linkage mapping, a genome-wide association study of plant height with a sorghum diversity panel pinpointed genomic regions underlying the trait variation, including Dw1, Dw2, Dw3, Dw4, and qHT7.1. Multilocus mixed model analysis confirmed the advantage of complex trait dissection using an integrated approach. Besides identifying a specific genetic example of heterosis, our research indicated that integrated molecular dissection of complex traits in different population types can enable plant breeders to fine tune the breeding process for crop production.

Genetic differentiation analysis for the identification of complementary parental pools for sorghum hybrid breeding in Ethiopia

The potential for exploiting heterosis for sorghum hybrid production in Ethiopia with improved local adaptation and farmers preferences has been investigated and populations suitable for initial hybrid development have been identified. Hybrids in sorghum have demonstrated increased productivity and stability of performance in the developed world. In Ethiopia, the uptake of hybrid sorghum has been limited to date, primarily due to poor adaptation and absence of farmer's preferred traits in existing hybrids. This study aimed to identify complementary parental pools to develop locally adapted hybrids, through an analysis of whole genome variability of 184 locally adapted genotypes and introduced hybrid parents (R and B). Genetic variability was assessed using genetic distance, model-based STRUCTURE analysis and pair-wise comparison of groups. We observed a high degree of genetic similarity between the Ethiopian improved inbred genotypes and a subset of landraces adapted to lowland agro-ecology with the introduced R lines. This coupled with the genetic differentiation from existing B lines, indicated that these locally adapted genotype groups are expected to have similar patterns of heterotic expression as observed between introduced R and B line pools. Additionally, the hybrids derived from these locally adapted genotypes will have the benefit of containing farmers preferred traits. The groups most divergent from introduced B lines were the Ethiopian landraces adapted to highland and intermediate agro-ecologies and a subset of lowland-adapted genotypes, indicating the potential for increased heterotic response of their hybrids. However, these groups were also differentiated from the R lines, and hence are different from the existing complementary heterotic pools. This suggests that although these groups could provide highly divergent parental pools, further research is required to investigate the extent of heterosis and their hybrid performance.

The regulatory network mediated by circadian clock genes is related to heterosis in rice

Exploitation of heterosis in rice (Oryza sativa L.) has contributed greatly to global food security. In this study, we generated three sets of reciprocal F1 hybrids of indica and japonica subspecies to evaluate the relationship between yield heterosis and the circadian clock. There were no differences in trait performance or heterosis between the reciprocal hybrids, indicating no maternal effects on heterosis. The indica-indica and indica-japonica reciprocal F1 hybrids exhibited pronounced heterosis for chlorophyll and starch content in leaves and for grain yield/biomass. In contrast, the japonica-japonica F1 hybrids showed low heterosis. The three circadian clock genes investigated expressed in an above-high-parent pattern (AHP) at seedling stage in all the hybrids. The five genes downstream of the circadian clock, and involved in chlorophyll and starch metabolic pathways, were expressed in AHP in hybrids with strong better-parent heterosis (BPH). Similarly, three of these five genes in the japonica-japonica F1 hybrids showing low BPH were expressed in positive overdominance, but the other two genes were expressed in additive or negative overdominance. These results indicated that the expression patterns of circadian clock genes and their downstream genes are associated with heterosis, which suggests that the circadian rhythm pathway may be related to heterosis in rice.

Consequences of hybridization and heterozygosity on plant vigor and phenotypic stability

The genomic makeup and phenotypes of plants are diversifying, in part due to artificial or natural selection in agricultural and natural environments. Utilization of these variations to enhance crop productivity requires an understanding of the relationships between genotype and phenotype in inbreds and hybrids derived from crosses between these populations. This review highlights recent studies on hybrid vigor (heterosis) and the related phenomenon of hybrid weakness - two types of non-additive inheritance. Heterosis is a phenomenon whereby the phenotype of first-generation hybrids is superior to that of their parents. Intralocus interactions between alleles, complementation of dominant alleles, or inter-loci epistatic interactions are genetic mechanisms that may cause non-additive phenotypic inheritance in hybrids. However, there are different views on what portion of the heterotic variation is modulated by each of these mechanisms. Another aspect of plant vigor is phenotypic stability or robustness in different environments and how this is influenced by gene heterozygosity. Hybrids are not necessarily more phenotypically stable than inbreds since local heterozygosity might be associated with negative effects on biochemical activities. This review integrates genetic and biochemical considerations to illustrate how these relationships may be tightly linked with breeding system and sequence divergence.

A Boolean gene regulatory model of heterosis and speciation

BACKGROUND

Modelling genetic phenomena affecting biological traits is important for the development of agriculture as it allows breeders to predict the potential of breeding for certain traits. One such phenomenon is heterosis or hybrid vigor: crossing individuals from genetically distinct populations often results in improvements in quantitative traits, such as growth rate, biomass production and stress resistance. Heterosis has become a very useful tool in global agriculture, but its genetic basis remains controversial and its effects hard to predict. We have taken a computational approach to studying heterosis, developing a simulation of evolution, independent reassortment of alleles and hybridization of Gene Regulatory Networks (GRNs) in a Boolean framework. These artificial regulatory networks exhibit topological properties that reflect those observed in biology, and fitness is measured as the ability of a network to respond to external inputs in a pre-defined way.

RESULTS

Our model reproduced common experimental observations on heterosis using only biologically justified parameters, such as mutation rates. Hybrid vigor was observed and its extent was seen to increase as parental populations diverged, up until a point of sudden collapse of hybrid fitness. Thus, the model also describes a process akin to speciation due to genetic incompatibility of the separated populations. We also reproduce, for the first time in a model, the fact that hybrid vigor cannot easily be fixed by within a breeding line, currently an important limitation of the use of hybrid crops. The simulation allowed us to study the effects of three standard models for the genetic basis of heterosis: dominance, over-dominance, and epistasis.

CONCLUSION

This study describes the most detailed simulation of heterosis using gene regulatory networks to date and reproduces several phenomena associated with heterosis for the first time in a model. The level of detail in our model allows us to suggest possible warning signs of the impending collapse of hybrid vigor in breeding. In addition, the simulation provides a framework that can be extended to study other aspects of heterosis and alternative evolutionary scenarios.

What is crop heterosis: new insights into an old topic

Heterosis (or hybrid vigor) refers to a natural phenomenon whereby hybrid offspring of genetically diverse individuals out-perform their parents in multiple traits including yield, adaptability and resistances to biotic and abiotic stressors. Innovations in technology and research continue to clarify the mechanisms underlying crop heterosis, however the intrinsic relationship between the biological basis of heterosis remain unclear. In this review, we aim to provide insight into the molecular genetic basis of heterosis by presenting recent advances in the 'omics' of heterosis and the role of non-coding regions, particularly in relation to energy-use efficiency. We propose that future research should focus on integrating the expanding datasets from different species and hybrid combinations, to mine key heterotic genes and unravel interactive 'omics' networks associated with heterosis. Improved understanding of heterosis and the biological basis for its manipulation in agriculture should help to streamline its use in enhancing crop productivity.

Mapping heterotic loci for yield and agronomic traits using chromosome segment introgression lines in cotton

In the present study, a set of chromosome segment introgression lines (CSILs) using Gossypium hirsutum L. TM-1 as the recipient parent and G. barbadense Hai7124 as the donor parent were used to explore the genetic basis of heterosis for interspecific hybrids. Two sets of F₁ populations individually derived from CSILs crossing with both parents were configured to investigate heterotic loci (HL) and substitution effect loci (SL). A total of 58 HL and 39 SL were identified in 3 years. One stable HL, hLP-A4-3, could be detected in all 3 years. Three HLs, hBS-A8-1, hLP-D6-1, and hSI-D7-11, could be detected in 2 years. Four SLs, sBS-D7-1, sLP-A8-1, sLP-D7-1, and sLP-D12-1, could be detected in 2 years. HL and SL tended to be distributed in some HL-rich chromosome segments with close positions. Compared with QTL detected in a former study, HL showed little overlap with QTL, indicating that trait phenotype and heterosis might be controlled by different sets of loci. All three forms of genetic effects (partial-, full-, over-dominant) were identified, while the over-dominant effect made the main contribution to heterosis. These results may help lay the foundation for clarifying the heredity mechanism of heterosis in cotton.