The Sister Chromatid Division of the Heteromorphic Sex Chromosomes in Silene Species and Their Transmissibility towards the Mitosis

Young sex chromosomes possess unique and ongoing dynamics that allow us to understand processes that have an impact on their evolution and divergence. The genus Silene includes species with evolutionarily young sex chromosomes, and two species of section Melandrium, namely Silene latifolia (24, XY) and Silene dioica (24, XY), are well-established models of sex chromosome evolution, Y chromosome degeneration, and sex determination. In both species, the X and Y chromosomes are strongly heteromorphic and differ in the genomic composition compared to the autosomes. It is generally accepted that for proper cell division, the longest chromosomal arm must not exceed half of the average length of the spindle axis at telophase. Yet, it is not clear what are the dynamics between males and females during mitosis and how the cell compensates for the presence of the large Y chromosome in one sex. Using hydroxyurea cell synchronization and 2D/3D microscopy, we determined the position of the sex chromosomes during the mitotic cell cycle and determined the upper limit for the expansion of sex chromosome non-recombining region. Using 3D specimen preparations, we found that the velocity of the large chromosomes is compensated by the distant positioning from the central interpolar axis, confirming previous mathematical modulations.     

Integration of Maps Enables a Cytogenomics Analysis of the Complete Karyotype in Solea senegalensis

The Pleuronectiformes order, which includes several commercially-important species, has undergone extensive chromosome evolution. One of these species is Solea senegalensis, a flatfish with 2n = 42 chromosomes. In this study, a cytogenomics approach and integration with previous maps was applied to characterize the karyotype of the species. Synteny analysis of S. senegalensis was carried out using two flatfish as a reference: Cynoglossus semilaevis and Scophthalmus maximus. Most S. senegalensis chromosomes (or chromosome arms for metacentrics and submetacentrics) showed a one-to-one macrosyntenic pattern with the other two species. In addition, we studied how repetitive sequences could have played a role in the evolution of S. senegalensis bi-armed (3, and 5–9) and acrocentric (11, 12 and 16) chromosomes, which showed the highest rearrangements compared with the reference species. A higher abundance of TEs (Transposable Elements) and other repeated elements was observed adjacent to telomeric regions on chromosomes 3, 7, 9 and 16. However, on chromosome 11, a greater abundance of DNA transposons was detected in interstitial BACs. This chromosome is syntenic with several chromosomes of the other two flatfish species, suggesting rearrangements during its evolution. A similar situation was also found on chromosome 16 (for microsatellites and low complexity sequences), but not for TEs (retroelements and DNA transposons). These differences in the distribution and abundance of repetitive elements in chromosomes that have undergone remodeling processes during the course of evolution also suggest a possible role for simple repeat sequences in rearranged regions.     

PGRP-LB: An Inside View into the Mechanism of the Amidase Reaction

Peptidoglycan recognition proteins (PGRPs) are ubiquitous among animals and play pivotal functions in insect immunity. Non-catalytic PGRPs are involved in the activation of immune pathways by binding to the peptidoglycan (PGN), whereas amidase PGRPs are capable of cleaving the PGN into non-immunogenic compounds. Drosophila PGRP-LB belongs to the amidase PGRPs and downregulates the immune deficiency (IMD) pathway by cleaving meso-2,6-diaminopimelic (meso-DAP or DAP)-type PGN. While the recognition process is well analyzed for the non-catalytic PGRPs, little is known about the enzymatic mechanism for the amidase PGRPs, despite their essential function in immune homeostasis. Here, we analyzed the specific activity of different isoforms of Drosophila PGRP-LB towards various PGN substrates to understand their specificity and role in Drosophila immunity. We show that these isoforms have similar activity towards the different compounds. To analyze the mechanism of the amidase activity, we performed site directed mutagenesis and solved the X-ray structures of wild-type Drosophila PGRP-LB and its mutants, with one of these structures presenting a protein complexed with the tracheal cytotoxin (TCT), a muropeptide derived from the PGN. Only the Y78F mutation abolished the PGN cleavage while other mutations reduced the activity solely. Together, our findings suggest the dynamic role of the residue Y78 in the amidase mechanism by nucleophilic attack through a water molecule to the carbonyl group of the amide function destabilized by Zn²⁺.     

Cytogenetically Elusive Sex Chromosomes in Scincoidean Lizards

The lizards of the species-rich clade Scincoidea including cordylids, gerrhosaurids, skinks, and xantusiids, show an almost cosmopolitan geographical distribution and a remarkable ecological and morphological divergence. However, previous studies revealed limited variability in cytogenetic traits. The sex determination mode was revealed only in a handful of gerrhosaurid, skink, and xantusiid species, which demonstrated either ZZ/ZW or XX/XY sex chromosomes. In this study, we explored the karyotypes of six species of skinks, two species of cordylids, and one gerrhosaurid. We applied conventional and molecular cytogenetic methods, including C-banding, fluorescence in situ hybridization with probes specific for telomeric motifs and rDNA loci, and comparative genomic hybridization. The diploid chromosome numbers are rather conserved among these species, but the chromosome morphology, the presence of interstitial telomeric sequences, and the topology of rDNA loci vary significantly. Notably, XX/XY sex chromosomes were identified only in Tiliqua scincoides, where, in contrast to the X chromosome, the Y chromosome lacks accumulations of rDNA loci. We confirm that within the lizards of the scincoidean clade, sex chromosomes remained in a generally poor stage of differentiation.     

Genome-Wide Analysis of the Peptidase M24 Superfamily in Triticum aestivum Demonstrates That TaM24-9 Is Involved in Abiotic Stress Response

The peptidase M24 (Metallopeptidase 24, M24) superfamily is essential for plant growth, stress response, and pathogen defense. At present, there are few systematic reports on the identification and classification of members of the peptidase M24 proteins superfamily in wheat. In this work, we identified 53 putative candidate TaM24 genes. According to the protein sequences characteristics, these members can be roughly divided into three subfamilies: I, II, III. Most TaM24 genes are complex with multiple exons, and the motifs are relatively conserved in each sub-group. Through chromosome mapping analysis, we found that the 53 genes were unevenly distributed on 19 wheat chromosomes (except 3A and 3D), of which 68% were in triads. Analysis of gene duplication events showed that 62% of TaM24 genes in wheat came from fragment duplication events, and there were no tandem duplication events to amplify genes. Analysis of the promoter sequences of TaM24 genes revealed that cis-acting elements were rich in response elements to drought, osmotic stress, ABA, and MeJA. We also studied the expression of TaM24 in wheat tissues at developmental stages and abiotic stress. Then we selected TaM24-9 as the target for further analysis. The results showed that TaM24-9 genes strengthened the drought and salt tolerance of plants. Overall, our analysis showed that members of the peptidase M24 genes may participate in the abiotic stress response and provided potential gene resources for improving wheat resistance.     

Genome-Wide Identification and Characterization of the CC-NBS-LRR Gene Family in Cucumber (Cucumis sativus L.)

The NBS-LRR (NLR) gene family plays a pivotal role in regulating disease defense response in plants. Cucumber is one of the most important vegetable crops in the world, and various plant diseases, including powdery mildew (PM), cause severe losses in both cucumber productivity and quality annually. To characterize and understand the role of the CC-NBS-LRR(CNL) family of genes in disease defense response in cucumber plants, we performed bioinformatical analysis to characterize these genes systematically. We identified 33 members of the CNL gene family in cucumber plants, and they are distributed on each chromosome with chromosome 4 harboring the largest cluster of five different genes. The corresponding CNL family member varies in the number of amino acids and exons, molecular weight, theoretical isoelectric point (pI) and subcellular localization. Cis-acting element analysis of the CNL genes reveals the presence of multiple phytohormone, abiotic and biotic responsive elements in their promoters, suggesting that these genes might be responsive to plant hormones and stress. Phylogenetic and synteny analysis indicated that the CNL proteins are conserved evolutionarily in different plant species, and they can be divided into four subfamilies based on their conserved domains. MEME analysis and multiple sequence alignment showed that conserved motifs exist in the sequence of CNLs. Further DNA sequence analysis suggests that CsCNL genes might be subject to the regulation of different miRNAs upon PM infection. By mining available RNA-seq data followed by real-time quantitative PCR (qRT-PCR) analysis, we characterized expression patterns of the CNL genes, and found that those genes exhibit a temporospatial expression pattern, and their expression is also responsive to PM infection, ethylene, salicylic acid, and methyl jasmonate treatment in cucumber plants. Finally, the CNL genes targeted by miRNAs were predicted in cucumber plants. Our results in this study provided some basic information for further study of the functions of the CNL gene family in cucumber plants.     

Genome-Wide Identification and Analysis of the Metallothionein Genes in Oryza Genus

Metallothionein (MT) proteins are low molecular mass, cysteine-rich, and metal-binding proteins that play an important role in maintaining metal homeostasis and stress response. However, the evolutionary relationships and functional differentiation of MT in the Oryza genus remain unclear. Here we identified 53 MT genes from six Oryza genera, including O. sativa ssp. japonica, O. rufipogon, O. sativa ssp. indica, O. nivara, O. glumaepatula, and O. barthii. The MT genes were clustered into four groups based on phylogenetic analysis. MT genes are unevenly distributed on chromosomes; almost half of the MT genes were clustered on chromosome 12, which may result from a fragment duplication containing the MT genes on chromosome 12. Five pairs of segmental duplication events and ten pairs of tandem duplication events were found in the rice MT family. The Ka/Ks values of the fifteen duplicated MT genes indicated that the duplicated MT genes were under a strong negative selection during evolution. Next, combining the promoter activity assay with gene expression analysis revealed different expression patterns of MT genes. In addition, the expression of OsMT genes was induced under different stresses, including NaCl, CdCl₂, ABA, and MeJ treatments. Additionally, we found that OsMT genes were mainly located in chloroplasts. These results imply that OsMT genes play different roles in response to these stresses. All results provide important insights into the evolution of the MT gene family in the Oryza genus, and will be helpful to further study the function of MT genes.     

Genes Involved in the Endoplasmic Reticulum N-Glycosylation Pathway of the Red Microalga Porphyridium sp.: A Bioinformatic Study

N-glycosylation is one of the most important post-translational modifications that influence protein polymorphism, including protein structures and their functions. Although this important biological process has been extensively studied in mammals, only limited knowledge exists regarding glycosylation in algae. The current research is focused on the red microalga Porphyridium sp., which is a potentially valuable source for various applications, such as skin therapy, food, and pharmaceuticals. The enzymes involved in the biosynthesis and processing of N-glycans remain undefined in this species, and the mechanism(s) of their genetic regulation is completely unknown. In this study, we describe our pioneering attempt to understand the endoplasmic reticulum N-Glycosylation pathway in Porphyridium sp., using a bioinformatic approach. Homology searches, based on sequence similarities with genes encoding proteins involved in the ER N-glycosylation pathway (including their conserved parts) were conducted using the TBLASTN function on the algae DNA scaffold contigs database. This approach led to the identification of 24 encoded-genes implicated with the ER N-glycosylation pathway in Porphyridium sp. Homologs were found for almost all known N-glycosylation protein sequences in the ER pathway of Porphyridium sp.; thus, suggesting that the ER-pathway is conserved; as it is in other organisms (animals, plants, yeasts, etc.).     

Genome-Wide Analysis of the Late Embryogenesis Abundant (LEA) and Abscisic Acid-, Stress-, and Ripening-Induced (ASR) Gene Superfamily from Canavalia rosea and Their Roles in Salinity/Alkaline and Drought Tolerance

Canavalia rosea (bay bean), distributing in coastal areas or islands in tropical and subtropical regions, is an extremophile halophyte with good adaptability to seawater and drought. Late embryogenesis abundant (LEA) proteins typically accumulate in response to various abiotic stresses, including dehydration, salinity, high temperature, and cold, or during the late stage of seed development. Abscisic acid-, stress-, and ripening-induced (ASR) genes are stress and developmentally regulated plant-specific genes. In this study, we reported the first comprehensive survey of the LEA and ASR gene superfamily in C. rosea. A total of 84 CrLEAs and three CrASRs were identified in C. rosea and classified into nine groups. All CrLEAs and CrASRs harbored the conserved motif for their family proteins. Our results revealed that the CrLEA genes were widely distributed in different chromosomes, and all of the CrLEA/CrASR genes showed wide expression features in different tissues in C. rosea plants. Additionally, we introduced 10 genes from different groups into yeast to assess the functions of the CrLEAs/CrASRs. These results contribute to our understanding of LEA/ASR genes from halophytes and provide robust candidate genes for functional investigations in plant species adapted to extreme environments.     

Molecular Karyotyping on Populus simonii × P. nigra and the Derived Doubled Haploid

The molecular karyotype could represent the basic genetic make-up in a cell nucleus of an organism or species. A doubled haploid (DH) is a genotype formed from the chromosome doubling of haploid cells. In the present study, molecular karyotype analysis of the poplar hybrid Populus simonii × P. nigra (P. xiaohei) and the derived doubled haploids was carried out with labeled telomeres, rDNA, and two newly repetitive sequences as probes by fluorescence in situ hybridization (FISH). The tandem repeats, pPC349_XHY and pPD284_XHY, with high-sequence homology were used, and the results showed that they presented the colocalized distribution signal in chromosomes. For P. xiaohei, pPD284_XHY produced hybridizations in chromosomes 1, 5, 8, and 9 in the hybrid. The combination of pPD284_XHY, 45S rDNA, and 5S rDNA distinctly distinguished six pairs of chromosomes, and the three pairs of chromosomes showed a significant difference in the hybridization between homologous chromosomes. The repeat probes used produced similar FISH hybridizations in the DH; nevertheless, pPD284_XHY generated an additional hybridization site in the telomere region of chromosome 14. Moreover, two pairs of chromosomes showed differential hybridization distributions between homologous chromosomes. Comparisons of the distinguished chromosomes between hybrid and DH poplar showed that three pairs of chromosomes in the DH presented hybridization patterns that varied from those of the hybrid. The No. 8 chromosome in DH and one of the homologous chromosomes in P. xiaohei shared highly similar FISH patterns, which suggested the possibility of intact or mostly partial transfer of the chromosome between the hybrid and DH. Our study will contribute to understanding the genetic mechanism of chromosomal variation in P. xiaohei and derived DH plants.     

Genome-Wide Comparative Analysis of the R2R3-MYB Gene Family in Five Solanaceae Species and Identification of Members Regulating Carotenoid Biosynthesis in Wolfberry

The R2R3-MYB is a large gene family involved in various plant functions, including carotenoid biosynthesis. However, this gene family lacks a comprehensive analysis in wolfberry (Lycium barbarum L.) and other Solanaceae species. The recent sequencing of the wolfberry genome provides an opportunity for investigating the organization and evolutionary characteristics of R2R3-MYB genes in wolfberry and other Solanaceae species. A total of 610 R2R3-MYB genes were identified in five Solanaceae species, including 137 in wolfberry. The LbaR2R3-MYB genes were grouped into 31 subgroups based on phylogenetic analysis, conserved gene structures, and motif composition. Five groups only of Solanaceae R2R3-MYB genes were functionally divergent during evolution. Dispersed and whole duplication events are critical for expanding the R2R3-MYB gene family. There were 287 orthologous gene pairs between wolfberry and the other four selected Solanaceae species. RNA-seq analysis identified the expression level of LbaR2R3-MYB differential gene expression (DEGs) and carotenoid biosynthesis genes (CBGs) in fruit development stages. The highly expressed LbaR2R3-MYB genes are co-expressed with CBGs during fruit development. A quantitative Real-Time (qRT)-PCR verified seven selected candidate genes. Thus, Lba11g0183 and Lba02g01219 are candidate genes regulating carotenoid biosynthesis in wolfberry. This study elucidates the evolution and function of R2R3-MYB genes in wolfberry and the four Solanaceae species.     

The Landscape of the Genomic Distribution and the Expression of the F-Box Genes Unveil Genome Plasticity in Hexaploid Wheat during Grain Development and in Response to Heat and Drought Stress

FBX proteins are subunits of the SCF complex (Skp1–cullin–FBX) belonging to the E3 ligase family, which is involved in the ubiquitin–proteasome 26S (UPS) pathway responsible for the post-translational protein turnover. By targeting, in a selective manner, key regulatory proteins for ubiquitination and 26S proteasome degradation, FBX proteins play a major role in plant responses to diverse developmental and stress conditions. Although studies on the genomic organization of the FBX gene family in various species have been reported, knowledge related to bread wheat (Triticum aestivum) is scarce and needs to be broadened. Using the latest assembly of the wheat genome, we identified 3670 TaFBX genes distributed non-homogeneously within the three subgenomes (A, B and D) and between the 21 chromosomes, establishing it as one of the richest gene families among plant species. Based on the presence of the five different chromosomal regions previously identified, the present study focused on the genomic distribution of the TaFBX family and the identification of differentially expressed genes during the embryogenesis stages and in response to heat and drought stress. Most of the time, when comparing the expected number of genes (taking into account the formal gene distribution on the entire wheat genome), the TaFBX family harbors a different pattern at the various stratum of observation (subgenome, chromosome, chromosomal regions). We report here that the local gene expansion of the TaFBX family must be the consequence of multiple and complex events, including tandem and small-scale duplications. Regarding the differentially expressed TaFBX genes, while the majority of the genes are localized in the distal chromosomal regions (R1 and R3), differentially expressed genes are more present in the interstitial regions (R2a and R2b) than expected, which could be an indication of the preservation of major genes in those specific chromosomal regions.     

A Comprehensive Identification and Function Analysis of Serine/Arginine-Rich (SR) Proteins in Cotton (Gossypium spp.)

As one of the most important factors in alternative splicing (AS) events, serine/arginine-rich (SR) proteins not only participate in the growth and development of plants but also play pivotal roles in abiotic stresses. However, the research about SR proteins in cotton is still lacking. In this study, we performed an extensive comparative analysis of SR proteins and determined their phylogeny in the plant lineage. A total of 169 SR family members were identified from four Gossypium species, and these genes could be divided into eight distinct subfamilies. The domain, motif distribution and gene structure of cotton SR proteins are conserved within each subfamily. The expansion of SR genes is mainly contributed by WGD and allopolyploidization events in cotton. The selection pressure analysis showed that all the paralogous gene pairs were under purifying selection pressure. Many cis-elements responding to abiotic stress and phytohormones were identified in the upstream sequences of the GhSR genes. Expression profiling suggested that some GhSR genes may involve in the pathways of plant resistance to abiotic stresses. The WGCNA analysis showed that GhSCL-8 co-expressed with many abiotic responding related genes in a salt-responding network. The Y2H assays showed that GhSCL-8 could interact with GhSRs in other subfamilies. The subcellular location analysis showed that GhSCL-8 is expressed in the nucleus. The further VIGS assays showed that the silencing of GhSCL-8 could decrease salt tolerance in cotton. These results expand our knowledge of the evolution of the SR gene family in plants, and they will also contribute to the elucidation of the biological functions of SR genes in the future.