Intron Positions Research Articles

Whole-Genome Sequence Resource of Fusarium oxysporum Strain TH15, a Plant Growth Promoting Endophytic Fungus Isolated from Tetrastigma hemsleyanum

Whole-Genome Sequence Resource of <i>Fusarium oxysporum</i> Strain TH15, a Plant Growth Promoting Endophytic Fungus Isolated from <i>Tetrastigma hemsleyanum</i>

Comparative mitochondrial genome analyses reveal conserved gene arrangement but massive expansion/contraction in two closely related Exserohilum pathogens

Exserohilum turcicum and E. rostratum, two closely related fungal species, are both economically important pathogens but have quite different target hosts (specific to plants and cross-kingdom infection, respectively). In the present study, complete circular mitochondrial genomes of the two Exserohilum species were sequenced and de novo assembled, which mainly comprised the same set of 13 core protein-coding genes (PCGs), two rRNAs, and a certain number of tRNAs and unidentified open reading frames (ORFs). Comparative analyses indicated that these two fungi had significant mitogenomic collinearity and consistent mitochondrial gene arrangement, yet with vastly different mitogenome sizes, 264,948 bp and 64,620 bp, respectively. By contrast with the 17 introns containing 17 intronic ORFs (one-to-one) in the E. rostratum mitogenome, E. turcicum involved far more introns (70) and intronic ORFs (126), which was considered as the main contributing factors of their mitogenome expansion/contraction. Within the generally intron-rich gene cox1, a total of 18 and 10 intron position classes (Pcls) were identified separately in the two mitogenomes. Moreover, 16.16% and 10.85% ratios of intra-mitogenomic repetitive regions were detected in E. turcicum and E. rostratum, respectively. Based on the combined mitochondrial gene dataset, we established a well-supported topology of phylogeny tree of 98 ascomycetes, implying that mitogenomes may act as an effective molecular marker for fungal phylogenetic reconstruction. Our results served as the first report on mitogenomes in the genus Exserohilum, and would have significant implications in understanding the origin, evolution and pathogenic mechanisms of this fungal lineage.

Novel Lentiviral Vectors for Gene Therapy of Sickle Cell Disease Combining Gene Addition and Gene Silencing Strategies

Sickle cell disease (SCD) is due to a mutation in the β-globin (HBB) gene causing the production of the sickle β S-globin chain. The sickle Hb (HbS, a 2β S2) polymerizes, leading to the formation of sickle-shaped red blood cells that cause vaso-occlusions and organ damage. Transplantation of autologous hematopoietic stem/progenitor cells (HSPCs) transduced with lentiviral vectors (LV) expressing an anti-sickling β-globin transgene (βAS LV) is a promising curative treatment; however, it is partially effective in SCD patients, who still present elevated HbS levels.Here, we aim to improve LVs to boost therapeutic β-like globin levels without increasing the mutagenic vector load in HSPCs. We developed 2 novel LVs expressing βAS together with an artificial microRNA (amiR) targeting either the fetal Hb (HbF) repressor BCL11A (βAS/amiRBCL11A) or the β S-globin (βAS/amiRHBB). By downregulating BCL11A, amiRBCL11A re-activates the expression of the endogenous anti-sickling fetal γ-globin, which, together with βAS, should improve the clinical course of SCD; β S-globin downregulation should favor βAS incorporation in Hb tetramers, increase therapeutic Hb levels and ameliorate the SCD phenotype.First, we developed βAS/amiRBCL11A LV by inserting the amiR in multiple position of the βAS intron 2 under the control of HBB promoter/enhancers to limit BCL11A downregulation to the erythroid lineage and reduce potential amiR-related cellular toxicity and off-target effects. We showed that amiR insertion site did not affect LV titer nor βAS expression in a human erythroid cell line (HUDEP2). BCL11A downregulation in HUDEP2 led to γ-globin gene de-repression and a high proportion of HbF + cells (RTqPCR, HPLC, flow cytometry). Importantly, the total amount of therapeutic β-like globins was substantially higher in βAS/amiRBCL11A LV- than in βAS LV-transduced cells, with no impairment in cell viability or erythroid differentiation.In parallel, we designed 17 amiRs targeting HBB and generated the corresponding βAS/amiRHBB LVs. We tested these LVs in HUDEP2 and selected 2 amiRs efficiently downregulating β-globin at mRNA and protein levels (RT-qPCR and Western Blot). Of note, we modified the βAS transgene by inserting silent mutations that prevent its recognition by the amiR (βASm).Finally, we tested βAS/amiRBCL11A and βAS/amiRHBB LVs in HSPCs from SCD patients. HSPC-derived erythroid cells transduced with βAS/amiRBCL11A LV showed increased HbF levels, although HbS levels remained high. To further reduce β S-globin levels, we targeted the β S-globin mRNA using the βAS/amiRHBB LV. Efficient HSPC transduction by βASm/amiRHBB LV led to a substantial decrease of β S-globin transcripts in HSPC-derived erythroid cells compared to the βAS LV-transduced cells (RTqPCR) at a VCN/cell of 2. Notably, the amiR specifically down-regulated β S-globin, without affecting βAS expression. In βASm/amiRHBB- vs βAS LV-transduced cells, HPLC analysis showed that β S-globin downregulation led to a significant decrease of HbS, which represented 58% and 71% of the total Hb, respectively). This was associated with a significant increase of the therapeutic Hb in βASm/amiRHBB LV- vs βAS LV-transduced erythroid cells (38% and 27% of the total Hb, respectively). Importantly, we observed a substantial reduction of the proportion of HbS-positive cells in βASm/amiRHBB- vs βAS LV-transduced samples (from 96% to 70%; Figure 1A). The increased incorporation of βAS in Hb tetramers and the decrease in β S-globin led to a better correction of the sickling phenotype in mature RBCs derived from HSPCs transduced with βASm/amiRHBB LV- compared to βAS LV (55% and 84% of sickling cells, respectively; Figure 1B). A clonal assay of hematopoietic progenitors showed no impairment in HSPC viability and differentiation towards the erythroid and myeloid lineages upon transduction with bifunctional LVs. βASm/amiRHBB LV showed a standard lentiviral integration profile. Finally, we performed RNAseq to further evaluate the safety of our therapeutic strategy.In conclusion, we created a LV able to concomitantly silence the β S-globin and express βAS, achieving clinically relevant levels of therapeutic Hb and efficient correction of the sickling phenotype. Therefore, the combination of gene addition and gene silencing strategies can improve the efficacy of current therapeutic approaches, representing a novel treatment for SCD. [Display omitted] DisclosuresCavazzana: Smart Immune: Other: co-founder.

Nucleotides in both donor and acceptor splice sites are responsible for choice in NAGNAG tandem splice sites.

Among alternative splicing events in the human transcriptome, tandem NAGNAG acceptor splice sites represent an appreciable proportion. Both proximal and distal NAG can be used to produce two splicing isoforms differing by three nucleotides. In some cases, the upstream exon can be alternatively spliced as well, which further increases the number of possible transcripts. In this study, we showed that NAG choice in tandem splice site depends considerably not only on the concerned acceptor, but also on the upstream donor splice site sequence. Using an extensive set of experiments with systematically modified two-exonic minigene systems of AFAP1L2 or CSTD gene, we recognized the third and fifth intronic upstream donor splice site position and the tandem acceptor splice site region spanning from -10 to +2, including NAGNAG itself, as the main drivers. In addition, competition between different branch points and their composition were also shown to play a significant role in NAG choice. All these nucleotide effects appeared almost additive, which explained the high variability in proximal versus distal NAG usage.

Modeling in yeast how rDNA introns slow growth and increase desiccation tolerance in lichens.

We connect ribosome biogenesis to desiccation tolerance in lichens, widespread symbioses between specialized fungi (mycobionts) and unicellular phototrophs. We test whether the introns present in the nuclear ribosomal DNA of lichen mycobionts contribute to their anhydrobiosis. Self-splicing introns are found in the rDNA of several eukaryotic microorganisms, but most introns populating lichen rDNA are unable to self-splice, being either catalytically impaired group I introns, or spliceosomal introns ectopically present in rDNA. Although the mycobiont’s splicing machinery removes all introns from rRNA, Northern analysis indicates delayed post-transcriptional removal during rRNA processing, suggesting interference with ribosome assembly. To study the effects of lichen introns in a model system, we used CRISPR to introduce a spliceosomal rDNA intron from the lichen fungus Cladonia grayi into all nuclear rDNA copies of Saccharomyces cerevisiae, which lacks rDNA introns. Three intron-bearing yeast mutants were constructed with the intron inserted either in the 18S rRNA genes, the 25S rRNA genes, or in both. The mutants removed the introns correctly but had half the rDNA genes of the wildtype, grew 4.4–6 times slower, and were 40–1700 times more desiccation tolerant depending on intron position and number. Intracellular trehalose, a disaccharide implicated in desiccation tolerance, was detected at low concentration. Our data suggest that the interference of the splicing machinery with ribosome assembly leads to fewer ribosomes and proteins and to slow growth and increased desiccation tolerance in the yeast mutants. The relevance of these findings for slow growth and desiccation tolerance in lichens is discussed.

IntSplice2: Prediction of the Splicing Effects of Intronic Single-Nucleotide Variants Using LightGBM Modeling.

Prediction of the effect of a single-nucleotide variant (SNV) in an intronic region on aberrant pre-mRNA splicing is challenging except for an SNV affecting the canonical GU/AG splice sites (ss). To predict pathogenicity of SNVs at intronic positions −50 (Int-50) to −3 (Int-3) close to the 3’ ss, we developed light gradient boosting machine (LightGBM)-based IntSplice2 models using pathogenic SNVs in the human gene mutation database (HGMD) and ClinVar and common SNVs in dbSNP with 0.01 ≤ minor allelic frequency (MAF) < 0.50. The LightGBM models were generated using features representing splicing cis-elements. The average recall/sensitivity and specificity of IntSplice2 by fivefold cross-validation (CV) of the training dataset were 0.764 and 0.884, respectively. The recall/sensitivity of IntSplice2 was lower than the average recall/sensitivity of 0.800 of IntSplice that we previously made with support vector machine (SVM) modeling for the same intronic positions. In contrast, the specificity of IntSplice2 was higher than the average specificity of 0.849 of IntSplice. For benchmarking (BM) of IntSplice2 with IntSplice, we made a test dataset that was not used to train IntSplice. After excluding the test dataset from the training dataset, we generated IntSplice2-BM and compared it with IntSplice using the test dataset. IntSplice2-BM was superior to IntSplice in all of the seven statistical measures of accuracy, precision, recall/sensitivity, specificity, F1 score, negative predictive value (NPV), and matthews correlation coefficient (MCC). We made the IntSplice2 web service at https://www.med.nagoya-u.ac.jp/neurogenetics/IntSplice2.

Abstract 803: Molecular diagnosis and immunotherapy of a rare lung carcinoma patient associated with PMS2 c.1144+1G&amp;gt;A mutation-driven Lynch syndrome

Abstract Background: Lynch syndrome (LS) is an autosomal dominant genetic disease caused by mutations in DNA mismatch repair (MMR) genes, leading to an increased risk of colorectal cancer and other malignancies. Tumors with MMR-deficiency often have microsatellite instability (MSI) and high mutational burden (TMB) which predict high response rate to immunotherapy. Therefore, accurate diagnosis of LS can create valuable opportunities for the clinical management of a broad spectrum of cancers. Methods: We conducted next-generation sequencing to assess both somatic and germline mutations of the proband. LS diagnosis was confirmed by mismatch repair protein immunochemistry and germline genetic testing results. Furthermore, we applied Sanger sequencing to screen the candidate variant in his family members. Results: Our proband is a 62-year-old Chinese diagnosed with invasive lung adenocarcinoma (pT1aN2M0, stage IIIA) in the right inferior lobe. Based on the near completely loss of PMS2 in cancer cells by MMR IHC, as well as instability of three MSI markers (NR21, BAT-26 and MONO-27), we suspected LS and performed a germline NGS panel analysis which revealed a G&amp;gt;A nucleotide substitution (c.1144+1G&amp;gt;A) located at the +1 position of intron 10 of the PMS2 gene. We also found that the proband has an oncogenic KRAS p.G13D somatic mutation and his young daughter with the same PMS2 variant has no cancer yet. To characterize the pathogenicity of this PMS2 germline variant, we conducted a comprehensive genetic analysis according to the ACMG/Sherloc guidelines. In silico analysis with three splice prediction programs suggested that the PMS2 c.1144+1G&amp;gt;A variant will generate aberrant splicing transcripts. The population data in the gnomAD web portal showed that the allele frequency of PMS2 c.1144+1G&amp;gt;A variant is 2/282496 without homozygotes, falling into the pathogenic range defined by the Sherloc guideline. Importantly, in a family with germline PMS2 c.1144+2T&amp;gt;A variant, which disrupts the same splice donor site as c.1144+1G&amp;gt;A, the cosegregation between this variant and LS-related cancers was clearly seen. Based on all these results, we classified the PMS2 c.1144+1G&amp;gt;A variant as likely pathogenic, and disruption of this specific PMS2 splicing site results in LS. Combination chemotherapy with cisplatin (75 mg/m2), pemetrexed (500 mg/m2), and nivolumab (200 mg) was administered after surgery every 3 weeks for a total of 4 cycles to the patient. He achieved 16 months DFS. We then provided an up-to-date genetic counseling to the patient's daughter with the same PMS2 variant. Conclusion: Our findings suggest that accurate LS diagnosis with modern pathology approaches including MSI/MMR-d testing and NGS sequencing offers extracolonic cancer patients the opportunity to reduce cancer-related morbidity and mortality through innovative treatments such as immunotherapy. Citation Format: Quanli Han, Xiaomo Li, Tonghui Ma, Qi Wang, Zhi Cui. Molecular diagnosis and immunotherapy of a rare lung carcinoma patient associated with PMS2 c.1144+1G&amp;gt;A mutation-driven Lynch syndrome [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 803.

Co-transcriptional splicing efficiencies differ within genes and between cell types.

Pre-mRNA splicing is carried out by the spliceosome and involves splice site recognition, removal of introns, and ligation of exons. Components of the spliceosome have been shown to interact with the elongating RNA polymerase II (RNAPII), which is thought to allow splicing to occur concurrently with transcription. However, little is known about the regulation and efficiency of cotranscriptional splicing in human cells. In this study, we used Bru-seq and BruChase-seq to determine the cotranscriptional splicing efficiencies of 17,000 introns expressed across six human cell lines. We found that less than half of all introns across these six cell lines were cotranscriptionally spliced. Splicing efficiencies for individual introns showed variations across cell lines, suggesting that splicing may be regulated in a cell type–specific manner. Moreover, the splicing efficiency of introns varied within genes. The efficiency of cotranscriptional splicing did not correlate with gene length, intron position, splice site strengths, or the intron/neighboring exons GC content. However, we identified binding signals from multiple RNA-binding proteins (RBPs) that correlated with splicing efficiency, including core spliceosomal machinery components—such as SF3B4, U2AF1, and U2AF2 showing higher binding signals in poorly spliced introns. In addition, multiple RBPs, such as BUD13, PUM1, and SND1, showed preferential binding in exons that flank introns with high splicing efficiencies. The nascent RNA splicing patterns presented here across multiple cell types add to our understanding of the complexity in RNA splicing, wherein RNA-binding proteins may play important roles in determining splicing outcomes in a cell type– and intron-specific manner.

Kindler epidermolysis bullosa-like skin phenotype and downregulated basement membrane zone gene expression in poikiloderma with neutropenia and a homozygous USB1 mutation

Epidermolysis bullosa (EB) is a genotypically heterogeneous group of disorders characterized by cutaneous blistering and erosions with a tremendous spectrum of severity. One of the distinct forms of EB, Kindler EB (KEB), manifests with blistering and poikiloderma; this subtype of EB is caused by mutations in the FERMT1 gene encoding kindlin-1. In this study, we investigated a patient clinically diagnosed as KEB with reduced FERMT1 gene expression and intensity of immunostaining for kindlin-1. Transmission electron microscopy showed lamina densa reduplication, frequently observed in KEB. However, no mutations were identified in FERMT1 in this patient with consanguineous parents, and this gene resided outside of genomic regions of homozygosity (ROH). Instead, whole-exome sequencing and homozygosity mapping identified a homozygous sequence variant at the +4 position of intron 2 in the USB1 gene, encoding an exoribonuclease required for processing of U6 snRNA, a critical component of spliceosomes. Examination of the patient's RNA by RNA-Seq confirmed the pathogenicity of this variant, causing aberrant splicing predicted to result in loss of function of USB1. Mutations in this gene have been reported in patients with poikiloderma and neutropenia, with a few reported cases in association with skin fragility, a condition distinct from the KEB phenotype. Transcriptome analysis revealed that several genes, expressed in the cutaneous basement membrane zone and previously associated with different subtypes of EB, were differentially downregulated at the mRNA level. EB-associated mRNA downregulation was confirmed at protein levels by skin immunofluorescence. These observations provide a novel mechanism for blistering and erosions in the skin as a result reduced presence of adhesion complexes critical for stable association of epidermis and dermis at the level of cutaneous basement membrane zone.

Glycosyltransferase B4GALNT2 as a Predictor of Good Prognosis in Colon Cancer: Lessons from Databases.

Background: glycosyltransferase B4GALNT2 and its cognate carbohydrate antigen Sda are highly expressed in normal colon but strongly downregulated in colorectal carcinoma (CRC). We previously showed that CRC patients expressing higher B4GALNT2 mRNA levels displayed longer survival. Forced B4GALNT2 expression reduced the malignancy and stemness of colon cancer cells. Methods: Kaplan–Meier survival curves were determined in “The Cancer Genome Atlas” (TCGA) COAD cohort for several glycosyltransferases, oncogenes, and tumor suppressor genes. Whole expression data of coding genes as well as miRNA and methylation data for B4GALNT2 were downloaded from TCGA. Results: the prognostic potential of B4GALNT2 was the best among the glycosyltransferases tested and better than that of many oncogenes and tumor suppressor genes; high B4GALNT2 expression was associated with a lower malignancy gene expression profile; differential methylation of an intronic B4GALNT2 gene position and miR-204-5p expression play major roles in B4GALNT2 regulation. Conclusions: high B4GALNT2 expression is a strong predictor of good prognosis in CRC as a part of a wider molecular signature that includes ZG16, ITLN1, BEST2, and GUCA2B. Differential DNA methylation and miRNA expression contribute to regulating B4GALNT2 expression during colorectal carcinogenesis.

Panorama of intron dynamics and gene rearrangements in the phylum Basidiomycota as revealed by the complete mitochondrial genome of Turbinellus floccosus.

In the present study, the complete mitogenome of Turbinellus floccosus was sequenced, assembled, and compared with other basidiomycete mitogenomes. The mitogenome of T. floccosus consists of a circular DNA molecule, with a size of 62,846 bp. Gene arrangement analysis indicated that large-scale gene rearrangements occurred in the levels of family and genus of basidiomycete species, and the mitogenome of T. floccosus contained a unique gene order. A significant correlation between the number of introns and the mitochondrial genome size of Basidiomycota were detected (P < 0.01). A total of 896 introns were detected in the core protein-coding genes (PCGs) of 74 basidiomycete species, and the cox1 gene was the largest host gene of basidiomycete introns. Intron position class (Pcls) P383 in the cox1 gene was the most common intron in Basidiomycota, which distributed in 40 of 74 basidiomycete species. In addition, frequent intron loss/gain events were detected in basidiomycete species. More than 50% of bases around insertion sites (- 15 bp to 15 bp) of Pcls from different species were conservative, indicating site preferences of intron insertions in Basidiomycota. Further analysis showed that 76.09% of introns tended to insert downstream to a T base in Basidiomycota. Phylogenetic analysis for 74 basidiomycetes indicated mitochondrial genes are effective molecular markers for phylogeny of basidiomycetes. The study served as the first report on the mitogenome from the family Gomphaceae, which will help to understand the intron origin and evolution in Basidiomycota. KEY POINTS: • The mitogenome of Turbinellus floccosus had a unique gene arrangement. • Intron loss/gain events were detected in the 74 basidiomycete species. • Introns tend to insert downstream of a T base in basidiomycete mitogenomes.

Origin and Evolution of the Human Bcl2-Associated Athanogene-1 (BAG-1).

Molecular chaperones, particularly the 70-kDa heat shock proteins (Hsp70s), are key orchestrators of the cellular stress response. To perform their critical functions, Hsp70s require the presence of specific co-chaperones, which include nucleotide exchange factors containing the BCL2-associated athanogene (BAG) domain. BAG-1 is one of these proteins that function in a wide range of cellular processes, including apoptosis, protein refolding, and degradation, as well as tumorigenesis. However, the origin of BAG-1 proteins and their evolution between and within species are mostly uncharacterized. This report investigated the macro- and micro-evolution of BAG-1 using orthologous sequences and single nucleotide polymorphisms (SNPs) to elucidate the evolution and understand how natural variation affects the cellular stress response. We first collected and analyzed several BAG-1 sequences across animals, plants, and fungi; mapped intron positions and phases; reconstructed phylogeny; and analyzed protein characteristics. These data indicated that BAG-1 originated before the animals, plants, and fungi split, yet most extant fungal species have lost BAG-1. Furthermore, although BAG-1’s structure has remained relatively conserved, kingdom-specific conserved differences exist at sites of known function, suggesting functional specialization within each kingdom. We then analyzed SNPs from the 1000 genomes database to determine the evolutionary patterns within humans. These analyses revealed that the SNP density is unequally distributed within the BAG1 gene, and the ratio of non-synonymous/synonymous SNPs is significantly higher than 1 in the BAG domain region, which is an indication of positive selection. To further explore this notion, we performed several biochemical assays and found that only one out of five mutations tested altered the major co-chaperone properties of BAG-1. These data collectively suggest that although the co-chaperone functions of BAG-1 are highly conserved and can probably tolerate several radical mutations, BAG-1 might have acquired specialized and potentially unexplored functions during the evolutionary process.

Low-Grade Appendiceal Mucinous Neoplasm (LAMN)-Case Report

The evolutionary origin and significance of spliceosomal introns have been the subject of many investigations. Two theories, "introns-early" theory and "introns-late" theory, have been proposed to explain the evolution of introns in eukaryotic genes. Intron position is generally conserved in paralogue and orthologue genes. Some introns occur at similar but not necessarily identical positions in homologous genes, which were separated by great evolutionary distances. This event can be explained by insertion, loss or movement of the intron over short distances. Intron loss and gain events are unique in evolution and can be useful as markers for phylogenetic analyses. The insertion of introns at an identical position suggests a common ancestor gene. Here we analyzed, using PCR and RT-PCR, the structure of the 1,3-beta-glucan synthase gene (FKS) in several clinical isolates of Paracoccidioides brasiliensis (Pb): isolates Pb 01, Pb 4940, Pb 8515, Pb 8311, Pb 8334, Pb 4268, Pb 1668, and Pb E. Our results showed that seven of the isolates examined showed identical structures concerning the position of introns in PbFKS1. PbFKS4940 showed the intron described at the 3' end and had lost that one at the 5' end. The presence of the PbFKS4940 transcript suggests that it could be a functional gene. These data suggest a divergent evolution for introns with regard to the 1,3-beta-glucan synthase gene in P. brasiliensis isolates.

Mitogenomes of Two Phallus Mushroom Species Reveal Gene Rearrangement, Intron Dynamics, and Basidiomycete Phylogeny.

Phallus indusiatus and Phallus echinovolvatus are edible bamboo mushrooms with pharmacological properties. We sequenced, assembled, annotated, and compared the mitogenomes of these species. Both mitogenomes were composed of circular DNA molecules, with sizes of 89,139 and 50,098 bp, respectively. Introns were the most important factor in mitogenome size variation within the genus Phallus. Phallus indusiatus, P. echinovolvatus, and Turbinellus floccosus in the subclass Phallomycetidae have conservative gene arrangements. Large-scale gene rearrangements were observed in species representing 42 different genera of Basidiomycetes. A variety of intron position classes were found in the 44 Basidiomycete species analyzed. A novel group II intron from the P. indusiatus mitogenome was compared with other fungus species containing the same intron, and we demonstrated that the insertion sites of the intron had a base preference. Phylogenetic analyses based on combined gene datasets yielded well-supported Bayesian posterior probability (BPP = 1) topologies. This indicated that mitochondrial genes are reliable molecular markers for analyzing the phylogenetic relationships of the Basidiomycetes. This is the first study of the mitogenome of the genus Phallus, and it increases our understanding of the population genetics and evolution of bamboo mushrooms and related species.

SMAD3 Hypomethylation as a Biomarker for Early Prediction of Colorectal Cancer.

The incidence and mortality rates of colorectal cancer (CRC) have been high in recent years. Prevention and early detection are crucial for decreasing the death rate. Therefore, this study aims to characterize the alteration patterns of mothers against decapentaplegic homolog 3 (SMAD3) in patients with CRC and its applications in early detection by using a genome-wide methylation array to identify an aberrant hypomethylation site in the intron position of the SMAD3 gene. Quantitative methylation-specific polymerase chain reaction showed that hypomethylated SMAD3 occurred in 91.4% (501/548) of Taiwanese CRC tissues and 66.6% of benign tubular adenoma polyps. In addition, SMAD3 hypomethylation was observed in 94.7% of patients with CRC from The Cancer Genome Atlas dataset. A decrease in circulating cell-free methylation SMAD3 was detected in 70% of CRC patients but in only 20% of healthy individuals. SMAD3 mRNA expression was low in 42.9% of Taiwanese CRC tumor tissues but high in 29.4% of tumors compared with paired adjacent normal tissues. Hypomethylated SMAD3 was found in cancers of the digestive system, such as liver cancer, gastric cancer, and colorectal cancer, but not in breast cancer, endometrial cancer, and lung cancer. In conclusion, SMAD3 hypomethylation is a potential diagnostic marker for CRC in Western and Asian populations.

The First Mitochondrial Genome for Geastrales (Sphaerobolus stellatus) Reveals Intron Dynamics and Large-Scale Gene Rearrangements of Basidiomycota.

In this study, the mitogenome of artillery fungus, Sphaerobolus stellatus, was assembled and compared with other Basidiomycota mitogenomes. The Sphaerobolus stellatus mitogenome was composed of circular DNA molecules, with a total size of 152,722 bp. Accumulation of intergenic and intronic sequences contributed to the Sphaerobolus stellatus mitogenome becoming the fourth largest mitogenome among Basidiomycota. We detected large-scale gene rearrangements in Basidiomycota mitogenomes, and the Sphaerobolus stellatus mitogenome contains a unique gene order. The quantity and position classes of intron varied between 75 Basidiomycota species we tested, indicating frequent intron loss/gain events occurred in the evolution of Basidiomycota. A novel intron position classes (P1281) was detected in the Sphaerobolus stellatus mitogenome, without any homologous introns from other Basidiomycota species. A pair of fragments with a total length of 9.12 kb in both the nuclear and mitochondrial genomes of Sphaerobolus stellatus was detected, indicating possible gene transferring events. Phylogenetic analysis based on the combined mitochondrial gene set obtained well-supported tree topologies (Bayesian posterior probabilities ≥ 0.99; bootstrap values ≥98). This study served as the first report on the mitogenome from the order Geastrales, which will promote the understanding of the phylogeny, population genetics, and evolution of the artillery fungus, Sphaerobolus stellatus.

Genome-Wide Identification and Expression Analysis of the bZIP Transcription Factors in the Mycoparasite Coniothyrium minitans.

The basic leucine zipper (bZIP) proteins family is one of the largest and most diverse transcription factors, widely distributed in eukaryotes. However, no information is available regarding the bZIP gene family in Coniothyrium minitans, an important biocontrol agent of the plant pathogen Sclerotinia sclerotiorum. In this study, we identified 34 bZIP genes from the C. minitans genome, which were classified into 8 groups based on their phylogenetic relationships. Intron analysis showed that 28 CmbZIP genes harbored a variable number of introns, and 15 of them shared a feature that intron inserted into the bZIP domain. The intron position in bZIP domain was highly conserved, which was related to recognize the arginine (R) and could be treated as a genomic imprinting. Expression analysis of the CmbZIP genes in response to abiotic stresses indicated that they might play distinct roles in abiotic stress responses. Results showed that 22 CmbZIP genes were upregulated during the later stage of conidial development. Furthermore, transcriptome analysis indicated that CmbZIP genes are involved in different stages of mycoparasitism. Among deletion mutants of four CmbZIPs (CmbZIP07, -09, -13, and -16), only ΔCmbZIP16 mutants significantly reduced its tolerance to the oxidative stress. The other mutants exhibited no significant effects on colony morphology, mycelial growth, conidiation, and mycoparasitism. Taken together, our results suggested that CmbZIP genes play important roles in the abiotic stress responses, conidial development, and mycoparasitism. These results provide comprehensive information of the CmbZIP gene family and lay the foundation for further research on the bZIP gene family regarding their biological functions and evolutionary history.

Cloning and functional characterization of a carotenoid cleavage dioxygenase 2 gene in safranal and crocin biosynthesis from Freesia hybrida

Safranal and crocin, commonly derived from the oxidative cleavage reaction of zeaxanthin in plants, are two kinds of apocarotenoids with versatile functions, which were only found in limited number of plant species. In this study, both metabolites were detected and varied concomitantly with the expression of carotenoid cleavage dioxygenase (CCD) genes in Freesia hybrida, Red River® and Ambiance cultivars. The newly isolated CCD, denoted here as FhCCD2, was phylogenetically clustered with other reported saffron CCD2s. Besides, ten introns were also observed in the genomic DNA sequence of FhCCD2 and the presence of N-terminal transporter peptide suggested its plastidial sub-localization. Biochemical analysis showed that the FhCCD2 cleaved zeaxanthin at the 7, 8 and 7′, 8’ double bonds to generate intermediates prerequisite for the biosynthesis of safranal and crocin. Further, gene transient expression analysis showed that the promoter of FhCCD2 was functional in Ambiance as well as Red River® cultivars, even with slight variation in their promoter sequence. At present, CCD2 proteins have only been found in Freesia and Crocus genus of Iridaceae family. Phylogenetic and intron position analysis infer that CCD2 perhaps emerged after the intron loss during evolutionary process of CCD1 or their shared ancestry.

Genome-Wide Identification, Evolution, and Expression of GDSL-Type Esterase/Lipase Gene Family in Soybean.

GDSL-type esterase/lipase proteins (GELPs) belong to the SGNH hydrolase superfamily and contain a conserved GDSL motif at their N-terminus. GELPs are widely distributed in nature, from microbes to plants, and play crucial roles in growth and development, stress responses and pathogen defense. However, the identification and functional analysis of GELP genes are hardly explored in soybean. This study describes the identification of 194 GELP genes in the soybean genome and their phylogenetic classification into 11 subfamilies (A–K). GmGELP genes are disproportionally distributed on 20 soybean chromosomes. Large-scale WGD/segmental duplication events contribute greatly to the expansion of the soybean GDSL gene family. The Ka/Ks ratios of more than 70% of duplicated gene pairs ranged from 0.1–0.3, indicating that most GmGELP genes were under purifying selection pressure. Gene structure analysis indicate that more than 74% of GmGELP genes are interrupted by 4 introns and composed of 5 exons in their coding regions, and closer homologous genes in the phylogenetic tree often have similar exon-intron organization. Further statistics revealed that approximately 56% of subfamily K members contain more than 4 introns, and about 28% of subfamily I members consist of less than 4 introns. For this reason, the two subfamilies were used to simulate intron gain and loss events, respectively. Furthermore, a new model of intron position distribution was established in current study to explore whether the evolution of multi-gene families resulted from the diversity of gene structure. Finally, RNA-seq data were used to investigate the expression profiles of GmGELP gene under different tissues and multiple abiotic stress treatments. Subsequently, 7 stress-responsive GmGELP genes were selected to verify their expression levels by RT-qPCR, the results were consistent with RNA-seq data. Among 7 GmGELP genes, GmGELP28 was selected for further study owing to clear responses to drought, salt and ABA treatments. Transgenic Arabidopsis thaliana and soybean plants showed drought and salt tolerant phenotype. Overexpression of GmGELP28 resulted in the changes of several physiological indicators, which allowed plants to adapt adverse conditions. In all, GmGELP28 is a potential candidate gene for improving the salinity and drought tolerance of soybean.

Comparative Mitogenome Analysis Reveals Mitochondrial Genome Differentiation in Ectomycorrhizal and Asymbiotic Amanita Species.

In this present study, we assembled and analyzed the mitogenomes of two asymbiotic and six ectomycorrhizal Amanita species based on next-generation sequencing data. The size of the eight Amanita mitogenomes ranged from 37,341 to 137,428 bp, and we considered introns to be one of the main factors contributing to the size variation of Amanita. The introns of the cox1 gene experienced frequent gain/loss events in Amanita; and the intron position class cox1P386 was lost in the six ectomycorrhizal Amanita species. In addition, ectomycorrhizal Amanita species had more repetitive sequences and fewer intergenic sequences than asymbiotic Amanita species in their mitogenomes. Large-scale gene rearrangements were detected in the Amanita species we tested, including gene displacements and inversions. On the basis of the combined mitochondrial gene set, we reconstructed the phylogenetic relationships of 66 Basidiomycetes. The six ectomycorrhizal Amanita species were of single origin, and the two saprophytic Amanita species formed two distinct clades. This study is the first to elucidate the functions of the mitogenome in the evolution and ecological adaptation of Amanita species.