Heterosis: revisiting the magic

Haplotype-resolved T2T genome assemblies and pangenome graph of pear reveal diverse patterns of allele-specific expression and the genomic basis of fruit quality traits

Hybrid crops often exhibit increased yield and greater resilience, yet the genomic mechanism(s) underlying hybrid vigor or heterosis remain unclear, hindering our ability to predict the expression of phenotypic traits in hybrid breeding. Here, we generated haplotype-resolved T2T genome assemblies of two pear hybrid varieties, 'Yuluxiang' (YLX) and 'Hongxiangsu' (HXS), which share the same maternal parent but differ in their paternal parents. We then used these assemblies to explore the genome-scale landscape of allele-specific expression (ASE) and create a pangenome graph for pear. ASE was observed for close to 6000 genes in both hybrid cultivars. A subset of ASE genes related to aspects of fruit quality such as sugars, organic acids, and cuticular wax were identified, suggesting their important contributions to heterosis. Specifically, Ma1, a gene regulating fruit acidity, is absent in the paternal haplotypes of HXS and YLX. A pangenome graph was built based on our assemblies and seven published pear genomes. Resequencing data for 139 cultivated pear genotypes (including 97 genotypes sequenced here) were subsequently aligned to the pangenome graph, revealing numerous structural variant hotspots and selective sweeps during pear diversification. As predicted, the Ma1 allele was found to be absent in varieties with low organic acid content, and this association was functionally validated by Ma1 overexpression in pear fruit and calli. Overall, these results reveal the contributions of ASE to fruit-quality heterosis and provide a robust pangenome reference for high-resolution allele discovery and association mapping.

Integrative multi-omics analysis reveals genetic and heterotic contributions to male fertility and yield in potato

The genetic analysis of potato is hampered by the complexity of tetrasomic inheritance. An ongoing effort aims to transform the clonally propagated tetraploid potato into a seed-propagated diploid crop, which would make genetic analyses much easier owing to disomic inheritance. Here, we construct and report the large-scale genetic and heterotic characteristics of a diploid F2 potato population derived from the cross of two highly homozygous inbred lines. We investigate 20,382 traits generated from multi-omics dataset and identify 25,770 quantitative trait loci (QTLs). Coupled with gene expression data, we construct a systems-genetics network for gene discovery in potatoes. Importantly, we explore the genetic basis of heterosis in this population, especially for yield and male fertility heterosis. We find that positive heterotic effects of yield-related QTLs and negative heterotic effects of metabolite QTLs (mQTLs) contribute to yield heterosis. Additionally, we identify a PME gene with a dominance heterotic effect that plays an important role in male fertility heterosis. This study provides genetic resources for the potato community and will facilitate the application of heterosis in diploid potato breeding.

Breed of origin analysis in genome-wide association studies: enhancing SNP-based insights into production traits in a commercial Brangus population

BACKGROUND

Carcass weight (HCW) and marbling (MARB) are critical for meat quality and market value in beef cattle. In composite breeds like Brangus, which meld the genetics of Angus and Brahman, SNP-based analyses have illuminated some genetic influences on these traits, but they fall short in fully capturing the nuanced effects of breed of origin alleles (BOA) on these traits. Focus on the impacts of BOA on phenotypic features within Brangus populations can result in a more profound understanding of the specific influences of Angus and Brahman genetics. Moreover, the consideration of BOA becomes particularly significant when evaluating dominance effects contributing to heterosis in crossbred populations. BOA provides a more comprehensive measure of heterosis due to its ability to differentiate the distinct genetic contributions originating from each parent breed. This detailed understanding of genetic effects is essential for making informed breeding decisions to optimize the benefits of heterosis in composite breeds like Brangus.

OBJECTIVE

This study aims to identify quantitative trait loci (QTL) influencing HCW and MARB by utilizing SNP and BOA information, incorporating additive, dominance, and overdominance effects within a multi-generational Brangus commercial herd.

METHODS

We analyzed phenotypic data from 1,066 genotyped Brangus steers. BOA inference was performed using LAMP-LD software using Angus and Brahman reference sets. SNP-based and BOA-based GWAS were then conducted considering additive, dominance, and overdominance models.

RESULTS

The study identified numerous QTLs for HCW and MARB. A notable QTL for HCW was associated to the SGCB gene, pivotal for muscle growth, and was identified solely in the BOA GWAS. Several BOA GWAS QTLs exhibited a dominance effect underscoring their importance in estimating heterosis.

CONCLUSIONS

Our findings demonstrate that SNP-based methods may not detect all genetic variation affecting economically important traits in composite breeds. BOA inclusion in genomic evaluations is crucial for identifying genetic regions contributing to trait variation and for understanding the dominance value underpinning heterosis. By considering BOA, we gain a deeper understanding of genetic interactions and heterosis, which is integral to advancing breeding programs. The incorporation of BOA is recommended for comprehensive genomic evaluations to optimize trait improvements in crossbred cattle populations.

Multi-tissue transcriptome profiling linked the association between tissue-specific circRNAs and the heterosis for feed intake and efficiency in chicken

Heterosis has been widely utilized in chickens. The nonadditive inheritance of genes contributes to this biological phenomenon. However, the role of circRNAs played in the heterosis is poorly determined. In this study, we observed divergent heterosis for residual feed intake (RFI) between 2 crossbreds derived from a reciprocal cross between White Leghorns and Beijing You chickens. Then, circRNA landscape for 120 samples covering the hypothalamus, liver, duodenum mucosa and ovary were profiled to elucidate the regulatory mechanisms of heterosis. We detected that a small proportion of circRNAs (7.83-20.35%) were additively and non-additively expressed, in which non-additivity was a major inheritance of circRNAs in the crossbreds. Tissue-specific expression of circRNAs was prevalent across 4 tissues. Weighted gene co-expression network analysis revealed circRNA-mRNA co-expression modules associated with feed intake and RFI in the hypothalamus and liver, and the co-expressed genes were enriched in oxidative phosphorylation pathway. We further identified 8 nonadditive circRNAs highly correlated with 16 nonadditive genes regulating negative heterosis for RFI in the 2 tissues. Circ-ITSN2 was validated in the liver tissue for its significantly positive correlation with PGPEP1L. Moreover, the bioinformatic analysis indicated that candidate circRNAs might be functioned by binding the microRNAs and interacting with the RNA binding proteins. The integration of multi-tissue transcriptome firstly linked the association between tissue-specific circRNAs and the heterosis for feed intake and efficiency in chicken, which provide novel insights into the molecular mechanism underlying heterosis for feed efficiency. The validated circRNAs can act as potential biomarkers for predicting RFI and its heterosis.

Schistosoma haematobium and Schistosoma bovis first generation hybrids undergo gene expressions changes consistent with species compatibility and heterosis

When two species hybridize, the two parental genomes are brought together and some alleles might interact for the first time. To date, the extent of the transcriptomic changes in first hybrid generations, along with their functional outcome constitute an important knowledge gap, especially in parasite species. Here we explored the molecular and functional outcomes of hybridization in first-generation hybrids between the blood fluke parasites Schistosoma haematobium and S. bovis. Through a transcriptomic approach, we measured gene expression in both parental species and hybrids. We described and quantified expression profiles encountered in hybrids along with the main biological processes impacted. Up to 7,100 genes fell into a particular hybrid expression profile (intermediate between the parental expression levels, over-expressed, under-expressed, or expressed like one of the parental lines). Most of these genes were different depending on the direction of the parental cross (S. bovis mother and S. haematobium father or the reverse) and depending on the sex. For a given sex and cross direction, the vast majority of genes were hence unassigned to a hybrid expression profile: either they were differentially expressed genes but not typical of any hybrid expression profiles or they were not differentially expressed neither between hybrids and parental lines nor between parental lines. The most prevalent profile of gene expression in hybrids was the intermediate one (24% of investigated genes). These results suggest that transcriptomic compatibility between S. haematobium and S. bovis remains quite high. We also found support for an over-dominance model (over- and under-expressed genes in hybrids compared to parental lines) potentially associated with heterosis. In females in particular, processes such as reproductive processes, metabolism and cell interactions as well as signaling pathways were indeed affected. Our study hence provides new insight on the biology of Schistosoma hybrids with evidences supporting compatibility and heterosis.

A higher-yield hybrid rice is achieved by assimilating a dominant heterotic gene in inbred parental lines

The exploitation of heterosis to integrate parental advantages is one of the fastest and most efficient ways of rice breeding. The genomic architecture of heterosis suggests that the grain yield is strongly correlated with the accumulation of numerous rare superior alleles with positive dominance. However, the improvements in yield of hybrid rice have shown a slowdown or even plateaued due to the limited availability of complementary superior alleles. In this study, we achieved a considerable increase in grain yield of restorer lines by inducing an alternative splicing event in a heterosis gene OsMADS1 through CRISPR-Cas9, which accounted for approximately 34.1%-47.5% of yield advantage over their corresponding inbred rice cultivars. To achieve a higher yield in hybrid rice, we crossed the gene-edited restorer parents harbouring OsMADS1GW3p6 with the sterile lines to develop new rice hybrids. In two-line hybrid rice Guang-liang-you 676 (GLY676), the yield of modified hybrids carrying the homozygous heterosis gene OsMADS1GW3p6 significantly exceeded that of the original hybrids with heterozygous OsMADS1. Similarly, the gene-modified F1 hybrids with heterozygous OsMADS1GW3p6 increased grain yield by over 3.4% compared to the three-line hybrid rice Quan-you-si-miao (QYSM) with the homozygous genotype of OsMADS1. Our study highlighted the great potential in increasing the grain yield of hybrid rice by pyramiding a single heterosis gene via CRISPR-Cas9. Furthermore, these results demonstrated that the incomplete dominance of heterosis genes played a major role in yield-related heterosis and provided a promising strategy for breeding higher-yielding rice varieties above what is currently achievable.

Imprinting but not cytonuclear interactions determines seed size heterosis in Arabidopsis hybrids

The parent-of-origin effect on seeds can result from imprinting (unequal expression of paternal and maternal alleles) or combinational effects between cytoplasmic and nuclear genomes, but their relative contributions remain unknown. To discern these confounding factors, we produced cytoplasmic-nuclear substitution (CNS) lines using recurrent backcrossing in Arabidopsis (Arabidopsis thaliana) ecotypes Col-0 and C24. These CNS lines differed only in the nuclear genome (imprinting) or cytoplasm. The CNS reciprocal hybrids with the same cytoplasm displayed ∼20% seed size difference, whereas the seed size was similar between the reciprocal hybrids with fixed imprinting. Transcriptome analyses in the endosperm of CNS hybrids using laser-capture microdissection identified 104 maternally expressed genes (MEGs) and 90 paternally expressed genes (PEGs). These imprinted genes were involved in pectin catabolism and cell wall modification in the endosperm. Homeodomain Glabrous9 (HDG9), an epiallele and one of 11 cross-specific imprinted genes, affected seed size. In the embryo, there were a handful of imprinted genes in the CNS hybrids but only 1 was expressed at higher levels than in the endosperm. AT4G13495 was found to encode a long-noncoding RNA (lncRNA), but no obvious seed phenotype was observed in lncRNA knockout lines. Nuclear RNA Polymerase D1 (NRPD1), encoding the largest subunit of RNA Pol IV, was involved in the biogenesis of small interfering RNAs. Seed size and embryos were larger in the cross using nrpd1 as the maternal parent than in the reciprocal cross, supporting a role of the maternal NRPD1 allele in seed development. Although limited ecotypes were tested, these results suggest that imprinting and the maternal NRPD1-mediated small RNA pathway play roles in seed size heterosis in plant hybrids.

Morphological variation of tail bone among two chicken breeds and their F1 progeny

Fancy breeds of Japanese indigenous chicken display extensive morphological diversity, particularly in tail feathers. Although marked differences in tail and bone traits have been reported between Tosa-jidori (wild type) and Minohikichabo (rich type) breeds, little is known about the pattern of genetic inheritance in cross experiments. Therefore, this study aimed to investigate the strain and sex effects, and inheritance patterns, in the morphometric variation of pygostyle bones among Tosa-jidori, Minohikichabo, and their F1 hybrids. Five morphological traits, angle of the apex of the pygostyle, pygostyle length, margo cranialis length, tail feather number, and body weight, were evaluated at the adult stage. A significant strain difference was detected in all traits, whereas significant sex differences were observed in only three traits, but not in the angle of the apex of the pygostyle and tail feather number. In F1 hybrids, the angle of the apex of the pygostyle was significantly different to that of Tosa-jidori but not that of Minohikichabo, whereas the pygostyle length and tail number of F1 hybrids were significantly different from those of Minohikichabo but not those of Tosa-jidori. A significant heterosis effect was found in the margo cranialis length and body weight. All five traits showed nonadditive inheritance patterns but varied in each trait between partial dominance (angle of the apex of pygostyle), full dominance (pygostyle length and tail feather number), and over-dominance (margo cranialis length and body weight). Interestingly, different patterns of genetic inheritance in the F1 hybrid were observed at different locations, even within the same pygostyle bone. Using the Japanese indigenous chicken model, these results provide a substantial step toward understanding the genetic architecture of morphology in chickens.

Editing of ORF138 restores fertility of Ogura cytoplasmic male sterile broccoli via mitoTALENs

Cytoplasmic male sterility (CMS), encoded by the mitochondrial open reading frames (ORFs), has long been used to economically produce crop hybrids. However, the utilization of CMS also hinders the exploitation of sterility and fertility variation in the absence of a restorer line, which in turn narrows the genetic background and reduces biodiversity. Here, we used a mitochondrial targeted transcription activator-like effector nuclease (mitoTALENs) to knock out ORF138 from the Ogura CMS broccoli hybrid. The knockout was confirmed by the amplification and re-sequencing read mapping to the mitochondrial genome. As a result, knockout of ORF138 restored the fertility of the CMS hybrid, and simultaneously manifested a cold-sensitive male sterility. ORF138 depletion is stably inherited to the next generation, allowing for direct use in the breeding process. In addition, we proposed a highly reliable and cost-effective toolkit to accelerate the life cycle of fertile lines from CMS-derived broccoli hybrids. By applying the k-mean clustering and interaction network analysis, we identified the central gene networks involved in the fertility restoration and cold-sensitive male sterility. Our study enables mitochondrial genome editing via mitoTALENs in Brassicaceae vegetable crops and provides evidence that the sex production machinery and its temperature-responsive ability are regulated by the mitochondria.

Self-incompatibility phenotypes of SRK mutants can be predicted with high accuracy

Only very limited information is available on why some non-synonymous variants severely alter gene function while others have no effect. To identify the characteristic features of mutations that strongly influence gene function, this study focused on S-locus receptor kinase, SRK, which encodes a highly polymorphic receptor kinase expressed in stigma papillary cells that underlies a female determinant of self-incompatibility in Brassicaceae. A set of 299 Arabidopsis thaliana transformants expressing mutated SRKb from A. lyrata was constructed and analyzed to determine the genotype and self-incompatibility phenotype of each transformant. Almost all the transformants showing the self-incompatibility defect contained mutations in AlSRKb that altered localization to the plasma membrane. The observed mutations occurred in amino acid residues that were highly conserved across S haplotypes and whose predicted locations were in the interior of the protein. These mutations were likely to underlie the self-incompatibility defect as they caused significant changes to amino acid properties. Such findings suggested that mutations causing the self-incompatibility defect were more likely to result from changes to AlSRKb biosynthesis than from loss of function. In addition, this study showed the RandomForest and Extreme Gradient Boosting methods could predict self-incompatibility phenotypes of SRK mutants with high accuracy.

A Case Study from the Overexpression of OsTZF5, Encoding a CCCH Tandem Zinc Finger Protein, in Rice Plants Across Nineteen Yield Trials

BACKGROUND

Development of transgenic rice overexpressing transcription factors involved in drought response has been previously reported to confer drought tolerance and therefore represents a means of crop improvement. We transformed lowland rice IR64 with OsTZF5, encoding a CCCH-tandem zinc finger protein, under the control of the rice LIP9 stress-inducible promoter and compared the drought response of transgenic lines and nulls to IR64 in successive screenhouse paddy and field trials up to the T6 generation.

RESULTS

Compared to the well-watered conditions, the level of drought stress across experiments varied from a minimum of - 25 to - 75 kPa at a soil depth of 30 cm which reduced biomass by 30-55% and grain yield by 1-92%, presenting a range of drought severities. OsTZF5 transgenic lines showed high yield advantage under drought over IR64 in early generations, which was related to shorter time to flowering, lower shoot biomass and higher harvest index. However, the increases in values for yield and related traits in the transgenics became smaller over successive generations despite continued detection of drought-induced transgene expression as conferred by the LIP9 promoter. The decreased advantage of the transgenics over generations tended to coincide with increased levels of homozygosity. Background cleaning of the transgenic lines as well as introgression of the transgene into an IR64 line containing major-effect drought yield QTLs, which were evaluated starting at the BC3F1 and BC2F3 generation, respectively, did not result in consistently increased yield under drought as compared to the respective checks.

CONCLUSIONS

Although we cannot conclusively explain the genetic factors behind the loss of yield advantage of the transgenics under drought across generations, our results help in distinguishing among potential drought tolerance mechanisms related to effectiveness of the transgenics, since early flowering and harvest index most closely reflected the levels of yield advantage in the transgenics across generations while reduced biomass did not.

Increased long-distance and homo-trans interactions related to H3K27me3 in Arabidopsis hybrids

In plants, the genome structure of hybrids changes compared with their parents, but the effects of these changes in hybrids remain elusive. Comparing reciprocal crosses between Col × C24 and C24 × Col in Arabidopsis using high-throughput chromosome conformation capture assay (Hi-C) analysis, we found that hybrid three-dimensional (3D) chromatin organization had more long-distance interactions relative to parents, and this was mainly located in promoter regions and enriched in genes with heterosis-related pathways. The interactions between euchromatin and heterochromatin were increased, and the compartment strength decreased in hybrids. In compartment domain (CD) boundaries, the distal interactions were more in hybrids than their parents. In the hybrids of CURLY LEAF (clf) mutants clfCol  × clfC24 and clfC24  × clfCol , the heterosis phenotype was damaged, and the long-distance interactions in hybrids were fewer than in their parents with lower H3K27me3. ChIP-seq data revealed higher levels of H3K27me3 in the region adjacent to the CD boundary and the same interactional homo-trans sites in the wild-type (WT) hybrids, which may have led to more long-distance interactions. In addition, the differentially expressed genes (DEGs) located in the boundaries of CDs and loop regions changed obviously in WT, and the functional enrichment for DEGs was different between WT and clf in the long-distance interactions and loop regions. Our findings may therefore propose a new epigenetic explanation of heterosis in the Arabidopsis hybrids and provide new insights into crop breeding and yield increase.

Gene expression bias between the subgenomes of allopolyploid hybrids is an emergent property of the kinetics of expression

Hybridization coupled to polyploidy, or allopolyploidy, has dramatically shaped the evolution of flowering plants, teleost fishes, and other lineages. Studies of recently formed allopolyploid plants have shown that the two subgenomes that merged to form that new allopolyploid do not generally express their genes equally. Instead, one of the two subgenomes expresses its paralogs more highly on average. Meanwhile, older allopolyploidy events tend to show biases in duplicate losses, with one of the two subgenomes retaining more genes than the other. Since reduced expression is a pathway to duplicate loss, understanding the origins of expression biases may help explain the origins of biased losses. Because we expect gene expression levels to experience stabilizing selection, our conceptual frameworks for how allopolyploid organisms form tend to assume that the new allopolyploid will show balanced expression between its subgenomes. It is then necessary to invoke phenomena such as differences in the suppression of repetitive elements to explain the observed expression imbalances. Here we show that, even for phenotypically identical diploid progenitors, the inherent kinetics of gene expression give rise to biases between the expression levels of the progenitor genes in the hybrid. Some of these biases are expected to be gene-specific and not give rise to global differences in progenitor gene expression. However, particularly in the case of allopolyploids formed from progenitors with different genome sizes, global expression biases favoring one subgenome are expected immediately on formation. Hence, expression biases are arguably the expectation upon allopolyploid formation rather than a phenomenon needing explanation. In the future, a deeper understanding of the kinetics of allopolyploidy may allow us to better understand both biases in duplicate losses and hybrid vigor.

Nonadditive and allele-specific expression of ghrelin in hybrid tilapia

Background

Interspecies hybridization is an important breeding method to generate fishes with heterosis in aquaculture. Using this method, hybrid Nile tilapia (Oreochromis niloticus, ♀) × blue tilapia (Oreochromis aureus, ♂) has been produced and widely farmed due to its growth and appetite superiorities. However, the genetic mechanism of these advanced traits is still not well understood. Ghrelin is a crucial gene that regulates growth and appetite in fishes. In the present study, we focused on the expression characteristics and its regulation of ghrelin in the hybrid.

Results

The tissue distribution analysis showed that ghrelin was predominantly expressed in the stomach in the hybrid. Ghrelin was more highly expressed in the stomach in the hybrid and Nile tilapia, compared to blue tilapia, showing a nonadditive pattern. Two single-nucleotide polymorphism (SNP) sites were identified including T/C and C/G from the second exon in the ghrelin gene from Nile tilapia and blue tilapia. By pyrosequencing based on the SNP sites, the allele-specific expression (ASE) of ghrelin in the hybrid was assayed. The result indicated that ghrelin in the hybrid showed higher maternal allelic transcript ratios. Fasting significantly increased ghrelin overall expression at 4, 8, 12, 24, and 48 h. In addition, higher maternal allelic transcript ratios were not changed in the fasting hybrids at 48 h. The cis and trans effects were determined by evaluating the overall expression and ASE values in the hybrid. The expression of ghrelin was mediated by compensating cis and trans effects in hybrid.

Conclusion

In summary, the present lines of evidence showed the nonadditive expression of ghrelin in the hybrid tilapia and its regulation by subgenomes, offering new insight into gene expression characteristics in hybrids.

Hybridization in agricultural weeds: A review from ecological, evolutionary, and management perspectives

Agricultural weeds frequently hybridize with each other or with related crop species. Some hybrid weeds exhibit heterosis (hybrid vigor), which may be stabilized through mechanisms like genome duplication or vegetative reproduction. Even when heterosis is not stabilized, hybridization events diversify weed gene pools and often enable adaptive introgression. Consequently, hybridization may promote weed evolution and exacerbate weed-crop competition. However, hybridization does not always increase weediness. Even when viable and fertile, hybrid weeds sometimes prove unsuccessful in crop fields. This review provides an overview of weed hybridization and its management implications. We describe intrinsic and extrinsic factors that influence hybrid fitness in agroecosystems. We also survey the rapidly growing literature on crop-weed hybridization and the link between hybridization and invasiveness. These topics are increasingly relevant in this era of genetic tools for crop improvement, intensive and simplified cropping systems, and globalized trade. The review concludes with suggested research priorities, including hybridization in the context of climate change, plant-insect interactions, and redesigned weed management programs. From a weed management perspective, hybridization is one of many reasons that researchers and land managers must diversify their weed control toolkits.

Introduction of barnase/barstar in soybean produces a rescuable male sterility system for hybrid breeding

Hybrid breeding for increased vigour has been used for over a century to boost agricultural outputs without requiring higher inputs. While this approach has led to some of the most substantial gains in crop productivity, breeding barriers have fundamentally limited soybean (Glycine max) from reaping the benefits of hybrid vigour. Soybean flowers self-pollinate prior to opening and thus are not readily amenable to outcrossing. In this study, we demonstrate that the barnase/barstar male sterility/rescue system can be used in soybean to produce hybrid seeds. By expressing the cytotoxic ribonuclease, barnase, under a tapetum-specific promoter in soybean anthers, we are able to completely block pollen maturation, creating male sterile plants. We show that fertility can be rescued in the F1 generation of these barnase-expressing lines when they are crossed with pollen from plants that express the barnase inhibitor, barstar. Importantly, we found that the successful rescue of male fertility is dependent on the relative dosage of barnase and barstar. When barnase and barstar were expressed under the same tapetum-specific promoter, the F1 offspring remained male sterile. When we expressed barstar under a relatively stronger promoter than barnase, we were able to achieve a successful rescue of male fertility in the F1 generation. This work demonstrates the successful implementation of a biotechnology approach to produce fertile hybrid offspring in soybean.

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.

The mRNA-lncRNA landscape of multiple tissues uncovers key regulators and molecular pathways that underlie heterosis for feed intake and efficiency in laying chickens

BACKGROUND

Heterosis is routinely exploited to improve animal performance. However, heterosis and its underlying molecular mechanism for feed intake and efficiency have been rarely explored in chickens. Feed efficiency continues to be an important breeding goal trait since feed accounts for 60 to 70% of the total production costs in poultry. Here, we profiled the mRNA-lncRNA landscape of 96 samples of the hypothalamus, liver and duodenum mucosa from White Leghorn (WL), Beijing-You chicken (YY), and their reciprocal crosses (WY and YW) to elucidate the regulatory mechanisms of heterosis.

RESULTS

We observed negative heterosis for both feed intake and residual feed intake (RFI) in YW during the laying period from 43 to 46 weeks of age. Analysis of the global expression pattern showed that non-additivity was a major component of the inheritance of gene expression in the three tissues for YW but not for WY. The YW-specific non-additively expressed genes (YWG) and lncRNA (YWL) dominated the total number of non-additively expressed genes and lncRNA in the hypothalamus and duodenum mucosa. Enrichment analysis of YWG showed that mitochondria components and oxidation phosphorylation (OXPHOS) pathways were shared among the three tissues. The OXPHOS pathway was enriched by target genes for YWL with non-additive inheritance of expression in the liver and duodenum mucosa. Weighted gene co-expression network analysis revealed divergent co-expression modules associated with feed intake and RFI in the three tissues from WL, YW, and YY. Among the negatively related modules, the OXPHOS pathway was enriched by hub genes in the three tissues, which supports the critical role of oxidative phosphorylation. Furthermore, protein quantification of ATP5I was highly consistent with ATP5I expression in the liver, which suggests that, in crossbred YW, non-additive gene expression is down-regulated and decreases ATP production through oxidative phosphorylation, resulting in negative heterosis for feed intake and efficiency.

CONCLUSIONS

Our results demonstrate that non-additively expressed genes and lncRNA involved in oxidative phosphorylation in the hypothalamus, liver, and duodenum mucosa are key regulators of the negative heterosis for feed intake and RFI in layer chickens. These findings should facilitate the rational choice of suitable parents for producing crossbred chickens.

Can heterosis and inbreeding depression explain the maintenance of outcrossing in a cleistogamous perennial?

PREMISE

What maintains mixed mating is an evolutionary enigma. Cleistogamy-the production of both potentially outcrossing chasmogamous and obligately selfing cleistogamous flowers on the same individual plant-is an excellent system to study the costs of selfing. Inbreeding depression can prevent the evolution of greater selfing within populations, and heterosis in crosses between populations may further tip the balance in favor of outcrossing. Few empirical estimates of inbreeding depression and heterosis in the same system exist for cleistogamous species.

METHODS

We investigate the potential costs of selfing by quantifying inbreeding depression and heterosis in three populations of the cleistogamous perennial Ruellia humilis Nutt (Acanthaceae). We performed three types of hand-pollinations-self, outcross-within, and outcross-between populations-and measured seed number, germination, total flower production, and estimated cumulative fitness for the resulting progeny in a greenhouse experiment.

RESULTS

We found moderate inbreeding depression for cumulative fitness (<30%) in two populations, but outbreeding depression for crosses within a third population (-26%). For between-population crosses, there was weak to modest heterosis (11-47%) in two of the population combinations, but modest to strong outbreeding depression (-21 to -71%) in the other four combinations.

CONCLUSIONS

Neither inbreeding depression nor heterosis was of sufficient magnitude to explain the continued production of chasmogamous flowers given the relative energetic advantage of cleistogamous flowers previously estimated for these populations. Outbreeding depression either within or between populations makes the maintenance of chasmogamous flowers even harder to explain. More information is needed on the genetic basis of cleistogamy to resolve this conundrum.

Imprinting but not cytonuclear interactions affects parent-of-origin effect on seed size in Arabidopsis hybrids

The parent-of-origin effect on seed size can result from imprinting or a combinational effect between cytoplasmic and nuclear genomes, but their relative contributions remain unknown. To discern these confounding effects, we generated cytoplasmic-nuclear substitution (CNS) lines using recurrent backcrossing in the Arabidopsis thaliana ecotypes Col-0 and C24. These CNS lines differ only in the nuclear genome (imprinting) or in the cytoplasm. The CNS reciprocal hybrids with the same cytoplasm display a ~20% seed size difference as observed in the conventional hybrids. However, seed size is similar between the reciprocal cybrids with fixed imprinting. Transcriptome analyses in the endosperm of CNS hybrids using laser-capture microdissection have identified 104 maternally expressed genes (MEGs) and 90 paternally-expressed genes (PEGs). These imprinted genes are involved in pectin catabolism and cell wall modification in the endosperm. HDG9, an epiallele and one of 11 cross-specific imprinted genes, controls seed size. In the embryo, a handful of imprinted genes is found in the CNS hybrids but only one is expressed higher in the embryo than endosperm. AT4G13495 encodes a long-noncoding RNA (lncRNA), but no obvious seed phenotype is observed in the lncRNA knockout lines. NRPD1, encoding the largest subunit of RNA Pol IV, is involved in the biogenesis of small interfering RNAs. Seed size and embryo is larger in the cross using nrpd1 as the maternal parent than in the reciprocal cross. In spite of limited ecotypes tested, these results suggest potential roles of imprinting and NRPD1-mediated small RNA pathway in seed size variation in hybrids.

Pervasive under-dominance in gene expression underlying emergent growth trajectories in Arabidopsis thaliana hybrids

BACKGROUND

Complex traits, such as growth and fitness, are typically controlled by a very large number of variants, which can interact in both additive and non-additive fashion. In an attempt to gauge the relative importance of both types of genetic interactions, we turn to hybrids, which provide a facile means for creating many novel allele combinations.

RESULTS

We focus on the interaction between alleles of the same locus, i.e., dominance, and perform a transcriptomic study involving 141 random crosses between different accessions of the plant model species Arabidopsis thaliana. Additivity is rare, consistently observed for only about 300 genes enriched for roles in stress response and cell death. Regulatory rare-allele burden affects the expression level of these genes but does not correlate with F1 rosette size. Non-additive, dominant gene expression in F1 hybrids is much more common, with the vast majority of genes (over 90%) being expressed below the parental average. Unlike in the additive genes, regulatory rare-allele burden in the dominant gene set is strongly correlated with F1 rosette size, even though it only mildly covaries with the expression level of these genes.

CONCLUSIONS

Our study underscores under-dominance as the predominant gene action associated with emergence of rosette growth trajectories in the A. thaliana hybrid model. Our work lays the foundation for understanding molecular mechanisms and evolutionary forces that lead to dominance complementation of rare regulatory alleles.

Unlikely heroes on the long and winding road to potato inbreeding

Conversion of potato from a tetraploid, heterozygous, vegetatively propagated crop to a diploid F1 hybrid, propagated via botanical seed, would constitute a considerable advance for global agriculture, but faces multiple challenges. One such challenge is the difficulty in inbreeding potato, which involves purging deleterious alleles from its genome. This commentary discusses possible reasons for this difficulty and highlights a recent sequence-based effort to classify SNP variation, in potato germplasm, according to its deleterious potential. Tools and strategies connected to this database may facilitate development of F1 hybrids.

Mating compatibility and offspring traits evaluation among different strains of Tenebrio molitor

The fast-growing sector of insects for food and feed stimulates researchers and the industry to explore uncharted territories, such as insect breeding, to improve economically important insect fitness traits. The yellow mealworm, Tenebrio molitor L., is one of the most thoroughly studied insect species as food and feed. However, data on mating compatibility and the effect of cross-breeding between different strains on the performance and fitness of the hybrids are scarce. In the present study, we comparatively evaluated the mating compatibility between two T. molitor inbred strains (Greek and Italian) and their outbred strains, i.e., Italian (♀)-Greek (♂) and Greek (♀)-Italian (♂), as well as the performance of their hybrid offspring. Based on the results, there is good mating compatibility among adults of the strains tested. Offspring performance, quantified as larval survival and final larval weight, followed a similar pattern among the crossings examined. Even though differences were insignificant, the outbreeding of females of the Italian strain resulted in a higher cumulative number of eggs and hatching rate and higher offspring weight. The present study sheds light on the mating compatibility of different T. molitor strains and their hybrids' economically important life table characteristics to take the mass production of insects one step further.

Plant biochemical genetics in the multiomics era

Our understanding of plant biology has been revolutionized by modern genetics and biochemistry. However, biochemical genetics can be traced back to the foundation of Mendelian genetics; indeed, one of Mendel's milestone discoveries of seven characteristics of pea plants later came to be ascribed to a mutation in a starch branching enzyme. Here, we review both current and historical strategies for the elucidation of plant metabolic pathways and the genes that encode their component enzymes and regulators. We use this historical review to discuss a range of classical genetic phenomena including epistasis, canalization, and heterosis as viewed through the lens of contemporary high-throughput data obtained via the array of approaches currently adopted in multiomics studies.

Maximization through optimization? On the relationship between hybrid performance and parental genetic distance

Heterosis is the improved performance of hybrids compared with their parental components and is widely exploited in agriculture. According to quantitative genetic theory, genetic distance between parents at heterotic quantitative trait loci is required for heterosis, but how heterosis varies with genetic distance has remained elusive, despite intensive research on the topic. Experimental studies have often found a positive association between heterosis and genetic distance that, however, varied in strength. Most importantly, it has remained unclear whether heterosis increases continuously with genetic distance or whether there is an optimum genetic distance after which heterosis declines again. Here, we revisit the relationship between heterosis and genetic distance and provide perspectives on how to maximize heterosis and hybrid performance in breeding, as well as the consequences for the design of heterotic groups and the utilization of more exotic material and genetic resources.

Environmentally independent selection for hybrids between divergent freshwater stickleback lineages in semi-natural ponds

Hybridization following secondary contact of genetically divergent populations can influence the range expansion of invasive species, though specific outcomes depend on the environmental dependence of hybrid fitness. Here, using two genetically and ecologically divergent threespine stickleback lineages that differ in their history of freshwater colonization, we estimate fitness variation of parental lineages and hybrids in semi-natural freshwater ponds with contrasting histories of nutrient loading. In our experiment, we found that fish from the older freshwater lineage (Lake Geneva) and hybrids outperformed fish from the younger freshwater lineage (Lake Constance) in terms of both growth and survival, regardless of the environmental context of our ponds. Across all ponds, hybrids exhibited the highest survival. Although wild-caught adult populations differed in their functional and defence morphology, it is unclear which of these traits underlie the fitness differences observed among juveniles in our experiment. Overall, our work suggests that when hybrid fitness is insensitive to environmental conditions, as observed here, introgression may promote population expansion into unoccupied habitats and accelerate invasion success.

Fine Mapping of fw6.3, a Major-Effect Quantitative Trait Locus That Controls Fruit Weight in Tomato

Tomato (Solanum lycopersicum) is a widely consumed vegetable, and the tomato fruit weight is a key yield component. Many quantitative trait loci (QTLs) controlling tomato fruit weight have been identified, and six of them have been fine-mapped and cloned. Here, four loci controlling tomato fruit weight were identified in an F₂ population through QTL seq.; fruit weight 6.3 (fw6.3) was a major-effect QTL and its percentage of variation explanation (R²) was 0.118. This QTL was fine-mapped to a 62.6 kb interval on chromosome 6. According to the annotated tomato genome (version SL4.0, annotation ITAG4.0), this interval contained seven genes, including Solyc06g074350 (the SELF-PRUNING gene), which was likely the candidate gene underlying variation in fruit weight. The SELF-PRUNING gene contained a single-nucleotide polymorphism that resulted in an amino acid substitution in the protein sequence. The large-fruit allele of fw6.3 (fw6.3_HG) was overdominant to the small-fruit allele fw6.3_RG. The soluble solids content was also increased by fw6.3_HG. These findings provide valuable information that will aid the cloning of the FW6.3 gene and ongoing efforts to breed tomato plants with higher yield and quality via molecular marker-assisted selection.

Gene expression variations and allele-specific expression of two rice and their hybrid in caryopses at single-nucleus resolution

Seeds are an indispensable part of the flowering plant life cycle and a critical determinant of agricultural production. Distinct differences in the anatomy and morphology of seeds separate monocots and dicots. Although some progress has been made with respect to understanding seed development in Arabidopsis, the transcriptomic features of monocotyledon seeds at the cellular level are much less understood. Since most important cereal crops, such as rice, maize, and wheat, are monocots, it is essential to study transcriptional differentiation and heterogeneity during seed development at a finer scale. Here, we present single-nucleus RNA sequencing (snRNA-seq) results of over three thousand nuclei from caryopses of the rice cultivars Nipponbare and 9311 and their intersubspecies F₁ hybrid. A transcriptomics atlas that covers most of the cell types present during the early developmental stage of rice caryopses was successfully constructed. Additionally, novel specific marker genes were identified for each nuclear cluster in the rice caryopsis. Moreover, with a focus on rice endosperm, the differentiation trajectory of endosperm subclusters was reconstructed to reveal the developmental process. Allele-specific expression (ASE) profiling in endosperm revealed 345 genes with ASE (ASEGs). Further pairwise comparisons of the differentially expressed genes (DEGs) in each endosperm cluster among the three rice samples demonstrated transcriptional divergence. Our research reveals differentiation in rice caryopsis from the single-nucleus perspective and provides valuable resources to facilitate clarification of the molecular mechanism underlying caryopsis development in rice and other monocots.

Performances of three-way cross hybrids over their respective single crosses and related heterosis of maize (Zea mays L.) hybrids evaluated in Ethiopia

Less attention had been given to the performances of three-way crosses and its comparative advantages of these hybrids over single crosses. This study was carried out to evaluate the performances of three-way crosses in comparison to single crosses for yield and related agronomic traits and to estimate the magnitude of heterosis. The trial was laid out in a simple alpha lattice design of 10 × 6 for lines, 6 × 5 for single crosses (SC), and 9 × 5 for three way-crosses and planted in adjacent plots in the 2019 cropping season in three locations namely Ambo, Abala-Farcha and Melkassa. Single cross hybrids showed a highly significant (P<1%) variation for grain yield, plant height, ear height, and ear length at three locations. These single cross hybrids had showed also a highly significant genotype by environment interaction (P < 1%) for grain yield, plant height, ear height and kernel per ear. Regarding three-way crosses, there was a significant variation (P<5%) on grain yield in Ambo and Melkassa but on ear height and rows per ear in Abala-Faracho. The genotype × environment interaction was significantly varied for grain yield, ear height and ear length. In the comparison, 80% crosses in Ambo, 73% in Abala-Faracho and 67% in Melkassa showed that three-way crosses were better in their performance than that of their respective single crosses. On the other hand, the single crosses that out-performed their respective three-way crosses were higher in Melkassa than Abala-Faracho and the least were reported from Ambo. Similarly, the maximum better and mid-parent heterosis was from single cross 1(769%) in Ambo and single cross 7 (104%) in Melkassa whereas TWC 14 (52%) and TWC 24 (78%) were the highest better and mid-parent heterosis, respectively in Ambo, TWC1 (56%), and TWC30 (25%) were the highest BPH, and MPH, respectively in Melkassa.

Correlation between Parental Transcriptome and Field Data for the Characterization of Heterosis in Chinese Cabbage

In Chinese cabbage breeding, hybrids have made a terrific contribution due to heterosis, the superior performance of offspring compared to their inbred parents. Since the development of new, top-performing hybrids requires a large scale of human and material resources, the prediction of hybrid performance is of utmost interest to plant breeders. In our research, leaf transcriptome data from eight parents were used to investigate if they might be employed as markers to predict hybrid performance and heterosis. In Chinese cabbage, heterosis of plant growth weight (PGW) and heterosis of head weight (HW) were more obvious than other traits. The number of differential expression genes (DEGs) between parents was related to the PGW, length of the biggest outer leaf (LOL), leaf head height (LHH), leaf head width (LHW), HW, leaf number of head (LNH) and plant height (PH) of hybrids, and up-regulated DEGs number was also associated with these traits. Euclidean and binary distances of parental gene expression levels were significantly correlated with the PGW, LOL, LHH, LHW, HW and PH of hybrids. Additionally, there was a significant correlation between the parental expression levels of multiple genes involved in the ribosomal metabolic pathway and hybrid observations and heterosis in PGW, with the BrRPL23A gene showing the highest correlation with the MPH of PGW(r = 0.75). Therefore, leaf transcriptome data can preliminarily predict the hybrid performance and select parents in Chinese cabbage.

SiMYBS3, Encoding a Setaria italica Heterosis-Related MYB Transcription Factor, Confers Drought Tolerance in Arabidopsis

Drought is a major limiting factor affecting grain production. Drought-tolerant crop varieties are required to ensure future grain production. Here, 5597 DEGs were identified using transcriptome data before and after drought stress in foxtail millet (Setaria italica) hybrid Zhangza 19 and its parents. A total of 607 drought-tolerant genes were screened through WGCNA, and 286 heterotic genes were screened according to the expression level. Among them, 18 genes overlapped. One gene, Seita.9G321800, encoded MYBS3 transcription factor and showed upregulated expression after drought stress. It is highly homologous with MYBS3 in maize, rice, and sorghum and was named SiMYBS3. Subcellular localization analysis showed that the SiMYBS3 protein was located in the nucleus and cytoplasm, and transactivation assay showed SiMYBS3 had transcriptional activation activity in yeast cells. Overexpression of SiMYBS3 in Arabidopsis thaliana conferred drought tolerance, insensitivity to ABA, and earlier flowering. Our results demonstrate that SiMYBS3 is a drought-related heterotic gene and it can be used for enhancing drought resistance in agricultural crop breeding.

Dynamic patterns of gene expression and regulatory variation in the maize seed coat

BACKGROUND

Seed size is an important factor contributing to maize yield, but its molecular mechanism remains unclear. The seed coat, which serves as one of the three components of the maize grain, determines seed size to a certain extent. The seed coat also shares the maternal genotype and is an ideal material for studying heterosis.

RESULTS

In this study, the self-pollinated seeds of the maize hybrid Yudan888 and its parental lines were continuously collected from 0 day after pollination (DAP) to 15 DAP for phenotyping, cytological observation and RNA-seq. The phenotypic data showed that 3 DAP and 8 DAP are the best time points to study maize seed coat heterosis. Cytological observations indicated that maize seed coat heterosis might be the result of the coordination between cell number and cell size. Furthermore, the RNA-seq results showed that the nonadditive genes changed significantly between 3 and 8 DAP. However, the number of genes expressed additively was not significantly different. Our findings suggest that seed coat heterosis in hybrid is the result of nonadditive expression caused by dynamic changes in genes at different time points during seed expansion and seed coat development. Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment indicated that genes related to DNA replication, cell cycle regulation, circadian rhythms and metabolite accumulation contributed significantly to hybrid seed coat heterosis.

CONCLUSION

Maize seed coat phenotyping allowed us to infer that 3 DAP and 8 DAP are important time points in the study of seed coat heterosis. Our findings provide evidence for genes involved in DNA replication, cell cycle regulation, circadian rhythms and metabolite accumulation in hybrid with high or low parental expression as major contributors to hybrid seed coat heterosis.

Combined transcriptome and metabolome analysis reveals the effects of light quality on maize hybrids

BACKGROUND

Heterosis, or hybrid vigor, refers to the phenotypic superiority of an F₁ hybrid relative to its parents in terms of growth rate, biomass production, grain yield, and stress tolerance. Light is an energy source and main environmental cue with marked impacts on heterosis in plants. Research into the production applications and mechanism of heterosis has been conducted for over a century and a half, but little is known about the effect of light on plant heterosis.

RESULTS

In this study, an integrated transcriptome and metabolome analysis was performed using maize (Zea mays L.) inbred parents, B73 and Mo17, and their hybrids, B73 × Mo17 (BM) and Mo17 × B73 (MB), grown in darkness or under far-red, red, or blue light. Most differentially expressed genes (73.72-92.50%) and differentially accumulated metabolites (84.74-94.32%) exhibited non-additive effects in BM and MB hybrids. Gene Ontology analysis revealed that differential genes and metabolites were involved in glutathione transfer, carbohydrate transport, terpenoid biosynthesis, and photosynthesis. The darkness, far-red, red, and blue light treatments were all associated with phenylpropanoid-flavonoid biosynthesis by Weighted Gene Co-expression Network Analysis and Kyoto Encyclopedia of Genes and Genomes enrichment analysis. Five genes and seven metabolites related to phenylpropanoid-flavonoid biosynthesis pathway were identified as potential contributors to the interactions between maize heterosis and light conditions. Consistent with the strong mid-parent heterosis observed for metabolites, significant increases in both fresh and dry weights were found in the MB and BM hybrids compared with their inbred parents. Unexpectedly, increasing light intensity resulted in higher biomass heterosis in MB, but lower biomass heterosis in BM.

CONCLUSIONS

The transcriptomic and metabolomic results provide unique insights into the effects of light quality on gene expression patterns and genotype-environment interactions, and have implications for gene mining of heterotic loci to improve maize production.

Loci Associated with Negative Heterosis for Viability and Meat Productivity in Interspecific Sheep Hybrids

Negative heterosis can occur on different economically important traits, but the exact biological mechanisms of this phenomenon are still unknown. The present study focuses on determining the genetic factors associated with negative heterosis in interspecific hybrids between domestic sheep (Ovis aries) and argali (Ovis ammon). One locus (rs417431015) associated with viability and two loci (rs413302370, rs402808951) associated with meat productivity were identified. One gene (ARAP2) was prioritized for viability and three for meat productivity (PDE2A, ARAP1, and PCDH15). The loci associated with meat productivity were demonstrated to fit the overdominant inheritance model and could potentially be involved int negative heterosis mechanisms.

Large-scale genomic and transcriptomic profiles of rice hybrids reveal a core mechanism underlying heterosis

BACKGROUND

Heterosis is widely used in agriculture. However, its molecular mechanisms are still unclear in plants. Here, we develop, sequence, and record the phenotypes of 418 hybrids from crosses between two testers and 265 rice varieties from a mini-core collection.

RESULTS

Phenotypic analysis shows that heterosis is dependent on genetic backgrounds and environments. By genome-wide association study of 418 hybrids and their parents, we find that nonadditive QTLs are the main genetic contributors to heterosis. We show that nonadditive QTLs are more sensitive to the genetic background and environment than additive ones. Further simulations and experimental analysis support a novel mechanism, homo-insufficiency under insufficient background (HoIIB), underlying heterosis. We propose heterosis in most cases is not due to heterozygote advantage but homozygote disadvantage under the insufficient genetic background.

CONCLUSION

The HoIIB model elucidates that genetic background insufficiency is the intrinsic mechanism of background dependence, and also the core mechanism of nonadditive effects and heterosis. This model can explain most known hypotheses and phenomena about heterosis, and thus provides a novel theory for hybrid rice breeding in future.

The effect of experimental hybridization on cognition and brain anatomy: Limited phenotypic variation and transgression in Poeciliidae

Hybridization can promote phenotypic variation and often produces trait combinations distinct from the parental species. This increase in available variation can lead to the manifestation of functional novelty when new phenotypes bear adaptive value under the environmental conditions in which they occur. Although the role of hybridization as a driver of variation and novelty in traits linked to fitness is well recognized, it remains largely unknown whether hybridization can fuel behavioral novelty by promoting phenotypic variation in brain morphology and/or cognitive traits. To address this question, we investigated the effect of hybridization on brain anatomy, learning ability, and cognitive flexibility in first- and second-generation hybrids of two closely related fish species (Poecilia reticulata and Poecilia wingei). Overall, we found that F1 and F2 hybrids showed intermediate brain morphology and cognitive traits compared to parental groups. Moreover, as phenotypic dispersion and transgression were low for both brain and cognitive traits, we suggest that hybridization is not a strong driver of brain anatomical and cognitive diversification in these Poeciliidae. To determine the generality of this conclusion, hybridization experiments with cognitive tests need to be repeated in other families.

Morphological and Physiological Responses of Hybrid Aspen (Populus tremuloides Michx. × Populus tremula L.) Clones to Light In Vitro

Micropropagation of fast-growing tree genotypes such as the hybrid aspen (Populus tremuloides Michx. × Populus tremula L.) is increasing. The efficiency of micropropagation depends on the luminaires, hence luminescent electric diodes (LED), which emit light of a narrow spectrum, are gaining popularity. Mostly, different LEDs are combined to increase the photosynthetic efficiency. However, light also acts as an environmental signal, which triggers specific responses in plants, which are genotype specific, and regarding hybrid aspen, are likely affected by heterosis. In this study, morphological and physiological responses of clones of hybrid aspen with contrasting field performance to the spectral composition of illumination were studied in vitro. Among the 15 variables measured, area of leaves and concentration and ratio of chlorophyll a and b explained most of the variance (58.6%), thereby linking a specific combination of traits to productivity. These traits and their responses to light were affected by heterosis, as indicated by the clone-treatment interaction, particularly for the clone's moderate productivity. The top-performing clones were little sensitive to illumination due to efficient photosystems. Nevertheless, illumination with wider spectral composition had generally positive effects on plantlet performance. Accordingly, clone-specific illumination protocols and luminaries capable of it are advantageous for the efficiency of micropropagation of hybrid aspen.

The repetitive DNA sequence landscape and DNA methylation in chromosomes of an apomictic tropical forage grass, Cenchrus ciliaris

Cenchrus ciliaris is an apomictic, allotetraploid pasture grass widely distributed in the tropical and subtropical regions of Africa and Asia. In this study, we aimed to investigate the genomic organization and characterize some of the repetitive DNA sequences in this species. Due to the apomictic propagation, various aneuploid genotypes are found, and here, we analyzed a 2n = 4x + 3 = 39 accession. The physical mapping of Ty1-copia and Ty3-gypsy retroelements through fluorescence in situ hybridization with a global assessment of 5-methylcytosine DNA methylation through immunostaining revealed the genome-wide distribution pattern of retroelements and their association with DNA methylation. Approximately one-third of Ty1-copia sites overlapped or spanned centromeric DAPI-positive heterochromatin, while the centromeric regions and arms of some chromosomes were labeled with Ty3-gypsy. Most of the retroelement sites overlapped with 5-methylcytosine signals, except for some Ty3-gypsy on the arms of chromosomes, which did not overlap with anti-5-mC signals. Universal retrotransposon probes did not distinguish genomes of C. ciliaris showing signals in pericentromeric regions of all 39 chromosomes, unlike highly abundant repetitive DNA motifs found in survey genome sequences of C. ciliaris using graph-based clustering. The probes developed from RepeatExplorer clusters gave strong in situ hybridization signals, mostly in pericentromeric regions of about half of the chromosomes, and we suggested that they differentiate the two ancestral genomes in the allotetraploid C. ciliaris, likely having different repeat sequence variants amplified before the genomes came together in the tetraploid.

Genetic diversity analysis and development of molecular markers for the identification of largemouth bass (Micropterus salmoides L.) based on whole-genome re-sequencing

Largemouth bass (Micropterus salmoides L.) is generally considered to comprise two subspecies, Florida bass (M. floridanus) and Northern Largemouth bass (M. salmoides), which have biological characteristic differences because of their geographical distribution. In this study, whole-genome re-sequencing was performed among 10 Florida and 10 Northern largemouth bass, respectively. In total, 999,793 SNPs and 227,797 InDels were finally identified, and 507,401 SNPs (50.75%) and 116,213 InDels (51.01%) were successfully mapped to annotated 18,629 genes and 14,060 genes, respectively. KEGG classification indicated that most of these genes were focused on the pathways including signal transduction, transport and catabolism, and endocrine system. Genetic diversity analysis indicated that Florida largemouth bass had higher genetic diversity than Northern largemouth bass, indicating that the germplasm quality of Northern largemouth bass had markedly reduced in China. To examine the accuracies of the identified markers, 23 SNPs and eight InDels (the insertions or deletions of more than 45 bp) were randomly selected and detected among Florida largemouth bass, Northern largemouth bass, and their F1 hybrids. The detection efficiencies of all the markers were higher than 95%; nineteen SNPs and three InDels could accurately distinguish the two subspecies and their F1 hybrids with 100% efficiencies. Moreover, the three InDel markers could clearly distinguish the two subspecies and their F1 hybrids with a PCR-based agarose gel electrophoresis. In conclusion, our study established a simple PCR-based method for the germplasm identification of largemouth bass, which will be useful in the germplasm protection, management, and hybridization breeding of largemouth bass.

Diversity and population structure of Nordic potato cultivars and breeding clones

BACKGROUND

The genetic diversity and population structure of breeding germplasm is central knowledge for crop improvement. To gain insight into the genetic potential of the germplasm used for potato breeding in a Nordic breeding program as well as all available accessions from the Nordic genebank (NordGen), 133 potato genotypes were genotyped using the Infinium Illumina 20 K SNP array. After SNP filtering, 11 610 polymorphic SNPs were included in the analysis. In addition, data from three important breeding traits - percent dry matter and uniformity of tuber shape and eye - were scored to measure the variation potato cultivars and breeding clones.

RESULTS

The genetic diversity among the genotypes was estimated using principal coordinate analysis based on the genetic distance between individuals, as well as by using the software STRUCTURE. Both methods suggest that the collected breeding material and the germplasm from the gene-bank are closely related, with a low degree of population structure between the groups. The phenotypic distribution among the genotypes revealed significant differences, especially between farmer's cultivars and released cultivars and breeding clones. The percent heterozygosity was similar between the groups, with a mean average of 58-60%. Overall, the breeding germplasm and the accessions from the Nordic genebank seems to be closely related with similar genetic background.

CONCLUSION

The genetic potential of available Nordic potato breeding germplasm is low, and for genetic hybridization purposes, genotypes from outside the Nordic region should be employed.

Comparative transcriptomic analysis of maize ear heterosis during the inflorescence meristem differentiation stage

BACKGROUND

Heterosis is widely used in many crops and is important for global food safety, and maize is one of the most successful crops to take advantage of heterosis. Gene expression patterns control the development of the maize ear, but the mechanisms by which heterosis affects transcriptional-level control are not fully understood.

RESULTS

In this study, we sampled ear inflorescence meristems (IMs) from the single-segment substitution maize (Zea mays) line lx9801^hlEW2b, which contains the heterotic locus hlEW2b associated with ear width, as well as the receptor parent lx9801, the test parent Zheng58, and their corresponding hybrids Zheng58 × lx9801^hlEW2b (HY) and Zheng58 × lx9801 (CK). After RNA sequencing and transcriptomic analysis, 2531 unique differentially expressed genes (DEGs) were identified between the two hybrids (HY vs. CK). Our results showed that approximately 64% and 48% of DEGs exhibited additive expression in HY and CK, whereas the other genes displayed a non-additive expression pattern. The DEGs were significantly enriched in GO functional categories of multiple metabolic processes, plant organ morphogenesis, and hormone regulation. These essential processes are potentially associated with heterosis performance during the maize ear developmental stage. In particular, 125 and 100 DEGs from hybrids with allele-specific expression (ASE) were specifically identified in HY and CK, respectively. Comparison between the two hybrids suggested that ASE genes were involved in different development-related processes that may lead to the hybrid vigor phenotype during maize ear development. In addition, several critical genes involved in auxin metabolism and IM development were differentially expressed between the hybrids and showed various expression patterns (additive, non-additive, and ASE). Changes in the expression levels of these genes may lead to differences in auxin homeostasis in the IM, affecting the transcription of core genes such as WUS that control IM development.

CONCLUSIONS

Our research suggests that additive, non-additive, and allele-specific expression patterns may fine-tune the expression of crucial DEGs that modulate carbohydrate and protein metabolic processes, nitrogen assimilation, and auxin metabolism to optimal levels, and these transcriptional changes may play important roles in maize ear heterosis. The results provide new information that increases our understanding of the relationship between transcriptional variation and heterosis during maize ear development, which may be helpful for clarifying the genetic and molecular mechanisms of heterosis.

Heterosis of fitness and phenotypic variance in the evolution of a diploid gene regulatory network

Heterosis describes the phenomenon, whereby a hybrid population has higher fitness than an inbred population, which has previously been explained by either Mendelian dominance or overdominance under the general assumption of a simple genotype-phenotype relationship. However, recent studies have demonstrated that genes interact through a complex gene regulatory network (GRN). Furthermore, phenotypic variance is reportedly lower for heterozygotes, and the origin of such variance-related heterosis remains elusive. Therefore, a theoretical analysis linking heterosis to GRN evolution and stochastic gene expression dynamics is required. Here, we investigated heterosis related to fitness and phenotypic variance in a system with interacting genes by numerically evolving diploid GRNs. According to the results, the heterozygote population exhibited higher fitness than the homozygote population, indicating fitness-related heterosis resulting from evolution. In addition, the heterozygote population exhibited lower noise-related phenotypic variance in expression levels than the homozygous population, implying that the heterozygote population is more robust to noise. Furthermore, the distribution of the ratio of heterozygote phenotypic variance to homozygote phenotypic variance exhibited quantitative similarity with previous experimental results. By applying dominance and differential gene expression rather than only a single gene expression model, we confirmed the correlation between heterosis and differential gene expression. We explain our results by proposing that the convex high-fitness region is evolutionarily shaped in the genetic space to gain noise robustness under genetic mixing through sexual reproduction. These results provide new insights into the effects of GRNs on variance-related heterosis and differential gene expression.

RNA-Seq Analysis Reveals Expression Regulatory Divergence of W-Linked Genes between Two Contrasting Chicken Breeds

Simple Summary

Understanding the mode of gene expression and regulation is essential for understanding the evolutionary process. Many previous studies tried to explain regulatory changes at the autosomal level, but little research has extended these explorations to the field of sex chromosomes due to their complex sex-limit features. Here, we first adopted an innovative method of identifying regulatory divergence of W-linked genes. Compared with cis-regulatory divergence, trans acting genes were more extensive in the W chromosome. We also found that divergent sex specific selection cannot strongly affect the expression evolution of the W chromosome. This insensitivity to selection may be one of the reasons why regulatory divergence is so small between autosomal and sex chromosomes.

Abstract

The regulation of gene expression is a complex process involving organism function and phenotypic diversity, and is caused by cis- and trans- regulation. While prior studies identified the regulatory pattern of the autosome rewiring in hybrids, the role of gene regulation in W sex chromosomes is not clear due to their degradation and sex-limit expression. Here, we developed reciprocal crosses of two chicken breeds, White Leghorn and Cornish Game, which exhibited broad differences in gender-related traits, and assessed the expression of the genes on the W chromosome to disentangle the contribution of cis- and trans-factors to expression divergence. We found that female-specific selection does not have a significant effect on W chromosome gene-expression patterns. For different tissues, there were most parental divergence expression genes in muscle, and also more heterosis compared with two other tissues. Notably, a broader pattern of trans regulation in the W chromosome was observed, which is consistent with autosomes. Taken together, this work describes the regulatory divergence of W-linked genes between two contrasting breeds and indicates sex chromosomes have a unique regulation and expression mechanism.

Mapping QTLs for spring green-up, plant vigor, and plant biomass in two lowland switchgrass populations

Switchgrass (Panicum virgatum L.) is an important perennial C4 species due to its large potential for cellulosic bioenergy feedstock production. Identification of quantitative trait loci (QTL) controlling important developmental traits is valuable to understanding the genetic basis and using marker-assisted selection (MAS) in switchgrass breeding. One F1 hybrid population derived from NL94 (♀) × SL93 (♂) and one S1 (first-generation selfed) population from NL94 were used in this study. Both the populations showed significant variations for genotype and genotype by environment interactions for three traits studied: plant vigor, spring green-up, and plant biomass. Plant vigor had strong and positive correlations with plant biomass in both populations. Broad-sense heritability estimates for plant vigor ranged from 0.46 to 0.74 and 0.45 to 0.74 in the hybrid and selfed population, respectively. Spring green-up had similar heritability estimates, 0.42-0.78 in the hybrid population, and 0.47-0.82 in the selfed population. Heritability of plant biomass was 0.54-0.64 in the hybrid population and 0.64-0.74 in the selfed population. Fifteen QTLs for spring green-up, 6 QTLs for plant vigor, and 3 QTLs for biomass yield were detected in the hybrid population, whereas 4 QTLs for spring green-up, 4 QTLs for plant vigor, and 1 QTL for biomass yield were detected in the selfed population. Markers associated with these QTLs can be used in MAS to accelerate switchgrass breeding program. This study provided new information in understanding the genetic control of biomass components and demonstrated substantial heterotic vigor that could be explored for breeding hybrid cultivars in switchgrass.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01296-7.

Quantitative Trait Loci Analysis Based on High-Density Mapping of Single-Nucleotide Polymorphisms by Genotyping-by-Sequencing Against Pine Wilt Disease in Japanese Black Pine (Pinus thunbergii)

Identifying genes/loci for resistance to pine wilt disease (PWD) caused by the pine wood nematode (PWN) is beneficial for improving resistance breeding in Pinus thunbergii, but to date, genetic information using molecular markers has been limited. Here, we constructed a high-density linkage map using genotyping-by-sequencing (GBS) and conducted quantitative trait loci (QTL) analysis for PWD resistance for the self-pollinated progeny of "Namikata 73," which is the most resistant variety among resistant varieties of P. thunbergii, following inoculation tests with PWN. An S1 mapping population consisting of the 116 progenies derived from self-pollination of the resistant variety, "Namikata 73" (resistance rank 5 to PWN), was inoculated with PWN isolate Ka-4 and evaluated for disease symptoms. To construct a high-density linkage map, we used single-nucleotide polymorphisms (SNPs) identified by GBS based on next-generation sequencing technology and some anchor DNA markers, expressed sequence tag (EST)-derived SNP markers and EST-derived simple sequence repeat (SSR) markers, and genomic SSR markers. The linkage map had 13 linkage groups (LGs) consisting of 2,365 markers including 2,243 GBS-SNP markers over a total map distance of 1968.4 centimorgans (cM). Results from QTL analysis using phenotype data and the linkage map indicated that PWD resistance is controlled by a single locus located on LG-3, as identified in a previous study. This locus showed overdominant genetic action in the present study. With the confirmation of PWD1 in two different mapping populations (present study and a previous study), the locus associated with this region is thought to be a good target for marker-assisted selection in P. thunbergii breeding programs in order to obtain high levels of resistance to PWD caused by PWN.

Identification of an anther-specific promoter from a male sterile AB line in Chinese cabbage (Brassica rapa L. ssp. pekinensis)

The promoter of the male sterile gene is important for studying male sterility. In this study, BraA08g014780.3C which differentially expressed between male sterile and fertile plants was identified from a genetic male sterile AB line of Chinese cabbage by RNA-seq. qRT-PCR revealed that BraA08g014780.3C was mainly expressed in the early stage of floral bud development in fertile plants, and preferentially expressed in their anthers. The promoter of BraA08g014780.3C was cloned and analyzed. Cis acting element analysis showed that the promoter of BraA08g014780.3C contains POLLEN1LELAT52 and GTGANTG10, which are both pollen-specific expression elements. The BraA08g014780.3Cp::GUS vector was constructed, then transformed to Arabidopsis thaliana Col-0. PCR analysis and sequencing of the transgenic Arabidopsis revealed that the GUS gene driven by the BraA08g014780.3C promoter was successfully transformed to the Arabidopsis. GUS staining indicated that the promoter of BraA08g014780.3C was an anther-specific promoter. Identifying the anther-specific promoter laid a foundation for revealing BraA08g014780.3C function in male sterility of Chinese cabbage.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03160-z.

Assembly of yield heterosis of an elite rice hybrid is promising by manipulating dominant quantitative trait loci

In the past, rice hybrids with strong heterosis have been obtained empirically, by developing and testing thousands of combinations. Here, we aimed to determine whether heterosis of an elite hybrid could be achieved by manipulating major quantitative trait loci. We used 202 chromosome segment substitution lines from the elite hybrid Shanyou 63 to evaluate single segment heterosis (SSH) of yield per plant and identify heterotic loci. All nine detected heterotic loci acted in a dominant fashion, and no SSH exhibited overdominance. Functional alleles of key yield-related genes Ghd7, Ghd7.1, Hd1, and GS3 were dispersed in both parents. No functional alleles of three investigated genes were expressed at higher levels in the hybrids than in the more desirable parents. A hybrid pyramiding eight heterotic loci in the female parent Zhenshan 97 background had a comparable yield to Shanyou 63 and much higher yield than Zhenshan 97. Five hybrids pyramiding eight or nine heterotic loci in the combined parental genome background showed similar yield performance to that of Shanyou 63. These results suggest that dominance underlying functional complementation is an important contributor to yield heterosis and that heterosis assembly might be successfully promised by manipulating several major dominant heterotic loci.

Advances in lentil production through heterosis: Evaluating generations and breeding systems

Heterosis is defined as increased performance of the F1 hybrid relative to its parents. In the current study, a cohort of populations and parents were created to evaluate and understand heterosis across generations (i.e., F1 to F3) in lentil, a self-pollinated annual diploid (2n = 2× = 14) crop species. Lentil plants were evaluated for heterotic traits in terms of plant height, biomass fresh weight, seed number, yield per plant and 100 grain weight. A total of 47 selected lentil genotypes were cross hybridized to generate 72 F1 hybrids. The F1 hybrids from the top five crosses exhibited between 31%-62% heterosis for seed number with reference to the better parent. The five best performing heterotic crosses were selected with a negative control for evaluation at the subsequent F2 generation and only the tails of the distribution taken forward to be assessed in the F3 generation as a sub selection. Overall, heterosis decreases across the subsequent generations for all traits studied. However, some individual genotypes were identified at the F2 and sub-selected F3 generations with higher levels of heterosis than the best F1 mean value (hybrid mimics). The phenotypic data for the selected F2 and sub selected F3 hybrids were analysed, and the study suggested that 100 grain weight was the biggest driver of yield followed by seed number. A genetic diversity analysis of all the F1 parents failed to correlate genetic distance and divergence among parents with heterotic F1's. Therefore, genetic distance was not a key factor to determine heterosis in lentil. The study highlights the challenges associated with different breeding systems for heterosis (i.e., F1 hybrid-based breeding systems and/or via hybrid mimics) but demonstrates the potential significant gains that could be achieved in lentil productivity.

Environment-Dependent Heterosis and Transgressive Gene Expression in Reciprocal Hybrids between the Channel Catfish Ictalurus punctatus and the Blue Catfish Ictalurus furcatus

The hybrid between female channel catfish (Ictalurus punctatus) and male blue catfish (Ictalurus furcatus) is superior in feed conversion, disease resistance, carcass yield, and harvestability compared to both parental species. However, heterosis and heterobeltiosis only occur in pond culture, and channel catfish grow much faster than the other genetic types in small culture units. This environment-dependent heterosis is intriguing, but the underlying genetic mechanisms are not well understood. In this study, phenotypic characterization and transcriptomic analyses were performed in the channel catfish, blue catfish, and their reciprocal F₁s reared in tanks. The results showed that the channel catfish is superior in growth-related morphometrics, presumably due to significantly lower innate immune function, as investigated by reduced lysozyme activity and alternative complement activity. RNA-seq analysis revealed that genes involved in fatty acid metabolism/transport are significantly upregulated in channel catfish compared to blue catfish and hybrids, which also contributes to the growth phenotype. Interestingly, hybrids have a 40-80% elevation in blood glucose than the parental species, which can be explained by a phenomenon called transgressive expression (overexpression/underexpression in F₁s than the parental species). A total of 1140 transgressive genes were identified in F₁ hybrids, indicating that 8.5% of the transcriptome displayed transgressive expression. Transgressive genes upregulated in F₁s are enriched for glycan degradation function, directly related to the increase in blood glucose level. This study is the first to explore molecular mechanisms of environment-dependent heterosis/heterobeltiosis in a vertebrate species and sheds light on the regulation and evolution of heterosis vs. hybrid incompatibility.