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.

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.

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.

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.

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.

Removal of the large inverted repeat from the plastid genome reveals gene dosage effects and leads to increased genome copy number.

The chloroplast genomes of most plants and algae contain a large inverted repeat (IR) region that separates two single-copy regions and harbours the ribosomal RNA operon. We have addressed the functional importance of the IR region by removing an entire copy of the 25.3-kb IR from the tobacco plastid genome. Using plastid transformation and subsequent selectable marker gene elimination, we precisely excised the IR, thus generating plants with a substantially reduced plastid genome size. We show that the lack of the IR results in a mildly reduced plastid ribosome number, suggesting a gene dosage benefit from the duplicated presence of the ribosomal RNA operon. Moreover, the IR deletion plants contain an increased number of plastid genomes, suggesting that genome copy number is regulated by measuring total plastid DNA content rather than by counting genomes. Together, our findings (1) demonstrate that the IR can enhance the translation capacity of the plastid, (2) reveal the relationship between genome size and genome copy number, and (3) provide a simplified plastid genome structure that will facilitate future synthetic biology applications.

Rapid response of fly populations to gene dosage across development and generations.

Although the effects of genetic and environmental perturbations on multicellular organisms are rarely restricted to single phenotypic layers, our current understanding of how developmental programs react to these challenges remains limited. Here, we have examined the phenotypic consequences of disturbing the bicoid regulatory network in early Drosophila embryos. We generated flies with two extra copies of bicoid, which causes a posterior shift of the network's regulatory outputs and a decrease in fitness. We subjected these flies to EMS mutagenesis, followed by experimental evolution. After only 8-15 generations, experimental populations have normalized patterns of gene expression and increased survival. Using a phenomics approach, we find that populations were normalized through rapid increases in embryo size driven by maternal changes in metabolism and ovariole development. We extend our results to additional populations of flies, demonstrating predictability. Together, our results necessitate a broader view of regulatory network evolution at the systems level.

Escherichia coli cells evade inducible parE toxin expression by reducing plasmid copy number.

Plasmids play important roles in microbial ecosystems, serving as carriers of antibiotic resistance and virulence. In the laboratory, they are essential tools for genetic manipulation and recombinant protein expression. We uncovered an intriguing survival phenotype in a fraction of the bacterial population while using plasmid-mediated arabinose-inducible gene expression to monitor the production of toxic ParE proteins. This phenotype was not correlated with changes to the plasmid sequence and could not be rescued by increasing arabinose uptake. Instead, survival correlates with a marked reduction in plasmid copy number (PCN). Reduced PCN is reproducible, not a function of the pre-existing population, and can be sequentially enriched by continual passage with induction. The reduction in PCN appears to allow mitigation of toxicity from the expression of ParE proteins while balancing the need to maintain a threshold PCN to withstand selection conditions. This indicates an adaptive cellular response to stressful conditions, likely by altering the regulation of plasmid replication. Furthermore, this survival mechanism appears to not be limited to a specific bacterial strain of Escherichia coli or ParE toxin family member, suggesting a generalized response. Finally, bacterial whole genome sequencing indicated an N845S residue substitution in DNA polymerase I, which correlates with the observed reduction in PCN and has been previously reported to impact plasmid replication. Further understanding this molecular mechanism has broader implications for this adaptive response of the dynamics of plasmid-mediated gene expression, microbial adaptation, and genetic engineering methodologies. IMPORTANCE: This research has increased our understanding of how bacteria respond to the pressure from plasmid-borne toxic genes, such as those found in toxin-antitoxin systems. Surprisingly, we found that bacteria survived toxic ParE protein expression by reducing the number of these plasmids in the cells. This discovery reveals another way in which bacteria can balance toxin expression with antibiotic selection to attenuate the effects of deleterious genes. This insight is not only valuable for understanding bacterial survival strategies but may also influence the development of better tools in biotechnology, where plasmids are often used to study the functional roles of genes.

Enhancing the expression of the unspecific peroxygenase in Komagataella phaffii through a combination strategy.

The unspecific peroxygenase (UPO) from Cyclocybe aegerita (AaeUPO) can selectively oxidize C-H bonds using hydrogen peroxide as an oxygen donor without cofactors, which has drawn significant industrial attention. Many studies have made efforts to enhance the overall activity of AaeUPO expressed in Komagataella phaffii by employing strategies such as enzyme-directed evolution, utilizing appropriate promoters, and screening secretion peptides. Building upon these previous studies, the objective of this study was to further enhance the expression of a mutant of AaeUPO with improved activity (PaDa-I) by increasing the gene copy number, co-expressing chaperones, and optimizing culture conditions. Our results demonstrated that a strain carrying approximately three copies of expression cassettes and co-expressing the protein disulfide isomerase showed an approximately 10.7-fold increase in volumetric enzyme activity, using the 2,2'-azino-bis(3-ethylbenzothiazoline-6-sulfonic acid) as the substrate. After optimizing the culture conditions, the volumetric enzyme activity of this strain further increased by approximately 48.7%, reaching 117.3 U/mL. Additionally, the purified catalytic domain of PaDa-I displayed regioselective hydroxylation of R-2-phenoxypropionic acid. The results of this study may facilitate the industrial application of UPOs. KEY POINTS: • The secretion of the catalytic domain of PaDa-I can be significantly enhanced through increasing gene copy numbers and co-expressing of protein disulfide isomerase. • After optimizing the culture conditions, the volumetric enzyme activity can reach 117.3 U/mL, using the 2,2'-azino-bis(3-ethylbenzothiazoline-6-sulfonic acid) as the substrate. • The R-2-phenoxypropionic acid can undergo the specific hydroxylation reaction catalyzed by catalytic domain of PaDa-I, resulting in the formation of R-2-(4-hydroxyphenoxy)propionic acid.

Tandem gene duplication selected by activation of horizontally transferred gene in bacteria.

Horizontal gene transfer occurs frequently in bacteria, but the mechanism driving activation and optimization of the expression of horizontally transferred genes (HTGs) in new recipient strains is not clear. Our previous study found that spontaneous tandem DNA duplication resulted in rapid activation of HTGs. Here, we took advantage of this finding to develop a novel technique for tandem gene duplication, named tandem gene duplication selected by activation of horizontally transferred gene in bacteria (TDAH), in which tandem duplication was selected by the activation of horizontally transferred selectable marker gene. TDAH construction does not contain any reported functional elements based on homologous or site-specific recombination and DNA amplification. TDAH only contains an essential selectable marker for copy number selection and 9-bp-microhomology border sequences for precise illegitimate recombination. One transformation and 3 days were enough to produce a high-copy strain, so its procedure is simple and fast. Without subsequent knockout of the endogenous recombination system, TDAH could also generate the relatively stable high-copy tandem duplication for plasmid-carried and genome-integrated DNA. TDAH also showed an excellent capacity for increase gene expression and worked well in different industrial bacteria. We also applied TDAH to select the optimal high copy number of ribA for vitamin B2 production in E. coli; the yield was improved by 3.5 times and remained stable even after 12 subcultures. TDAH is a useful tool for recombinant protein production and expression optimization of biosynthetic pathways. KEY POINTS: • We develop a novel and efficient technique (TDAH) for tandem gene duplication in bacterium. TDAH is based on the mechanism of HTG rapid activation. TDAH does not contain any reported functional elements based on homologous recombination and DNA amplification. TDAH only contains an essential selectable marker for copy number selection, so its construction and procedure are very simple and fast. • TDAH is the first reported selected and stable tandem-gene-duplication technique in which the selected high-copy plasmid-carried and genome-integrated DNA could remain stable without the subsequent knockout of recombination system. • TDAH showed an excellent capacity for regulating gene expression and worked well in different industrial bacteria, indicating it is a useful tool for recombinant protein production and expression optimization of biosynthetic pathways. • TDAH was applied to select the optimal high copy number of ribA for vitamin B2 production in E. coli; the yield was improved by 3.5-fold and remained stable even after 12 subcultures.

Porin deficiency or plasmid copy number increase mediated carbapenem-resistant Escherichia coli resistance evolution.

This study investigated resistance evolution mechanisms of conjugated plasmids and bacterial hosts under different concentrations of antibiotic pressure. Ancestral strain ECNX52 was constructed by introducing the blaNDM-5-carrying IncX3 plasmid into E. coli C600, and was subjected to laboratory evolution under different concentrations of meropenem pressure. Minimal inhibitory concentrations and conjugation frequency were determined. Fitness of these strains was assessed. Whole genome sequencing and transcriptional changes were performed. Ancestral host or plasmids were recombined with evolved hosts or plasmids to verify plasmid or host factors in resistance evolution. Role of the repA mutation on plasmid copy number was determined. Two out of the four clones (EM2N1 and EM2N3) exhibited four-fold increase in MIC when exposed to a continuous pressure of 2 µg/mL MEM (1/32 MIC), by down regulating expression of outer membrane protein ompF. Besides, all four clones displayed four-fold increase in MIC and higher conjugation frequency when subjected to a continuous pressure of 4 µg/mL MEM (1/16 MIC), attributing to increasing plasmid copy number generated by repA D140Y (GAT→TAT) mutation. Bacterial hosts and conjugative plasmids can undergo resistance evolution under certain concentrations of antimicrobial pressure by reducing the expression of outer membrane proteins or increasing plasmid copy numbers.

Three concurrent mechanisms generate gene copy number variation and transient antibiotic heteroresistance.

Heteroresistance is a medically relevant phenotype where small antibiotic-resistant subpopulations coexist within predominantly susceptible bacterial populations. Heteroresistance reduces treatment efficacy across diverse bacterial species and antibiotic classes, yet its genetic and physiological mechanisms remain poorly understood. Here, we investigated a multi-resistant Klebsiella pneumoniae isolate and identified three primary drivers of gene dosage-dependent heteroresistance for several antibiotic classes: tandem amplification, increased plasmid copy number, and transposition of resistance genes onto cryptic plasmids. All three mechanisms imposed fitness costs and were genetically unstable, leading to fast reversion to susceptibility in the absence of antibiotics. We used a mouse gut colonization model to show that heteroresistance due to elevated resistance-gene dosage can result in antibiotic treatment failures. Importantly, we observed that the three mechanisms are prevalent among Escherichia coli bloodstream isolates. Our findings underscore the necessity for treatment strategies that address the complex interplay between plasmids, resistance cassettes, and transposons in bacterial populations.

Differences in 5'untranslated regions highlight the importance of translational regulation of dosage sensitive genes.

BACKGROUND: Untranslated regions (UTRs) are important mediators of post-transcriptional regulation. The length of UTRs and the composition of regulatory elements within them are known to vary substantially across genes, but little is known about the reasons for this variation in humans. Here, we set out to determine whether this variation, specifically in 5'UTRs, correlates with gene dosage sensitivity. RESULTS: We investigate 5'UTR length, the number of alternative transcription start sites, the potential for alternative splicing, the number and type of upstream open reading frames (uORFs) and the propensity of 5'UTRs to form secondary structures. We explore how these elements vary by gene tolerance to loss-of-function (LoF; using the LOEUF metric), and in genes where changes in dosage are known to cause disease. We show that LOEUF correlates with 5'UTR length and complexity. Genes that are most intolerant to LoF have longer 5'UTRs, greater TSS diversity, and more upstream regulatory elements than their LoF tolerant counterparts. We show that these differences are evident in disease gene-sets, but not in recessive developmental disorder genes where LoF of a single allele is tolerated. CONCLUSIONS: Our results confirm the importance of post-transcriptional regulation through 5'UTRs in tight regulation of mRNA and protein levels, particularly for genes where changes in dosage are deleterious and lead to disease. Finally, to support gene-based investigation we release a web-based browser tool, VuTR, that supports exploration of the composition of individual 5'UTRs and the impact of genetic variation within them.

Gene Dosage of F5 c.3481C>T Stop-Codon (p.R1161Ter) Switches the Clinical Phenotype from Severe Thrombosis to Recurrent Haemorrhage: Novel Hypotheses for Readthrough Strategy.

Inherited defects in the genes of blood coagulation essentially express the severity of the clinical phenotype that is directly correlated to the number of mutated alleles of the candidate leader gene (e.g., heterozygote vs. homozygote) and of possible additional coinherited traits. The F5 gene, which codes for coagulation factor V (FV), plays a two-faced role in the coagulation cascade, exhibiting both procoagulant and anticoagulant functions. Thus, defects in this gene can be predisposed to either bleeding or thrombosis. A Sanger sequence analysis detected a premature stop-codon in exon 13 of the F5 gene (c.3481C>T; p.R1161Ter) in several members of a family characterised by low circulating FV levels and contrasting clinical phenotypes. The propositus, a 29 y.o. male affected by recurrent haemorrhages, was homozygous for the F5 stop-codon and for the F5 c.1691G>A (p.R506Q; FV-Leiden) inherited from the heterozygous parents, which is suggestive of combined cis-segregation. The homozygous condition of the stop-codon completely abolished the F5 gene expression in the propositus (FV:Ag < 1%; FV:C < 1%; assessed by ELISA and PT-based one-stage clotting assay respectively), removing, in turn, any chance for FV-Leiden to act as a prothrombotic molecule. His father (57 y.o.), characterised by severe recurrent venous thromboses, underwent a complete molecular thrombophilic screening, revealing a heterozygous F2 G20210A defect, while his mother (56 y.o.), who was negative for further common coagulation defects, reported fully asymptomatic anamnesis. To dissect these conflicting phenotypes, we performed the ProC®Global (Siemens Helthineers) coagulation test aimed at assessing the global pro- and anticoagulant balance of each family member, investigating the responses to the activated protein C (APC) by means of an APC-sensitivity ratio (APC-sr). The propositus had an unexpectedly poor response to APC (APC-sr: 1.09; n.v. > 2.25), and his father and mother had an APC-sr of 1.5 and 2.0, respectively. Although ProC®Global prevalently detects the anticoagulant side of FV, the exceptionally low APC-sr of the propositus and his discordant severe-moderate haemorrhagic phenotype could suggest a residual expression of mutated FV p.506QQ through a natural readthrough or possible alternative splicing mechanisms. The coagulation pathway may be physiologically rebalanced through natural and induced strategies, and the described insights might be able to track the design of novel treatment approaches and rebalancing molecules.

Duplicate Gene Expression and Possible Mechanisms of Paralog Retention During Bacterial Genome Expansion.

Gene duplication contributes to the evolution of expression and the origin of new genes, but the relative importance of different patterns of duplicate gene expression and mechanisms of retention remains debated and particularly poorly understood in bacteria. Here, we investigated gene expression patterns for two lab strains of the cyanobacterium Acaryochloris marina with expanding genomes that contain about 10-fold more gene duplicates compared with most bacteria. Strikingly, we observed a generally stoichiometric pattern of greater combined duplicate transcript dosage with increased gene copy number, in contrast to the prevalence of expression reduction reported for many eukaryotes. We conclude that increased transcript dosage is likely an important mechanism of initial duplicate retention in these bacteria and may persist over long periods of evolutionary time. However, we also observed that paralog expression can diverge rapidly, including possible functional partitioning, for which different copies were respectively more highly expressed in at least one condition. Divergence may be promoted by the physical separation of most Acaryochloris duplicates on different genetic elements. In addition, expression pattern for ancestrally shared duplicates could differ between strains, emphasizing that duplicate expression fate need not be deterministic. We further observed evidence for context-dependent transcript dosage, where the aggregate expression of duplicates was either greater or lower than their single-copy homolog depending on physiological state. Finally, we illustrate how these different expression patterns of duplicated genes impact Acaryochloris biology for the innovation of a novel light-harvesting apparatus and for the regulation of recA paralogs in response to environmental change.

DEFA1A3 DNA gene-dosage regulates the kidney innate immune response during upper urinary tract infection.

Antimicrobial peptides (AMPs) are host defense effectors with potent neutralizing and immunomodulatory functions against invasive pathogens. The AMPs α-Defensin 1-3/DEFA1A3 participate in innate immune responses and influence patient outcomes in various diseases. DNA copy-number variations in DEFA1A3 have been associated with severity and outcomes in infectious diseases including urinary tract infections (UTIs). Specifically, children with lower DNA copy numbers were more susceptible to UTIs. The mechanism of action by which α-Defensin 1-3/DEFA1A3 copy-number variations lead to UTI susceptibility remains to be explored. In this study, we use a previously characterized transgenic knock-in of the human DEFA1A3 gene mouse to dissect α-Defensin 1-3 gene dose-dependent antimicrobial and immunomodulatory roles during uropathogenic Escherichia coli (UPEC) UTI. We elucidate the relationship between kidney neutrophil- and collecting duct intercalated cell-derived α-Defensin 1-3/DEFA1A3 expression and UTI. We further describe cooperative effects between α-Defensin 1-3 and other AMPs that potentiate the neutralizing activity against UPEC. Cumulatively, we demonstrate that DEFA1A3 directly protects against UPEC meanwhile impacting pro-inflammatory innate immune responses in a gene dosage-dependent manner.

Using digital PCR to predict ciliate abundance from ribosomal RNA gene copy numbers.

Ciliates play a key role in most ecosystems. Their abundance in natural samples is crucial for answering many ecological questions. Traditional methods of quantifying individual species, which rely on microscopy, are often labour-intensive, time-consuming and can be highly biassed. As a result, we investigated the potential of digital polymerase chain reaction (dPCR) for quantifying ciliates. A significant challenge in this process is the high variation in the copy number of the taxonomic marker gene (ribosomal RNA [rRNA]). We first quantified the rRNA gene copy numbers (GCN) of the model ciliate, Paramecium tetraurelia, during different stages of the cell cycle and growth phases. The per-cell rRNA GCN varied between approximately 11,000 and 130,000, averaging around 50,000 copies per cell. Despite these variations in per-cell rRNA GCN, we found a highly significant correlation between GCN and cell numbers. This is likely due to the coexistence of different cellular stages in an uncontrolled (environmental) ciliate population. Thanks to the high sensitivity of dPCR, we were able to detect the target gene in a sample that contained only a single cell. The dPCR approach presented here is a valuable addition to the molecular toolbox in protistan ecology. It may guide future studies in quantifying and monitoring the abundance of targeted (even rare) ciliates in natural samples.