Integration of eQTL and machine learning to dissect causal genes with pleiotropic effects in genetic regulation networks of seed cotton yield.

The dissection of a gene regulatory network (GRN) that complements the genome-wide association study (GWAS) locus and the crosstalk underlying multiple agronomical traits remains a major challenge. In this study, we generate 558 transcriptional profiles of lint-bearing ovules at one day post-anthesis from a selective core cotton germplasm, from which 12,207 expression quantitative trait loci (eQTLs) are identified. Sixty-six known phenotypic GWAS loci are colocalized with 1,090 eQTLs, forming 38 functional GRNs associated predominantly with seed yield. Of the eGenes, 34 exhibit pleiotropic effects. Combining the eQTLs within the seed yield GRNs significantly increases the portion of narrow-sense heritability. The extreme gradient boosting (XGBoost) machine learning approach is applied to predict seed cotton yield phenotypes on the basis of gene expression. Top-ranking eGenes (NF-YB3, FLA2, and GRDP1) derived with pleiotropic effects on yield traits are validated, along with their potential roles by correlation analysis, domestication selection analysis, and transgenic plants.

Global identification of natural antisense transcripts in Gossypium hirsutum and Gossypium barbadense under chilling stress.

Natural antisense transcripts (NATs) in model plants have been recognized as important regulators of gene expression under abiotic stresses. However, the functional roles of NATs in crops under low temperature are still unclear. Here, we identified 815 and 689 NATs from leaves of Gossypium hirsutum and G. barbadense under chilling stress. Among those, 224 NATs were identified as interspecific homologs between the two species. The correlation coefficients for expression of NATs and their cognate sense genes (CSG) were 0.43 and 0.37 in G. hirsutum and G. barbadense, respectively. Furthermore, expression of interspecific NATs and CSGs alike was highly consistent under chilling stress with correlation coefficients of 0.90-0.91. Four cold-associated NATs were selected for functional validation using virus-induced gene silencing (VIGS). Our results suggest that CAN1 engage in the molecular regulation of chilling stress by regulating SnRK2.8 expression. This highly conserved NAT have valuable potential for applications in breeding cold-tolerant cotton.

Establishment of qRT-PCR for absolute quantitative detection of Chikungunya virus / 中国热带医学

@#Abstract: Objective　To develop a real-time fluorescent quantitative RT-PCR (qRT-PCR) method for qualitative and quantitative Chikungunya virus (CHIKV) analysis. Methods Based on the systematic analysis of the genomic sequences of Chikungunya and its related arboviruses, the specific nucleic acid sequences for Chikungunya virus were screened and identified, and then the primers and TaqMan probe were designed. Meanwhile, the human GAPDH gene was used as an internal reference. The reaction system for qRT-PCR was systematically optimized by L9(34) orthogonal design, and a rapid detection method for Chikungunya by qRT-PCR based on TaqMan probe methods was established. The sensitivity, specificity, reproducibility, and coverage of the established method were analyzed in detail. The standard curve was made, and the absolute quantitative method was established using the cloned nucleic acid fragments as positive samples. Results　A real-time fluorescent quantitative RT-PCR assay was developed for the qualitative and quantitative analysis of Chikungunya virus. The reaction system included Chikungunya virus and reference internal gene specific primers and probe, RT/Taq enzyme mixture, reaction buffer, and negative and positive reference. The established method obtained positive results with the ROSS strain of ECSA subtype, LR2006 strain of IOL branch, 181/25 strain of Asian type and Dongguan 2010 epidemic strains of Chikungunya virus, but there was no cross-reaction with other 18 arboviruses belonging to Flaviviruses, Alphaviruses and Bunyavirus. The minimum detection limit of the established method was 5.80 copies/mL, and a linear relationship was observed between the amount of input plasmid DNA and fluorescence signal value over a range of 5.80×102 copies/mL to 5.80×1010 copies/mL, and the correlation coefficient was 0.999 5. The qRT-PCR amplification efficiency was 91%, and the intra-assay variations and inter-assay variations were 0.01-0.07 and 0.03-0.11, respectively. Conclusions　The TaqMan qRT-PCR method developed in this study can qualitatively and quantitatively detect Chikungunya virus rapidly with specificity and sensitivity, providing a technical method for the prevention and control of this viral disease.

Profiling of H3K4me3 Modification in Plants using Cleavage under Targets and Tagmentation.

Epigenomic regulation at the chromatic level, including DNA and histone modifications, behaviors of transcription factors, and non-coding RNAs with their recruited proteins, lead to temporal and spatial control of gene expression. Cleavage under targets and tagmentation (CUT&Tag) is an enzyme-tethering method in which the specific chromatin protein is firstly recognized by its specific antibody, and then the antibody tethers a protein A-transposase (pA-Tn5) fusion protein, which cleaves the targeted chromatin in situ by the activation of magnesium ions. Here, we provide our previously published CUT&Tag protocol using intact nuclei isolated from allortetraploid cotton leaves with modification. This step-by-step protocol can be used for epigenomic research in plants. In addition, substantial modifications for plant nuclei isolation are provided with critical comments.

Functional examination of lncRNAs in allotetraploid Gossypium hirsutum.

BACKGROUND: An evolutionary model using diploid and allotetraploid cotton species identified 80 % of non-coding transcripts in allotetraploid cotton as being uniquely activated in comparison with its diploid ancestors. The function of the lncRNAs activated in allotetraploid cotton remain largely unknown. RESULTS: We employed transcriptome analysis to examine the relationship between the lncRNAs and mRNAs of protein coding genes (PCGs) in cotton leaf tissue under abiotic stresses. LncRNA expression was preferentially associated with that of the flanking PCGs. Selected highly-expressed lncRNA candidates (n = 111) were subjected to a functional screening pilot test in which virus-induced gene silencing was integrated with abiotic stress treatment. From this low-throughput screen, we obtained candidate lncRNAs relating to plant height and tolerance to drought and other abiotic stresses. CONCLUSIONS: Low-throughput screen is an effective method to find functional lncRNA for further study. LncRNAs were more active in abiotic stresses than PCG expression, especially temperature stress. LncRNA XLOC107738 may take a cis-regulatory role in response to environmental stimuli. The degree to which lncRNAs are constitutively expressed may impact expression patterns and functions on the individual gene level rather than in genome-wide aggregate.

Neofunctionalization of a polyploidization-activated cotton long intergenic non-coding RNA DAN1 during drought stress regulation.

The genomic shock of whole-genome duplication (WGD) and hybridization introduces great variation into transcriptomes, for both coding and noncoding genes. An altered transcriptome provides a molecular basis for improving adaptation during the evolution of new species. The allotetraploid cotton, together with the putative diploid ancestor species compose a fine model for study the rapid gene neofunctionalization over the genome shock. Here we report on Drought-Associated Non-coding gene 1 (DAN1), a long intergenic noncoding RNA (lincRNA) that arose from the cotton progenitor A-diploid genome after hybridization and WGD events during cotton evolution. DAN1 in allotetraploid upland cotton (Gossypium hirsutum) is a drought-responsive lincRNA predominantly expressed in the nucleoplasm. Chromatin isolation by RNA purification profiling and electrophoretic mobility shift assay analysis demonstrated that GhDAN1 RNA can bind with DNA fragments containing AAAG motifs, similar to DNA binding with one zinc finger transcription factor binding sequences. The suppression of GhDAN1 mainly regulates genes with AAAG motifs in auxin-response pathways, which are associated with drought stress regulation. As a result, GhDAN1-silenced plants exhibit improved tolerance to drought stress. This phenotype resembles the drought-tolerant phenotype of the A-diploid cotton ancestor species, which has an undetectable expression of DAN1. The role of DAN1 in cotton evolution and drought tolerance regulation suggests that the genomic shock of interspecific hybridization and WGD stimulated neofunctionalization of non-coding genes during the natural evolutionary process.

Functional analysis of the cotton CLE polypeptide signaling gene family in plant growth and development.

The CLAVATA3 (CLV3)/EMBRYO SURROUNDING REGION (ESR)-RELATED (CLE) gene family encodes a large number of polypeptide signaling molecules involved in the regulation of shoot apical meristem division and root and vascular bundle development in a variety of plants. CLE family genes encode important short peptide hormones; however, the functions of these signaling polypeptides in cotton remain largely unknown. In the current work, we studied the effects of the CLE family genes on growth and development in cotton. Based on the presence of a conserved CLE motif of 13 amino acids, 93 genes were characterized as GhCLE gene family members, and these were subcategorized into 7 groups. A preliminary analysis of the cotton CLE gene family indicated that the activity of its members tends to be conserved in terms of both the 13-residue conserved domain at the C-terminus and their subcellular localization pattern. Among the 14 tested genes, the ectopic overexpression of GhCLE5::GFP partially mimicked the phenotype of the clv3 mutant in Arabidopsis. GhCLE5 could affect the endogenous CLV3 in binding to the receptor complex, comprised of CLV1, CLV2, and CRN, in the yeast two-hybrid assay and split-luciferase assay. Silencing GhCLE5 in cotton caused a short seedling phenotype. Therefore, we concluded that the cotton GhCLE gene family is functionally conserved in apical shoot development regulation. These results indicate that CLE also plays roles in cotton development as a short peptide hormone.

Pseudomonas species isolated via high-throughput screening significantly protect cotton plants against verticillium wilt.

Verticillium wilt (VW) caused by Verticillium dahliae is a devastating soil-borne disease that causes severe yield losses in cotton and other major crops worldwide. Here we conducted a high-throughput screening of isolates recovered from 886 plant rhizosphere samples taken from the three main cotton-producing areas of China. Fifteen isolates distributed in different genera of bacteria that showed inhibitory activity against V. dahliae were screened out. Of these, two Pseudomonas strains, P. protegens XY2F4 and P. donghuensis 22G5, showed significant inhibitory action against V. dahliae. Additional comparative genomic analyses and phenotypical assays confirmed that P. protegens XY2F4 and P. donghuensis 22G5 were the strains most efficient at protecting cotton plants against VW due to specific biological control products they produced. Importantly, we identified a significant efficacy of the natural tropolone compound 7-hydroxytropolone (7-HT) against VW. By phenotypical assay using the wild-type 22G5 and its mutant strain in 7-HT production, we revealed that the 7-HT produced by P. donghuensis is the major substance protecting cotton against VW. This study reveals that Pseudomonas specifically has gene clusters that allow the production of effective antipathogenic metabolites that can now be used as new agents in the biocontrol of VW.

Role of phasiRNAs from two distinct phasing frames of GhMYB2 loci in cis- gene regulation in the cotton genome.

BACKGROUND: Phased small interfering RNA (phasiRNA) is primarily derived from the 22-nt miRNA targeting loci. GhMYB2, a gene with potential roles in cotton fiber cell fate determination, is a target gene of miR828 and miR858 in the generation of phasiRNAs. RESULTS: In the presented work, through the evaluation of phasing scores and phasiRNA distribution pattern, we found that phasiRNAs from GhMYB2 were derived from the 3' cleavage fragments of 22-nt miR828 and 21-nt miR858 respectively. These two miRNA targeting sites initiated two phasing frames on transcripts of one locus. By means of RNA ligase-mediated rapid amplification of cDNA ends (RLM-RACE), we further demonstrated that phasiRNAs derived from the two phasing frames played a role in cis-regulation of GhMYB2. The phasiRNAs derived from GhMYB2 were expressed in the somatic tissues, especially in anther and hypocotyl. We further employed our previous small RNA sequencing data as well as the degradome data of cotton fiber bearing ovules, anthers, hypocotyls and embryogenic calli tissues published in public databases, to validate the expression, phasing pattern and functions of phasiRNAs. CONCLUSIONS: The presenting research provide insights of the molecular mechanism of phasiRNAs in regulation of GhMYB2 loci.

Root Morphological Traits and Spatial Distribution under Different Nitrogen Treatments and Their Relationship with Grain Yield in Super Hybrid Rice.

Succeeding in breeding super hybrid rice has been considered as a great progress in rice production in China. This on-farm study was conducted with Minirhizotron techniques to identify dynamic root morphological traits and distribution (0-30 cm) under different nitrogen treatments. Five elite super hybrid rice cultivars, Liangyoupeijiu (LYPJ), Yliangyou 1(YLY1), Yliangyou 2(YLY2), Yliangyou 900(YLY900) and Super 1000(S1000), were grown at four N levels: 0 kg ha-1 (N1), 210 kg ha-1 (N2), 300 kg ha-1 (N3) and 390 kg ha-1 (N4) in 2015 and 2016. Results showed these cultivars had greater root traits and higher grain yield under N3. Total root number (TRN) and total root length (TRL) of these cultivars reached maximum at 55 days after transplanting (DAT). The new released cultivars YLY900 and S1000 were featured with an improved root system among these cultivars. The percentage of root number on 10-20 cm soil was over 50% compared with other soil layer. A significant positive correlation was found between grain yield and both TRN and TRL at 10-20 cm soil layer (P < 0.01). Given this situation, the grain yield of super rice cultivars could be further improved by increasing the proportion of roots at 10-20 cm soil layer under suitable nitrogen management.