Genomic insights into glume pubescence in durum wheat: GWAS and haplotype analysis implicates TdELD1-1A as a candidate gene

Glume pubescence is an important morphological trait for the characterization of wheat cultivars. It shows tolerance to biotic and abiotic stresses to some extent. Hg1 (formerly named Hg) locus on chromosome 1AS controls glume pubescence in wheat. Its genetic analysis, fine-mapping and candidate gene analysis have been widely studied recently, however, the cloning of Hg1 has not yet been reported. Here, we conducted a GWAS between a dense panel of 171,103 SNPs and glume pubescence (Gp) in a durum wheat population of 145 lines, and further analyzed the candidate genes of Hg1 combined with the gene expression, functional annotation, and haplotype analysis. As a results, TRITD0Uv1G104670 (TdELD1-1A), encoding glycosyltransferase-like ELD1/KOBITO 1, was detected as the most promising candidate gene of Hg1 for glume pubescence in durum wheat. Our findings not only contribute to a deeper understanding of its cloning and functional validation but also underscore the significance of accurate genome sequences and annotations. Additionally, our study highlights the relevance of unanchored sequences in chrUn and the application of bioinformatics analysis for gene discovery in durum wheat.

Genome-Wide Identification of Calmodulin-Binding Protein 60 Gene Family and the Function of GhCBP60B in Cotton Growth and Development and Abiotic Stress Response

The calmodulin-binding protein 60 (CBP60) family is a gene family unique to plants, and its members play a crucial role in plant defense responses to pathogens and growth and development. Considering that cotton is the primary source of natural cotton textile fiber, the functional study of its CBP60 gene family members is critical. In this research, we successfully identified 162 CBP60 members from the genomes of 21 species. Of these, 72 members were found in four cotton species, divided into four clades. To understand the function of GhCBP60B in cotton in depth, we conducted a detailed analysis of its sequence, structure, cis-acting elements, and expression patterns. Research results show that GhCBP60B is located in the nucleus and plays a crucial role in cotton growth and development and response to salt and drought stress. After using VIGS (virus-induced gene silencing) technology to conduct gene silencing experiments, we found that the plants silenced by GhCBP60B showed dwarf plants and shortened stem nodes, and the expression of related immune genes also changed. In further abiotic stress treatment experiments, we found that GhCBP60B-silenced plants were more sensitive to drought and salt stress, and their POD (peroxidase) activity was also significantly reduced. These results imply the vital role of GhCBP60B in cotton, especially in regulating plant responses to drought and salt stress. This study systematically analyzed CBP60 gene family members through bioinformatics methods and explored in depth the biological function of GhCBP60B in cotton. These research results lay a solid foundation for the future use of the GhCBP60B gene to improve cotton plant type and its drought and salt resistance.

Characterization and expression analysis of GLABRA3 (GL3) genes in cotton: insights into trichome development and hormonal regulation

BACKGROUND

GLABRA3 (GL3) and ENHANCER OF GLABRA3 (EGL3) genes encode a typical helix-loop-helix (bHLH) transcription factors that primarily regulate trichome branching and root hair development, DNA endoreduplication, trichoblast size, and stomatal formation. The functions of GL3 genes in cotton crop have been poorly characterized. In this study, we performed comprehensive genome-wide scans for GL3 and EGL3 homologs to enhance our comprehension of their potential roles in trichome and fiber development in cotton crop.

METHODS AND RESULTS

Our findings paraded that Gossypium hirsutum and G. barbadense have 6 GL3s each, unevenly distributed on 4 chromosomes whereas, G. arboreum, and G. raimondii have 3 GL3s each, unevenly distributed on 2 chromosomes. Gh_A08G2088 and Gb_A09G2187, despite having the same bHLH domain as the other GL3 genes, were excluded due to remarkable short sequences and limited number of motifs, indicating a lack of potential functional activity. The phylogenetic analysis categorized remaining 16 GL3s into three subfamilies (Group I-III) closely related to A. thaliana. The 16 GL3s have complete bHLH domain, encompassing 590-631 amino acids, with molecular weights (MWs) ranging from 65.92 to 71.36 kDa. Within each subfamily GL3s depicted shared similar gene structures and motifs, indicating conserved characteristics within respective groups. Promoter region analysis revealed 27 cis-acting elements, these elements were responsive to salicylic acid, abscisic acid (ABA), methyl jasmonate (MeJA), and gibberellin. The expression of GL3 genes was analyzed across 12 tissues in both G. barbadense and G. hirsutum using the publicly available RNA-seq data. Among GL3s, Gb_D11G0219, Gb_D11G0214, and Gb_D08G2182, were identified as relatively highly expressed across different tissues, consequently selected for hormone treatment and expression validation in G. barbadense. RT-qPCR results demonstrated significant alterations in the expression levels of Gb_D11G0219 and Gb_D11G0214 following MeJA, GA, and ABA treatment. Subcellular localization prediction revealed that most GL3 proteins were predominantly expressed in the nucleus, while a few were localized in the cytoplasm and chloroplasts.

CONCLUSIONS

In summary, this study lays the foundation for subsequent functional validation of GL3 genes by identifying hormonal regulation patterns and probable sites of action in cotton trichome formation and fiber development. The results stipulate a rationale to elucidate the roles and regulatory mechanisms of GL3 genes in the intricate process of cotton fibre and trichome development.

The GaKAN2, a KANADI transcription factor, modulates stem trichomes in Gossypium arboreum

KEY MESSAGE

GaKAN2, a member of the KANADI family, was found to be widely expressed in the cotton tissues and regulates trichome development through complex pathways. Cotton trichomes are believed to be the defense barrier against insect pests. Cotton fiber and trichomes are single-cell epidermal extensions with shared regulatory mechanisms. Despite several studies underlying mechanism of trichome development remains elusive. The KANADI is one of the key transcription factors (TFs) family, regulating Arabidopsis trichomes growth. However, the function of KANADI genes in cotton remains unknown. In the current study genome-wide scanning, transcriptomic analysis, gene silencing, subcellular localization, and yeast two-hybrid techniques were employed to decipher the function of KANADI TFs family genes in cotton crop. A total of 7 GaKAN genes were found in the Gossypium arboreum. Transcriptomic data revealed that these genes were significantly expressed in stem and root. Moreover, GaKAN2 was widely expressed in other tissues also. Subsequently, we selected GaKAN2 to validate the function of KANADI genes. Silencing of GaKAN2 resulted in a 24.99% decrease in single-cell trichomes and an 11.33% reduction in internodal distance, indicating its potential role in regulating trichomes and plant growth. RNA-Seq analysis elucidated that GaSuS and GaERS were the downstream genes of GaKAN2. The transcriptional activation and similarity in silencing phenotype between GaKAN2 and GaERS suggested that GaKAN2 regulates trichomes development through GaERS. Moreover, KEGG analysis revealed that a significant number of genes were enriched in the biosynthesis of secondary metabolites and plant hormone signal transduction pathways, thereby suggesting that GaKAN2 regulates the stem trichomes and plant growth. The GFP subcellular localization and yeast transcriptional activation analysis elucidated that GaKAN2 was located in the nucleus and capable of regulating the transcription of downstream genes. This study elucidated the function and characteristics of the KANADI gene family in cotton, providing a fundamental basis for further research on GaKAN2 gene in cotton plant trichomes and plant developmental processes.

Matrix-bound nanovesicle-associated IL-33 supports functional recovery after skeletal muscle injury by initiating a pro-regenerative macrophage phenotypic transition

Injuries to skeletal muscle are among the most common injuries in civilian and military populations, accounting for nearly 60% of extremity injuries. The standard of care for severe extremity injury has been focused upon limb salvage procedures and the utilization of tissue grafts or orthotics in conjunction with rehabilitation to avoid amputation. Nonetheless, many patients have persistent strength and functional deficits that permanently impact their quality of life. Preclinical and clinical studies have shown that partial restoration of functional skeletal muscle tissue following injury can be achieved by the implantation of a biologic scaffold composed of extracellular matrix (ECM). These favorable outcomes are mediated, at least in part, through local immunomodulation. The mechanisms underlying this immunomodulatory effect, however, are poorly understood. The present study investigates a potential mechanistic driver of the immunomodulatory effects; specifically, the effect of selected ECM components upon inflammation resolution and repair. Results show that the host response to skeletal muscle injury is profoundly altered and functional recovery decreased in il33-/- mice compared to age- and sex-matched wildtype counterparts by 14 days post-injury. Results also show that IL-33, contained within matrix-bound nanovesicles (MBV), supports skeletal muscle regeneration by regulating local macrophage activation toward a pro-remodeling phenotype via canonical and non-canonical pathways to improve functional recovery from injury compared to untreated il33-/- counterparts. Taken together, these data suggest that MBV and their associated IL-33 cargo represent a novel homeostatic signaling mechanism that contributes to skeletal muscle repair.

Alternative Splicing during Fiber Development in G. hirsutum

Cotton is a valuable cash crop in many countries. Cotton fiber is a trichome that develops from a single epidermal cell and serves as an excellent model for understanding cell differentiation and other life processes. Alternative splicing (AS) of genes is a common post-transcriptional regulatory process in plants that is essential for plant growth and development. The process of AS during cotton fiber formation, on the other hand, is mainly unknown. A substantial number of multi-exon genes were discovered to be alternatively spliced during cotton fiber formation in this study, accounting for 23.31% of the total number of genes in Gossypium hirsutum. Retention intron (RI) is not necessarily the most common AS type, indicating that AS genes and processes during fiber development are very temporal and tissue-specific. When compared to fiber samples, AS is more prevalent at the fiber initiation stages and in the ovule, indicating that development stages and tissues use different AS strategies. Genes involved in fiber development have gone through stage-specific AS, demonstrating that AS regulates cotton fiber development. Furthermore, AS can be regulated by trans-regulation elements such as splicing factor and cis-regulation elements such as gene length, exon numbers, and GC content, particularly at exon-intron junction sites. Our findings also suggest that increased DNA methylation may aid in the efficiency of AS, and that gene body methylation is key in AS control. Finally, our research will provide useful information about the roles of AS during the cotton fiber development process.

Status and prospects of genome-wide association studies in cotton

Over the last two decades, the use of high-density SNP arrays and DNA sequencing have allowed scientists to uncover the majority of the genotypic space for various crops, including cotton. Genome-wide association study (GWAS) links the dots between a phenotype and its underlying genetics across the genomes of populations. It was first developed and applied in the field of human disease genetics. Many areas of crop research have incorporated GWAS in plants and considerable literature has been published in the recent decade. Here we will provide a comprehensive review of GWAS studies in cotton crop, which includes case studies on biotic resistance, abiotic tolerance, fiber yield and quality traits, current status, prospects, bottlenecks of GWAS and finally, thought-provoking question. This review will serve as a catalog of GWAS in cotton and suggest new frontiers of the cotton crop to be studied with this important tool.

Expression Patterns Divergence of Reciprocal F1 Hybrids Between Gossypium hirsutum and Gossypium barbadense Reveals Overdominance Mediating Interspecific Biomass Heterosis

Hybrid breeding has provided an impetus to the process and achievement of a higher yield and quality of crops. Interspecific hybridization is critical for resolving parental genetic diversity bottleneck problems. The reciprocal interspecific hybrids and their parents (Gossypium hirsutum and Gossypium barbadense) have been applied in this study to elucidate the transcription regulatory mechanism of early biomass heterosis. Phenotypically, the seed biomass, plant height over parent heterosis, leaf area over parent heterosis, and fresh and dry biomass were found to be significantly higher in hybrids than in parents. Analysis of leaf areas revealed that the one-leaf stage exhibits the most significant performance in initial vegetative growth vigor and larger leaves in hybrids, increasing the synthesis of photosynthesis compounds and enhancing photosynthesis compound synthesis. Comparative transcriptome analysis showed that transgressive down-regulation (TDR) is the main gene expression pattern in the hybrids (G. hirsutum × G. barbadense, HB), and it was found that the genes of photosystem I and Adenosine triphosphate (ATP)-binding may promote early growth vigor. Transgressive up-regulation (TUR) is the major primary gene expression pattern in the hybrids (G. barbadense × G. hirsutum, BH), and photosystem II-related genes mediated the performance of early biomass heterosis. The above results demonstrated that overdominance mediates biomass heterosis in interspecific hybrid cotton and the supervisory mechanism divergence of hybrids with different females. Photosynthesis and other metabolic process are jointly involved in controlling early biomass heterosis in interspecific hybrid cotton. The expression pattern data of transcriptome sequencing were supported using the qRT-PCR analysis. Our findings could be useful in theoretical and practical studies of early interspecific biomass heterosis, and the results provide potential resources for the theoretical and applied research on early interspecific biomass heterosis.

IL-1β induces changes in expression of core circadian clock components PER2 and BMAL1 in primary human chondrocytes through the NMDA receptor/CREB and NF-κB signalling pathways

The circadian clock is a specialised cell signalling circuit present in almost all cells. It controls the timing of key cell activities such as proliferation and differentiation. In osteoarthritis, expression of two components of the circadian clock, BMAL1 and PER2 is altered in chondrocytes and this change has been causally linked with the increase in proliferation and altered chondrocyte differentiation in disease. IL-1β, an inflammatory cytokine abundant in OA joints, has previously been shown to induce changes in BMAL1 and PER2 expression in chondrocytes. The purpose of this study is to identify the mechanism involved. We found IL-1β treatment of primary human chondrocytes led to activation of NMDA receptors as evidenced by an increase in phosphorylation of GluN1 and an increase in intracellular calcium which was blocked by the NMDAR antagonist MK801. Levels of phosphorylated CREB were also elevated in IL-1β treated cells and this effect was blocked by co-treatment of cells with IL-1β and the NMDAR antagonist MK-801. Knockdown of CREB or inhibition of CREB activity prevented the IL-1β induced increase in PER2 expression in chondrocytes but had no effect on BMAL1. Phosphorylated p65 levels were elevated in IL-1β treated chondrocytes indicating increased NF-κB activation. Inhibition of NF-κB activity prevented the IL-1β induced reduction in BMAL1 expression and partially mitigated the IL-1β induced increase in PER2 expression in chondrocytes. These data indicate that the NMDAR/CREB and NF-κB signalling pathways regulate the core circadian clock components PER2 and BMAL1 in chondrocytes. Given that changes in expression of these clock components have been observed in a wide range of diseases, these findings may be broadly relevant for understanding the mechanism leading to circadian clock changes in pathology.

Genome-Wide Analysis of SRNF Genes in Gossypium hirsutum Reveals the Role of GhSRNF18 in Primary Root Growth

Root systems are instrumental for water and nutrient uptake and the anchorage of plants in the soil. Root regulating GL2-interacting repressors (GIRs) contain a Short RING-like Zinc-Finger (SRNF) domain, but there has been no comprehensive characterization about this gene family in any plant species. Here, we renamed the GIR-like proteins as SRNF proteins due to their conserved domain and identified 140 SRNF genes from 16 plant species including 24 GhSRNF genes in Gossypium hirsutum. Phylogenetic analysis of the SRNFs revealed both similarities and divergences between five subfamilies. Notably, synteny analysis revealed that polyploidization and whole-genome duplication contribute to the expansion of the GhSRNF gene family. Various cis-acting regulatory elements were shown to be pertinent to light, phytohormone, defense responsive, and meristem regulation. Furthermore, GhSRNF2/15 were predominantly expressed in root, whereas the expression of GhSRNF18 is positively correlated with the primary root (PR) length in G. hirsutum, quantified by quantitative real-time PCR (qRT-PCR). Over-expression of GhSRNF18 in Arabidopsis and virus-induced gene silencing (VIGS) of GhSRNF18 in G. hirsutum has revealed the role of GhSRNF18 in PR growth. The over-expression of GhSRNF18 in Arabidopsis resulted in an increase of meristematic activities and auxin accumulations in PRs, which were consistent with the transcriptomic data. Our results suggested that GhSRNF18 positively regulates PR growth. This study increased our understanding of the SRNF gene family in plants and provided a novel rationale for the further investigation of cotton root morphogenesis regulated by the GhSRNFs.

Differentiation in the genetic basis of stem trichome development between cultivated tetraploid cotton species

BACKGROUND

Cotton stem trichomes and seed fibers are each single celled structures formed by protrusions of epidermal cells, and were found sharing the overlapping molecular mechanism. Compared with fibers, cotton stem trichomes are more easily observed, but the molecular mechanisms underlying their development are still poorly understood.

RESULTS

In this study, Gossypium hirsutum (Gh) and G. barbadense (Gb) were found to differ greatly in percentages of varieties/accessions with glabrous stems and in trichome density, length, and number per trichopore. Gh varieties normally had long singular and clustered trichomes, while Gb varieties had short clustered trichomes. Genetic mapping using five F₂ populations from crosses between glabrous varieties and those with different types of stem trichomes revealed that much variation among stem trichome phenotypes could be accounted for by different combinations of genes/alleles on Chr. 06 and Chr. 24. The twenty- six F₁ generations from crosses between varieties with different types of trichomes had varied phenotypes, further suggesting that the trichomes of tetraploid cotton were controlled by different genes/alleles. Compared to modern varieties, a greater proportion of Gh wild accessions were glabrous or had shorter and denser trichomes; whereas a smaller proportion of Gb primitive accessions had glabrous stems. A close correlation between fuzz fiber number and stem trichome density was observed in both Gh and Gb primitive accessions and modern varieties.

CONCLUSION

Based on these findings, we hypothesize that stem trichomes evolved in parallel with seed fibers during the domestication of cultivated tetraploid cotton. In addition, the current results illustrated that stem trichome can be used as a morphological index of fiber quality in cotton conventional breeding.

GaHD1, a candidate gene for the Gossypium arboreum SMA-4 mutant, promotes trichome and fiber initiation by cellular H2O2 and Ca2+ signals

Cotton fibers are initiated from the epidermal cells of the ovule before or on the day of anthesis. Gossypium arboreum SMA-4 mutant contains recessive mutation (sma-4(ha)) and has the phenotypes of fibreless seeds and glabrous stems. In this study, fine mapping and alternative splicing analysis indicated a nucleotide substitution (AG → AC) at splicing site in a homeodomain-leucine zipper IV family gene (GaHD1) might cause gene A3S (Alternative 3' splicing) mistake, suggested that GaHD1 was the candidate gene of sma-4(ha). Many genes related to the fiber initiation are identified to be differentially expressed in the mutant which could result in the blocked fiber initiation signals such as H₂O_2, or Ca in the mutant. Further comparative physiological analysis of H₂O₂ production and Ca²⁺ flux in the SMA-4 and wide type cotton confirmed that H₂O₂ and Ca were important fiber initiation signals and regulated by GaHD1. The in vitro ovule culture of the mutant with hormones recovered the fibered phenotype coupled with the restoration of these signals. Overexpressing of GaHD1 in Arabidopsis increased trichome densities on the sepal, leaf, and stem tissues while transient silencing of the GaHD1 gene in G. arboreum reduced the trichome densities. These phenotypes indicated that GaHD1 is the candidate gene of SMA-4 with a crucial role in acting upstream molecular switch of signal transductions for cotton trichome and fiber initiations.

The Ligon lintless -2 Short Fiber Mutation Is Located within a Terminal Deletion of Chromosome 18 in Cotton

Extreme elongation distinguishes about one-fourth of cotton (Gossypium sp.) seed epidermal cells as "lint" fibers, useful for the textile industry, from "fuzz" fibers (<5 mm). Ligon lintless-2 (Li ₂ ), a dominant mutation that results in no lint fiber but normal fuzz fiber, offers insight into pathways and mechanisms that differentiate spinnable cotton from its progenitors. A genetic map developed using 1,545 F2 plants showed that marker CISP15 was 0.4 cM from Li ₂ , and "dominant" simple sequence repeat (SSR) markers (i.e. with null alleles in the Li ₂ genotype) SSR7 and SSR18 showed complete linkage with Li ₂ Nonrandom distribution of markers with null alleles suggests that the Li ₂ phenotype results from a 176- to 221-kb deletion of the terminal region of chromosome 18 that may have been masked in prior pooled-sample mapping strategies. The deletion includes 10 genes with putative roles in fiber development. Two Glycosyltransferase Family 1 genes showed striking expression differences during elongation of wild-type versus Li ₂ fiber, and virus-induced silencing of these genes in the wild type induced Li ₂ -like phenotypes. Further, at least 7 of the 10 putative fiber development genes in the deletion region showed higher expression in the wild type than in Li ₂ mutants during fiber development stages, suggesting coordinated regulation of processes in cell wall development and cell elongation, consistent with the hypothesis that some fiber-related quantitative trait loci comprise closely spaced groups of functionally diverse but coordinately regulated genes.

Preferential insertion of a Ty1 LTR-retrotransposon into the A sub-genome's HD1 gene significantly correlated with the reduction in stem trichomes of tetraploid cotton

Stem trichomes and seed fibers originate from epidermal cells and partially share a regulatory pathway at the molecular level. In Gossypium barbadense, two insertions of a Ty1 long-terminal repeat-retrotransposon [transposable element TE1 and TE2] in a homeodomain-leucine zipper gene (HD1) result in glabrous stems. The primers used to identify the TE insertions in G. barbadense were applied to screen for the same events in 81 modern G. hirsutum varieties and 31 wild races. Three wild races were found carrying the same TEs as G. barbadense. However, the TE insertions in two of these wild races occurred at different sites (4th exon), therefore, named TE3, while the TE in the other wild race occurred at the same site as TE2. An RNA sequencing and qRT-PCR analysis indicated that the loss of HD1 function was caused by the TE insertion. Genetic mapping revealed a strong association between glabrous stems and TE3 insertions, confirming that HD1 is a critical gene for stem trichome initiation in G. hirsutum, as in G. barbadense. Using the long-terminal repeat sequence as a query to search against the Texas Marker-1 reference genome sequence, we found that the TE occurred after tetraploid cotton formation and evolved at different rates in G. hirsutum and G. barbadense. Interestingly, at least three independent insertion events of the same retrotransposon occurred preferentially in the A sub-genome's HD1 gene, but not in the D sub-genome of G. hirsutum or G. barbadense, suggesting that an unknown TE insertion mechanism and resultant gene function changes may have hastened cotton speciation.

Identification and Expression Profiling of the Regulator of Chromosome Condensation 1 (RCC1) Gene Family in Gossypium Hirsutum L. under Abiotic Stress and Hormone Treatments

The regulator of chromosome condensation 1 (RCC1) is the nucleotide exchange factor for a GTPase called the Ras-related nuclear protein, and it is important for nucleo-plasmic transport, mitosis, nuclear membrane assembly, and control of chromatin agglutination during the S phase of mitosis in animals. In plants, RCC1 molecules act mainly as regulating factors for a series of downstream genes during biological processes such as the ultraviolet-B radiation (UV-B) response and cold tolerance. In this study, 56 genes were identified in upland cotton by searching the associated reference genomes. The genes were found to be unevenly distributed on 26 chromosomes, except A06, A12, D03, and D12. Phylogenetic analysis by maximum-likelihood revealed that the genes were divided into five subgroups. The RCC1 genes within the same group shared similar exon/intron patterns and conserved motifs in their encoded proteins. Most genes of the RCC1 family are expressed differently under various hormone treatments and are negatively controlled by salt stress. Gh_A05G3028 and Gh_D10G2310, which encode two proteins located in the nucleus, were strongly induced under salt treatment, while mutants of their homoeologous gene (UVR8) in Arabidopsis and VIGS (virus induced gene silencing) lines of the two genes above in G. hirsutum exhibited a salt-sensitive phenotype indicating their potential role in salt resistance in cotton. These results provide valuable reference data for further study of RCC1 genes in cotton.

A genetic system on chromosome arm 1BL of wild emmer causes distorted segregation in common wheat

Nonrandom segregation ratios of alleles 'segregation distortion' can have a striking impact on transmission genetics, and with widespread availability of genetic markers has been shown to be a frequent phenomenon. To investigate the possible effect of genetic interaction on segregation distortion and genetic map construction, the segregation and mapping of genetic markers locatedon wheat chromosomes 1A and 1B were followed in four recombinant substitution line (RSL) populations, produced using four chromosome-arm substitution lines (CASLs 1AS, 1AL, 1BS and 1BL) of wild emmer (Triticum turgidum var. dicoccoides, accession TTD140) in the background of the common wheat (T. aestivum) cultivar Bethlehem (BLH), each crossed to BLH itself. Using these four RSL populations, four genetic maps of chromosome 1 arms were constructed. A total of 22 genetic markers representing 19 loci were assigned to chromosome 1A, and 32 markers representing 30 loci were assigned to 1B. For chromosome 1B, two linkage maps were also constructed using RFLP data of an F₂ population derived from the same cross combination as the RSLs. The RSL and F₂ maps varied in genetic distances, but showed the same linear order of DNA markers. Segregation analysis revealed strong selection against BLH alleles on chromosome 1B, skewing the allelic frequency distribution in favour of TTD in both F₂ and RSL populations at all marker loci. On the contrary, strong selection against TTD alleles on chromosome 1A was detected for some loci in the BLH × CASL1AL RSLs, and their distribution was significantly skewed to BLH. F₂ populations always showed more segregation distortion than the corresponding RSLs. More markers near the region of chromosome 1B shared by both CASL1BS and 1BL (∼55 cM on chromosome 1B across the centromere) showed significantly distorted segregation in the BLH × CASL1BL population than in thecorresponding BLH × CASL1BS populations. Six markers located on chromosome 1A region shared by CASL1AS and 1AL showed significantly distorted segregation in 1AL-RSL, while no marker showed distorted segregation in 1AS-RSL. These results indicated that genetic factor(s) in the centromere region cause the distorted segregation of genetic markers on wheat chromosome 1B.

Altered N-methyl D-aspartate receptor subunit expression causes changes to the circadian clock and cell phenotype in osteoarthritic chondrocytes

The chondrocyte circadian clock is altered in osteoarthritis. This change is implicated in the disease-associated changes in chondrocyte phenotype and cartilage loss. Why the clock is changed is unknown. N-methyl-D-aspartate receptors (NMDAR) are critical for regulating the hypothalamic clock. Chondrocytes also express NMDAR and the type of NMDAR subunits expressed changes in osteoarthritis.

OBJECTIVE

To determine if NMDAR regulate the chondrocyte clock and phenotype.

DESIGN

Chondrocytes isolated from macroscopically-normal (MN) and osteoarthritic human cartilage were treated with NMDAR antagonists or transfected with GRIN2A or GRIN2B-targetting siRNA. H5 chondrocytes were transfected with GluN2B-expression plasmids. Clock genes and chondrocyte phenotypic markers were measured by RT-qPCR.

RESULTS

PER2 amplitude was higher and BMAL1 amplitude lower in osteoarthritic compared to MN chondrocytes. In osteoarthritic chondrocytes, NMDAR inhibition restored PER2 and BMAL1 expression to levels similar to MN chondrocytes, and resulted in reduced MMP13 and COL10A1. Paradoxically, NMDAR inhibition in MN chondrocytes resulted in increased PER2, decreased BMAL1 and increased MMP13 and COL10A1. Osteoarthritic, but not MN chondrocytes expressed GluN2B NMDAR subunits. GluN2B knockdown in osteoarthritic chondrocytes restored expression of circadian clock components and phenotypic markers to levels similar to MN chondrocytes. Ectopic expression of GluN2B resulted in reduced BMAL1, increased PER2 and altered SOX9, RUNX2 and MMP13 expression. Knockdown of PER2 mitigated the effects of GluN2B on SOX9 and MMP13.

CONCLUSIONS

NMDAR regulate the chondrocyte clock and phenotype suggesting NMDAR may also regulate clocks in other peripheral tissues. GluN2B expression in osteoarthritis may contribute to pathology by altering the chondrocyte clock.

Obesity, osteoarthritis and genetic risk: The rs182052 polymorphism in the ADIPOQ gene is potentially associated with risk of knee osteoarthritis

Objectives

Given the function of adiponectin (ADIPOQ) on the inflammatory condition of obesity and osteoarthritis (OA), we hypothesized that the ADIPOQ gene might be a candidate gene for a marker of susceptibility to OA.

Methods

We systematically screened three tagging polymorphisms (rs182052, rs2082940 and rs6773957) in the ADIPOQ gene, and evaluated the association between the genetic variants and OA risk in a case-controlled study that included 196 OA patients and 442 controls in a northern Chinese population. Genotyping was performed using the Sequenom MassARRAY iPLEX platform.

Results

The single nucleotide polymorphism (SNP) rs182052 was found to be potentially associated with knee OA risk (additive model: odds ratio = 1.38; 95% confidence interval 1.07 to 1.76; p = 0.012). Furthermore, a non-significant association was observed for rs182052 and body mass index with regard to OA risk in interaction analyses (p = 0.063). Similarly, no significant interaction was detected for rs182052 and age with regard to OA risk (p = 0.614).

Conclusion

These findings suggest that the SNP rs182052 in the ADIPOQ gene may potentially modify individual susceptibility to knee OA in the Chinese population. Further studies are warranted to investigate our findings in more depth.

Cite this article: L. Jiang, X. Zhu, J. Rong, B. Xing, S. Wang, A. Liu, M. Chu, G. Huang. Obesity, osteoarthritis and genetic risk: The rs182052 polymorphism in the ADIPOQ gene is potentially associated with risk of knee osteoarthritis. Bone Joint Res 2018;7:494–500. DOI: 10.1302/2046-3758.77.BJR-2017-0274.R1.

Divergence and evolution of cotton bHLH proteins from diploid to allotetraploid

BACKGROUND

Polyploidy is considered a major driving force in genome expansion, yielding duplicated genes whose expression may be conserved or divergence as a consequence of polyploidization.

RESULTS

We compared the genome sequences of tetraploid cotton (Gossypium hirsutum) and its two diploid progenitors, G. arboreum and G. raimondii, and found that the bHLH genes were conserved over the polyploidization. Oppositely, the expression of the homeolgous gene pairs was diversified. The biased homeologous proportion for bHLH family is significantly higher (64.6%) than the genome wide homeologous expression bias (40%). Compared with cacao (T. cacao), orthologous genes only accounted for a small proportion (41.7%) of whole cotton bHLHs family. The further Ks analysis indicated that bHLH genes underwent at least two distinct episodes of whole genome duplication: a recent duplication (1.0-60.0 million years ago, MYA, 0.005 < Ks < 0.312) and an old duplication (> 60.0 MYA, 0.312 < Ks < 3.0). The old duplication event might have played a key role in the expansion of the bHLH family. Both recent and old duplicated pairs (68.8%) showed a divergent expression profile, indicating specialized functions. The expression diversification of the duplicated genes suggested it might be a universal feature of the long-term evolution of cotton.

CONCLUSIONS

Overview of cotton bHLH proteins indicated a conserved and divergent evolution from diploids to allotetraploid. Our results provided an excellent example for studying the long-term evolution of polyploidy.

Skin resistance to UVB-induced oxidative stress and hyperpigmentation by the topical use of Lactobacillus helveticus NS8-fermented milk supernatant

AIMS

In this study, we investigated the preventive properties of the supernatant of Lactobacillus helveticus NS8-fermented milk (NS8-FS) against UV light-induced skin oxidative damage and hyperpigmentation.

METHODS AND RESULTS

NS8-FS exhibited significant radical scavenging activity in tests with ABST⁺ and DPPH scavenging methods, and as well strongly inhibited 3-morpholinosydnonimine (Sin-1)-induced ROS generation in HaCat keratinocytes. Unexpectedly, NS8-FS was found to inhibit melanin production in B16F10 melanoma cells and to exhibit inhibitory effects both to the enzymatic activity of tyrosinase (TYR) and the expression of proteins required for melanin synthesis. In SKH-1 hairless mice, topical application of NS8-FS alleviated UVB-induced skin photodamage, including the improvement of the appearance of epidermal thickness, transepidennal water loss and lipid peroxidation levels. In the tanning guinea pig model, the whitening effect of NS8-FS was demonstrated using Masson-Fontana staining and TYR staining. Furthermore, NS8-FS was shown to stimulate the nuclear translocation and activation of the Nrf2 protein, along with recovery of antioxidant enzyme activities.

CONCLUSION

NS8-FS exhibits the protective capacities against UV light-induced skin oxidative damage and hyperpigmentation.

SIGNIFICANCE AND IMPACT OF THE STUDY

Our findings indicate the potential of cell-free fermented products of lactic acid bacteria in topical photoprotection.

The rs4238326 polymorphism in ALDH1A2 gene potentially associated with non-post traumatic knee osteoarthritis susceptibility: a two-stage population-based study

OBJECTIVE

A recent genome-wide association study reported significant associations of genetic variants within the ALDH1A2 gene with osteoarthritis (OA) of the hand in European populations. However, these findings have not been well generalized to other joints, or to other populations.

METHODS

We performed a two-stage population-based case-control study including 196 non-post traumatic knee OA cases and 442 controls in the first stage and independent 143 non-post traumatic knee OA cases and 238 controls in the second stage in a Chinese population by genotyping eight tagging polymorphisms in ALDH1A2.

RESULTS

In the first stage, the single nucleotide polymorphism (SNP) rs4238326 was found to be potentially associated with knee OA risk (additive model: odds ratio [OR] = 0.70; 95% confidence interval [95% CI] = 0.49-1.01; P = 0.055), which was further confirmed in the second stage with similar effect (additive model: OR = 0.60; 95% CI = 0.38-0.95; P = 0.029). After combining the two stages, we found that the variant C allele of rs4238326 was probably associated with decreased risk of knee OA (additive model: OR = 0.65; 95% CI = 0.49-0.86; P = 0.003). Furthermore, interaction analyses showed that rs4238326 interacted multiplicatively with age to contribute to knee OA risk (interaction P = 0.041).

CONCLUSIONS

These findings indicate that the SNP rs4238326 in ALDH1A2 gene may potentially modify individual susceptibility to knee OA in the Chinese population. Beyond that, further studies are warranted to validate and extend our findings, and future functional studies are required to clarify the possible mechanisms.

Clinical interrogation and application of super-selective intracranial artery infusion chemotherapy for lung cancer patients with brain metastases

OBJECTIVE

The purpose of this study was to evaluate the clinical efficacy of super-selective intracranial artery infusion chemotherapy and to determine correlated prognostic parameters for advanced lung cancer patients with brain metastases.

PATIENTS AND METHODS

Fifty-four lung cancer patients with brain metastasis who had no previous treatment were enrolled for the study. These patients received super-selective intracranial artery infusion chemotherapy, as well as arterial infusion chemotherapy for primary and metastatic lesions. The procedure was performed once every 4 weeks. Patients were monitored to evaluate short-term clinical outcomes 4 weeks after the first 2 treatments, and follow-up visits performed every 4 weeks after the first 4 treatments until the appearance of disease progression or intolerable toxicity.

RESULTS

All 54 cases were treated at least 4 times. The overall response rate was 55.56% (30/54), and the disease control rate was 85.19% (46/54). The median overall survival was 7 months, with a 95% confidence interval (CI) of 5.87-8.13 months, and the median progression-free survival was 4 months, with a 95% CI of 3.20-4.80 months. The 6-month survival rate and 1-year survival rate were 81.48% (44/54) and 18.52% (10/54), respectively.

CONCLUSION

Super-selective intracranial artery infusion chemotherapy provides a clinically efficacious avenue of treatment for lung cancer patients with brain metastases. Pathological classification, Karnofsky performance status, and extracranial metastases may serve as reliable prognostic parameters in determining the clinical outcomes for lung cancer patients with brain metastases.

No association of the single nucleotide polymorphism rs8044769 in the fat mass and obesity-associated gene with knee osteoarthritis risk and body mass index: A population-based study in China

OBJECTIVES

Previous genome-wide association studies (GWAS) have reported significant association of the single nucleotide polymorphism (SNP) rs8044769 in the fat mass and obesity-associated gene (FTO) with osteoarthritis (OA) risk in European populations. However, these findings have not been confirmed in Chinese populations.

METHODS

We systematically genotyped rs8044769 and evaluated the association between the genetic variants and OA risk in a case-controlled study including 196 OA cases and 442 controls in a northern Chinese population. Genotyping was performed using the Sequenom MassARRAY iPLEX platform.

RESULTS

We found that the variant T allele of rs8044769 showed no significant association of OA risk (p = 0.791), or association with body mass index (BMI) (pmeta = 0.786) in an additive genetic model. However, we detected a significant interaction between rs8044769 genotypes and BMI on OA risk (p = 0.037), as well as a borderline interaction between rs8044769 genotypes and age on OA risk (p = 0.062).

CONCLUSIONS

Our findings indicate that rs8044769 in the FTO gene may not modify individual susceptibility to OA or increased BMI in the Chinese population. Further studies are warranted to validate and extend our findings.Cite this article: Prof L. Jiang. No association of the single nucleotide polymorphism rs8044769 in the fat mass and obesity-associated gene with knee osteoarthritis risk and body mass index: A population-based study in China. Bone Joint Res 2016;5:169-174. DOI: 10.1302/2046-3758.55.2000589.

Comparative transmission genetics of introgressed chromatin in Gossypium (cotton) polyploids

PREMISE OF THE STUDY

Introgression is widely acknowledged as a potential source of valuable genetic variation, and growing effort is being invested in analysis of interspecific crosses conferring transgressive variation. Experimental backcross populations provide an opportunity to study transmission genetics following interspecific hybridization, identifying opportunities and constraints to introgressive crop improvement. The evolutionary consequences of introgression have been addressed at the theoretical level, however, issues related to levels and patterns of introgression among (plant) species remain inadequately explored, including such factors as polyploidization, subgenome interaction inhabiting a common nucleus, and the genomic distribution and linkage relationships of introgressant alleles.

METHODS

We analyze introgression into the polyploid Gossypium hirsutum (upland cotton) from its sister G. tomentosum and compare the level and pattern with that of G. barbadense representing a different clade tracing to the same polyploidization.

KEY RESULTS

Across the genome, recurrent backcrossing to Gossypium hirsutum yielded only one-third of the expected average frequency of the G. tomentosum allele, although one unusual region showed preferential introgression. Although a similar rate of introgression is found in the two subgenomes of polyploid (AtDt) G. hirsutum, a preponderance of multilocus interactions were largely within the Dt subgenome.

CONCLUSIONS

Skewed G. tomentosum chromatin transmission is polymorphic among two elite G. hirsutum genotypes, which suggests that genetic background may profoundly affect introgression of particular chromosomal regions. Only limited correspondence is found between G. hirsutum chromosomal regions that are intolerant to introgression from the two species, G. barbadense and G. tomentosum, concentrated near possible inversion polymorphisms. Complex transmission of introgressed chromatin highlights the challenges to utilization of exotic germplasm in crop improvement.

Mapping of Ppd-B1, a Major Candidate Gene for Late Heading on Wild Emmer Chromosome Arm 2BS and Assessment of Its Interactions with Early Heading QTLs on 3AL

Wheat heading date is an important agronomic trait determining maturation time and yield. A set of common wheat (Triticum aestivum var. Chinese Spring; CS)-wild emmer (T. turgidum L. subsp. dicoccoides (TDIC)) chromosome arm substitution lines (CASLs) was used to identify and allocate QTLs conferring late or early spike emergence by examining heading date. Genetic loci accelerating heading were found on TDIC chromosome arms 3AL and 7BS, while loci delaying heading were located on 4AL and 2BS. To map QTLs conferring late heading on 2BS, F2 populations derived from two cross combinations of CASL2BS × CS and CASL3AL × CASL2BS were developed and each planted at two times, constituting four F2 mapping populations. Heading date varied continuously among individuals of these four populations, suggesting quantitative characteristics. A genetic map of 2BS, consisting of 23 SSR and one single-stranded conformation polymorphism (SSCP) marker(s), was constructed using these F2 populations. This map spanned a genetic length of 53.2 cM with average marker density of 2.3 cM. The photoperiod-sensitivity gene Ppd-B1 was mapped to chromosome arm 2BS as a SSCP molecular marker, and was validated as tightly linked to a major QTL governing late heading of CASL2BS in all mapping populations. A significant dominance by additive effect of Ppd-B1 with the LUX gene located on 3AL was also detected. CS had more copies of Ppd-B1 than CASL2BS, implying that increased copy number could elevate the expression of Ppd-1 in CS, also increasing expression of LUX and FT genes and causing CS to have an earlier heading date than CASL2BS in long days.

Genetic fine mapping and candidate gene analysis of the Gossypium hirsutum Ligon lintless-1 (Li1) mutant on chromosome 22(D)

Ligon lintless-1 (Li1) is a Gossypium hirsutum mutant that is controlled by a dominant gene that arrests the development of cotton fiber after anthesis. Two F2 mapping populations were developed from mutant (Li1 × H7124) F1 plants in 2012 and 2013; each was composed of 142 and 1024 plants, respectively. Using these populations, Li1 was mapped to a 0.3-cM region in which nine single-strand conformation polymorphism markers co-segregated with the Li1 locus. In the published G. raimondii genome, these markers were mapped to a region of about 1.2 Mb (the Li1 region) and were separated by markers that flanked the Li1 locus in the genetic map, dividing the Li1 region into three segments. Thirty-six genes were annotated by the gene prediction software FGENESH (Softberry) in the Li1 region. Twelve genes were candidates of Li1, while the remaining 24 genes were identified as transposable elements, DNA/RNA polymerase superfamily or unknown function genes. Among the 12 candidate genes, those encoding ribosomal protein s10, actin protein, ATP synthase, and beta-tubulin 5 were the most-promising candidates of the Li1 mutant because the function of these genes is closely related to fiber development. High-throughput RNA sequencing and quantitative PCR revealed that these candidate genes had obvious differential gene expression between mutant and wild-type plants at the fiber elongation stage, strengthening the inference that they could be the most likely candidate gene of the Li1 mutant phenotype.

Molecular characterization of a transcriptionally active Ty1/copia-like retrotransposon in Gossypium

A transcriptionally active Ty1/copia -like retrotransposon was identified in the genome of Gossypium barbadense. The different heat activation of this element was observed in two tetraploid cotton species. Most retrotransposons from plants are transcriptionally silent, or activated under certain conditions. Only a small portion of elements are transcriptionally active under regular condition. A long terminal repeat (LTR) retrotransposon was isolated from the cultivated Sea Island cotton (H7124) genome during the investigation of the function of a homeodomain leucine zipper gene (HD1) in trichome growth. Insertion of this element in HD1 gene of At sub-genome was related to the trichomeless stem in Gossypium barbadense. The element, named as GBRE-1, had all features of a typical Ty1/copia retrotransposon and possessed high similarity to the members of ONSEN retrotransposon family. It was 4997 bp long, comprising a single 4110 bp open reading frame, which encoded 1369 amino acids including the conserved domains of gag and pol. The expression of GBRE-1 was detected under regular condition in G. barbadense and G. hirsutum, and its expression level was increased under heat-stress condition in G. hirsutum. Besides, its expression pattern was similar to that of the ONSEN retrotransposon. Abundant cis-regulatory motifs related to stress-response and transcriptional regulation were found in the LTR sequence. These results suggested that GBRE-1 was a transcriptionally active retrotransposon in Gossypium. To our knowledge, this is the first report of the isolation of a complete Ty1/copia-type retrotransposon with present-day transcriptional activity in cotton.

Natural antisense RNAs are involved in the regulation of CD45 expression in autoimmune diseases

CD45 is a transmembrane protein tyrosine phosphatase that is specifically expressed in hematopoietic cells and can initiate signal transduction via the dephosphorylation of tyrosine. Alternatively spliced transcript variants of this gene encode distinct isoforms, which indicate different functional states of CD45. Among these variants, CD45RO, which contains neither exon 4, 5, or 6, is over-expressed in lymphocytes in autoimmune diseases, including systemic lupus erythematosus, rheumatoid arthritis, multiple sclerosis, and type I diabetes. The CD45 RO serves as a marker of the immune response activity and lymphocyte development. Previous studies have indicated that exon splicing is generally correlated with local hypermethylated DNA and acetylated histone modification, while autoimmune diseases are commonly associated with global hypomethylation and histone deacetylation in lymphocytes. Thus, the question arises of how exons 4, 5, and 6 of CD45RO are excluded under the status of global DNA hypomethylation and histone deacetylation in these autoimmune diseases. On the basis of the analyses of the context sequence of CD45 and its natural antisense RNA in GenBank, we proposed that the long noncoding RNA encoded by the natural antisense gene of CD45 contributes to the expressional regulation of the CD45RO splicing variant via recruitment of DNA methyltransferase and histone modification modulators specific to the sense gene CD45; thus, it is associated with the over-expression of CD45RO and the functional regulation of lymphocytes in the pathogenic development of autoimmune diseases.

Gene expression of the p16(INK4a)-Rb and p19(Arf)-p53-p21(Cip/Waf1) signaling pathways in the regulation of hematopoietic stem cell aging by ginsenoside Rg1

The elucidation of the molecular mechanisms underlying the effects of traditional Chinese medicines in clinical practice is a key step toward their worldwide application, and this topic is currently a subject of intense research interest. Rg1, a component of ginsenoside, has recently been shown to perform several pharmacological functions; however, the underlying mechanisms of these effects remain unclear. In the present study, we investigated whether Rg1 has an anti-senescence effect on hematopoietic stem cells (HSCs) and the possible molecular mechanisms driving any effects. The results showed that Rg1 could effectively delay tert-butyl hydroperoxide (t-BHP)-induced senescence and inhibit gene expression in the p16(INK4a)-Rb and p19(Arf)-p53-p21(Cip/Waf1) signaling pathways in HSCs. Our study suggested that these two signaling pathways might be potential targets for elucidating the molecular mechanisms of the Rg1 anti-senescence effect.

Epidemiological comparison of clinical manifestations according to HLA-B*27 carrier status of Chinese ankylosing spondylitis patients

The aim of the study was to investigate and compare the clinical manifestations between HLA-B27(+) and HLA-B27(-) ankylosing spondylitis (AS) patients in order to obtain knowledge of the impact of HLA-B27 status on AS, and to inform clinical treatment. A nationwide epidemiological investigation was performed from November 2008 to October 2010. The demographic data and clinical characteristics, and the status of HLA-B27 were collected using questionnaires and laboratory assay, respectively. A total of 2144 patients (78.5% males and 78.4% HLA-B27(+) AS patients) participated in this study. The percentages of males, patients with family history, and involvement of lumbar spine, thoracic spine and hip joints, were observed to be significantly higher in the HLA-B27(+) AS patients than in their HLA-B27(-) AS peers.

Synthesis of new derivatives of copper complexes of Josiphos family ligands for applications in asymmetric catalysis

New derivatives of copper complexes of Josiphos family ligands have been prepared and studied in asymmetric catalytic addition of Grignard reagents to enones, enoates and aromatic ketones.

Gene expression profile analysis of Ligon lintless-1 (Li1) mutant reveals important genes and pathways in cotton leaf and fiber development

Ligon lintless-1 (Li1) is a monogenic dominant mutant of Gossypium hirsutum (upland cotton) with a phenotype of impaired vegetative growth and short lint fibers. Despite years of research involving genetic mapping and gene expression profile analysis of Li1 mutant ovule tissues, the gene remains uncloned and the underlying pathway of cotton fiber elongation is still unclear. In this study, we report the whole genome-level deep-sequencing analysis of leaf tissues of the Li1 mutant. Differentially expressed genes in leaf tissues of mutant versus wild-type (WT) plants are identified, and the underlying pathways and potential genes that control leaf and fiber development are inferred. The results show that transcription factors AS2, YABBY5, and KANDI-like are significantly differentially expressed in mutant tissues compared with WT ones. Interestingly, several fiber development-related genes are found in the downregulated gene list of the mutant leaf transcriptome. These genes include heat shock protein family, cytoskeleton arrangement, cell wall synthesis, energy, H2O2 metabolism-related genes, and WRKY transcription factors. This finding suggests that the genes are involved in leaf morphology determination and fiber elongation. The expression data are also compared with the previously published microarray data of Li1 ovule tissues. Comparative analysis of the ovule transcriptomes of Li1 and WT reveals that a number of pathways important for fiber elongation are enriched in the downregulated gene list at different fiber development stages (0, 6, 9, 12, 15, 18dpa). Differentially expressed genes identified in both leaf and fiber samples are aligned with cotton whole genome sequences and combined with the genetic fine mapping results to identify a list of candidate genes for Li1.