Multiple routes to fungicide resistance: Interaction of Cyp51 gene sequences, copy number and expression.

We examined the molecular basis of triazole resistance in Blumeria graminis f. sp. tritici (wheat mildew, Bgt), a model organism among powdery mildews. Four genetic models for responses to triazole fungicides were identified among US and UK isolates, involving multiple genetic mechanisms. Firstly, only two amino acid substitutions in CYP51B lanosterol demethylase, the target of triazoles, were associated with resistance, Y136F and S509T (homologous to Y137F and S524T in the reference fungus Zymoseptoria tritici). As sequence variation did not explain the wide range of resistance, we also investigated Cyp51B copy number and expression, the latter using both reverse transcription-quantitative PCR and RNA-seq. The second model for resistance involved higher copy number and expression in isolates with a resistance allele; thirdly, however, moderate resistance was associated with higher copy number of wild-type Cyp51B in some US isolates. A fourth mechanism was heteroallelism with multiple alleles of Cyp51B. UK isolates, with significantly higher mean resistance than their US counterparts, had higher mean copy number, a high frequency of the S509T substitution, which was absent from the United States, and in the most resistant isolates, heteroallelism involving both sensitivity residues Y136+S509 and resistance residues F136+T509. Some US isolates were heteroallelic for Y136+S509 and F136+S509, but this was not associated with higher resistance. The obligate biotrophy of Bgt may constrain the tertiary structure and thus the sequence of CYP51B, so other variation that increases resistance may have a selective advantage. We describe a process by which heteroallelism may be adaptive when Bgt is intermittently exposed to triazoles.

Enhanced heterologous gene expression in Trichoderma reesei by promoting multicopy integration.

Trichoderma reesei displays a high capability to produce extracellular proteins and therefore is used as a platform for the expression of heterologous genes. In a previous study, an expression cassette with the constitutive tef1 promoter and the cbh1 terminator compatible with flow cytometry analysis was developed. Independent transformants obtained by a random integration into the genome of a circular plasmid containing the expression cassette showed a wide range of fluorescence levels. Whole genome sequencing was conducted on eight of the transformed strains using two next-generation sequencing (NGS) platforms: Illumina paired-end sequencing and Oxford Nanopore. In all strains, the expression plasmid was inserted at the same position in the genome, i.e., upstream of the tef1 gene, indicating an integration by homologous recombination. The different levels of fluorescence observed correspond to different copy numbers of the plasmid. Overall, the integration of a circular plasmid with the green fluorescence protein (egfp) transgene under the control of tef1 promoter favors multicopy integration and allows over-production of this heterologous protein on glucose. In conclusion, an expression system based on using the tef1 promotor could be one of the building blocks for improving high-value heterologous protein production by increasing the copy number of the encoding genes into the genome of the platform strain. KEY POINTS: • Varied eGFP levels from tef1 promoter and cbh1 terminator expression. • Whole genome sequencing on short and long reads platforms reveals various plasmid copy numbers in strains. • Plasmids integrate at the same genomic site by homologous recombination in all strains.

Comprehensive Analyses of Somatic Copy Number Alterations and Mutations Based on the Adenoma-Carcinoma Sequence.

AIMS: Identifying molecular alterations in the adenoma and carcinoma components within the same tumor would greatly contribute to understanding the neoplastic progression of early colorectal cancer. METHODS AND RESULTS: We examined somatic copy number alterations (SCNAs) and mutations involved in the adenoma and carcinoma components obtained from the same tumor in 46 cases of microsatellite-stable carcinoma in adenoma, using a genome-wide SNP array and gene mutation panel. In addition, we also performed hierarchical clustering to determine the SCNA frequencies in the tumors, resulting in stratification of the samples into two subgroups according to SCNA frequency. Subgroup 1 was characterized by multiple SCNAs and carcinoma components exclusively, while Subgroup 2 was characterized by a low frequency of SCNAs and both the adenoma and carcinoma components. The numbers of total genes and genes with gains were higher in the carcinoma than adenoma components. The three most frequent gains in both components were located at 1p36.33-1q44, 2p25.3-2q37.3, and 3p26.3-3q29. However, no candidate genes mapped to these regions. APC and KRAS mutations were common in both components, whereas the frequency of TP53 mutations was statistically higher in the carcinoma than adenoma component. However, TP53 mutations were not correlated with SCNA frequency. CONCLUSIONS: We suggest that considerable SCNAs and TP53 mutations are required for progression from adenoma to carcinoma within the same intramucosal neoplastic lesion.

Prevalent Fast Evolution of Genes Involved in Heterochromatin Functions.

Heterochromatin is a gene-poor and repeat-rich genomic compartment universally found in eukaryotes. Despite its low transcriptional activity, heterochromatin plays important roles in maintaining genome stability, organizing chromosomes, and suppressing transposable elements. Given the importance of these functions, it is expected that genes involved in heterochromatin regulation would be highly conserved. Yet, a handful of these genes were found to evolve rapidly. To investigate whether these previous findings are anecdotal or general to genes modulating heterochromatin, we compile an exhaustive list of 106 candidate genes involved in heterochromatin functions and investigate their evolution over short and long evolutionary time scales in Drosophila. Our analyses find that these genes exhibit significantly more frequent evolutionary changes, both in the forms of amino acid substitutions and gene copy number change, when compared to genes involved in Polycomb-based repressive chromatin. While positive selection drives amino acid changes within both structured domains with diverse functions and intrinsically disordered regions, purifying selection may have maintained the proportions of intrinsically disordered regions of these proteins. Together with the observed negative associations between the evolutionary rate of these genes and the genomic abundance of transposable elements, we propose an evolutionary model where the fast evolution of genes involved in heterochromatin functions is an inevitable outcome of the unique functional roles of heterochromatin, while the rapid evolution of transposable elements may be an effect rather than cause. Our study provides an important global view of the evolution of genes involved in this critical cellular domain and provides insights into the factors driving the distinctive evolution of heterochromatin.

Compensation of gene dosage on the mammalian X.

Changes in gene dosage can have tremendous evolutionary potential (e.g. whole-genome duplications), but without compensatory mechanisms, they can also lead to gene dysregulation and pathologies. Sex chromosomes are a paradigmatic example of naturally occurring gene dosage differences and their compensation. In species with chromosome-based sex determination, individuals within the same population necessarily show 'natural' differences in gene dosage for the sex chromosomes. In this Review, we focus on the mammalian X chromosome and discuss recent new insights into the dosage-compensation mechanisms that evolved along with the emergence of sex chromosomes, namely X-inactivation and X-upregulation. We also discuss the evolution of the genetic loci and molecular players involved, as well as the regulatory diversity and potentially different requirements for dosage compensation across mammalian species.

Evolutionary variation in gene conversion at the avian MHC is explained by fluctuating selection, gene copy numbers and life history.

The major histocompatibility complex (MHC) multigene family encodes key pathogen-recognition molecules of the vertebrate adaptive immune system. Hyper-polymorphism of MHC genes is de novo generated by point mutations, but new haplotypes may also arise by re-shuffling of existing variation through intra- and inter-locus gene conversion. Although the occurrence of gene conversion at the MHC has been known for decades, we still have limited understanding of its functional importance. Here, I took advantage of extensive genetic resources (~9000 sequences) to investigate broad scale macroevolutionary patterns in gene conversion processes at the MHC across nearly 200 avian species. Gene conversion was found to constitute a universal mechanism in birds, as 83% of species showed footprints of gene conversion at either MHC class and 25% of all allelic variants were attributed to gene conversion. Gene conversion processes were stronger at MHC-II than MHC-I, but inter-specific variation at both MHC classes was explained by similar evolutionary scenarios, reflecting fluctuating selection towards different optima and drift. Gene conversion showed uneven phylogenetic distribution across birds and was driven by gene copy number variation, supporting significant role of inter-locus gene conversion processes in the evolution of the avian MHC. Finally, MHC gene conversion was stronger in species with fast life histories (high fecundity) and in long-distance migrants, likely reflecting variation in population sizes and host-pathogen coevolutionary dynamics. The results provide a robust comparative framework for understanding macroevolutionary variation in gene conversion at the avian MHC and reinforce important contribution of this mechanism to functional MHC diversity.

X chromosome dosage drives statin-induced dysglycemia and mitochondrial dysfunction.

Statin drugs lower blood cholesterol levels for cardiovascular disease prevention. Women are more likely than men to experience adverse statin effects, particularly new-onset diabetes (NOD) and muscle weakness. Here we find that impaired glucose homeostasis and muscle weakness in statin-treated female mice are associated with reduced levels of the omega-3 fatty acid, docosahexaenoic acid (DHA), impaired redox tone, and reduced mitochondrial respiration. Statin adverse effects are prevented in females by administering fish oil as a source of DHA, by reducing dosage of the X chromosome or the Kdm5c gene, which escapes X chromosome inactivation and is normally expressed at higher levels in females than males. As seen in female mice, we find that women experience more severe reductions than men in DHA levels after statin administration, and that DHA levels are inversely correlated with glucose levels. Furthermore, induced pluripotent stem cells from women who developed NOD exhibit impaired mitochondrial function when treated with statin, whereas cells from men do not. These studies identify X chromosome dosage as a genetic risk factor for statin adverse effects and suggest DHA supplementation as a preventive co-therapy.

Therapeutic targeting of immunometabolism reveals a critical reliance on hexokinase 2 dosage for microglial activation and Alzheimer's progression.

Neuroinflammation is a prominent feature of Alzheimer's disease (AD). Activated microglia undergo a reprogramming of cellular metabolism necessary to power their cellular activities during disease. Thus, selective targeting of microglial immunometabolism might be of therapeutic benefit for treating AD. In the AD brain, the levels of microglial hexokinase 2 (HK2), an enzyme that supports inflammatory responses by promoting glycolysis, are significantly increased. In addition, HK2 displays non-metabolic activities that extend its inflammatory role beyond glycolysis. The antagonism of HK2 affects microglial phenotypes and disease progression in a gene-dose-dependent manner. HK2 complete loss fails to improve pathology by exacerbating inflammation, while its haploinsufficiency reduces pathology in 5xFAD mice. We propose that the partial antagonism of HK2 is effective in slowing disease progression by modulating NF-κB signaling through its cytosolic target, IKBα. The complete loss of HK2 affects additional inflammatory mechanisms related to mitochondrial dysfunction.

Relative rDNA copy number is not associated with resistance training-induced skeletal muscle hypertrophy and does not affect myotube anabolism in vitro.

Ribosomal DNA (rDNA) copies exist across multiple chromosomes, and interindividual variation in copy number is speculated to influence the hypertrophic response to resistance training. Thus, we examined if rDNA copy number was associated with resistance training-induced skeletal muscle hypertrophy. Participants (n = 53 male, 21 ± 1 yr old; n = 29 female, 21 ± 2 yr old) performed 10-12 wk of full-body resistance training. Hypertrophy outcomes were determined, as was relative rDNA copy number from preintervention vastus lateralis (VL) biopsies. Pre- and postintervention VL biopsy total RNA was assayed in all participants, and mRNA/rRNA markers of ribosome content and biogenesis were also assayed in the 29 female participants before training, 24 h following training bout 1, and in the basal state after 10 wk of training. Across all participants, no significant associations were evident between relative rDNA copy number and training-induced changes in whole body lean mass (r = -0.034, P = 0.764), vastus lateralis thickness (r = 0.093, P = 0.408), mean myofiber cross-sectional area (r = -0.128, P = 0.259), or changes in muscle RNA concentrations (r = 0.026, P = 0.818), and these trends were similar when examining each gender. However, all Pol-I regulon mRNAs as well as 45S pre-rRNA, 28S rRNA, and 18S rRNA increased 24 h following the first training bout in female participants. Follow-up studies using LHCN-M2 myotubes demonstrated that a reduction in relative rDNA copy number induced by bisphenol A did not significantly affect insulin-like-growth factor-induced myotube hypertrophy. These findings suggest that relative rDNA copy number is not associated with myofiber hypertrophy.NEW & NOTEWORTHY We examined ribosomal DNA (rDNA) copy numbers in men and women who resistance trained for 10-12 wk and found no significant associations with skeletal muscle hypertrophy outcomes. These data, along with in vitro data in immortalized human myotubes whereby rDNA copy number was reduced, provide strong evidence that relative rDNA copy number is not associated with anabolism.

ZBTB21 suppresses CRE-mediated transcription to impair synaptic function in Down syndrome.

Down syndrome (DS) is the most common chromosomal disorder and a major cause of intellectual disability. The genetic etiology of DS is the extra copy of chromosome 21 (HSA21)-encoded genes; however, the contribution of specific HSA21 genes to DS pathogenesis remains largely unknown. Here, we identified ZBTB21, an HSA21-encoded zinc-finger protein, as a transcriptional repressor in the regulation of synaptic function. We found that normalization of the Zbtb21 gene copy number in DS mice corrected deficits in cognitive performance, synaptic function, and gene expression. Moreover, we demonstrated that ZBTB21 binds to canonical cAMP-response element (CRE) DNA and that its binding to CRE could be competitive with CRE-binding factors such as CREB. ZBTB21 represses CRE-dependent gene expression and results in the negative regulation of synaptic plasticity, learning and memory. Together, our results identify ZBTB21 as a CRE-binding protein and repressor in cAMP-dependent gene regulation, contributing to cognitive defects in DS.

Cohesin composition and dosage independently affect early development in zebrafish.

Cohesin, a chromatin-associated protein complex with four core subunits (Smc1a, Smc3, Rad21 and either Stag1 or 2), has a central role in cell proliferation and gene expression in metazoans. Human developmental disorders termed 'cohesinopathies' are characterized by germline variants of cohesin or its regulators that do not entirely eliminate cohesin function. However, it is not clear whether mutations in individual cohesin subunits have independent developmental consequences. Here, we show that zebrafish rad21 or stag2b mutants independently influence embryonic tailbud development. Both mutants have altered mesoderm induction, but only homozygous or heterozygous rad21 mutation affects cell cycle gene expression. stag2b mutants have narrower notochords and reduced Wnt signaling in neuromesodermal progenitors as revealed by single-cell RNA sequencing. Stimulation of Wnt signaling rescues transcription and morphology in stag2b, but not rad21, mutants. Our results suggest that mutations altering the quantity versus composition of cohesin have independent developmental consequences, with implications for the understanding and management of cohesinopathies.

A cell competition system with one gene expression from a single-copy gene in one cell.

Even with advanced plasmid and viral vectors, attaining copy numbers of multiple genes among different transfected cells is challenging. We achieved one gene expression from a single-copy gene in one cell using a transgene competition system, a combination of the Kazusa cDNA clones and our dual recombinase-mediated cassette exchange system. All 48 nuclear receptors were simultaneously expressed in one dish at the same expression level in HEK293 using this system, and the cell proliferation rate was compared. Significant differences were observed between cells transfected with CMV- or EF1 promoter-driven expression of the 48 nuclear receptors after 8 weeks. The EF1-NR1I2 cell line, which exhibited the highest increase from 2 to 8 weeks, showed 1.13-fold higher proliferation than the EF1-DsRed line. On the other hand, the EF1-NR4A1 cell line, which showed the maximum decrease at 8 weeks, showed 0.88-fold lower proliferation than the EF1-DsRed line. The results were confirmed in both our transgene competition system and long-term growth experiments. Our transgene competition system offers a wide-range, simple, and accurate cell competition method.

Intermolecular Gene Conversion for the Equalization of Genome Copies in the Polyploid Haloarchaeon Haloferax volcanii: Identification of Important Proteins.

The model haloarchaeon Haloferax volcanii is polyploid with about 20 copies of its major chromosome. Recently it has been described that highly efficient intermolecular gene conversion operates in H. volcanii to equalize the chromosomal copies. In the current study, 24 genes were selected that encode proteins with orthologs involved in gene conversion or homologous recombination in archaea, bacteria, or eukaryotes. Single gene deletion strains of 22 genes and a control gene were constructed in two parent strains for a gene conversion assay; only radA and radB were shown to be essential. Protoplast fusions were used to generate strains that were heterozygous for the gene HVO_2528, encoding an enzyme for carotinoid biosynthesis. It was revealed that a lack of six of the proteins did not influence the efficiency of gene conversion, while sixteen mutants had severe gene conversion defects. Notably, lack of paralogous proteins of gene families had very different effects, e.g., mutant Δrad25b had no phenotype, while mutants Δrad25a, Δrad25c, and Δrad25d were highly compromised. Generation of a quadruple rad25 and a triple sph deletion strain also indicated that the paralogs have different functions, in contrast to sph2 and sph4, which cannot be deleted simultaneously. There was no correlation between the severity of the phenotypes and the respective transcript levels under non-stressed conditions, indicating that gene expression has to be induced at the onset of gene conversion. Phylogenetic trees of the protein families Rad3/25, MutL/S, and Sph/SMC/Rad50 were generated to unravel the history of the paralogous proteins of H. volcanii. Taken together, unselected intermolecular gene conversion in H. volcanii involves at least 16 different proteins, the molecular roles of which can be studied in detail in future projects.

Siderophore-harboring gut bacteria and fecal siderophore genes for predicting the responsiveness of fecal microbiota transplantation for active ulcerative colitis.

BACKGROUND: Predictive markers for fecal microbiota transplantation (FMT) outcomes in patients with active ulcerative colitis (UC) are poorly defined. We aimed to investigate changes in gut microbiota pre- and post-FMT and to assess the potential value in determining the total copy number of fecal bacterial siderophore genes in predicting FMT responsiveness. METHODS: Patients with active UC (Mayo score ≥ 3) who had undergone two FMT procedures were enrolled. Fecal samples were collected before and 8 weeks after each FMT session. Patients were classified into clinical response and non-response groups, based on their Mayo scores. The fecal microbiota profile was accessed using metagenomic sequencing, and the total siderophore genes copy number via quantitative real-time polymerase chain reaction. Additionally, we examined the association between the total siderophore genes copy number and FMT efficacy. RESULTS: Seventy patients with UC had undergone FMT. The clinical response and remission rates were 50% and 10% after the first FMT procedure, increasing to 72.41% and 27.59% after the second FMT. The cumulative clinical response and clinical remission rates were 72.86% and 25.71%. Compared with baseline, the response group showed a significant increase in Faecalibacterium, and decrease in Enterobacteriaceae, consisted with the changes of the total bacterial siderophore genes copy number after the second FMT (1889.14 vs. 98.73 copies/ng, P < 0.01). Virulence factor analysis showed an enriched iron uptake system, especially bacterial siderophores, in the pre-FMT response group, with a greater contribution from Escherichia coli. The total baseline copy number was significantly higher in the response group than non-response group (1889.14 vs. 94.86 copies/ng, P < 0.01). A total baseline copy number cutoff value of 755.88 copies/ng showed 94.7% specificity and 72.5% sensitivity in predicting FMT responsiveness. CONCLUSIONS: A significant increase in Faecalibacterium, and decrease in Enterobacteriaceae and the total fecal siderophore genes copy number were observed in responders after FMT. The siderophore genes and its encoding bacteria may be of predictive value for the clinical responsiveness of FMT to active ulcerative colitis.

Deep learning for predicting 16S rRNA gene copy number.

Culture-independent 16S rRNA gene metabarcoding is a commonly used method for microbiome profiling. To achieve more quantitative cell fraction estimates, it is important to account for the 16S rRNA gene copy number (hereafter 16S GCN) of different community members. Currently, there are several bioinformatic tools available to estimate the 16S GCN values, either based on taxonomy assignment or phylogeny. Here we present a novel approach ANNA16, Artificial Neural Network Approximator for 16S rRNA gene copy number, a deep learning-based method that estimates the 16S GCN values directly from the 16S gene sequence strings. Based on 27,579 16S rRNA gene sequences and gene copy number data from the rrnDB database, we show that ANNA16 outperforms the commonly used 16S GCN prediction algorithms. Interestingly, Shapley Additive exPlanations (SHAP) shows that ANNA16 can identify unexpected informative positions in 16S rRNA gene sequences without any prior phylogenetic knowledge, which suggests potential applications beyond 16S GCN prediction.

Chromosomal copy number and mutational status are required to authenticate ovarian cancer cell lines as appropriate cell models.

BACKGROUND: The mutational status of ovarian cancer cell line IGROV-1 is inconsistent across the literature, suggestive of multiple clonal populations of the cell line. IGROV-1 has previously been categorised as an inappropriate model for high-grade serous ovarian cancer. METHODS: IGROV-1 cells were obtained from the Netherlands Cancer Institute (IGROV-1-NKI) and the MD Anderson Cancer Centre (IGROV-1-MDA). Cell lines were STR fingerprinted and had their chromosomal copy number analysed and BRCA1/2 genes sequenced. Mutation status of ovarian cancer-related genes were extracted from the literature. RESULTS: The IGROV-1-NKI cell line has a tetraploid chromosomal profile. In contrast, the IGROV-1-MDA cell line has pseudo-normal chromosomes. The IGROV-1-NKI and IGROV-MDA are both STR matches (80.7% and 84.6%) to the original IGROV-1 cells isolated in 1985. However, IGROV-1-NKI and IGROV-1-MDA are not an STR match to each other (78.1%) indicating genetic drift. The BRCA1 and BRCA2 gene sequences are 100% identical between IGROV-1-MDA and IGROV-1-NKI, including a BRCA1 heterozygous deleterious mutation. The IGROV-1-MDA cells are more resistant to cisplatin and olaparib than IGROV-1-NKI. IGROV-1 has a mutational profile consistent with both Type I (PTEN, PIK3CA and ARID1A) and Type II ovarian cancer (BRCA1, TP53) and is likely to be a Type II high-grade serous carcinoma of the SET (Solid, pseudo-Endometroid and Transitional cell carcinoma-like morphology) subtype. CONCLUSIONS: Routine testing of chromosomal copy number as well as the mutational status of ovarian cancer related genes should become the new standard alongside STR fingerprinting to ensure that ovarian cancer cell lines are appropriate models.

A dynamic in vitro model of Down syndrome neurogenesis with trisomy 21 gene dosage correction.

Excess gene dosage from chromosome 21 (chr21) causes Down syndrome (DS), spanning developmental and acute phenotypes in terminal cell types. Which phenotypes remain amenable to intervention after development is unknown. To address this question in a model of DS neurogenesis, we derived trisomy 21 (T21) human induced pluripotent stem cells (iPSCs) alongside, otherwise, isogenic euploid controls from mosaic DS fibroblasts and equipped one chr21 copy with an inducible XIST transgene. Monoallelic chr21 silencing by XIST is near-complete and irreversible in iPSCs. Differential expression reveals that T21 neural lineages and iPSCs share suppressed translation and mitochondrial pathways and activate cellular stress responses. When XIST is induced before the neural progenitor stage, T21 dosage correction suppresses a pronounced skew toward astrogenesis in neural differentiation. Because our transgene remains inducible in postmitotic T21 neurons and astrocytes, we demonstrate that XIST efficiently represses genes even after terminal differentiation, which will empower exploration of cell type-specific T21 phenotypes that remain responsive to chr21 dosage.

Implementation of a High-Accuracy Targeted Gene Expression Panel for Clinical Care.

This study describes the validation of a clinical RNA expression panel with evaluation of concordance between gene copy gain by a next-generation sequencing (NGS) assay and high gene expression by an RNA expression panel. The RNA Salah Targeted Expression Panel (RNA STEP) was designed with input from oncologists to include 204 genes with utility for clinical trial prescreening and therapy selection. RNA STEP was validated with the nanoString platform using remnant formalin-fixed, paraffin-embedded-derived RNA from 102 patients previously tested with a validated clinical NGS panel. The repeatability, reproducibility, and concordance of RNA STEP results with NGS results were evaluated. RNA STEP demonstrated high repeatability and reproducibility, with excellent correlation (r > 0.97, P < 0.0001) for all comparisons. Comparison of RNA STEP high gene expression (log2 ratio ≥ 2) versus NGS DNA-based gene copy number gain (copies ≥ 5) for 38 mutually covered genes revealed an accuracy of 93.0% with a positive percentage agreement of 69.4% and negative percentage agreement of 93.8%. Moderate correlation was observed between platforms (r = 0.53, P < 0.0001). Concordance between high gene expression and gene copy number gain varied by specific gene, and some genes had higher accuracy between assays. Clinical implementation of RNA STEP provides gene expression data complementary to NGS and offers a tool for prescreening patients for clinical trials.