[DNA methylation in cotton hybrids and their parents]

Exploitation of epigenetic variation of crop wild relatives for crop improvement and agrobiodiversity preservation

Crop wild relatives (CWRs) are recognized as the best potential source of traits for crop improvement. However, successful crop improvement using CWR relies on identifying variation in genes controlling desired traits in plant germplasms and subsequently incorporating them into cultivars. Epigenetic diversity may provide an additional layer of variation within CWR and can contribute novel epialleles for key traits for crop improvement. There is emerging evidence that epigenetic variants of functional and/or agronomic importance exist in CWR gene pools. This provides a rationale for the conservation of epigenotypes of interest, thus contributing to agrobiodiversity preservation through conservation and (epi)genetic monitoring. Concepts and techniques of classical and modern breeding should consider integrating recent progress in epigenetics, initially by identifying their association with phenotypic variations and then by assessing their heritability and stability in subsequent generations. New tools available for epigenomic analysis offer the opportunity to capture epigenetic variation and integrate it into advanced (epi)breeding programmes. Advances in -omics have provided new insights into the sources and inheritance of epigenetic variation and enabled the efficient introduction of epi-traits from CWR into crops using epigenetic molecular markers, such as epiQTLs.

Transcriptome shock in interspecific F1 allotriploid hybrids between Brassica species

Interspecific hybridization drives the evolution of angiosperms and can be used to introduce novel alleles for important traits or to activate heterosis in crop breeding. Hybridization brings together gene expression networks from two different species, potentially causing global alterations of gene expression in the F1 plants which is called 'transcriptome shock'. Here, we explored such a transcriptome shock in allotriploid Brassica hybrids. We generated interspecific F1 allotriploid hybrids between the allotetraploid species Brassica napus and three accessions of the diploid species Brassica rapa. RNA-seq of the F1 hybrids and the parental plants revealed that 26.34-30.89% of genes were differentially expressed between the parents. We also analyzed expression level dominance and homoeolog expression bias between the parents and the F1 hybrids. The expression-level dominance biases of the Ar, An, and Cn subgenomes was genotype and stage dependent, whereas significant homoeolog expression bias was observed among three subgenomes from different parents. Furthermore, more genes were involved in trans regulation than in cis regulation in allotriploid F1 hybrids. Our findings provide new insights into the transcriptomic responses of cross-species hybrids and hybrids showing heterosis, as well as a new method for promoting the breeding of desirable traits in polyploid Brassica species.

Deltapine 15 contributes to the genomic architecture of modern upland cotton cultivars

Foundation parents play a critical role in the genetic constituents of the derived genotypes. Deltapine-15 (DLP-15), introduced to China in 1950, is one of the most commonly used parents for early breeding programs in China. However, the formation and inheritance patterns of genomic constituents have not been studied. Therefore, this study aimed at understanding and exploring the genomic architecture of 146 DLP-15 derived cultivars with a common foundation parent DLP-15. Population structure based on sequencing data clustered genotypes into two groups (G1 and G2) supported by principal component analysis. Further exploration led to the identification of Chr-A08 with significantly differentiated regions between two groups. Moreover, we identified genome-wide identity by descent (IBD) segments (840 segments) to understand the genomic inheritance pattern in DLP-15 derived cultivars, spanning the 20-95 Mb region on Chr-A08. Interestingly, Chr-A08 depicted a unique inheritance pattern from DLP-15 to its derived cultivars. IBD-segment-based haplotype analysis suggested significant differences among the two groups. Phenotypic trait association with DLP-derived haplotypes concerning Chr-A08 suggested a significant increase in yield and fiber quality. Furthermore, distinguished IBD segments overlapped with previously reported QTLs concerning fiber yield and quality. Our results systematically identified genomic signatures transmitted from the foundation parent DLP-15 to its derived cultivars and provided a basis for further exploiting excellent haplotypes associated with DLP-15.

An enhanced photosynthesis and carbohydrate metabolic capability contributes to heterosis of the cotton (Gossypium hirsutum) hybrid 'Huaza Mian H318', as revealed by genome-wide gene expression analysis

Background

Heterosis has been exploited for decades in different crops due to resulting in dramatic increases in yield, but relatively little molecular evidence on this topic was reported in cotton.

Results

The elite cotton hybrid variety ‘Huaza Mian H318’ (H318) and its parental lines were used to explore the source of its yield heterosis. A four-year investigation of yield-related traits showed that the boll number of H318 showed higher stability than that of its two parents, both in suitable and unsuitable climate years. In addition, the hybrid H318 grew faster and showed higher fresh and dry weights than its parental lines at the seedling stage. Transcriptome analysis of seedlings identified 17,308 differentially expressed genes (DEGs) between H318 and its parental lines, and 3490 extremely changed DEGs were screened out for later analysis. Most DEGs (3472/3490) were gathered between H318 and its paternal line (4–5), and only 64 DEGs were found between H318 and its maternal line (B0011), which implied that H318 displays more similar transcriptional patterns to its maternal parent at the seedling stage. GO and KEGG analyses showed that these DEGs were highly enriched in photosynthesis, lipid metabolic, carbohydrate metabolic and oxidation-reduction processes, and the expression level of these DEGs was significantly higher in H318 relative to its parental lines, which implied that photosynthesis, metabolism and stress resistances were enhanced in H318.

Conclusion

The enhanced photosynthesis, lipid and carbohydrate metabolic capabilities contribute to the heterosis of H318 at the seedling stage, and establishes a material foundation for subsequent higher boll-setting rates in complex field environments.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07580-8.

Comparative transcriptome analysis between inbred and hybrids reveals molecular insights into yield heterosis of upland cotton

BACKGROUND

Utilization of heterosis has greatly improved the productivity of many crops worldwide. Understanding the potential molecular mechanism about how hybridization produces superior yield in upland cotton is critical for efficient breeding programs.

RESULTS

In this study, high, medium, and low hybrids varying in the level of yield heterosis were screened based on field experimentation of different years and locations. Phenotypically, high hybrid produced a mean of 14% more seed cotton yield than its better parent. Whole-genome RNA sequencing of these hybrids and their four inbred parents was performed using different tissues of the squaring stage. Comparative transcriptomic differences in each hybrid parent triad revealed a higher percentage of differentially expressed genes (DEGs) in each tissue. Expression level dominance analysis identified majority of hybrids DEGs were biased towards parent like expressions. An array of DEGs involved in ATP and protein binding, membrane, cell wall, mitochondrion, and protein phosphorylation had more functional annotations in hybrids. Sugar metabolic and plant hormone signal transduction pathways were most enriched in each hybrid. Further, these two pathways had most mapped DEGs on known seed cotton yield QTLs. Integration of transcriptome, QTLs, and gene co-expression network analysis discovered genes Gh_A03G1024, Gh_D08G1440, Gh_A08G2210, Gh_A12G2183, Gh_D07G1312, Gh_D08G1467, Gh_A03G0889, Gh_A08G2199, and Gh_D05G0202 displayed a complex regulatory network of many interconnected genes. qRT-PCR of these DEGs was performed to ensure the accuracy of RNA-Seq data.

CONCLUSIONS

Through genome-wide comparative transcriptome analysis, the current study identified nine key genes and pathways associated with biological process of yield heterosis in upland cotton. Our results and data resources provide novel insights and will be useful for dissecting the molecular mechanism of yield heterosis in cotton.

A novel discovery of a long terminal repeat retrotransposon-induced hybrid weakness in rice

Hybrid weakness is a post-zygotic hybridization barrier frequently observed in plants, including rice. In this study, we describe the genomic variation among three temperate japonica rice (Oryza sativa ssp. japonica) varieties 'Aranghyangchalbyeo' ('CH7'), 'Sanghaehyangheolua' ('CH8') and 'Shinseonchalbyeo' ('CH9'), carrying different hybrid weakness genes. The reciprocal progeny obtained from crossing any two varieties displayed characteristic hybrid weakness traits. We mapped and cloned a new locus, Hwc3 (hybrid weakness 3), on chromosome 4. Sequence analysis identified that a long terminal repeat (LTR) retrotransposon was inserted into the promoter region of the Hwc3 gene in 'CH7'. A 4-kb DNA fragment from 'CH7' containing the Hwc3 gene with the inserted LTR retrotransposon was able to induce hybrid weakness in hybrids with 'CH8' plants carrying the Hwc1 gene by genetic complementation. We investigated the differential gene expression profile of F1 plants exhibiting hybrid weakness and detected that the genes associated with energy metabolism were significantly down-regulated compared with the parents. Based on our results, we propose that LTR retrotransposons could be a potential cause of hybrid weakness in intrasubspecific hybrids in japonica rice. Understanding the molecular mechanisms underlying intrasubspecific hybrid weakness is important for increasing our knowledge on reproductive isolation and could have significant implications for rice improvement and hybrid breeding.

Physical mapping and candidate gene prediction of fertility restorer gene of cytoplasmic male sterility in cotton

Background

Cytoplasmic male sterility (CMS) is a maternally inherited trait failing to produce functional pollen. It plays a pivotal role in the exploitation of crop heterosis. The specific locus amplified fragment sequencing (SLAF-seq) as a high-resolution strategy for the identification of new SNPs on a large-scale is gradually applied for functional gene mining. The current study combined the bulked segregant analysis (BSA) with SLAF-seq to identify the candidate genes associated with fertility restorer gene (Rf) in CMS cotton.

Methods

Illumina sequencing systematically investigated the parents. A segregating population comprising of 30 + 30 F₂ individuals was developed using 3096A (female parent) as sterile and 866R (male parent) as a restorer. The original data obtained by dual-index sequencing were analyzed to obtain the reads of each sample that were compared to the reference genome in order to identify the SLAF tag with a polymorphism in parent lines and the SNP with read-associated coverage. Based on SLAF tags, SNP-index analysis, Euclidean distance (ED) correlation analysis, and whole genome resequencing, the hot regions were annotated.

Results

A total of 165,007 high-quality SLAF tags, with an average depth of 47.90× in the parents and 50.78× in F₂ individuals, were sequenced. In addition, a total of 137,741 SNPs were detected: 113,311 and 98,861 SNPs in the male and female parent, respectively. A correlation analysis by SNP-index and ED initially located the candidate gene on 1.35 Mb of chrD05, and 20 candidate genes were identified. These genes were involved in genetic variations, single base mutations, insertions, and deletions. Moreover, 42 InDel markers of the whole genome resequencing were also detected.

Conclusions

In this study, associated markers identified by super-BSA could accelerate the study of CMS in cotton, and as well as in other crops. Some of the 20 genes’ preliminary characteristics provided useful information for further studies on CMS crops.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4406-y) contains supplementary material, which is available to authorized users.

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.