Acyl-CoA N-acyltransferase influences fertility by regulating lipid metabolism and jasmonic acid biogenesis in cotton

The functions of laccase gene GhLAC15 in fiber colouration and development in brown-colored cotton

The monotonicity of color type in naturally colored cottons (NCCs) has become the main limiting factor to their widespread use, simultaneously coexisting with poor fiber quality. The synchronous improvement of fiber quality and color become more urgent and crucial as the demand for sustainable development increases. The homologous gene of wild cotton Gossypium stocksii LAC15 in G. hirsutum, GhLAC15, was also dominantly expressed in the developing fibers of brown cotton XC20 from 5 DPA (day post anthesis) to 25 DPA, especially at the secondary cell wall thickening stage (20 DPA and 25 DPA). In XC20 plants with downregulated GhLAC15 (GhLAC15i), a remarkable reduction in proanthocyanidins (PAs) and lignin contents was observed. Some of the key genes in the phenylpropane and flavonoid biosynthesis pathway were down-regulated in GhLAC15i plants. Notably, the fiber length of GhLAC15i plants showed an obvious increase and the fiber color was lightened. Moreover, we found that the thickness of cotton fiber cell wall was decreased in GhLAC15i plants and the fiber surface became smoother compared to that of WT. Taken together, this study revealed that GhLAC15 played an important role in PAs and lignin biosynthesis in naturally colored cotton fibers. It might mediate fiber color and fiber quality by catalyzing PAs oxidation and lignin polymerization, ultimately regulating fiber colouration and development.

The fertility research of "Huajin 6", a new variety of honeysuckle

The aim of this study was to investigate the fertility of "Huajin 6" and the effect of exogenous methyl jasmonate on its fertility. In this study, "Huajin 6", "Huajin 6" treated with methyl jasmonate and "Damaohua" were used as the research objects, the stamen phenotypes and the shape of pollen grains were observed, pollen viability and stigma receptivity were measured. The results showed that the pistil structure and function were normal, and although the stamen anthers did not dehisce, they were still capable of producing pollen with a certain amount of vigor. Methyl jasmonate could promote the opening of the flowers of "Huajin 6" and improve the development of pollen grains to a certain extent, but it could not promote anthers dehiscence of "Huajin 6". This study can provide theoretical guidance for the cultivation of new honeysuckle varieties using "Huajin 6".

Applications of Virus-Induced Gene Silencing in Cotton

Virus-induced gene silencing (VIGS) is an RNA-mediated reverse genetics technique that has become an effective tool to investigate gene function in plants. Cotton is one of the most important economic crops globally. In the past decade, VIGS has been successfully applied in cotton functional genomic studies, including those examining abiotic and biotic stress responses and vegetative and reproductive development. This article summarizes the traditional vectors used in the cotton VIGS system, the visible markers used for endogenous gene silencing, the applications of VIGS in cotton functional genomics, and the limitations of VIGS and how they can be addressed in cotton.

Heat-responsive microRNAs participate in regulating the pollen fertility stability of CMS-D2 restorer line under high-temperature stress

Anther development and pollen fertility of cytoplasmic male sterility (CMS) conditioned by Gossypium harknessii cytoplasm (CMS-D2) restorer lines are susceptible to continuous high-temperature (HT) stress in summer, which seriously hinders the large-scale application of "three-line" hybrids in production. Here, integrated small RNA, transcriptome, degradome, and hormone profiling was performed to explore the roles of microRNAs (miRNAs) in regulating fertility stability in mature pollens of isonuclear alloplasmic near-isogenic restorer lines NH and SH under HT stress at two environments. A total of 211 known and 248 novel miRNAs were identified, of which 159 were differentially expressed miRNAs (DEMs). Additionally, 45 DEMs in 39 miRNA clusters (PmCs) were also identified, and most highly expressed miRNAs were significantly induced in SH under extreme HT, especially four MIR482 and six MIR6300 family miRNAs. PmC28 was located in the fine-mapped interval of the Rf1 gene and contained two DEMs, gra-miR482_L-2R + 2 and gma-miR2118a-3p_R + 1_1ss18TG. Transcriptome sequencing identified 6281 differentially expressed genes, of which heat shock protein (HSP)-related genes, such as HSP70, HSP22, HSP18.5-C, HSP18.2 and HSP17.3-B, presented significantly reduced expression levels in SH under HT stress. Through integrating multi-omics data, we constructed a comprehensive molecular network of miRNA-mRNA-gene-KEGG containing 35 pairs of miRNA/target genes involved in regulating the pollen development in response to HT, among which the mtr-miR167a_R + 1, tcc-miR167c and ghr-miR390a, tcc-miR396c_L-1 and ghr-MIR169b-p3_1ss6AG regulated the pollen fertility by influencing ARF8 responsible for the auxin signal transduction, ascorbate and aldarate metabolism, and the sugar and lipid metabolism and transport pathways, respectively. Further combination with hormone analysis revealed that HT-induced jasmonic acid signaling could activate the expression of downstream auxin synthesis-related genes and cause excessive auxin accumulation, followed by a cascade of auxin signal transduction, ultimately resulting in pollen abortion. The results provide a new understanding of how heat-responsive miRNAs regulate the stability of fertility restoration for CMS-D2 cotton under heat stress.

Dose effects of restorer gene modulate pollen fertility in cotton CMS-D2 restorer lines via auxin signaling and flavonoid biosynthesis

KEY MESSAGE

Dose effects of Rf1 gene regulated retrieval mechanism of pollen fertility for CMS-D2 cotton. Cytoplasmic male sterility conditioned by Gossypium harknessii cytoplasm (CMS-D2) is an economical pollination control system for producing hybrid cotton seeds compared to artificial and chemical emasculation methods. However, the unstable restoring ability of restorer lines is a main barrier in the large-scale application of "three-line" hybrid cotton in China. Our phenotypic investigation determined that the homozygous Rf1Rf1 allelic genotype had a stronger ability to generate fertile pollen than the heterozygous Rf1rf1 allelic genotype. To decipher the genetic mechanisms that control the differential levels of pollen fertility, an integrated metabolomic and transcriptomic analysis was performed at two environments using pollen grains of four cotton genotypes differing in Rf1 alleles or cytoplasm. Totally 5,391 differential metabolite features were detected, and 369 specific differential metabolites (DMs) were identified between homozygous and heterozygous Rf1 allelic genotypes with CMS-D2 cytoplasm. In addition, transcriptome analysis identified 2,490 differentially expressed genes (DEGs) and 96 unique hub DEGs with dynamic regulation in this comparative combination. Further integrated analyses revealed that several key DEGs and DMs involved in indole biosynthesis, flavonoid biosynthesis, and sugar metabolism had strong network linkage with fertility restoration. In vitro application of auxin analogue NAA and inhibitor Auxinole confirmed that over-activated auxin signaling might inhibit pollen development, whereas suppressing auxin signaling partially promoted pollen development in CMS-D2 cotton. Our results provide new insight into how the dosage effects of the Rf1 gene regulate pollen fertility of CMS-D2 cotton.

Unveiling the hidden role of aquatic viruses in hydrocarbon pollution bioremediation

Hydrocarbon pollution poses substantial environmental risks to water and soil. Bioremediation, which utilizes microorganisms to manage pollutants, offers a cost-effective solution. However, the role of viruses, particularly bacteriophages (phages), in bioremediation remains unexplored. This study examines the diversity and activity of hydrocarbon-degradation genes encoded by environmental viruses, focusing on phages, within public databases. We identified 57 high-quality phage-encoded auxiliary metabolic genes (AMGs) related to hydrocarbon degradation, which we refer to as virus-encoded hydrocarbon degradation genes (vHYDEGs). These genes are encoded by taxonomically diverse aquatic phages and highlight the under-characterized global virosphere. Six protein families involved in the initial alkane hydroxylation steps were identified. Phylogenetic analyses revealed the diverse evolutionary trajectories of vHYDEGs across habitats, revealing previously unknown biodegraders linked evolutionarily with vHYDEGs. Our findings suggest phage AMGs may contribute to alkane and aromatic hydrocarbon degradation, participating in the initial, rate-limiting hydroxylation steps, thereby aiding hydrocarbon pollution bioremediation and promoting their propagation. To support future research, we developed vHyDeg, a database containing identified vHYDEGs with comprehensive annotations, facilitating the screening of hydrocarbon degradation AMGs and encouraging their bioremediation applications.

Deciphering the Atlas of Post-Translational Modification in Sugarcane

In plants, lysine acetylation (Kac), 2-hydroxyisobutyrylation (Khib), and lysine lactylation (Kla), the three new types of post-translational modification (PTM), play very important roles in growth, development, and resistance to adverse environmental stresses. Herein, we report the first global acetylome, 2-hydroxyisobutyrylome, and lactylome in sugarcane. A total of 8573 Kac, 4637 Khib, and 215 Kla sites across 3903, 1507, and 139 modified proteins were identified. Besides, homology analyses revealed the Kac, Khib, and Kla sites on histones were conserved between sugarcane and rice or poplar. Functional annotations demonstrated that the Kac, Khib, and Kla proteins were mainly involved in energy metabolism. In addition, a number of modified transcription factors and stress-related proteins, which were constitutively expressed in different tissues of sugarcane and induced by drought, cold or Sporisorium scitamineum stress, were identified. Finally, a proposed working mode on how PTM functions in sugarcane was depicted. We thus concluded that PTM should play a role in sugarcane growth, development, and response to biotic and abiotic stresses, but the mechanisms require further investigation. The present study provided the all-new comprehensive profile of proteins Kac, Khib, and Kla and a new perspective to understand the molecular mechanisms of protein PTMs in sugarcane.

The cotton pectin methyl esterase gene GhPME21 functions in microspore development and fertility in Gossypium hirsutum L

Pectin widely exists in higher plants' cell walls and intercellular space of higher plants and plays an indispensable role in plant growth and development. We identified 55 differentially expressed genes related to pectin degradation by transcriptomic analysis in the male sterile mutant, ms1. A gene encoding pectin methylesterase (GhPME21) was found to be predominantly expressed in the developing stamens of cotton but was significantly down-regulated in ms1 stamens. The tapetal layer of GhPME21 interfered lines (GhPME21i) was significantly thickened compared to that of WT at the early stage; anther compartment morphology of GhPME21i lines was abnormal, and the microspore wall was broken at the middle stage; Alexander staining showed that the pollen grains of GhPME21i lines differed greatly in volume at the late stage. The mature pollen surfaces of GhPME21i lines were deposited with discontinuous and broken sheets and prickles viewed under SEM. Fewer pollen tubes were observed to germinate in vitro in GhPME21i lines, while tiny of those in vivo were found to elongate to the ovary. The seeds harvested from GhPME21i lines as pollination donors were dry and hollow. The changes of phenotypes in GhPME21i lines at various stages illustrated that the GhPME21 gene played a vital role in the development of cotton stamens and controlled plant fertility by affecting stamen development, pollen germination, and pollen tube elongation. The findings of this study laid the groundwork for further research into the molecular mechanisms of PMEs involved in microspore formation and the creation of cotton male sterility materials.

Cloning and functional analysis of GhDFR1, a key gene of flavonoid synthesis pathway in naturally colored cotton

BACKGROUND

The naturally colored brown cotton fiber is the most widely used environmentally friendly textile material, which primarily contains proanthocyanidins and their derivatives. Many structural genes in the flavonoid synthesis pathway are known to improve the genetic resources of naturally colored cotton. Among them, DFR is a crucial late enzyme to synthesis both anthocyanins and proanthocyanidins in the plant flavonoid pathway.

METHODS

The protein sequences of GhDFRs were analyzed using bioinformatic tools. The expression levels of GhDFRs in various tissues and organs of upland cotton Zongxu1 (ZX1), were analyzed by quantitative real-time PCR, and the expression pattern of GhDFR1 during fiber development of white cotton and brown cotton was analyzed further. The function of GhDFR1 in NCC ZX1 was preliminarily analyzed by virus induced gene silencing (VIGS) technology.

RESULTS

Bioinformatic analysis revealed that GhDFRs sequences in upland cotton genome were extremely conserved. Furthermore, evolutionary tree analysis revealed that the functions of GhDFR1 and GhDFR2, and GhDFR3 and GhDFR4, presented different and shared some similarities. Our study showed GhDFR1 and GhDFR2 were specifically expressed in fibers, while GhDFR3 and GhDFR4 were specifically expressed in petals. GhDFR1 was exclusively expressed in brown cotton fiber at various stages of development and progressively increased with the growth of fiber, but the trend of expression in white cotton was quite the opposite. We silenced GhDFR1 expression in brown cotton fiber using VIGS technology, and observed the VIGS-interference plants. After reducing the expression level of GhDFR1, the period for significant GhDFR1 expression in the developing fibers changed, reducing the content of anthocyanins, and lightening the color of mature cotton fibers.

CONCLUSION

GhDFR1 was preferentially expressed in brown cotton during fiber development. The timing of GhDFR1 expression for flavonoid synthesis altered, resulting in anthocyanin contents reduced and the fiber color of the GhDFR1i lines lightened. These findings showed the role of GhDFR1 in fiber coloration of NCC and provided a new candidate for NCC genetic improvement.

High-temperature stress in crops: male sterility, yield loss and potential remedy approaches

Global food security is one of the utmost essential challenges in the 21st century in providing enough food for the growing population while coping with the already stressed environment. High temperature (HT) is one of the main factors affecting plant growth, development and reproduction and causes male sterility in plants. In male reproductive tissues, metabolic changes induced by HT involve carbohydrates, lipids, hormones, epigenetics and reactive oxygen species, leading to male sterility and ultimately reducing yield. Understanding the mechanism and genes involved in these pathways during the HT stress response will provide a new path to improve crops by using molecular breeding and biotechnological approaches. Moreover, this review provides insight into male sterility and integrates this with suggested strategies to enhance crop tolerance under HT stress conditions at the reproductive stage.

Responses of differential metabolites and pathways to high temperature in cucumber anther

Cucumber is one of the most important vegetable crops, which is widely planted all over the world. Cucumber always suffers from high-temperature stress in South China in summer. In this study, liquid chromatography-mass spectrometry (LC-MS) analysis was used to study the differential metabolites of cucumber anther between high-temperature (HT) stress and normal condition (CK). After HT, the pollen fertility was significantly reduced, and abnormal anther structures were observed by the paraffin section. In addition, the metabolomics analysis results showed that a total of 125 differential metabolites were identified after HT, consisting of 99 significantly upregulated and 26 significantly downregulated metabolites. Among these differential metabolites, a total of 26 related metabolic pathways were found, and four pathways showed significant differences, namely, porphyrin and chlorophyll metabolism; plant hormone signal transduction; amino sugar and nucleotide sugar metabolism; and glycine, serine, and threonine metabolism. In addition, pollen fertility was decreased by altering the metabolites of plant hormone signal transduction and amino acid and sugar metabolism pathway under HT. These results provide a comprehensive understanding of the metabolic changes in cucumber anther under HT.

Integrated analysis of metabolome and transcriptome reveals the cytoplasmic effects of CMS-D2 on pollen fertility resulting from disrupted lipid metabolism

Using cytoplasmic male sterility of Gossypium harknesii (CMS-D2) is an economical and effective method to produce cotton hybrids. However, the detrimental cytoplasmic effects of CMS-D2 on pollen fertility and fiber yields greatly limit the further development of three-line hybrid cotton in China. In this study, an integrated non-targeted metabolomics and transcriptome analysis was performed on mature pollens of maintainer line NB, isonuclear alloplasmic near-isogenic restorer lines NH and SH under two environments. A total of 820 metabolites were obtained, of which lipids and lipid-like molecules were the most, followed by organic acids derivatives, phenylpropanoids, and polyketides. Transcriptome analysis revealed significantly more differentially expressed genes (DEGs) in SH versus NH both in Anyang and Jiujiang, and most of the DEGs were significantly upregulated. Further KEGG analysis showed that most DEGs were enriched in the biosynthesis of unsaturated fatty acids, phenylalanine metabolism, and phagosome. Based on the weighted gene co-expression network analysis, totally 74 hub genes were also identified, of which three transcription factors, i.e., WRKY22, WRKY53, and ARF18 were significantly upregulated in SH and may play a negative regulatory role in pollen development by directly or indirectly regulating the jasmonic acid synthesis and signal transduction. Moreover, we found that the negative effects of CMS-D2 cytoplasm on pollen fertility were mainly due to disturbed lipid metabolism, especially the metabolic balance of unsaturated fatty acids, ultimately resulting in the decline of pollen fertility. Meanwhile, in the presence of CMS-D2 sterile cytoplasm, the cytoplasmic-nucleus interaction effects generated a substantial quantity of flavonoids involved in the fertility restoration process. This study preliminarily clarified some of the reasons for the negative effects of CMS-D2 cytoplasm on pollen fertility, and our results will provide an important theoretical reference for further breeding and improvement of three-line hybrid cotton in the future.

Genomic selection and genetic architecture of agronomic traits during modern rapeseed breeding

Rapeseed (Brassica napus L.) is an important oil-producing crop for the world. Its adaptation, yield and quality have been considerably improved in recent decades, but the genomic basis underlying successful breeding selection remains unclear. Hence, we conducted a comprehensive genomic assessment of rapeseed in the breeding process based on the whole-genome resequencing of 418 diverse rapeseed accessions. We unraveled the genomic basis for the selection of adaptation and agronomic traits. Genome-wide association studies identified 628 associated loci-related causative candidate genes for 56 agronomically important traits, including plant architecture and yield traits. Furthermore, we uncovered nonsynonymous mutations in plausible candidate genes for agronomic traits with significant differences in allele frequency distributions across the improvement process, including the ribosome recycling factor (BnRRF) gene for seed weight. This study provides insights into the genomic basis for improving rapeseed varieties and a valuable genomic resource for genome-assisted rapeseed breeding.

Insights into the Genomic Architecture of Seed and Pod Quality Traits in the U.S. Peanut Mini-Core Diversity Panel

Traits such as seed weight, shelling percent, percent sound mature kernels, and seed dormancy determines the quality of peanut seed. Few QTL (quantitative trait loci) studies using biparental mapping populations have identified QTL for seed dormancy and seed grade traits. Here, we report a genome-wide association study (GWAS) to detect marker–trait associations for seed germination, dormancy, and seed grading traits in peanut. A total of 120 accessions from the U.S. peanut mini-core collection were evaluated for seed quality traits and genotyped using Axiom SNP (single nucleotide polymorphism) array for peanut. We observed significant variation in seed quality traits in different accessions and different botanical varieties. Through GWAS, we were able to identify multiple regions associated with sound mature kernels, seed weight, shelling percent, seed germination, and dormancy. Some of the genomic regions that were SNP associated with these traits aligned with previously known QTLs. For instance, QTL for seed dormancy has been reported on chromosome A05, and we also found SNP on the same chromosome associated with seed dormancy, explaining around 20% of phenotypic variation. In addition, we found novel genomic regions associated with seed grading, seed germination, and dormancy traits. SNP markers associated with seed quality and dormancy identified here can accelerate the selection process. Further, exploring the function of candidate genes identified in the vicinity of the associated marker will assist in understanding the complex genetic network that governs seed quality.

GhPLP2 Positively Regulates Cotton Resistance to Verticillium Wilt by Modulating Fatty Acid Accumulation and Jasmonic Acid Signaling Pathway

Patatin-like proteins (PLPs) have non-specific lipid acyl hydrolysis (LAH) activity, which can hydrolyze membrane lipids into fatty acids and lysophospholipids. The vital role of PLPs in plant growth and abiotic stress has been well documented. However, the function of PLPs in plant defense responses against pathogens is still poorly understood. Here, we isolated and identified a novel cotton (Gossypium hirsutum) PLP gene GhPLP2. The expression of GhPLP2 was induced upon treatment with Verticillium dahliae, the signaling molecules jasmonic acid (JA) and ethylene (ETH) in cotton plants. Subcellular localization revealed that GhPLP2 was localized to the plasma membrane. GhPLP2-silenced cotton plants were more susceptible to infection by V. dahliae, while the overexpression of GhPLP2 in Arabidopsis enhanced its resistance to V. dahliae, which was apparent as mild symptoms, and a decrease in the disease index and fungal biomass. The hypersensitive response, deposition of callose, and H₂O₂ accumulation triggered by V. dahliae elicitor were reduced in GhPLP2-silenced cotton plants. The overexpression of GhPLP2 in Arabidopsis resulted in the accumulation of linoleic acid (LA, 18:2) and α-linolenic acid (ALA, 18:3) and facilitated the biosynthesis of JA and JA-mediated defensive responses. GhPLP2 silencing in cotton plants consistently reduced the accumulation of linoleic acid (LA, 18:2) and α-linolenic acid (ALA, 18:3) and suppressed the biosynthesis of JA and the defensive responses mediated by JA. These results indicate that GhPLP2 is involved in the resistance of cotton to V. dahliae by maintaining fatty acid metabolism pools for JA biosynthesis and activating the JA signaling pathway.

GhPLP2 Positively Regulates Cotton Resistance to Verticillium Wilt by Modulating Fatty Acid Accumulation and Jasmonic Acid Signaling Pathway

Patatin-like proteins (PLPs) have non-specific lipid acyl hydrolysis (LAH) activity, which can hydrolyze membrane lipids into fatty acids and lysophospholipids. The vital role of PLPs in plant growth and abiotic stress has been well documented. However, the function of PLPs in plant defense responses against pathogens is still poorly understood. Here, we isolated and identified a novel cotton (Gossypium hirsutum) PLP gene GhPLP2. The expression of GhPLP2 was induced upon treatment with Verticillium dahliae, the signaling molecules jasmonic acid (JA) and ethylene (ETH) in cotton plants. Subcellular localization revealed that GhPLP2 was localized to the plasma membrane. GhPLP2-silenced cotton plants were more susceptible to infection by V. dahliae, while the overexpression of GhPLP2 in Arabidopsis enhanced its resistance to V. dahliae, which was apparent as mild symptoms, and a decrease in the disease index and fungal biomass. The hypersensitive response, deposition of callose, and H₂O₂ accumulation triggered by V. dahliae elicitor were reduced in GhPLP2-silenced cotton plants. The overexpression of GhPLP2 in Arabidopsis resulted in the accumulation of linoleic acid (LA, 18:2) and α-linolenic acid (ALA, 18:3) and facilitated the biosynthesis of JA and JA-mediated defensive responses. GhPLP2 silencing in cotton plants consistently reduced the accumulation of linoleic acid (LA, 18:2) and α-linolenic acid (ALA, 18:3) and suppressed the biosynthesis of JA and the defensive responses mediated by JA. These results indicate that GhPLP2 is involved in the resistance of cotton to V. dahliae by maintaining fatty acid metabolism pools for JA biosynthesis and activating the JA signaling pathway.

Slow anion channel GhSLAC1 is essential for stomatal closure in response to drought stress in cotton

Drought is one of the abiotic stresses which affects the growth and development of plants, including cotton. The role of stomatal anion channel SLAC1 has been well established in regulating stomatal closure in response to drought stress in several plant species. However, the gene encoding for the main S-type anion channel SLAC1 in cotton has not been identified hence its role in drought stress response remains uncharacterized. In this study, we identified Gh_A08G1582 as the gene encoding for GhSLAC1 in cotton. The gene exhibited abundant expression in leaves and was localized in cell membrane. Furthermore, the expression of GhSLAC1 in Arabidopsis slac1-3 mutants rescued the defective stomatal movement phenotypes of the mutants, pointing to its role in stomata regulation. GhSLAC1 channel was activated by AtOST1 in Xenopus laevis oocytes and showed greater permeability for nitrate than chloride. Further data demonstrated that transgenic cotton lines with silenced GhSLAC1 exhibited obvious leaf wilting phenotype and strong stomatal closure insensitivity under drought stress. Taken together, these results demonstrate that GhSLAC1 is an essential element for stomatal closure in response to drought in cotton.

High day and night temperatures distinctively disrupt fatty acid and jasmonic acid metabolism, inducing male sterility in cotton

High temperature stress is an inevitable environmental factor in certain geographical regions. To study the effect of day and night high temperature stress on male reproduction, the heat-sensitive cotton line H05 was subjected to high temperature stress. High day/normal night (HN) and normal day/high night (NH) temperature treatments were compared with normal day/normal night (NN) temperature as a control. At the anther dehiscence stage, significant differences were observed, with a reduction in flower size and filament length, and sterility in pollen, seen in NH more than in HN. A total of 36 806 differentially expressed genes were screened, which were mainly associated with fatty acid and jasmonic acid (JA) metabolic pathways. Fatty acid and JA contents were reduced more in NH than HN. Under NH, ACYL-COA OXIDASE 2 (ACO2), a JA biosynthesis gene, was down-regulated. Interestingly, aco2 CRISPR-Cas9 mutants showed male sterility under the NN condition. The exogenous application of methyl jasmonate to early-stage buds of mutants rescued the sterile pollen and indehiscent anther phenotypes at the late stage. These data show that high temperature at night may affect fatty acid and JA metabolism in anthers by suppressing GhACO2 and generate male sterility more strongly than high day temperature.

Identification of miRNA-eQTLs in maize mature leaf by GWAS

Background

MiRNAs play essential roles in plant development and response to biotic and abiotic stresses through interaction with their target genes. The expression level of miRNAs shows great variations among different plant accessions, developmental stages, and tissues. Little is known about the content within the plant genome contributing to the variations in plants. This study aims to identify miRNA expression-related quantitative trait loci (miR-QTLs) in the maize genome.

Results

The miRNA expression level from next generation sequencing (NGS) small RNA libraries derived from mature leaf samples of the maize panel (200 maize lines) was estimated as phenotypes, and maize Hapmap v3.2.1 was chosen as the genotype for the genome-wide association study (GWAS). A total of four significant miR-eQTLs were identified contributing to miR156k-5p, miR159a-3p, miR390a-5p and miR396e-5p, and all of them are trans-eQTLs. In addition, a strong positive coexpression of miRNA was found among five miRNA families. Investigation of the effects of these miRNAs on the expression levels and target genes provided evidence that miRNAs control the expression of their targets by suppression and enhancement.

Conclusions

These identified significant miR-eQTLs contribute to the diversity of miRNA expression in the maize penal at the developmental stages of mature leaves in maize, and the positive and negative regulation between miRNA and its target genes has also been uncovered.

Genetic variants and underlying mechanisms influencing variance heterogeneity in maize

Traditional genetic studies focus on identifying genetic variants associated with the mean difference in a quantitative trait. Because genetic variants also influence phenotypic variation via heterogeneity, we conducted a variance-heterogeneity genome-wide association study to examine the contribution of variance heterogeneity to oil-related quantitative traits. We identified 79 unique variance-controlling single nucleotide polymorphisms (vSNPs) from the sequences of 77 candidate variance-heterogeneity genes for 21 oil-related traits using the Levene test (P < 1.0 × 10^-5 ). About 30% of the candidate genes encode enzymes that work in lipid metabolic pathways, most of which define clear expression variance quantitative trait loci. Of the vSNPs specifically associated with the genetic variance heterogeneity of oil concentration, 89% can be explained by additional linked mean-effects genetic variants. Furthermore, we demonstrated that gene × gene interactions play important roles in the formation of variance heterogeneity for fatty acid compositional traits. The interaction pattern was validated for one gene pair (GRMZM2G035341 and GRMZM2G152328) using yeast two-hybrid and bimolecular fluorescent complementation analyses. Our findings have implications for uncovering the genetic basis of hidden additive genetic effects and epistatic interaction effects, and we indicate opportunities to stabilize efficient breeding and selection of high-oil maize (Zea mays L.).

Functional analysis of GhCHS, GhANR and GhLAR in colored fiber formation of Gossypium hirsutum L

BACKGROUND

The formation of natural colored fibers mainly results from the accumulation of different anthocyanidins and their derivatives in the fibers of Gossypium hirsutum L. Chalcone synthase (CHS) is the first committed enzyme of flavonoid biosynthesis, and anthocyanidins are transported into fiber cells after biosynthesis mainly by Anthocyanidin reductase (ANR) and Leucoanthocyanidin reductase (LAR) to present diverse colors with distinct stability. The biochemical and molecular mechanism of pigment formation in natural colored cotton fiber is not clear.

RESULTS

The three key genes of GhCHS, GhANR and GhLAR were predominantly expressed in the developing fibers of colored cotton. In the GhCHSi, GhANRi and GhLARi transgenic cottons, the expression levels of GhCHS, GhANR and GhLAR significantly decreased in the developing cotton fiber, negatively correlated with the content of anthocyanidins and the color depth of cotton fiber. In colored cotton Zongxu1 (ZX1) and the GhCHSi, GhANRi and GhLARi transgenic lines of ZX1, HZ and ZH, the anthocyanidin contents of the leaves, cotton kernels, the mixture of fiber and seedcoat were all changed and positively correlated with the fiber color.

CONCLUSION

The three genes of GhCHS, GhANR and GhLAR were predominantly expressed early in developing colored cotton fibers and identified to be a key genes of cotton fiber color formation. The expression levels of the three genes affected the anthocyanidin contents and fiber color depth. So the three genes played a crucial part in cotton fiber color formation and has important significant to improve natural colored cotton quality and create new colored cotton germplasm resources by genetic engineering.

Comparative Transcriptome Analysis of Developing Seeds and Silique Wall Reveals Dynamic Transcription Networks for Effective Oil Production in Brassica napus L

Vegetable oil is an essential constituent of the human diet and renewable raw material for industrial applications. Enhancing oil production by increasing seed oil content in oil crops is the most viable, environmentally friendly, and sustainable approach to meet the continuous demand for the supply of vegetable oil globally. An in-depth understanding of the gene networks involved in oil biosynthesis during seed development is a prerequisite for breeding high-oil-content varieties. Rapeseed (Brassica napus) is one of the most important oil crops cultivated on multiple continents, contributing more than 15% of the world’s edible oil supply. To understand the phasic nature of oil biosynthesis and the dynamic regulation of key pathways for effective oil accumulation in B. napus, comparative transcriptomic profiling was performed with developing seeds and silique wall (SW) tissues of two contrasting inbred lines with ~13% difference in seed oil content. Differentially expressed genes (DEGs) between high- and low-oil content lines were identified across six key developmental stages, and gene enrichment analysis revealed that genes related to photosynthesis, metabolism, carbohydrates, lipids, phytohormones, transporters, and triacylglycerol and fatty acid synthesis tended to be upregulated in the high-oil-content line. Differentially regulated DEG patterns were revealed for the control of metabolite and photosynthate production in SW and oil biosynthesis and accumulation in seeds. Quantitative assays of carbohydrates and hormones during seed development together with gene expression profiling of relevant pathways revealed their fundamental effects on effective oil accumulation. Our results thus provide insights into the molecular basis of high seed oil content (SOC) and a new direction for developing high-SOC rapeseed and other oil crops.

Imbalance of tyrosine by modulating TyrA arogenate dehydrogenases impacts growth and development of Arabidopsis thaliana

l-Tyrosine is an essential aromatic amino acid required for the synthesis of proteins and a diverse array of plant natural products; however, little is known on how the levels of tyrosine are controlled in planta and linked to overall growth and development. Most plants synthesize tyrosine by TyrA arogenate dehydrogenases, which are strongly feedback-inhibited by tyrosine and encoded by TyrA1 and TyrA2 genes in Arabidopsis thaliana. While TyrA enzymes have been extensively characterized at biochemical levels, their in planta functions remain uncertain. Here we found that TyrA1 suppression reduces seed yield due to impaired anther dehiscence, whereas TyrA2 knockout leads to slow growth with reticulate leaves. The tyra2 mutant phenotypes were exacerbated by TyrA1 suppression and rescued by the expression of TyrA2, TyrA1 or tyrosine feeding. Low-light conditions synchronized the tyra2 and wild-type growth, and ameliorated the tyra2 leaf reticulation. After shifting to normal light, tyra2 transiently decreased tyrosine and subsequently increased aspartate before the appearance of the leaf phenotypes. Overexpression of the deregulated TyrA enzymes led to hyper-accumulation of tyrosine, which was also accompanied by elevated aspartate and reticulate leaves. These results revealed that TyrA1 and TyrA2 have distinct and overlapping functions in flower and leaf development, respectively, and that imbalance of tyrosine, caused by altered TyrA activity and regulation, impacts growth and development of Arabidopsis. The findings provide critical bases for improving the production of tyrosine and its derived natural products, and further elucidating the coordinated metabolic and physiological processes to maintain tyrosine levels in plants.

Improvement of Verticillium Wilt Resistance by Applying Arbuscular Mycorrhizal Fungi to a Cotton Variety with High Symbiotic Efficiency under Field Conditions

Arbuscular mycorrhizal fungi (AMF) play an important role in nutrient cycling processes and plant stress resistance. To evaluate the effect of Rhizophagus irregularis CD1 on plant growth promotion (PGP) and Verticillium wilt disease, the symbiotic efficiency of AMF (SEA) was first investigated over a range of 3% to 94% in 17 cotton varieties. The high-SEA subgroup had significant PGP effects in a greenhouse. From these results, the highest-SEA variety of Lumian 1 was selected for a two-year field assay. Consistent with the performance from the greenhouse, the AMF-mediated PGP of Lumian 1 also produced significant results, including an increased plant height, stem diameter, number of petioles, and phosphorus content. Compared with the mock treatment, AMF colonization obviously inhibited the symptom development of Verticillium dahliae and more strongly elevated the expression of pathogenesis-related genes and lignin synthesis-related genes. These results suggest that AMF colonization could lead to the mycorrhiza-induced resistance (MIR) of Lumian 1 to V. dahliae. Interestingly, our results indicated that the AMF endosymbiont could directly inhibit the growth of phytopathogenic fungi including V. dahliae by releasing undefined volatiles. In summary, our results suggest that stronger effects of AMF application result from the high-SEA.

Genome-wide association study reveals candidate genes influencing lipids and diterpenes contents in Coffea arabica L

Lipids, including the diterpenes cafestol and kahweol, are key compounds that contribute to the quality of coffee beverages. We determined total lipid content and cafestol and kahweol concentrations in green beans and genotyped 107 Coffea arabica accessions, including wild genotypes from the historical FAO collection from Ethiopia. A genome-wide association study was performed to identify genomic regions associated with lipid, cafestol and kahweol contents and cafestol/kahweol ratio. Using the diploid Coffea canephora genome as a reference, we identified 6,696 SNPs. Population structure analyses suggested the presence of two to three groups (K = 2 and K = 3) corresponding to the east and west sides of the Great Rift Valley and an additional group formed by wild accessions collected in western forests. We identified 5 SNPs associated with lipid content, 4 with cafestol, 3 with kahweol and 9 with cafestol/kahweol ratio. Most of these SNPs are located inside or near candidate genes related to metabolic pathways of these chemical compounds in coffee beans. In addition, three trait-associated SNPs showed evidence of directional selection among cultivated and wild coffee accessions. Our results also confirm a great allelic richness in wild accessions from Ethiopia, especially in accessions originating from forests in the west side of the Great Rift Valley.

Proteome-wide Analysis of Lysine 2-hydroxyisobutyrylation in Developing Rice (Oryza sativa) Seeds

Lysine 2-hydroxyisobutyrylation is a recently identified protein post-translational modification that is known to affect the association between histone and DNA. However, non-histone protein lysine 2-hydroxyisobutyrylation remains largely unexplored. Utilizing antibody-based affinity enrichment and nano-HPLC/MS/MS analyses of 2-hydroxyisobutyrylation peptides, we efficaciously identified 9,916 2-hydroxyisobutyryl lysine sites on 2,512 proteins in developing rice seeds, representing the first lysine 2-hydroxyisobutyrylome dataset in plants. Functional annotation analyses indicated that a wide variety of vital biological processes were preferably targeted by lysine 2-hydroxyisobutyrylation, including glycolysis/gluconeogenesis, TCA cycle, starch biosynthesis, lipid metabolism, protein biosynthesis and processing. Our finding showed that 2-hydroxyisobutyrylated histone sites were conserved across plants, human, and mouse. A number of 2-hydroxyisobutyryl sites were shared with other lysine acylations in both histone and non-histone proteins. Comprehensive analysis of the lysine 2-hydroxyisobutyrylation sites illustrated that the modification sites were highly sequence specific with distinct motifs, and they had less surface accessibility than other lysine residues in the protein. Overall, our study provides the first systematic analysis of lysine 2-hydroxyisobutyrylation proteome in plants, and it serves as an important resource for future investigations of the regulatory mechanisms and functions of lysine 2-hydroxyisobutyrylation.

Proteome-wide lysine acetylation identification in developing rice (Oryza sativa) seeds and protein co-modification by acetylation, succinylation, ubiquitination, and phosphorylation

Protein lysine acetylation is a highly conserved post-translational modification with various biological functions. However, only a limited number of acetylation sites have been reported in plants, especially in cereals, and the function of non-histone protein acetylation is still largely unknown. In this report, we identified 1003 lysine acetylation sites in 692 proteins of developing rice seeds, which greatly extended the number of known acetylated sites in plants. Seven distinguished motifs were detected flanking acetylated lysines. Functional annotation analyses indicated diverse biological processes and pathways engaged in lysine acetylation. Remarkably, we found that several key enzymes in storage starch synthesis pathway and the main storage proteins were heavily acetylated. A comprehensive comparison of the rice acetylome, succinylome, ubiquitome and phosphorylome with available published data was conducted. A large number of proteins carrying multiple kinds of modifications were identified and many of these proteins are known to be key enzymes of vital metabolic pathways. Our study provides extending knowledge of protein acetylation. It will have critical reference value for understanding the mechanisms underlying PTM mediated multiple signal integration in the regulation of metabolism and development in plants.

Evolution, gene expression profiling and 3D modeling of CSLD proteins in cotton

BACKGROUND

Among CESA-like gene superfamily, the cellulose synthase-like D (CSLD) genes are most similar to cellulose synthase genes and have been reported to be involved in tip-growing cell and stem development. However, there has been no genome-wide characterization of this gene subfamily in cotton. We thus sought to analyze the evolution and functional characterization of CSLD proteins in cotton based on fully sequenced cotton genomes.

RESULTS

A total of 23 full-length CSLD proteins were identified in Gossypium raimondii, Gossypium arboreum and Gossypium hirsutum. The phylogenetic tree divided the CSLD proteins into five clades with strong support: CSLD1, CSLD2/3, CSLD4, CSLD5 and CSLD6. The total expression of GhCSLD genes was the highest in androecium & gynoecium (mostly contributed by CSLD1 and CSLD4) compared with other CSL genes. CSLD1 and CSLD4 were only highly expressed in androecium & gynoecium (A&G), and showed tissue-specific expression. The total expression of CSLD2/3, 5 and 6 was highest in the specific tissues. These results suggest that CSLD genes showed the different pattern of expression. Cotton CSLD proteins were subjected to different evolutionary pressures, and the CSLD1 and CSLD4 proteins exhibited episodic and long-term shift positive selection. The predicted three-dimensional structure of GrCSLD1 suggested that GrCSLD1 belongs to glycosyltransferase family 2. The amino acid residues under positive selection in the CSLD1 lineage are positioned in a region adjacent to the class-specific region (CSR), β1-strand and transmembrane helices (TMHs) in the GrCSLD1structure.

CONCLUSION

Our results characterized the CSLD proteins by an integrated approach containing phylogeny, transcriptional profiling and 3D modeling. The study added to the understanding about the importance of the CSLD family and provide a useful reference for selecting candidate genes and their associations with the biosynthesis of the cell wall in cotton.

Differential transcript profiling alters regulatory gene expression during the development of Gossypium arboreum, G.stocksii and somatic hybrids

Polyploidy or genome doubling (i.e., the presence of two or more diploid parental genome sets within an organism) are very important in higher plants. Of particular interest are the mechanisms in the new microenvironment of the common nucleus, where doubled regulatory networks interact to generate a viable genetic system capable of regulating growth, development and responses to the environment. To determine the effects of whole genome merging and doubling on the global gene expression architecture of a new polyploid, derived from protoplast fusion of the A₁A₁ genome of Gossypium arboreum and the E₁E₁ genome of Gossypium stocksii, we monitored gene expression through cDNA-AFLP in the somatic hybrids (G. arboreum + G. stocksii). The genomic expression patterns of the somatic hybrids revealed that changes in expression levels mainly involved regulatory genes (31.8% of the gene expression profiles), and the AA and EE genomes contributed equally to genome-wide expression in the newly formed AAEE genome from additivity and dominance effects. These results provide a novel perspective on polyploid gene regulation and hint at the underlying genetic basis of allopolyploid adaption in the new microenvironmental nucleus.

Two Lysin-Motif Receptor Kinases, Gh-LYK1 and Gh-LYK2, Contribute to Resistance against Verticillium wilt in Upland Cotton

Lysin-motif (LysM) receptor kinases (LYKs) play essential roles in recognition of chitin and activation of defense responses against pathogenic fungi in the model plants Arabidopsis and rice. The function of LYKs in non-model plants, however, remains elusive. In the present work, we found that the transcription of two LYK-encoding genes from cotton, Gh-LYK1 and Gh-LYK2, was induced after Verticillium dahliae infection. Virus-induced gene silencing (VIGS) of Gh-LYK1 and Gh-LYK2 in cotton plants compromises resistance to V. dahliae. As putative pattern recognition receptors (PRRs), both Gh-LYK1 and Gh-LYK2 are membrane-localized, and all three LysM domains of Gh-LYK1 and Gh-LYK2 are required for their chitin-binding ability. However, since Gh-LYK2, but not Gh-LYK1, is a pseudo-kinase and, on the other hand, the ectodomain (ED) of Gh-LYK2 can induce reactive oxygen species (ROS) burst in planta, Gh-LYK2 and Gh-LYK1 may contribute differently to cotton defense. Taken together, our results establish that both Gh-LYK1 and Gh-LYK12 are required for defense against V. dahliae in cotton, possibly through different mechanisms.

Comparative Transcriptomic Analysis of Two Brassica napus Near-Isogenic Lines Reveals a Network of Genes That Influences Seed Oil Accumulation

Rapeseed (Brassica napus) is an important oil seed crop, providing more than 13% of the world's supply of edible oils. An in-depth knowledge of the gene network involved in biosynthesis and accumulation of seed oil is critical for the improvement of B. napus. Using available genomic and transcriptomic resources, we identified 1,750 acyl-lipid metabolism (ALM) genes that are distributed over 19 chromosomes in the B. napus genome. B. rapa and B. oleracea, two diploid progenitors of B. napus, contributed almost equally to the ALM genes. Genome collinearity analysis demonstrated that the majority of the ALM genes have arisen due to genome duplication or segmental duplication events. In addition, we profiled the expression patterns of the ALM genes in four different developmental stages. Furthermore, we developed two B. napus near isogenic lines (NILs). The high oil NIL, YC13-559, accumulates significantly higher (∼10%) seed oil compared to the other, YC13-554. Comparative gene expression analysis revealed upregulation of lipid biosynthesis-related regulatory genes in YC13-559, including SHOOTMERISTEMLESS, LEAFY COTYLEDON 1 (LEC1), LEC2, FUSCA3, ABSCISIC ACID INSENSITIVE 3 (ABI3), ABI4, ABI5, and WRINKLED1, as well as structural genes, such as ACETYL-CoA CARBOXYLASE, ACYL-CoA DIACYLGLYCEROL ACYLTRANSFERASE, and LONG-CHAIN ACYL-CoA SYNTHETASES. We observed that several genes related to the phytohormones, gibberellins, jasmonate, and indole acetic acid, were differentially expressed in the NILs. Our findings provide a broad account of the numbers, distribution, and expression profiles of acyl-lipid metabolism genes, as well as gene networks that potentially control oil accumulation in B. napus seeds. The upregulation of key regulatory and structural genes related to lipid biosynthesis likely plays a major role for the increased seed oil in YC13-559.

RNA Interference for Functional Genomics and Improvement of Cotton (Gossypium sp.)

RNA interference (RNAi), is a powerful new technology in the discovery of genetic sequence functions, and has become a valuable tool for functional genomics of cotton (Gossypium sp.). The rapid adoption of RNAi has replaced previous antisense technology. RNAi has aided in the discovery of function and biological roles of many key cotton genes involved in fiber development, fertility and somatic embryogenesis, resistance to important biotic and abiotic stresses, and oil and seed quality improvements as well as the key agronomic traits including yield and maturity. Here, we have comparatively reviewed seminal research efforts in previously used antisense approaches and currently applied breakthrough RNAi studies in cotton, analyzing developed RNAi methodologies, achievements, limitations, and future needs in functional characterizations of cotton genes. We also highlighted needed efforts in the development of RNAi-based cotton cultivars, and their safety and risk assessment, small and large-scale field trials, and commercialization.

Cotton Leaf Curl Multan Virus-Derived Viral Small RNAs Can Target Cotton Genes to Promote Viral Infection

RNA silencing is a conserved mechanism in plants that targets viruses. Viral small RNAs (vsiRNAs) can be generated from viral double-stranded RNA replicative intermediates within the infected host, or from host RNA-dependent RNA polymerases activity on viral templates. The abundance and profile of vsiRNAs in viral infections have been reported previously. However, the involvement of vsiRNAs during infection of the Geminiviridae family member cotton leaf curl virus (CLCuD), which causes significant economic losses in cotton growing regions, remains largely uncharacterized. Cotton leaf curl Multan virus (CLCuMuV) associated with a betasatellite called Cotton leaf curl Multan betasatellite (CLCuMuB) is a major constraint to cotton production in South Asia and is now established in Southern China. In this study, we obtained the profiles of vsiRNAs from CLCuMV and CLCuMB in infected upland cotton (Gossypium hirsutum) plants by deep sequencing. Our data showed that vsiRNA that were derived almost equally from sense and antisense CLCuD DNA strands accumulated preferentially as 21- and 22-nucleotide (nt) small RNA population and had a cytosine bias at the 5'-terminus. Polarity distribution revealed that vsiRNAs were almost continuously present along the CLCuD genome and hotspots of sense and antisense strands were mainly distributed in the Rep proteins region of CLCuMuV and in the C1 protein of CLCuMuB. In addition, hundreds of host transcripts targeted by vsiRNAs were predicted, many of which encode transcription factors associated with biotic and abiotic stresses. Quantitative real-time polymerase chain reaction analysis of selected potential vsiRNA targets showed that some targets were significantly down-regulated in CLCuD-infected cotton plants. We also verified the potential function of vsiRNA targets that may be involved in CLCuD infection by virus-induced gene silencing (VIGS) and 5'-rapid amplification of cDNA end (5'-RACE). Here, we provide the first report on vsiRNAs responses to CLCuD infection in cotton.