Segmental Duplications Research Articles

Genome-Wide Identification and Expression Analysis of Sucrose Transporter Gene Family in Wheat Lines under Heat Stress

Sucrose transporters (SUTs) play vital roles in phloem sucrose unloading and transportation in wheat grains. However, the genomic information regarding the SUT gene family and their expression patterns in response to heat stress in grains of male-sterile wheat (Triticum aestivum L.) lines has not been systematically studied. In this study, a thorough examination of the wheat SUT gene family was conducted, focusing on their expression patterns in male-sterile lines under heat stress conditions in grain tissues. A total of 19 SUT genes were identified, with phylogenetic analysis indicating their classification into five distinct groups. Polyploidization was identified as a substantial factor in the expansion of SUT genes, with segmental duplication being the predominant mechanism driving the evolutionary expansion of the SUT gene family in wheat. Transcriptome data indicate that the expression levels of TaSUT1 and TaSUT2 were higher than other SUT genes in grains of male-sterile lines. The TaSUT1 expression showed a gradual decreasing trend, while TaSUT2 showed a reverse trend with the process of grain filling. After heat stress, the TaSUT1 expression in grains of male-sterile lines was first significantly increased and then significantly decreased with the filling stage extension, aligning with the observed trend of sucrose levels, indicating that heat stress may decrease the grain weight by reducing sucrose unloading and transportation process in grains. These results provide a systematic analysis of the SUT gene family and lay a theoretical foundation for us to study the grain filling of male-sterile lines in response to abiotic stress.

Genome-Wide Identification and Expression Analysis of ADK Gene Family Members in Cotton under Abiotic Stress.

Adenosine kinase (ADK) is a key enzyme widely distributed in plants, playing an important role in maintaining cellular energy homeostasis and regulating plant growth, development, and responses to environmental stresses. However, research on ADK genes in cotton (Gossypium hirsutum), an economically significant crop, has been limited. This study identified 92 ADK genes from four cotton species (G. arboreum, G. raimondii, G. hirsutum, and G. barbadense) using HMMER and Local BLASTP methods and classified them into six groups. Chromosomal localization revealed a random distribution of ADK genes in G. hirsutum, with 13 genes located on the At subgenome and 14 genes on the Dt subgenome. Gene structure analysis showed consistency in exon-intron organization within subgroups, while conserved motif analysis identified subgroup-specific motifs, indicating functional diversity. Synteny and collinearity mapping analysis revealed that the primary expansion mechanisms of the ADK gene family in cotton are polyploidy and segmental duplication. Cis-regulatory elements in GhADK promoters were classified into light response, hormone response, developmental regulation, and stress response. We also analyzed the expression patterns of GhADK genes under a low temperature (4 °C) and drought conditions. Most GhADK genes responded to cold stress with different expression patterns, indicating their roles in rapid response and long-term cold adaptation. Under drought stress, expression patterns varied, with some genes showing sustained high expression levels. The qRT-PCR validation of transcriptomic data confirmed the stress-induced expression patterns of selected GhADK genes. Functional analysis through the VIGS silencing of GhADK25 demonstrated its importance in cold and drought stress responses, with silencing resulting in poor growth under stress, highlighting its significance in stress tolerance. This study provides a basis for further understanding the evolutionary relationships and functions of the cotton ADK gene family.

IS1-mediated chromosomal amplification of the arn operon leads to polymyxin B resistance in Escherichia coli B strains.

Polymyxins [colistin and polymyxin B (PMB)] comprise an important class of natural product lipopeptide antibiotics used to treat multidrug-resistant Gram-negative bacterial infections. These positively charged lipopeptides interact with lipopolysaccharide (LPS) located in the outer membrane and disrupt the permeability barrier, leading to increased uptake and bacterial cell death. Many bacteria counter polymyxins by upregulating genes involved in the biosynthesis and transfer of amine-containing moieties to increase positively charged residues on LPS. Although 4-deoxy-l-aminoarabinose (Ara4N) and phosphoethanolamine (PEtN) are highly conserved LPS modifications in Escherichia coli, different lineages exhibit variable PMB susceptibilities and frequencies of resistance for reasons that are poorly understood. Herein, we describe a mechanism prevalent in E. coli B strains that depends on specific insertion sequence 1 (IS1) elements that flank genes involved in the biosynthesis and transfer of Ara4N to LPS. Spontaneous and transient chromosomal amplifications mediated by IS1 raise the frequency of PMB resistance by 10- to 100-fold in comparison to strains where a single IS1 element located 90 kb away from the end of the arn operon has been deleted. Amplification involving IS1 becomes the dominant resistance mechanism in the absence of PEtN modification. Isolates with amplified arn operons gradually lose their PMB-resistant phenotype with passaging, consistent with classical PMB heteroresistance behavior. Analysis of the whole genome transcriptome profile showed altered expression of genes residing both within and outside of the duplicated chromosomal segment, suggesting complex phenotypes including PMB resistance can result from tandem amplification events.IMPORTANCEPhenotypic variation in susceptibility and the emergence of resistant subpopulations are major challenges to the clinical use of polymyxins. While a large database of genes and alleles that can confer polymyxin resistance has been compiled, this report demonstrates that the chromosomal insertion sequence (IS) content and distribution warrant consideration as well. Amplification of large chromosomal segments containing the arn operon by IS1 increases the Ara4N content of the lipopolysaccharide layer in Escherichia coli B lineages using a mechanism that is orthogonal to transcriptional upregulation through two-component regulatory systems. Altogether, our work highlights the importance of IS elements in modulating gene expression and generating diverse subpopulations that can contribute to phenotypic polymyxin B heteroresistance.

Genome-Wide Characterization of Fructose 1,6-Bisphosphate Aldolase Genes and Expression Profile Reveals Their Regulatory Role in Abiotic Stress in Cucumber.

The fructose-1,6-bisphosphate aldolase (FBA) gene family exists in higher plants, with the genes of this family playing significant roles in plant growth and development, as well as response to abiotic stresses. However, systematic reports on the FBA gene family and its functions in cucumber are lacking. In this study, we identified five cucumber FBA genes, named CsFBA1-5, that are distributed randomly across chromosomes. Phylogenetic analyses involving these cucumber FBAs, alongside eight Arabidopsis FBA proteins and eight tomato FBA proteins, were conducted to assess their homology. The CsFBAs were grouped into two clades. We also analyzed the physicochemical properties, motif composition, and gene structure of the cucumber FBAs. This analysis highlighted differences in the physicochemical properties and revealed highly conserved domains within the CsFBA family. Additionally, to explore the evolutionary relationships of the CsFBA family further, we constructed comparative syntenic maps with Arabidopsis and tomato, which showed high homology but only one segmental duplication event within the cucumber genome. Expression profiles indicated that the CsFBA gene family is responsive to various abiotic stresses, including low temperature, heat, and salt. Taken together, the results of this study provide a theoretical foundation for understanding the evolution of and future research into the functional characterization of cucumber FBA genes during plant growth and development.

Genome-Wide Characterization and Expression Profiling of the AP2/ERF Gene Family in Fragaria vesca L.

The wild strawberry (Fragaria vesca L.; F. vesca) represents a resilient and extensively studied model organism. While the AP2/ERF gene family plays a pivotal role in plant development, its exploration within F. vesca remains limited. In this study, we characterized the AP2/ERF gene family in wild strawberries using the recently released genomic data (F. vesca V6.0). We conducted an analysis of the gene family expansion pattern, we examined gene expression in stem segments and leaves under cold conditions, and we explored its functional attributes. Our investigation revealed that the FvAP2/ERF family comprises 86 genes distributed among four subfamilies: AP2 (17), RAV (6), ERF (62), and Soloist (1). Tandem and segmental duplications significantly contributed to the growth of this gene family. Furthermore, predictive analysis identified several cis-acting elements in the promoter region associated with meristematic tissue expression, hormone regulation, and resistance modulation. Transcriptomic analysis under cold stress unveiled diverse responses among multiple FvAP2/ERFs in stem segments and leaves. Real-time fluorescence quantitative reverse transcription PCR (RT-qPCR) results confirmed elevated expression levels of select genes following the cold treatment. Additionally, overexpression of FvERF23 in Arabidopsis enhanced cold tolerance, resulting in significantly increased fresh weight and root length compared to the wild-type control. These findings lay the foundation for further exploration into the functional roles of FvAP2/ERF genes.

Genome-Wide Analysis of the Auxin/Indoleacetic Acid (Aux/IAA) Gene Family in Autopolyploid Sugarcane (Saccharum spontaneum).

The auxin/indoleacetic acid (Aux/IAA) family plays a central role in regulating gene expression during auxin signal transduction. Nonetheless, there is limited knowledge regarding this gene family in sugarcane. In this study, 92 members of the IAA family were identified in Saccharum spontaneum, distributed on 32 chromosomes, and classified into three clusters based on phylogeny and motif compositions. Segmental duplication and recombination events contributed largely to the expansion of this superfamily. Additionally, cis-acting elements in the promoters of SsIAAs involved in plant hormone regulation and stress responsiveness were predicted. Transcriptomics data revealed that most SsIAA expressions were significantly higher in stems and basal parts of leaves, and at nighttime, suggesting that these genes might be involved in sugar transport. QRT-PCR assays confirmed that cold and salt stress significantly induced four and five SsIAAs, respectively. GFP-subcellular localization showed that SsIAA23 and SsIAA12a were localized in the nucleus, consistent with the results of bioinformatics analysis. In conclusion, to a certain extent, the functional redundancy of family members caused by the expansion of the sugarcane IAA gene family is related to stress resistance and regeneration of sugarcane as a perennial crop. This study reveals the gene evolution and function of the SsIAA gene family in sugarcane, laying the foundation for further research on its mode of action.

Genome-Wide Identification, Evolution, and Expression Analysis of the DIR Gene Family in Schima superba.

Schima superba, commonly known as the Chinese guger tree, is highly adaptable and tolerant of poor soil conditions. It is one of the primary species forming the evergreen broad-leaved forests in southern China. Dirigent proteins (DIRs) play crucial roles in the synthesis of plant lignin and lignans, secondary metabolism, and response to adversity stress. However, research on the DIR gene family in S. superba is currently limited. This study identified 24 SsDIR genes, categorizing them into three subfamilies. These genes are unevenly distributed across 13 chromosomes, with 83% being intronless. Collinearity analysis indicated that tandem duplication played a more significant role in the expansion of the gene family compared to segmental duplication. Additionally, we analyzed the expression patterns of SsDIRs in different tissues of S. superba. The SsDIR genes exhibited distinct expression patterns across various tissues, with most being specifically expressed in the roots. Further screening identified SsDIR genes that may regulate drought stress, with many showing differential expression under drought stress conditions. In the promoter regions of SsDIRs, various cis-regulatory elements involved in developmental regulation, hormone response, and stress response were identified, which may be closely related to their diverse regulatory functions. This study will contribute to the further functional identification of SsDIR genes, providing insights into the biosynthetic pathways of lignin and lignans and the mechanisms of plant stress resistance.

Genome-wide and functional analysis of late embryogenesis abundant (LEA) genes during dormancy and sprouting periods of kernel consumption apricots (P. armeniaca L. × P. sibirica L.)

Late embryogenesis abundant (LEA) proteins play a crucial role in protecting cells from stress, making them potential contributors to abiotic stress tolerance. This study focuses on apricot (P. armeniaca L. × P. sibirica L.), where a comprehensive genome-wide analysis identified 54 LEA genes, categorized into eight subgroups based on phylogenetic relationships. Synteny analysis revealed 14 collinear blocks containing LEA genes between P. armeniaca × P. sibirica and Arabidopsis thaliana, with an additional 9 collinear blocks identified between P. armeniaca × P. sibirica and poplar. Examination of gene structure and conserved motifs indicated that these subgroups exhibit consistent exon-intron patterns and shared motifs. The expansion and duplication of LEA genes in P. armeniaca × P. sibirica were driven by whole-genome duplication (WGD), segmental duplication, and tandem duplication events. Expression analysis, utilizing RNA-seq data and quantitative real-time RT-PCR (qRT-PCR), indicated induction of PasLEA2-20, PasLEA3-2, PasLEA6-1, Pasdehydrin-3, and Pasdehydrin-5 in flower buds during dormancy and sprouting phases. Coexpression network analysis linked LEA genes with 15 cold-resistance genes. Remarkably, during the four developmental stages of flower buds in P. armeniaca × P. sibirica - physiological dormancy, ecological dormancy, sprouting period, and germination stage - the expression patterns of all PasLEAs coexpressed with cold stress-related genes remained consistent. Protein-protein interaction networks, established using Arabidopsis orthologs, emphasized connections between PasLEA proteins and cold resistance pathways. Overexpression of certain LEA genes in yeast and Arabidopsis conferred advantages under cold stress, including increased pod length, reduced bolting time and flowering time, improved survival and seed setting rates, elevated proline accumulation, and enhanced antioxidative enzymatic activities. Furthermore, these overexpressed plants exhibited upregulation of genes related to flower development and cold resistance. The Y1H assay confirmed that PasGBF4 and PasDOF3.5 act as upstream regulatory factors by binding to the promoter region of PasLEA3-2. PasDOF2.4, PasDnaJ2, and PasAP2 were also found to bind to the promoter of Pasdehydrin-3, regulating the expression levels of downstream genes. This comprehensive study explores the evolutionary relationships among PasLEA genes, protein interactions, and functional analyses during various stages of dormancy and sprouting in P. armeniaca × P. sibirica. It offers potential targets for enhancing cold resistance and manipulating flower bud dormancy in this apricot hybrid.

Genome-wide analysis of wheat xyloglucan endotransglucosylase/hydrolase (XTH) gene family revealed TaXTH17 involved in abiotic stress responses

BackgroundEnvironmental stresses, including high salinity and drought, severely diminish wheat yield and quality globally. The xyloglucan endotransglucosylase/hydrolase (XTH) family represents a class of cell wall-modifying enzymes and plays important roles in plants growth, development and stress adaptation. However, systematic analyses of XTH family genes and their functions under salt and drought stresses have not been undertaken in wheat.ResultsIn this study, we identified a total of 135 XTH genes in wheat, which were clustered into three evolutionary groups. These TaXTHs were unevenly distributed on 21 chromosomes of wheat with a majority of TaXTHs located on homelogous groups 2, 3 and 7. Gene duplication analysis revealed that segmental and tandem duplication were the main reasons for the expansion of XTH family in wheat. Interaction network predictions indicated that TaXTHs could interact with multiple proteins, including three kinases, one methyltransferase and one gibberellin-regulated protein. The promoters of the TaXTH genes harbored various cis-acting elements related to stress and hormone responses. RNA-seq data analyses showed that some TaXTH genes were induced by salt and drought stresses. Furthermore, we verified that TaXTH17 was induced by abiotic stresses and phytohormone treatments, and demonstrated that TaXTH17 was localized in the secretory pathway and cell wall. Functional analyses conducted in heterologous expression systems and in wheat established that TaXTH17 plays a negative role in plant resistance to salt and drought.ConclusionsWe identified 135 XTH genes in wheat and conducted comprehensive analyses of their phylogenetic relationships, gene structures, conserved motifs, gene duplication events, chromosome locations, interaction networks, cis-acting elements and gene expression patterns. Furthermore, we provided solid evidence supporting the notion that TaXTH17 plays a negative role in plant resistance to salt and drought stresses. Collectively, our results provide valuable insights into understanding wheat XTHs, particularly their involvement in plant stress responses, and establish a foundation for further functional and mechanistic studies of TaXTHs.

Genome-wide analysis of R2R3-MYB transcription factors in poplar and functional validation of PagMYB147 in defense against Melampsora magnusiana

Main conclusionTranscription of PagMYB147 was induced in poplar infected by Melampsora magnusiana, and a decline in its expression levels increases the host's susceptibility, whereas its overexpression promotes resistance to rust disease.Poplars are valuable tree species with diverse industrial and silvicultural applications. The R2R3-MYB subfamily of transcription factors plays a crucial role in response to biotic stresses. However, the functional studies on poplar R2R3-MYB genes in resistance to leaf rust disease are still insufficient. We identified 191 putative R2R3-MYB genes in the Populus trichocarpa genome. A phylogenetic analysis grouped poplar R2R3-MYBs and Arabidopsis R2R3-MYBs into 33 subgroups. We detected 12 tandem duplication events and 148 segmental duplication events, with the latter likely being the main contributor to the expansion of poplar R2R3-MYB genes. The promoter regions of these genes contained numerous cis-acting regulatory elements associated with response to stress and phytohormones. Analyses of RNA-Seq data identified a multiple R2R3-MYB genes response to Melampsora magnusiana (Mmag). Among them, PagMYB147 was significantly up-regulated under Mmag inoculation, salicylic acid (SA) and methyl jasmonate (MeJA) treatment, and its encoded product was primarily localized to the cell nucleus. Silencing of PagMYB147 exacerbated the severity of Mmag infection, likely because of decreased reactive oxygen species (ROS) production and phenylalanine ammonia-lyase (PAL) enzyme activity, and up-regulation of genes related to ROS scavenging and down-regulation of genes related to PAL, SA and JA signaling pathway. In contrast, plants overexpressing PagMYB147 showed the opposite ROS accumulation, PAL enzyme activity, SA and JA-related gene expressions, and improved Mmag resistance. Our findings suggest that PagMYB147 acts as a positive regulatory factor, affecting resistance in poplar to Mmag by its involvement in the regulation of ROS homeostasis, SA and JA signaling pathway.

Functional and evolutionary comparative analysis of the DIR gene family in Nicotiana tabacum L. and Solanum tuberosum L.

BackgroundThe dirigent (DIR) genes encode proteins that act as crucial regulators of plant lignin biosynthesis. In Solanaceae species, members of the DIR gene family are intricately related to plant growth and development, playing a key role in responding to various biotic and abiotic stresses. It will be of great application significance to analyze the DIR gene family and expression profile under various pathogen stresses in Solanaceae species.ResultsA total of 57 tobacco NtDIRs and 33 potato StDIRs were identified based on their respective genome sequences. Phylogenetic analysis of DIR genes in tobacco, potato, eggplant and Arabidopsis thaliana revealed three distinct subgroups (DIR-a, DIR-b/d and DIR-e). Gene structure and conserved motif analysis showed that a high degree of conservation in both exon/intron organization and protein motifs among tobacco and potato DIR genes, especially within members of the same subfamily. Total 8 pairs of tandem duplication genes (3 pairs in tobacco, 5 pairs in potato) and 13 pairs of segmental duplication genes (6 pairs in tobacco, 7 pairs in potato) were identified based on the analysis of gene duplication events. Cis-regulatory elements of the DIR promoters participated in hormone response, stress responses, circadian control, endosperm expression, and meristem expression. Transcriptomic data analysis under biotic stress revealed diverse response patterns among DIR gene family members to pathogens, indicating their functional divergence. After 96 h post-inoculation with Ralstonia solanacearum L. (Ras), tobacco seedlings exhibited typical symptoms of tobacco bacterial wilt. The qRT-PCR analysis of 11 selected NtDIR genes displayed differential expression pattern in response to the bacterial pathogen Ras infection. Using line 392278 of potato as material, typical symptoms of potato late blight manifested on the seedling leaves under Phytophthora infestans infection. The qRT-PCR analysis of 5 selected StDIR genes showed up-regulation in response to pathogen infection. Notably, three clustered genes (NtDIR2, NtDIR4, StDIR3) exhibited a robust response to pathogen infection, highlighting their essential roles in disease resistance.ConclusionThe genome-wide identification, evolutionary analysis, and expression profiling of DIR genes in response to various pathogen infection in tobacco and potato have provided valuable insights into the roles of these genes under various stress conditions. Our results could provide a basis for further functional analysis of the DIR gene family under pathogen infection conditions.

Genome-wide identification and expression analysis of calmodulin and calmodulin-like genes in passion fruit (Passiflora edulis) and their involvement in flower and fruit development

BackgroundThe calmodulin (CaM) and calmodulin-like (CML) proteins play regulatory roles in plant growth and development, responses to biotic and abiotic stresses, and other biological processes. As a popular fruit and ornamental crop, it is important to explore the regulatory mechanism of flower and fruit development of passion fruit.ResultsIn this study, 32 PeCaM/PeCML genes were identified from passion fruit genome and were divided into 9 groups based on phylogenetic analysis. The structural analysis, including conserved motifs, gene structure and homologous modeling, illustrates that the PeCaM/PeCML in the same subgroup have relative conserved structural features. Collinearity analysis suggested that the expansion of the CaM/CML gene family likely took place mainly by segmental duplication, and the whole genome replication events were closely related with the rapid expansion of the gene group. PeCaM/PeCMLs were potentially required for different floral tissues development. Significantly, PeCML26 had extremely high expression levels during ovule and fruit development compared with other PeCML genes, suggesting that PeCML26 had potential functions involved in the development of passion fruit flowers and fruits. The co-presence of various cis-elements associated with growth and development, hormone responsiveness, and stress responsiveness in the promoter regions of these PeCaM/PeCMLs might contribute to their diverse regulatory roles. Furthermore, PeCaM/PeCMLs were also induced by various abiotic stresses. This work provides a comprehensive understanding of the CaM/CML gene family and valuable clues for future studies on the function and evolution of CaM/CML genes in passion fruit.ConclusionA total of 32 PeCaM/PeCML genes were divided into 9 groups. The PeCaM/PeCML genes showed differential expression patterns in floral tissues at different development stages. It is worth noting that PeCML26, which is highly homologous to AtCaM2, not only interacts with multiple BBR-BPC TFs, but also has high expression levels during ovule and fruit development, suggesting that PeCML26 had potential functions involved in the development of passion fruit flowers and fruits. This research lays the foundation for future investigations and validation of the potential function of PeCaM/PeCML genes in the growth and development of passion fruit.

Genome-wide characterization and expression analysis of the cultivated peanut AhPR10 gene family mediating resistance to Aspergillusflavus L.

The pathogenesis-related protein PR10 is essential for plant growth, development, and stress responses. In this study, PR10 genes in cultivated peanut (Arachis hypogaea L.) were systematically identified, after which their phylogenetic relationships, conserved motifs, gene structures, and syntenic relationships were analyzed. A total of 54 AhPR10 genes were identified. They were then divided into eight groups according to their phylogenetic relationships, which were supported by the characterization of gene structures and conserved motifs. Analyses of chromosomal distribution and synteny revealed that segmental duplications were critical for the expansion of the AhPR10 gene family. In addition, the identified AhPR10 genes had constitutive and inducible expression patterns. Notably, AhPR10-7, AhPR10-33, and AhPR10-41 may have crucial functions affecting the resistance of peanut to A. flavus. In vitro fungistatic experiments indicated that recombinant AhPR10-33 can effectively inhibit A. flavus mycelial growth. The study results provide useful insights for future research on AhPR10 functions that protect peanut from the detrimental effects of A. flavus.

Detection of an 8p23.1 Inversion Using High-Resolution Optical Genome Mapping

Abstract Objective To evaluate the performance of optical genome mapping (OGM) in identifying an inversion located in the short arm of chromosome 8 (8p, 8p23.1), flanked by regions of complex segmental duplication (SD), using the GRCh38 and telomere-to-telomere (T2T) genome references. Methods We investigated a couple suspected of carrying the 8p23.1 inversion due to a terminal deletion combined with an interstitial duplication of 8p found in their abortus. OGM was performed on both individuals. The data were mapped to the current GRCh38 and the updated T2T genome references, respectively. Results The 8p23.1 inversion was observed in the female when mapping OGM data to the T2T assembly. In contrast, under the GRCh38 reference, the orientation between the suspected breakpoints within the SD regions could not be distinguished. Additional variants of uncertain significance were also identified in both individuals. Conclusion Our findings highlight the superiority of the T2T reference in recognizing structural variations involving SD regions. The enhanced SV detection using the T2T reference may contribute to a better understanding of genome instability and human diseases.

Pathogenic variants in KMT2C result in a neurodevelopmental disorder distinct from Kleefstra and Kabuki syndromes

Trithorax-related H3K4 methyltransferases, KMT2C and KMT2D, are critical epigenetic modifiers. Haploinsufficiency of KMT2C was only recently recognized as a cause of neurodevelopmental disorder (NDD), so the clinical and molecular spectrums of the KMT2C-related NDD (now designated as Kleefstra syndrome 2) are largely unknown. We ascertained 98 individuals with rare KMT2C variants, including 75 with protein-truncating variants (PTVs). Notably, ∼15% of KMT2C PTVs were inherited. Although the most highly expressed KMT2C transcript consists of only the last four exons, pathogenic PTVs were found in almost all the exons of this large gene. KMT2C variant interpretation can be challenging due to segmental duplications and clonal hematopoesis-induced artifacts. Using samples from 27 affected individuals, divided into discovery and validation cohorts, we generated a moderate strength disorder-specific KMT2C DNA methylation (DNAm) signature and demonstrate its utility in classifying non-truncating variants. Based on 81 individuals with pathogenic/likely pathogenic variants, we demonstrate that the KMT2C-related NDD is characterized by developmental delay, intellectual disability, behavioral and psychiatric problems, hypotonia, seizures, short stature, and other comorbidities. The facial module of PhenoScore, applied to photographs of 34 affected individuals, reveals that the KMT2C-related facial gestalt is significantly different from the general NDD population. Finally, using PhenoScore and DNAm signatures, we demonstrate that the KMT2C-related NDD is clinically and epigenetically distinct from Kleefstra and Kabuki syndromes. Overall, we define the clinical features, molecular spectrum, and DNAm signature of the KMT2C-related NDD and demonstrate they are distinct from Kleefstra and Kabuki syndromes highlighting the need to rename this condition.

Diagnosing missed cases of spinal muscular atrophy in genome, exome, and panel sequencing datasets.

Spinal muscular atrophy (SMA) is a genetic disorder that causes progressive degeneration of lower motor neurons and the subsequent loss of muscle function throughout the body. It is the second most common recessive disorder in individuals of European descent and is present in all populations. Accurate tools exist for diagnosing SMA from genome sequencing data. However, there are no publicly available tools for GRCh38-aligned data from panel or exome sequencing assays which continue to be used as first line tests for neuromuscular disorders. This deficiency creates a critical gap in our ability to diagnose SMA in large existing rare disease cohorts, as well as newly sequenced exome and panel datasets. We therefore developed and extensively validated a new tool - SMA Finder - that can diagnose SMA not only in genome, but also exome and panel sequencing samples aligned to GRCh37, GRCh38, or T2T-CHM13. It works by evaluating aligned reads that overlap the c.840 position of SMN1 and SMN2 in order to detect the most common molecular causes of SMA. We applied SMA Finder to 16,626 exomes and 3,911 genomes from heterogeneous rare disease cohorts sequenced at the Broad Institute Center for Mendelian Genomics as well as 1,157 exomes and 8,762 panel sequencing samples from Tartu University Hospital. SMA Finder correctly identified all 16 known SMA cases and reported nine novel diagnoses which have since been confirmed by clinical testing, with another four novel diagnoses undergoing validation. Notably, out of the 29 total SMA positive cases, 23 had an initial clinical diagnosis of muscular dystrophy, congenital myasthenic syndrome, or myopathy. This underscored the frequency with which SMA can be misdiagnosed as other neuromuscular disorders and confirmed the utility of using SMA Finder to reanalyze phenotypically diverse neuromuscular disease cohorts. Finally, we evaluated SMA Finder on 198,868 individuals that had both exome and genome sequencing data within the UK Biobank (UKBB) and found that SMA Finder's overall false positive rate was less than 1 / 200,000 exome samples, and its positive predictive value (PPV) was 97%. We also observed 100% concordance between UKBB exome and genome calls. This analysis showed that, even though it is located within a segmental duplication, the most common causal variant for SMA can be detected with comparable accuracy to monogenic disease variants in non-repetitive regions. Additionally, the high PPV demonstrated by SMA Finder, the existence of treatment options for SMA in which early diagnosis is imperative for therapeutic benefit, as well as widespread availability of clinical confirmatory testing for SMA, warrants the addition of SMN1 to the ACMG list of genes with reportable secondary findings after genome and exome sequencing.

The TIR1/AFB Family in Solanum melongena: Genome-Wide Identification and Expression Profiling under Stresses and Picloram Treatment

TIR1/AFB proteins are a class of auxin receptors with key roles in plant development and biotic and abiotic stress responses; several have been identified as targets of the auxin-mimicking herbicide picloram. In this study, we identified five putative TIR1/AFB gene family members in the important vegetable crop Solanum melongena (eggplant) and characterized them using bioinformatics tools and gene expression analyses. Phylogenetic analysis of the TIR1/AFBs classified them into three subgroups based on their Arabidopsis and Solanum lycopersicum homologs. AFB6 homologs were present only in S. melongena and S. lycopersicum, whereas AFB2/3 homologs were found only in Arabidopsis. One pair of S. melongena TIR1 homologs were located in syntenic regions in the genome and appeared to have arisen by segmental duplication. Promoter analysis revealed 898 cis-elements in the TIR1/AFB promoters, 125 of which were related to hormones, stress, light, or growth responses, but only SmAFB5 had a cis-acting regulatory element involved in auxin responsiveness (AuxRR-core). RNA sequencing and expression profiling showed that the TIR1/AFB genes were differentially expressed at different growth stages and in response to light, temperature, and drought. Only SmTIR1A expression was significantly induced by picloram treatment and different growth stages. TIR1/AFB expression is regulated by microRNAs (miRNAs) in other plant species, and we identified 6 or 29 miRNAs that potentially targeted the five TIR1/AFB genes on the basis of comparisons with S. lycopersicum and S. tuberosum miRNAs, respectively. Three-dimensional protein structure predictions revealed that all the TIR1/AFB proteins were very similar in structure, differing only in the numbers of alpha helices and in one angle linking an α helix and a β sheet. For measuring the function of TIR1/AFB genes in response to drought, SmAFB5 was selected, and knockdown by virus-induced gene silence (VIGS) 35S::SmAFB5 lines showed resistance to drought compared to controls. These analyses provide insight into the potential functions of TIR1/AFBs during growth and in response to stress; they highlight differences among the SmTIR1/AFBs that may be useful for eggplant breeding.

Genome-wide identification, expression analysis of WRKY transcription factors in Citrus ichangensis and functional validation of CiWRKY31 in response to cold stress

BackgroundIchang papeda (Citrus ichangensis), a wild perennial plant of the Rutaceae family, is a cold-hardy plant. WRKY transcription factors are crucial regulators of plant growth and development as well as abiotic stress responses. However, the WRKY genes in C. ichangensis (CiWRKY) and their expression patterns under cold stress have not been thoroughly investigated, hindering our understanding of their role in cold tolerance.ResultsIn this study, a total of 52 CiWRKY genes identified in the genome of C. ichangensis were classified into three main groups and five subgroups based on phylogenetic analysis. Comprehensive analyses of motif features, conserved domains, and gene structures were performed. Segmental duplication plays a significant role in the CiWRKY gene family expansion. Cis-acting element analysis revealed the presence of various stress-responsive elements in the promoters of the majority of CiWRKYs. Gene ontology (GO) analysis and protein-protein interaction predictions indicate that the CiWRKYs exhibit crucial roles in regulation of both development and stress response. Expression profiling analysis demonstrates that 14 CiWRKYs were substantially induced under cold stress. Virus-induced gene silencing (VIGS) assay confirmed that CiWRKY31, one of the cold-induced WRKYs, functions positively in regulation of cold tolerance.ConclusionSequence and protein properties of CiWRKYs were systematically analyzed. Among the 52 CiWRKY genes 14 members exhibited cold-responsive expression patterns, and CiWRKY31 was verified to be a positive regulator of cold tolerance. These findings pave way for future investigations to understand the molecular functions of CiWRKYs in cold tolerance and contribute to unravelling WRKYs that may be used for engineering cold tolerance in citrus.

Genome-wide identification of the Q-type C2H2 zinc finger protein gene family and expression analysis under abiotic stress in lotus (Nelumbo nucifera G.)

BackgroundLotus (Nelumbo nucifera G.) is an important aquatic plant with high ornamental, economic, cultural and ecological values, but abiotic stresses seriously affect its growth and distribution. Q-type C2H2 zinc finger proteins (ZFPs) play an important role in plant growth development and environmental stress responses. Although the Q-type C2H2 gene family has been identified in some plants, limited reports has been carried out it in lotus.ResultsIn this study, we identified 45 Q-type NnZFP members in lotus. Based on the phylogenetic tree, these Q-type NnZFP gene family members were divided into 4 groups, including C1-1i, C1-2i, C1-3i and C1-4i. Promoter cis-acting elements analysis indicated that most Q-type NnZFP gene family members in lotus were associated with response to abiotic stresses. Through collinearity analyses, no tandem duplication gene pairs and 14 segmental duplication gene pairs were identified, which showed that duplication events might play a key role in the expansion of the Q-type NnZFP gene family. The synteny results suggested that 54 and 28 Q-type NnZFP genes were orthologous to Arabidopsis and rice, respectively. The expression patterns of these Q-type NnZFP genes revealed that 30 Q-type NnZFP genes were expressed in at least one lotus tissue. Nn5g30550 showed relatively higher expression levels in all tested tissues. 12 genes were randomly selected with at least one gene from each phylogenetic clade, and the expression of these selected genes were confirmed by qRT-PCR (quantitative real-time polymerase chain reaction). The results indicated that Q-type NnZFP genes were extensively involved in cadmium, drought, salt and cold stresses responses. Among them, 11 genes responded to at least three different stress treatments, especially Nn2g12894, which induced by all four treatments.ConclusionsThese results could increase our understanding of the characterization of the Q-type NnZFP gene family and provide relevant information for further functional analysis of Q-type NnZFP genes in plant development, and abiotic stress tolerance in lotus.

Perforation in duplicated intestinal segment: a rare presentation

Perforation in duplicated intestinal segment: a rare presentation