Identification and molecular evolution of the GLX genes in 21 plant species: a focus on the Gossypium hirsutum.

BACKGROUND: The glyoxalase system includes glyoxalase I (GLXI), glyoxalase II (GLXII) and glyoxalase III (GLXIII), which are responsible for methylglyoxal (MG) detoxification and involved in abiotic stress responses such as drought, salinity and heavy metal. RESULTS: In this study, a total of 620 GLX family genes were identified from 21 different plant species. The results of evolutionary analysis showed that GLX genes exist in all species from lower plants to higher plants, inferring that GLX genes might be important for plants, and GLXI and GLXII account for the majority. In addition, motif showed an expanding trend in the process of evolution. The analysis of cis-acting elements in 21 different plant species showed that the promoter region of the GLX genes were rich in phytohormones and biotic and abiotic stress-related elements, indicating that GLX genes can participate in a variety of life processes. In cotton, GLXs could be divided into two groups and most GLXIs distributed in group I, GLXIIs and GLXIIIs mainly belonged to group II, indicating that there are more similarities between GLXII and GLXIII in cotton evolution. The transcriptome data analysis and quantitative real-time PCR analysis (qRT-PCR) show that some members of GLX family would respond to high temperature treatment in G.hirsutum. The protein interaction network of GLXs in G.hirsutum implied that most members can participate in various life processes through protein interactions. CONCLUSIONS: The results elucidated the evolutionary history of GLX family genes in plants and lay the foundation for their functions analysis in cotton.

Identification of MYB gene family and functional analysis of GhMYB4 in cotton (Gossypium spp.).

Myeloblastosis (MYB) transcription factors (TFs) form a large gene family involved in a variety of biological processes in plants. Little is known about their roles in the development of cotton pigment glands. In this study, 646 MYB members were identified in Gossypium hirsutum genome and phylogenetic classification was analyzed. Evolution analysis revealed assymetric evolution of GhMYBs during polyploidization and sequence divergence of MYBs in G. hirustum was preferentially happend in D sub-genome. WGCNA (weighted gene co-expression network analysis) showed that four modules had potential relationship with gland development or gossypol biosynthesis in cotton. Eight differentially expressed GhMYB genes were identified by screening transcriptome data of three pairs of glanded and glandless cotton lines. Of these, four were selected as candidate genes for cotton pigment gland formation or gossypol biosynthesis by qRT-PCR assay. Silencing of GH_A11G1361 (GhMYB4) downregulated expression of multiple genes in gossypol biosynthesis pathway, indicating it could be involved in gossypol biosynthesis. The potential protein interaction network suggests that several MYBs may have indirect interaction with GhMYC2-like, a key regulator of pigment gland formation. Our study was the systematic analysis of MYB genes in cotton pigment gland development, providing candidate genes for further study on the roles of cotton MYB genes in pigment gland formation, gossypol biosynthesis and future crop plant improvement.

Interpretation of convolutional neural networks reveals crucial sequence features involving in transcription during fiber development.

BACKGROUND: Upland cotton provides the most natural fiber in the world. During fiber development, the quality and yield of fiber were influenced by gene transcription. Revealing sequence features related to transcription has a profound impact on cotton molecular breeding. We applied convolutional neural networks to predict gene expression status based on the sequences of gene transcription start regions. After that, a gradient-based interpretation and an N-adjusted kernel transformation were implemented to extract sequence features contributing to transcription. RESULTS: Our models had approximate 80% accuracies, and the area under the receiver operating characteristic curve reached over 0.85. Gradient-based interpretation revealed 5' untranslated region contributed to gene transcription. Furthermore, 6 DOF binding motifs and 4 transcription activator binding motifs were obtained by N-adjusted kernel-motif transformation from models in three developmental stages. Apart from 10 general motifs, 3 DOF5.1 genes were also detected. In silico analysis about these motifs' binding proteins implied their potential functions in fiber formation. Besides, we also found some novel motifs in plants as important sequence features for transcription. CONCLUSIONS: In conclusion, the N-adjusted kernel transformation method could interpret convolutional neural networks and reveal important sequence features related to transcription during fiber development. Potential functions of motifs interpreted from convolutional neural networks could be validated by further wet-lab experiments and applied in cotton molecular breeding.

Assessment of the predictive ability of the Johns Hopkins Fall Risk Assessment Tool (Chinese Version) in inpatient settings.

AIMS: This study was to assess the predictive ability of the Johns Hopkins Fall Risk Assessment Tool (Chinese Version) in inpatient settings. DESIGN: A case-control study. METHODS: This study was conducted in a tertiary hospital based on 2019 data. With a case-control design in a 1:2 ratio, the predictive ability of the Johns Hopkins Fall Risk Assessment Tool (Chinese Version) was determined by ROC curve. The best cut point was identified based on sensitivity, specificity, positive predict value and negative predict value. Conditional logistical regression analysis was conducted to test the predictive ability of each indicator. RESULTS: The study included 309 patients, with 103 in the case group and 206 in the control groups. Generally, the predictive ability was acceptable with the area under ROC curve value at 0.73 (95% CI: 0.67-0.79). Positive predict value and negative predict value performed best at the cut point of 13. Sensitivity at cut point 6 was much higher than that at cut point 13, though specificity was lower. Except for age, all indicators in the Johns Hopkins Fall Risk Assessment Tool (Chinese Version) demonstrated significant predictive ability as to occurrence of fall. CONCLUSION: The Johns Hopkins Fall Risk Assessment Tool (Chinese Version) is a reliable assessment instrument in the inpatient settings. IMPACT: This is the first study that evaluated the predictive ability of the Johns Hopkins Fall Risk Assessment Tool (Chinese version) in the inpatient settings, and proved that the instrument is reliable for assessing inpatient fall risks. Further studies could be carried out to assess the predict ability of Johns Hopkins Fall Risk Assessment Tool (Chinese version) among specific populations.

Molecular cloning and characterization of GhERF105, a gene contributing to the regulation of gland formation in upland cotton (Gossypium hirsutum L.).

BACKGROUND: Gossypium hirsutum L. (cotton) is one of the most economically important crops in the world due to its significant source of fiber, feed, foodstuff, oil and biofuel products. However, the utilization of cottonseed was limited due to the presence of small and darkly pigmented glands that contain large amounts of gossypol, which is toxic to human beings and non-ruminant animals. To date, some progress has been made in the pigment gland formation, but the underlying molecular mechanism of its formation was still unclear. RESULTS: In this study, we identified an AP2/ERF transcription factor named GhERF105 (GH_A12G2166), which was involved in the regulation of gland pigmentation by the comparative transcriptome analysis of the leaf of glanded and glandless plants. It encoded an ERF protein containing a converved AP2 domain which was localized in the nucleus with transcriptional activity, and showed the high expression in glanded cotton accessions that contained much gossypol. Virus-induced gene silencing (VIGS) against GhERF105 caused the dramatic reduction in the number of glands and significantly lowered levels of gossypol in cotton leaves. GhERF105 showed the patterns of spatiotemporal and inducible expression in the glanded plants. CONCLUSIONS: These results suggest that GhERF105 contributes to the pigment gland formation and gossypol biosynthesis in partial organs of glanded plant. It also provides a potential molecular basis to generate 'glandless-seed' and 'glanded-plant' cotton cultivar.

Genome-wide identification and expression analysis of GL2-interacting-repressor (GIR) genes during cotton fiber and fuzz development.

MAIN CONCLUSION: GL2-interacting-repressor (GIR) family members may contribute to fiber/fuzz formation via a newly discovered unique pathway in Gossypium arboreum. There are similarities between cotton fiber development and the formation of trichomes and root hairs. The GL2-interacting-repressors (GIRs) are crucial regulators of root hair and trichome formation. The GaFzl gene, annotated as GaGIR1, is negatively associated with trichome development and fuzz initiation. However, there is relatively little available information regarding the other GIR genes in cotton, especially regarding their effects on cotton fiber development. In this study, 21 GIR family genes were identified in the diploid cotton species Gossypium arboreum; these genes were divided into three groups. The GIR genes were characterized in terms of their phylogenetic relationships, structures, chromosomal distribution and evolutionary dynamics. These GIR genes were revealed to be unequally distributed on 12 chromosomes in the diploid cotton genome, with no GIR gene detected on Ga06. The cis-acting elements in the promoter regions were predicted to be responsive to light, phytohormones, defense activities and stress. The transcriptomic data and qRT-PCR results revealed that most GIR genes were not differentially expressed between the wild-type control and the fuzzless mutant line. Moreover, 14 of 21 family genes were expressed at high levels, indicating these genes may play important roles during fiber development and fuzz formation. Furthermore, Ga01G0231 was predominantly expressed in root samples, suggestive of a role in root hair formation rather than in fuzz initiation and development. The results of this study have enhanced our understanding of the GIR genes and their potential utility for improving cotton fiber through breeding.

Gene Expression Correlation Analysis Reveals MYC-NAC Regulatory Network in Cotton Pigment Gland Development.

Plant NAC (NAM, ATAF1/2, and CUC2) family is involved in various development processes including Programmed Cell Death (PCD) associated development. However, the relationship between NAC family and PCD-associated cotton pigment gland development is largely unknown. In this study, we identified 150, 153 and 299 NAC genes in newly updated genome sequences of G. arboreum, G. raimondii and G. hirsutum, respectively. All NAC genes were divided into 8 groups by the phylogenetic analysis and most of them were conserved during cotton evolution. Using the vital regulator of gland formation GhMYC2-like as bait, expression correlation analysis screened out 6 NAC genes which were low-expressed in glandless cotton and high-expressed in glanded cotton. These 6 NAC genes acted downstream of GhMYC2-like and were induced by MeJA. Silencing CGF1(Cotton Gland Formation1), another MYC-coding gene, caused almost glandless phenotype and down-regulated expression of GhMYC2-like and the 6 NAC genes, indicating a MYC-NAC regulatory network in gland development. In addition, predicted regulatory mechanism showed that the 6 NAC genes were possibly regulated by light, various phytohormones and transcription factors as well as miRNAs. The interaction network and DNA binding sites of the 6 NAC transcription factors were also predicted. These results laid the foundation for further study of gland-related genes and gland development regulatory network.

Genome-wide identification of the expansin gene family reveals that expansin genes are involved in fibre cell growth in cotton.

BACKGROUND: Expansins (EXPs), a group of proteins that loosen plant cell walls and cellulosic materials, are involved in regulating cell growth and diverse developmental processes in plants. However, the biological functions of this gene family in cotton are still unknown. RESULTS: In this paper, we identified a total of 93 expansin genes in Gossypium hirsutum. These genes were classified into four subfamilies, including 67 GhEXPAs, 8 GhEXPBs, 6 GhEXLAs, and 12 GhEXLBs, and divided into 15 subgroups. The 93 expansin genes are distributed over 24 chromosomes, excluding Ghir_A02 and Ghir_D06. All GhEXP genes contain multiple exons, and each GhEXP protein has multiple conserved motifs. Transcript profiling and qPCR analysis revealed that the expansin genes have distinct expression patterns among different stages of cotton fibre development. Among them, 3 genes (GhEXPA4o, GhEXPA1A, and GhEXPA8h) were highly expressed in the initiation stage, 9 genes (GhEXPA4a, GhEXPA13a, GhEXPA4f, GhEXPA4q, GhEXPA8f, GhEXPA2, GhEXPA8g, GhEXPA8a, and GhEXPA4n) had high expression during the fast elongation stage, and GhEXLA1c and GhEXLA1f were preferentially expressed in the transition stage of fibre development. CONCLUSIONS: Our results provide a solid basis for further elucidation of the biological functions of expansin genes in relation to cotton fibre development and valuable genetic resources for future crop improvement.

Nasal Pressure Injuries Due to Nasal Continuous Positive Airway Pressure Treatment in Newborns: A Prospective Observational Study.

PURPOSE: The purpose of this study was to calculate incidence, severity, and risk factors of nasal pressure injuries due to nasal continuous positive airway pressure (NCPAP) treatment in newborns. DESIGN: A prospective observational study. SUBJECTS AND SETTING: Newborns admitted between March 2017 and February 2018 to the neonatal intensive care unit of the First Affiliated Hospital of Xiamen University, Xiamen, China. METHODS: All newborns' noses were examined during NCPAP application. Every NCPAP-related nasal pressure injury including occurrence date, injury severity, outcomes, and pressure injury treatment methods was recorded. These data were collected twice a week by a research nurse. Nasal pressure injuries were classified using the National Pressure Ulcer Advisory Panel/European Pressure Ulcer Advisory Panel pressure injury classification system. RESULTS: During the study period, 429 newborns received NCPAP treatment via nasal prongs. Nasal pressure injuries were observed in 149 (34.7%); 99 (66.44%) were classified as Stage 1, 48 (32.31%) were Stage 2, and 2 (1.25%) cases were classified as deep tissue injury. The risk of nasal pressure injury was significantly higher when gestational age was less than 32 weeks (odds ratio [OR], 3.728; 95% confidence interval [CI], 1.18-11.77; P ≤ .025) and in those who received NCPAP treatment for more than 6 days (OR, 0.262; 95% CI, 0.087-0.787; P ≤ .017). The mean interval between the application of NCPAP and onset of nasal pressure injury was 4.72 days (SD, 4.78; range, 0-30 days). CONCLUSIONS: Nasal pressure injuries are a prevalent complication of NCPAP use, especially in preterm newborns. Our results identified a gestational age of less than 32 weeks and longer use of NCPAP are important factors associated with nasal pressure injuries. Methods to prevent the development of injuries such as the use of a prophylactic dressing along and replacement of binasal prongs with nasal masks are advocated.

[The effects of transgenic Cry1Ac+Cry2Ab cotton on cotton bollworm control and functional response of predators on whitefly].

In this study, we detected and clarified the roles of transgenic Cry1Ac+Cry2Ab cotton "639020" in controlling cotton bollworm (Helicoverpa armigera) during critical periods of bud stage (second generation of bollworm), flowering stage (third generation of bollworm) and bolling stage (fourth generation of bollworm) as well as the influences of 639020 cotton on functional response of the main predators (Chrysopa sinica larvae, Propylaea japonica, Orius and Erigonidium graminicola ) on whitefly using transgenic Cry1Ac cotton "CCRI41" and conventional cotton "CCRI49" as the control. Our results showed that the 639020 cotton well controlled the second and third generation of bollworm, and the level of insect resistance increased by 52.85% and 16.22% separately compared with that of CCRI41, with a significant effect on the second generation of bollworm. Moreover, the number of bollworm eggs in 639020 cotton field was lower than that in CCRI41 and CCRI49 cotton fields (except the second generation of bollworm) during the cotton bud, flowering and bolling stages. Although the number of bollworm larvae in 639020 cotton field was significantly lower than that in CCRI49 field, and both under the controlling index, it has no significant difference compared with that in CCRI41 cotton field. There were also no obvious changes in predator functions of Chrysopa sinica, Propylaea japonica, Orius and Erigonidium graminicola on bemisia tabaci between 639020, CCRI41 and CCRI49 cotton filed. This study evaluated the safety of new transgenic cotton on environment, anti-insect activity of exogenous gene and the safety of production and application prospect.

Equatorin is not essential for acrosome biogenesis but is required for the acrosome reaction.

The acrosome is a specialized organelle that covers the anterior part of the sperm nucleus and plays an essential role in mammalian fertilization. However, the regulatory mechanisms controlling acrosome biogenesis and acrosome exocytosis during fertilization are largely unknown. Equatorin (Eqtn) is a membrane protein that is specifically localized to the acrosomal membrane. In the present study, the physiological functions of Eqtn were investigated using a gene knockout mouse model. We found that Eqtn(-/-) males were subfertile. Only approximately 50% of plugged females were pregnant after mating with Eqtn(-/-) males, whereas more than 90% of plugged females were pregnant after mating with control males. Sperm and acrosomes from Eqtn(-/-) mice presented normal motility and morphology. However, the fertilization and induced acrosome exocytosis rates of Eqtn-deficient sperm were dramatically reduced. Further studies revealed that the Eqtn protein might interact with Syntaxin1a and SNAP25, but loss of Eqtn did not affect the protein levels of these genes. Therefore, our study demonstrates that Eqtn is not essential for acrosome biogenesis but is required for the acrosome reaction. Eqtn is involved in the fusion of the outer acrosomal membrane and the sperm plasma membrane during the acrosome reaction, most likely via an interaction with the SNARE complex.

Wt1 is involved in leydig cell steroid hormone biosynthesis by regulating paracrine factor expression in mice.

Wilms tumor gene, Wt1, is abundantly expressed in testis Sertoli cells. Our recent study demonstrated that Wt1 is involved in spermatogenesis by regulating Sertoli cell polarity. In the present study, we found that Wt1 is also required for steroidogenesis in Leydig cells and that deletion of the Wt1 gene resulted in defects in testosterone biosynthesis and downregulation of steroidogenic gene expression, including cytochrome P450 side-chain cleavage (P450scc), steroidogenic acute regulatory protein (StAR), 3beta-hydroxysteroid dehydrogenase I (3beta-HSD), and cytochrome P450 17A1 (Cyp17a1). The expression of LHR was significantly decreased in Wt1(-/flox); Cre-ER(TM) testes after tamoxifen induction, whereas the luteinizing hormone level in serum was unchanged. Further studies revealed that desert hedgehog (Dhh) expression was regulated by Wt1 in Sertoli cells and that its expression was significantly reduced in Wt1-deficient testes. In vitro study demonstrated that the defect in testosterone production and decreased expression of several steroidogenic genes in Wt1-deficient testis explants was partially rescued by smoothened agonist (SAG), a hedgehog pathway agonist. These results indicate that Wt1 is most likely involved in Leydig cell steroidogenesis by regulating the expression of paracrine factors in seminiferous tubules. Dhh probably had important roles in this process, but we could not exclude the possibility that other factors were also required for Leydig cell steroidogenesis. Loss of Wt1 leads to downregulation of paracrine factors, which in turn causes a decrease in steroidogenic enzyme expression and reduces testosterone production in Leydig cells. The results of this study further confirm that the cross talk between Sertoli cells and Leydig cells has important roles in Leydig cell steroidogenesis.

Cotton transposon-related variome reveals roles of transposon-related variations in modern cotton cultivation.

INTRODUCTION: Transposon plays a vital role in cotton genome evolution, contributing to the expansion and divergence of genomes within the Gossypium genus. However, knowledge of transposon activity in modern cotton cultivation is limited. OBJECTIVES: In this study, we aimed to construct transposon-related variome within Gossypium genus and reveal role of transposon-related variations during cotton cultivation. In addition, we try to identify valuable transposon-related variations for cotton breeding. METHODS: We utilized graphical genome construction to build up the graphical transposon-related variome. Based on the graphical variome, we integrated t-test, eQTL analysis and Mendelian Randomization (MR) to identify valuable transposon activities and elite genes. In addition, a convolutional neural network (CNN) model was constructed to evaluate epigenomic effects of transposon-related variations. RESULTS: We identified 35,980 transposon activities among 10 cotton genomes, and the diversity of genomic and epigenomic features was observed among 21 transposon categories. The graphical cotton transposon-related variome was constructed, and 9,614 transposon-related variations with plasticity in the modern cotton cohort were used for eQTL, phenotypic t-test and Mendelian Randomization. 128 genes were identified as gene resources improving fiber length and strength simultaneously. 4 genes were selected from 128 genes to construct the elite gene panel whose utility has been validated in a natural cotton cohort and 2 accessions with phenotypic divergence. Based on the eQTL analysis results, we identified transposon-related variations involved in cotton's environmental adaption and human domestication, providing evidence of their role in cotton's adaption-domestication cooperation. CONCLUSIONS: The cotton transposon-related variome revealed the role of transposon-related variations in modern cotton cultivation, providing genomic resources for cotton molecular breeding.

Fingerprint Finder: Identifying Genomic Fingerprint Sites in Cotton Cohorts for Genetic Analysis and Breeding Advancement.

Genomic data in Gossypium provide numerous data resources for the cotton genomics community. However, to fill the gap between genomic analysis and breeding field work, detecting the featured genomic items of a subset cohort is essential for geneticists. We developed FPFinder v1.0 software to identify a subset of the cohort's fingerprint genomic sites. The FPFinder was developed based on the term frequency-inverse document frequency algorithm. With the short-read sequencing of an elite cotton pedigree, we identified 453 pedigree fingerprint genomic sites and found that these pedigree-featured sites had a role in cotton development. In addition, we applied FPFinder to evaluate the geographical bias of fiber-length-related genomic sites from a modern cotton cohort consisting of 410 accessions. Enriching elite sites in cultivars from the Yangtze River region resulted in the longer fiber length of Yangze River-sourced accessions. Apart from characterizing functional sites, we also identified 12,536 region-specific genomic sites. Combining the transcriptome data of multiple tissues and samples under various abiotic stresses, we found that several region-specific sites contributed to environmental adaptation. In this research, FPFinder revealed the role of the cotton pedigree fingerprint and region-specific sites in cotton development and environmental adaptation, respectively. The FPFinder can be applied broadly in other crops and contribute to genetic breeding in the future.

GhHDZ76, a cotton HD-Zip transcription factor, involved in regulating the initiation and early elongation of cotton fiber development in G. hirsutum.

In this study, the whole HD-Zip family members of G. hirsutum were identified, and GhHDZ76 was classified into the HD-Zip IV subgroup. GhHDZ76 was predominantly expressed in the 0-5 DPA of fiber development stage and localized in the nucleus. Overexpression of GhHDZ76 significantly increased the length and density of trichomes in Arabidopsis thaliana. The fiber length of GhHDZ76 knockout lines by CRISPR/Cas9 was significantly shorter than WT at the early elongation and mature stage, indicating that GhHDZ76 positively regulate the fiber elongation. Scanning electron microscopy showed that the number of ovule surface protrusion of 0 DPA of GhHDZ76 knockout lines was significantly lower than WT, suggesting that GhHDZ76 can also promote the initiation of fiber development. The transcript level of GhWRKY16, GhRDL1, GhEXPA1 and GhMYB25 genes related to fiber initiation and elongation in GhHDZ76 knockout lines were significantly decreased. Yeast two-hybrid and Luciferase complementation imaging (LCI) assays showed that GhHDZ76 can interact with GhWRKY16 directly. As a transcription factor, GhHDZ76 has transcriptional activation activity, which could bind to L1-box elements of the promoters of GhRDL1 and GhEXPA1. Double luciferase reporter assay showed that the GhWRKY16 could enhance the transcriptional activity of GhHDZ76 to pGhRDL1, but it did not promote the transcriptional activity of GhHDZ76 to pGhEXPA1. GhHDZ76 protein may also promote the transcriptional activity of GhWRKY16 to the downstream target gene GhMYB25. Our results provided a new gene resource for fiber development and a theoretical basis for the genetic improvement of cotton fiber quality.

Genome-wide analysis elucidates the roles of GhHMA genes in different abiotic stresses and fiber development in upland cotton.

The heavy metal-binding domain is involved in heavy metal transporting and plays a significant role in plant detoxification. However, the functions of HMAs are less well known in cotton. In this study, a total of 143 GhHMAs (heavy metal-binding domain) were detected by genome-wide identification in G. hirsutum L. All the GhHMAs were classified into four groups via phylogenetic analysis. The exon/intron structure and protein motifs indicated that each branch of the GhHMA genes was highly conserved. 212 paralogous GhHMA gene pairs were identified, and the segmental duplications were the main role to the expansion of GhHMAs. The Ka/Ks values suggested that the GhHMA gene family has undergone purifying selection during the long-term evolutionary process. GhHMA3 and GhHMA75 were located in the plasma membrane, while GhHMA26, GhHMA117 and GhHMA121 were located in the nucleus, respectively. Transcriptomic data and qRT-PCR showed that GhHMA26 exhibited different expression patterns in each tissue and during fiber development or under different abiotic stresses. Overexpressing GhHMA26 significantly promoted the elongation of leaf trichomes and also improved the tolerance to salt stress. Therefore, GhHMA26 may positively regulate fiber elongation and abiotic stress. Yeast two-hybrid assays indicated that GhHMA26 and GhHMA75 participated in multiple biological functions. Our results suggest some genes in the GhHMAs might be associated with fiber development and the abiotic stress response, which could promote further research involving functional analysis of GhHMA genes in cotton.

Infection of SARS-CoV-2 causes severe pathological changes in mouse testis.

Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has affected more than 600 million people worldwide. Several organs including lung, intestine, and brain are infected by SARS-CoV-2. It has been reported that SARS-CoV-2 receptor angiotensin-converting enzyme-2 (ACE2) is expressed in human testis. However, whether testis is also affected by SARS-CoV-2 is still unclear. In this study, we generate a human ACE2 (hACE2) transgenic mouse model in which the expression of hACE2 gene is regulated by hACE2 promoter. Sertoli and Leydig cells from hACE2 transgenic mice can be infected by SARS-CoV-2 pseudovirus in vitro, and severe pathological changes are observed after injecting the SARS-CoV-2 pseudovirus into the seminiferous tubules. Further studies reveal that Sertoli and Leydig cells from hACE2 transgenic mice are also infected by authentic SARS-CoV-2 virus in vitro. After testis interstitium injection, authentic SARS-CoV-2 viruses are first disseminated to the interstitial cells, and then detected inside the seminiferous tubules which in turn cause germ cell loss and disruption of seminiferous tubules. Our study demonstrates that testis is most likely a target of SARS-CoV-2 virus. Attention should be paid to the reproductive function in SARS-CoV-2 patients.

Comparative Transcriptome Analysis Revealed Key Genes Regulating Gossypol Synthesis in Tetraploid Cultivated Cotton.

Tetraploid cultivated cotton (Gossypium spp.) produces cottonseeds rich in protein and oil. Gossypol and related terpenoids, stored in the pigment glands of cottonseeds, are toxic to human beings and monogastric animals. However, a comprehensive understanding of the genetic basis of gossypol and gland formation is still lacking. We performed a comprehensive transcriptome analysis of four glanded versus two glandless tetraploid cultivars distributed in Gossypium hirsutum and Gossypium barbadense. A weighted gene co-expression network analysis (WGCNA) based on 431 common differentially expressed genes (DEGs) uncovered a candidate module that was strongly associated with the reduction in or disappearance of gossypol and pigment glands. Further, the co-expression network helped us to focus on 29 hub genes, which played key roles in the regulation of related genes in the candidate module. The present study contributes to our understanding of the genetic basis of gossypol and gland formation and serves as a rich potential source for breeding cotton cultivars with gossypol-rich plants and gossypol-free cottonseed, which is beneficial for improving food safety, environmental protection, and economic gains of tetraploid cultivated cotton.

Functional characterization of the GhNRT2.1e gene reveals its significant role in improving nitrogen use efficiency in Gossypium hirsutum.

Background: Nitrate is the primary type of nitrogen available to plants, which is absorbed and transported by nitrate transporter 2 (NRT2) at low nitrate conditions. Methods: Genome-wide identification of NRT2 genes in G. hirsutum was performed. Gene expression patterns were revealed using RNA-seq and qRT-PCR. Gene functions were characterized using overexpression in A. thaliana and silencing in G. hirsutum. Protein interactions were verified by yeast two-hybrid and luciferase complementation imaging (LCI) assays. Results: We identified 14, 14, seven, and seven NRT2 proteins in G. hirsutum, G. barbadense, G. raimondii, and G. arboreum. Most NRT2 proteins were predicted in the plasma membrane. The NRT2 genes were classified into four distinct groups through evolutionary relationships, with members of the same group similar in conserved motifs and gene structure. The promoter regions of NRT2 genes included many elements related to growth regulation, phytohormones, and abiotic stresses. Tissue expression pattern results revealed that most GhNRT2 genes were specifically expressed in roots. Under low nitrate conditions, GhNRT2 genes exhibited different expression levels, with GhNRT2.1e being the most up-regulated. Arabidopsis plants overexpressing GhNRT2.1e exhibited increased biomass, nitrogen and nitrate accumulation, nitrogen uptake and utilization efficiency, nitrogen-metabolizing enzyme activity, and amino acid content under low nitrate conditions. In addition, GhNRT2.1e-silenced plants exhibited suppressed nitrate uptake and accumulation, hampered plant growth, affected nitrogen metabolism processes, and reduced tolerance to low nitrate. The results showed that GhNRT2.1e could promote nitrate uptake and transport under low nitrate conditions, thus effectively increasing nitrogen use efficiency (NUE). We found that GhNRT2.1e interacts with GhNAR2.1 by yeast two-hybrid and LCI assays. Discussion: Our research lays the foundation to increase NUE and cultivate new cotton varieties with efficient nitrogen use.

Cotton pedigree genome reveals restriction of cultivar-driven strategy in cotton breeding.

BACKGROUND: Many elite genes have been identified from the available cotton genomic data, providing various genetic resources for gene-driven breeding. However, backbone cultivar-driven breeding is the most widely applied strategy. Revealing the genetic basis of cultivar-driven strategy's restriction is crucial for transition of cotton breeding strategy. RESULT: CRI12 is a backbone cultivar in cultivar-driven breeding. Here we sequence the pedigree of CRI12 using Nanopore long-read sequencing. We construct a graphical pedigree genome using the high-quality CRI12 genome and 13,138 structural variations within 20 different pedigree members. We find that low hereditary stability of elite segments in backbone cultivars is a drawback of cultivar-driven strategy. We also identify 623 functional segments in CRI12 for multiple agronomic traits in presence and absence variation-based genome-wide association study on three cohorts. We demonstrate that 25 deleterious segments are responsible for the geographical divergence of cotton in pathogen resistance. We also characterize an elite pathogen-resistant gene (GhKHCP) utilized in modern cotton breeding. In addition, we identify 386 pedigree fingerprint segments by comparing the segments of the CRI12 pedigree with those of a large cotton population. CONCLUSION: We characterize the genetic patterns of functional segments in the pedigree of CRI12 using graphical genome method, revealing restrictions of cultivar-driven strategies in cotton breeding. These findings provide theoretical support for transitioning from cultivar-driven to gene-driven strategy in cotton breeding.