Centromere sequence-independent but biased loading of subgenome-specific CENH3 variants in allopolyploid Arabidopsis suecica.

Centromeric nucleosomes are determined by the replacement of the canonical histone H3 with the centromere-specific histone H3 (CENH3) variant. Little is known about the centromere organization in allopolyploid species where different subgenome-specific CENH3s and subgenome-specific centromeric sequences coexist. Here, we analyzed the transcription and centromeric localization of subgenome-specific CENH3 variants in the allopolyploid species Arabidopsis suecica. Synthetic A. thaliana x A. arenosa hybrids were generated and analyzed to mimic the early evolution of A. suecica. Our expression analyses indicated that CENH3 has generally higher expression levels in A. arenosa compared to A. thaliana, and this pattern persists in the hybrids. We also demonstrated that despite a different centromere DNA composition, the centromeres of both subgenomes incorporate CENH3 encoded by both subgenomes, but with a positive bias towards the A. arenosa-type CENH3. The intermingled arrangement of both CENH3 variants demonstrates centromere plasticity and may be an evolutionary adaption to handle more than one CENH3 variant in the process of allopolyploidization.

Rational design engineering of a more thermostable Sulfurihydrogenibium yellowstonense carbonic anhydrase for potential application in carbon dioxide capture technologies.

Implementing hyperthermostable carbonic anhydrases into CO2 capture and storage technologies in order to increase the rate of CO2 absorption from the industrial flue gases is of great importance from technical and economical points of view. The present study employed a combination of in silico tools to further improve thermostability of a known thermostable carbonic anhydrase from Sulfurihydrogenibium yellowstonense. Experimental results showed that our rationally engineered K100G mutant not only retained the overall structure and catalytic efficiency but also showed a 3 °C increase in the melting temperature and a two-fold improvement in the enzyme half-life at 85 °C. Based on the molecular dynamics simulation results, rearrangement of salt bridges and hydrogen interactions network causes a reduction in local flexibility of the K100G variant. In conclusion, our study demonstrated that thermostability can be improved through imposing local structural rigidity by engineering a single-point mutation on the surface of the enzyme.

A holistic view of muscle metabolic reprogramming through personalized metabolic modeling in newly diagnosed diabetic patients.

Type 2 diabetes mellitus (T2DM) is a challenging and progressive metabolic disease caused by insulin resistance. Skeletal muscle is the major insulin-sensitive tissue that plays a pivotal role in blood sugar homeostasis. Dysfunction of muscle metabolism is implicated in the disturbance of glucose homeostasis, the development of insulin resistance, and T2DM. Understanding metabolism reprogramming in newly diagnosed patients provides opportunities for early diagnosis and treatment of T2DM as a challenging disease to manage. Here, we applied a system biology approach to investigate metabolic dysregulations associated with the early stage of T2DM. We first reconstructed a human muscle-specific metabolic model. The model was applied for personalized metabolic modeling and analyses in newly diagnosed patients. We found that several pathways and metabolites, mainly implicating in amino acids and lipids metabolisms, were dysregulated. Our results indicated the significance of perturbation of pathways implicated in building membrane and extracellular matrix (ECM). Dysfunctional metabolism in these pathways possibly interrupts the signaling process and develops insulin resistance. We also applied a machine learning method to predict potential metabolite markers of insulin resistance in skeletal muscle. 13 exchange metabolites were predicted as the potential markers. The efficiency of these markers in discriminating insulin-resistant muscle was successfully validated.

HDAC3 Knockdown Dysregulates Juvenile Hormone and Apoptosis-Related Genes in Helicoverpa armigera.

Insect development requires genes to be expressed in strict spatiotemporal order. The dynamic regulation of genes involved in insect development is partly orchestrated by the histone acetylation-deacetylation via histone acetyltransferases (HATs) and histone deacetylases (HDACs). Although histone deacetylase 3 (HDAC3) is required for mice during early embryonic development, its functions in Helicoverpa armigera (H. armigera) and its potential to be used as a target of insecticides remain unclear. We treated H. armigera with HDAC3 siRNA and RGFP966, a specific inhibitor, examining how the HDAC3 loss-of-function affects growth and development. HDAC3 siRNA and RGFP966 treatment increased mortality at each growth stage and altered metamorphosis, hampering pupation and causing abnormal wing development, reduced egg production, and reduced hatching rate. We believe that the misregulation of key hormone-related genes leads to abnormal pupa development in HDAC3 knockout insects. RNA-seq analysis identified 2788 differentially expressed genes (≥two-fold change; p ≤ 0.05) between siHDAC3- and siNC-treated larvae. Krüppel homolog 1 (Kr-h1), was differentially expressed in HDAC3 knockdown larvae. Pathway-enrichment analysis revealed the significant enrichment of genes involved in the Hippo, MAPK, and Wnt signaling pathways following HDAC3 knockdown. Histone H3K9 acetylation was increased in H. armigera after siHDAC3 treatment. In conclusion, HDAC3 knockdown dysregulated juvenile hormone (JH)-related and apoptosis-related genes in H. armigera. The results showed that the HDAC3 gene is a potential target for fighting H. armigera.

EARN: an ensemble machine learning algorithm to predict driver genes in metastatic breast cancer.

BACKGROUND: Today, there are a lot of markers on the prognosis and diagnosis of complex diseases such as primary breast cancer. However, our understanding of the drivers that influence cancer aggression is limited. METHODS: In this work, we study somatic mutation data consists of 450 metastatic breast tumor samples from cBio Cancer Genomics Portal. We use four software tools to extract features from this data. Then, an ensemble classifier (EC) learning algorithm called EARN (Ensemble of Artificial Neural Network, Random Forest, and non-linear Support Vector Machine) is proposed to evaluate plausible driver genes for metastatic breast cancer (MBCA). The decision-making strategy for the proposed ensemble machine is based on the aggregation of the predicted scores obtained from individual learning classifiers to be prioritized homo sapiens genes annotated as protein-coding from NCBI. RESULTS: This study is an attempt to focus on the findings in several aspects of MBCA prognosis and diagnosis. First, drivers and passengers predicted by SVM, ANN, RF, and EARN are introduced. Second, biological inferences of predictions are discussed based on gene set enrichment analysis. Third, statistical validation and comparison of all learning methods are performed by some evaluation metrics. Finally, the pathway enrichment analysis (PEA) using ReactomeFIVIz tool (FDR < 0.03) for the top 100 genes predicted by EARN leads us to propose a new gene set panel for MBCA. It includes HDAC3, ABAT, GRIN1, PLCB1, and KPNA2 as well as NCOR1, TBL1XR1, SIRT4, KRAS, CACNA1E, PRKCG, GPS2, SIN3A, ACTB, KDM6B, and PRMT1. Furthermore, we compare results for MBCA to other outputs regarding 983 primary tumor samples of breast invasive carcinoma (BRCA) obtained from the Cancer Genome Atlas (TCGA). The comparison between outputs shows that ROC-AUC reaches 99.24% using EARN for MBCA and 99.79% for BRCA. This statistical result is better than three individual classifiers in each case. CONCLUSIONS: This research using an integrative approach assists precision oncologists to design compact targeted panels that eliminate the need for whole-genome/exome sequencing. The schematic representation of the proposed model is presented as the Graphic abstract.

Unraveling the molecular heterogeneity in type 2 diabetes: a potential subtype discovery followed by metabolic modeling.

BACKGROUND: Type 2 diabetes mellitus (T2DM) is a complex multifactorial disease with a high prevalence worldwide. Insulin resistance and impaired insulin secretion are the two major abnormalities in the pathogenesis of T2DM. Skeletal muscle is responsible for over 75% of the glucose uptake and plays a critical role in T2DM. Here, we sought to provide a better understanding of the abnormalities in this tissue. METHODS: The muscle gene expression patterns were explored in healthy and newly diagnosed T2DM individuals using supervised and unsupervised classification approaches. Moreover, the potential of subtyping T2DM patients was evaluated based on the gene expression patterns. RESULTS: A machine-learning technique was applied to identify a set of genes whose expression patterns could discriminate diabetic subjects from healthy ones. A gene set comprising of 26 genes was found that was able to distinguish healthy from diabetic individuals with 94% accuracy. In addition, three distinct clusters of diabetic patients with different dysregulated genes and metabolic pathways were identified. CONCLUSIONS: This study indicates that T2DM is triggered by different cellular/molecular mechanisms, and it can be categorized into different subtypes. Subtyping of T2DM patients in combination with their real clinical profiles will provide a better understanding of the abnormalities in each group and more effective therapeutic approaches in the future.

Beyond Taxonomic Analysis of Microbiomes: A Functional Approach for Revisiting Microbiome Changes in Colorectal Cancer.

Colorectal cancer (CRC) is one of the most prevalent cancers in the world, especially in developed countries. In different studies, the association between CRC and dysbiosis of gut microbiome has been reported. However, most of these works focus on the taxonomic variation of the microbiome, which presents little, if any, functional insight about the reason behind and/or consequences of microbiome dysbiosis. In this study, we used a previously reported metagenome dataset which is obtained by sequencing 156 microbiome samples of healthy individuals as the control group (Co), as well as microbiome samples of patients with advanced colorectal adenoma (Ad) and colorectal carcinoma (Ca). Features of the microbiome samples have been analyzed at the level of species, as well as four functional levels, i.e., gene, KEGG orthology (KO) group, Enzyme Commission (EC) number, and reaction. It was shown that, at each of these levels, certain features exist which show significant changing trends during cancer progression. In the next step, a list of these features were extracted, which were shown to be able to predict the category of Co, Ad, and Ca samples with an accuracy of >85%. When only one group of features (species, gene, KO group, EC number, reaction) was used, KO-related features were found to be the most successful features for classifying the three categories of samples. Notably, species-related features showed the least success in sample classification. Furthermore, by applying an independent test set, we showed that these performance trends are not limited to our original dataset. We determined the most important classification features at each of the four functional levels. We propose that these features can be considered as biomarkers of CRC progression. Finally, we show that the intra-diversity of each sample at the levels of bacterial species and genes is much more than those of the KO groups, EC numbers, and reactions of that sample. Therefore, we conclude that the microbiome diversity at the species level, or gene level, is not necessarily associated with the diversity at the functional level, which again indicates the importance of KO-, EC-, and reaction-based features in metagenome analysis. The source code of proposed method is freely available from https://www.bioinformatics.org/mamed.

Cancerouspdomains: comprehensive analysis of cancer type-specific recurrent somatic mutations in proteins and domains.

BACKGROUND: Discriminating driver mutations from the ones that play no role in cancer is a severe bottleneck in elucidating molecular mechanisms underlying cancer development. Since protein domains are representatives of functional regions within proteins, mutations on them may disturb the protein functionality. Therefore, studying mutations at domain level may point researchers to more accurate assessment of the functional impact of the mutations. RESULTS: This article presents a comprehensive study to map mutations from 29 cancer types to both sequence- and structure-based domains. Statistical analysis was performed to identify candidate domains in which mutations occur with high statistical significance. For each cancer type, the corresponding type-specific domains were distinguished among all candidate domains. Subsequently, cancer type-specific domains facilitated the identification of specific proteins for each cancer type. Besides, performing interactome analysis on specific proteins of each cancer type showed high levels of interconnectivity among them, which implies their functional relationship. To evaluate the role of mitochondrial genes, stem cell-specific genes and DNA repair genes in cancer development, their mutation frequency was determined via further analysis. CONCLUSIONS: This study has provided researchers with a publicly available data repository for studying both CATH and Pfam domain regions on protein-coding genes. Moreover, the associations between different groups of genes/domains and various cancer types have been clarified. The work is available at http://www.cancerouspdomains.ir .

Rye B chromosomes encode a functional Argonaute-like protein with in vitro slicer activities similar to its A chromosome paralog.

B chromosomes (Bs) are supernumerary, dispensable parts of the nuclear genome, which appear in many different species of eukaryote. So far, Bs have been considered to be genetically inert elements without any functional genes. Our comparative transcriptome analysis and the detection of active RNA polymerase II (RNAPII) in the proximity of B chromatin demonstrate that the Bs of rye (Secale cereale) contribute to the transcriptome. In total, 1954 and 1218 B-derived transcripts with an open reading frame were expressed in generative and vegetative tissues, respectively. In addition to B-derived transposable element transcripts, a high percentage of short transcripts without detectable similarity to known proteins and gene fragments from A chromosomes (As) were found, suggesting an ongoing gene erosion process. In vitro analysis of the A- and B-encoded AGO4B protein variants demonstrated that both possess RNA slicer activity. These data demonstrate unambiguously the presence of a functional AGO4B gene on Bs and that these Bs carry both functional protein coding genes and pseudogene copies. Thus, B-encoded genes may provide an additional level of gene control and complexity in combination with their related A-located genes. Hence, physiological effects, associated with the presence of Bs, may partly be explained by the activity of B-located (pseudo)genes.

Point mutation impairs centromeric CENH3 loading and induces haploid plants.

The chromosomal position of the centromere-specific histone H3 variant CENH3 (also called "CENP-A") is the assembly site for the kinetochore complex of active centromeres. Any error in transcription, translation, modification, or incorporation can affect the ability to assemble intact CENH3 chromatin and can cause centromere inactivation [Allshire RC, Karpen GH (2008) Nat Rev Genet 9 (12):923-937]. Here we show that a single-point amino acid exchange in the centromere-targeting domain of CENH3 leads to reduced centromere loading of CENH3 in barley, sugar beet, and Arabidopsis thaliana. Haploids were obtained after cenh3 L130F-complemented cenh3-null mutant plants were crossed with wild-type A. thaliana. In contrast, in a noncompeting situation (i.e., centromeres possessing only mutated or only wild-type CENH3), no uniparental chromosome elimination occurs during early embryogenesis. The high degree of evolutionary conservation of the identified mutation site offers promising opportunities for application in a wide range of crop species in which haploid technology is of interest.

The differential loading of two barley CENH3 variants into distinct centromeric substructures is cell type- and development-specific.

The organization of centromeric chromatin of diploid barley (Hordeum vulgare) encoding two (α and ß) CENH3 variants was analysed by super-resolution microscopy. Antibody staining revealed that both CENH3 variants are organized in distinct but intermingled subdomains in interphase, mitotic and meiotic centromeres. Artificially extended chromatin fibres illustrate that these subdomains are formed by polynucleosome clusters. Thus, a CENH3 variant-specific loading followed by the arrangement into specific intermingling subdomains forming the centromere region appears. The CENH3 composition and transcription vary among different tissues. In young embryos, most interphase centromeres are composed of both CENH3 variants, while in meristematic root cells, a high number of nuclei contain ßCENH3 mainly dispersed within the nucleoplasm. A similar distribution and no preferential arrangement of the two CENH3 variants in relationship to the spindle poles suggest that both homologs meet the same function in metaphase cells.

Evolution and biology of supernumerary B chromosomes.

B chromosomes (Bs) are dispensable components of the genome exhibiting non-Mendelian inheritance and have been widely reported on over several thousand eukaryotes, but still remain an evolutionary mystery ever since their first discovery over a century ago [1]. Recent advances in genome analysis have significantly improved our knowledge on the origin and composition of Bs in the last few years. In contrast to the prevalent view that Bs do not harbor genes, recent analysis revealed that Bs of sequenced species are rich in gene-derived sequences. We summarize the latest findings on supernumerary chromosomes with a special focus on the origin, DNA composition, and the non-Mendelian accumulation mechanism of Bs.

Formation and expression of pseudogenes on the B chromosome of rye.

B chromosomes (Bs) are dispensable components of the genomes of numerous species. In contrast with the prevalent view that Bs do not harbor genes, our recent sequence analysis revealed that Bs of rye (Secale cereale) are rich in gene-derived sequences. We compared these gene-like fragments of the rye B with their ancestral A-located counterparts and confirmed an A chromosomal origin and the pseudogenization of B-located gene-like fragments. About 15% of the pseudogene-like fragments on Bs are transcribed in a tissue-type and genotype-specific manner. In addition, B-located sequences can cause in trans down- or upregulation of A chromosome-encoded genic fragments. Phenotypes and effects associated with the presence of Bs might be explained by the activity of B-located pseudogenes. We propose a model for the evolution of B-located pseudogenes.

High-copy sequences reveal distinct evolution of the rye B chromosome.

B chromosomes (Bs) are supernumerary chromosomes that vary in number among individuals of the same species. Because of their dispensable nature, their non-Mendelian inheritance and their origin from A chromosomes (As), one might assume that Bs followed a different evolutionary pathway from As, this being reflected in differences in their high-copy DNA constitution. We provide detailed insight into the composition and distribution of rye (Secale cereale) B-located high-copy sequences. A- and B-specific high-copy sequences were identified in silico. Mobile elements and satellite sequences were verified by fluorescence in situ hybridization (FISH). Replication was analyzed via EdU incorporation. Although most repeats are similarly distributed along As and Bs, several transposons are either amplified or depleted on the B. An accumulation of B-enriched satellites was found mostly in the nondisjunction control region of the B, which is transcriptionally active and late-replicating. All B-enriched sequences are not unique to the B but are also present in other Secale species, suggesting the origin of the B from As of the same genus. Our findings highlight the differences between As and Bs. Although Bs originated from As, they have since taken a separate evolutionary pathway.

Selfish supernumerary chromosome reveals its origin as a mosaic of host genome and organellar sequences.

Supernumerary B chromosomes are optional additions to the basic set of A chromosomes, and occur in all eukaryotic groups. They differ from the basic complement in morphology, pairing behavior, and inheritance and are not required for normal growth and development. The current view is that B chromosomes are parasitic elements comparable to selfish DNA, like transposons. In contrast to transposons, they are autonomously inherited independent of the host genome and have their own mechanisms of mitotic or meiotic drive. Although B chromosomes were first described a century ago, little is known about their origin and molecular makeup. The widely accepted view is that they are derived from fragments of A chromosomes and/or generated in response to interspecific hybridization. Through next-generation sequencing of sorted A and B chromosomes, we show that B chromosomes of rye are rich in gene-derived sequences, allowing us to trace their origin to fragments of A chromosomes, with the largest parts corresponding to rye chromosomes 3R and 7R. Compared with A chromosomes, B chromosomes were also found to accumulate large amounts of specific repeats and insertions of organellar DNA. The origin of rye B chromosomes occurred an estimated ∼1.1-1.3 Mya, overlapping in time with the onset of the genus Secale (1.7 Mya). We propose a comprehensive model of B chromosome evolution, including its origin by recombination of several A chromosomes followed by capturing of additional A-derived and organellar sequences and amplification of B-specific repeats.

Intraspecific hybrids of Arabidopsis thaliana revealed no gross alterations in endopolyploidy, DNA methylation, histone modifications and transcript levels.

Arabidopsis accessions Col-0 and C24 and their reciprocal hybrids were employed as a model system to investigate the potential relationship between changes in DNA methylation, chromatin structure, endopolyploidization and gene expression in heterotic genotypes. Nucleolus size, endopolyploidization level and distribution of DNA and histone H3 methylation at the microscopic level does not differ between parents and their hybrids. Methylation sensitive amplified polymorphism revealed a largely constant pattern of DNA methylation (97% of signals analyzed) after intraspecific crosses. The parental expression profile of selected genes was maintained in hybrid offspring. No correlation was found between expression pattern and DNA methylation levels at restriction sites within 5' regulatory regions. Thus, the results revealed only minor changes of chromatin properties and other nuclear features in response to intraspecific hybridization in Arabidopsis thaliana.