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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Analysis of the relationship between MYCN gene and serum NSE and urinary VMA levels and neuroblastoma pathological features and prognosis.

The purpose of this study was to explore the relationship between the MYCN gene, serum neuron-specific enolase (NSE), urinary vanillylmandelic acid (VMA) levels, and neuroblastoma pathological features and prognosis. Ninety-four children with neuroblastoma treated in the hospital were selected to compare the differences in MYCN gene amplification, serum NSE, and urinary VMA levels in children with different clinicopathological features and prognoses. The proportion of children with MYCN gene copy number ≥10 in INSS stage 3-4 was higher than that of children with INSS stage 1-2 (P < 0.05); the proportion of children with MYCN gene copy number ≥10 in high-risk children in the COG risk stratification was higher than that of children with intermediate and low risk (P < 0.05); the serum NSE of children aged >12 months higher than that of children aged ≤12 months (P < 0.05); serum NSE of children with tumors >500 cm3 higher than that of children with tumors ≤500 cm3 (P < 0.05); serum NSE and urinary VMA of children with INSS staging of stages 3-4 were higher than that of children with stages 1 to 2 (P < 0.05); serum NSE and urinary VMA in children with lymph node metastasis were higher than that of children without lymph node metastasis (P < 0.05); serum NSE of children with MYCN gene copy number ≥10 was higher than that of children without lymph node metastasis (P < 0.05); the proportion of children with MYCN gene copy number ≥10 who died, and the percentages of serum NSE and urinary VMA were higher than those of the surviving children (P < 0.05). MYCN gene amplification and serum NSE and urinary VMA levels were related to the age, tumor size, INSS stage, COG stage, lymph node metastasis, and prognosis of the children with neuroblastoma.

CDX2 dose-dependently influences the gene regulatory network underlying human extraembryonic mesoderm development.

Loss of Cdx2 in vivo leads to stunted development of the allantois, an extraembryonic mesoderm-derived structure critical for nutrient delivery and waste removal in the early embryo. Here, we investigate how CDX2 dose-dependently influences the gene regulatory network underlying extraembryonic mesoderm development. By engineering human induced pluripotent stem cells (hiPSCs) consisting of wild-type (WT), heterozygous (CDX2-Het), and homozygous null CDX2 (CDX2-KO) genotypes, differentiating these cells in a 2D gastruloid model, and subjecting these cells to single-nucleus RNA and ATAC sequencing, we identify several pathways that are dose-dependently regulated by CDX2 including VEGF and non-canonical WNT. snATAC-seq reveals that CDX2-Het cells retain a WT-like chromatin accessibility profile, suggesting accessibility alone is not sufficient to drive this variability in gene expression. Because the loss of CDX2 or TBXT phenocopy one another in vivo, we compared differentially expressed genes in our CDX2-KO to those from TBXT-KO hiPSCs differentiated in an analogous experiment. This comparison identifies several communally misregulated genes that are critical for cytoskeletal integrity and tissue permeability. Together, these results clarify how CDX2 dose-dependently regulates gene expression in the extraembryonic mesoderm and reveal pathways that may underlie the defects in vascular development and allantoic elongation seen in vivo.

Integrating high-throughput analysis to create an atlas of replication origins in Trypanosoma cruzi in the context of genome structure and variability.

Trypanosoma cruzi is the etiologic agent of the most prevalent human parasitic disease in Latin America, Chagas disease. Its genome is rich in multigenic families that code for virulent antigens and are present in the rapidly evolving genomic compartment named Disruptive. DNA replication is a meticulous biological process in which flaws can generate mutations and changes in chromosomal and gene copy numbers. Here, integrating high-throughput and single-molecule analyses, we were able to identify Predominant, Flexible, and Dormant Orc1Cdc6-dependent origins as well as Orc1Cdc6-independent origins. Orc1Cdc6-dependent origins were found in multigenic family loci, while independent origins were found in the Core compartment that contains conserved and hypothetical protein-coding genes, in addition to multigenic families. In addition, we found that Orc1Cdc6 density is related to the firing of origins and that Orc1Cdc6-binding sites within fired origins are depleted of a specific class of nucleosomes that we previously categorized as dynamic. Together, these data suggest that Orc1Cdc6-dependent origins may contribute to the rapid evolution of the Disruptive compartment and, therefore, to the success of T. cruzi infection and that the local epigenome landscape is also involved in this process.IMPORTANCETrypanosoma cruzi, responsible for Chagas disease, affects millions globally, particularly in Latin America. Lack of vaccine or treatment underscores the need for research. Parasite's genome, with virulent antigen-coding multigenic families, resides in the rapidly evolving Disruptive compartment. Study sheds light on the parasite's dynamic DNA replication, discussing the evolution of the Disruptive compartment. Therefore, the findings represent a significant stride in comprehending T. cruzi's biology and the molecular bases that contribute to the success of infection caused by this parasite.

Gene dosage of independent dynein arm motor preassembly factors influences cilia assembly in Chlamydomonas reinhardtii.

Motile cilia assembly utilizes over 800 structural and cytoplasmic proteins. Variants in approximately 58 genes cause primary ciliary dyskinesia (PCD) in humans, including the dynein arm (pre)assembly factor (DNAAF) gene DNAAF4. In humans, outer dynein arms (ODAs) and inner dynein arms (IDAs) fail to assemble motile cilia when DNAAF4 function is disrupted. In Chlamydomonas reinhardtii, a ciliated unicellular alga, the DNAAF4 ortholog is called PF23. The pf23-1 mutant assembles short cilia and lacks IDAs, but partially retains ODAs. The cilia of a new null allele (pf23-4) completely lack ODAs and IDAs and are even shorter than cilia from pf23-1. In addition, PF23 plays a role in the cytoplasmic modification of IC138, a protein of the two-headed IDA (I1/f). As most PCD variants in humans are recessive, we sought to test if heterozygosity at two genes affects ciliary function using a second-site non-complementation (SSNC) screening approach. We asked if phenotypes were observed in diploids with pairwise heterozygous combinations of 21 well-characterized ciliary mutant Chlamydomonas strains. Vegetative cultures of single and double heterozygous diploid cells did not show SSNC for motility phenotypes. When protein synthesis is inhibited, wild-type Chlamydomonas cells utilize the pool of cytoplasmic proteins to assemble half-length cilia. In this sensitized assay, 8 double heterozygous diploids with pf23 and other DNAAF mutations show SSNC; they assemble shorter cilia than wild-type. In contrast, double heterozygosity of the other 203 strains showed no effect on ciliary assembly. Immunoblots of diploids heterozygous for pf23 and wdr92 or oda8 show that PF23 is reduced by half in these strains, and that PF23 dosage affects phenotype severity. Reductions in PF23 and another DNAAF in diploids affect the ability to assemble ODAs and IDAs and impedes ciliary assembly. Thus, dosage of multiple DNAAFs is an important factor in cilia assembly and regeneration.

Plasmid-mediated phenotypic noise leads to transient antibiotic resistance in bacteria.

The rise of antibiotic resistance is a critical public health concern, requiring an understanding of mechanisms that enable bacteria to tolerate antimicrobial agents. Bacteria use diverse strategies, including the amplification of drug-resistance genes. In this paper, we showed that multicopy plasmids, often carrying antibiotic resistance genes in clinical bacteria, can rapidly amplify genes, leading to plasmid-mediated phenotypic noise and transient antibiotic resistance. By combining stochastic simulations of a computational model with high-throughput single-cell measurements of blaTEM-1 expression in Escherichia coli MG1655, we showed that plasmid copy number variability stably maintains populations composed of cells with both low and high plasmid copy numbers. This diversity in plasmid copy number enhances the probability of bacterial survival in the presence of antibiotics, while also rapidly reducing the burden of carrying multiple plasmids in drug-free environments. Our results further support the tenet that multicopy plasmids not only act as vehicles for the horizontal transfer of genetic information between cells but also as drivers of bacterial adaptation, enabling rapid modulation of gene copy numbers. Understanding the role of multicopy plasmids in antibiotic resistance is critical, and our study provides insights into how bacteria can transiently survive lethal concentrations of antibiotics.

Fusobacterium nucleatum subsp. polymorphum recovered from malignant and potentially malignant oral disease exhibit heterogeneity in adhesion phenotypes and adhesin gene copy number, shaped by inter-subspecies horizontal gene transfer and recombination-derived mosaicism.

Fusobacterium nucleatum is an anaerobic commensal of the oral cavity associated with periodontitis and extra-oral diseases, including colorectal cancer. Previous studies have shown an increased relative abundance of this bacterium associated with oral dysplasia or within oral tumours. Using direct culture, we found that 75 % of Fusobacterium species isolated from malignant or potentially malignant oral mucosa were F. nucleatum subsp. polymorphum. Whole genome sequencing and pangenome analysis with Panaroo was carried out on 76 F. nucleatum subsp. polymorphum genomes. F. nucleatum subsp. polymorphum was shown to possesses a relatively small core genome of 1604 genes in a pangenome of 7363 genes. Phylogenetic analysis based on the core genome shows the isolates can be separated into three main clades with no obvious genotypic associations with disease. Isolates recovered from healthy and diseased sites in the same patient are generally highly related. A large repertoire of adhesins belonging to the type V secretion system (TVSS) could be identified with major variation in repertoire and copy number between strains. Analysis of intergenic recombination using fastGEAR showed that adhesin complement is shaped by horizontal gene transfer and recombination. Recombination events at TVSS adhesin genes were not only common between lineages of subspecies polymorphum, but also between different subspecies of F. nucleatum. Strains of subspecies polymorphum with low copy numbers of TVSS adhesin encoding genes tended to have the weakest adhesion to oral keratinocytes. This study highlights the genetic heterogeneity of F. nucleatum subsp. polymorphum and provides a new framework for defining virulence in this organism.

The copy number of the eukaryotic rRNA gene can be counted comprehensively.

Gene sequence has been widely used in molecular ecology. For instance, the ribosomal RNA (rRNA) gene has been widely used as a biological marker to understand microbial communities. The variety of the detected rRNA gene sequences reflects the diversity of the microorganisms existing in the analyzed sample. Their biomass can also be estimated by applying quantitative sequencing with information on rRNA gene copy numbers in genomes; however, information on rRNA gene copy numbers is still limited. Especially, the copy number in microbial eukaryotes is much less understood than that of prokaryotes, possibly because of the large and complex structure of eukaryotic genomes. In this study, we report an alternative approach that is more appropriate than the existing method of quantitative sequencing and demonstrate that the copy number of eukaryotic rRNA can be measured efficiently and comprehensively. By applying this approach widely, information on the eukaryotic rRNA copy number can be determined, and their community structures can be depicted and compared more efficiently.

Droplet digital polymerase chain reaction-based quantitation of therapeutic lentiviral vector copies in transduced hematopoietic stem cells.

BACKGROUND AIMS: Gene therapy using lentiviral vectors (LVs) that harbor a functional ß-globin gene provides a curative treatment for hemoglobinopathies including beta-thalassemia and sickle cell disease. Accurate quantification of the vector copy number (VCN) and/or the proportion of transduced cells is critical to evaluate the efficacy of transduction and stability of the transgene during treatment. Moreover, commonly used techniques for LV quantification, including real-time quantitative polymerase chain reaction (PCR) or fluorescence-activated cell sorting, require either a standard curve or expression of a reporter protein for the detection of transduced cells. In the present study, we describe a digital droplet PCR (ddPCR) technique to measure the lentiviral VCN in transduced hematopoietic stem and progenitor cells (HSPCs). METHODS: After HSPCs were transduced with an LV encoding the therapeutic ß-globin (ßA-T87Q) gene, the integrated lentiviral sequence in the host genome was amplified with primers that targeted a sequence within the vector and the human RPP30 gene. The dynamic range of ddPCR was between 5 × 10-3 ng and 5 × 10-6 ng of target copy per reaction. RESULTS: We found that the ddPCR-based approach was able to estimate VCN with high sensitivity and a low standard deviation. Furthermore, ddPCR-mediated quantitation of lentiviral copy numbers in differentiated erythroblasts correlated with the level of ßA-T87Q protein detected by reverse-phase high-performance liquid chromatography. CONCLUSIONS: Taken together, the ddPCR technique has the potential to precisely detect LV copy numbers in the host genome, which can be used for VCN estimation, calculation of infectious titer and multiplicity of infection for HSPC transduction in a clinical setting.

Effects of gene dosage and development on subcortical nuclei volumes in individuals with 22q11.2 copy number variations.

The 22q11.2 locus contains genes critical for brain development. Reciprocal Copy Number Variations (CNVs) at this locus impact risk for neurodevelopmental and psychiatric disorders. Both 22q11.2 deletions (22qDel) and duplications (22qDup) are associated with autism, but 22qDel uniquely elevates schizophrenia risk. Understanding brain phenotypes associated with these highly penetrant CNVs can provide insights into genetic pathways underlying neuropsychiatric disorders. Human neuroimaging and animal models indicate subcortical brain alterations in 22qDel, yet little is known about developmental differences across specific nuclei between reciprocal 22q11.2 CNV carriers and typically developing (TD) controls. We conducted a longitudinal MRI study in a total of 385 scans from 22qDel (n = 96, scans = 191, 53.1% female), 22qDup (n = 37, scans = 64, 45.9% female), and TD controls (n = 80, scans = 130, 51.2% female), across a wide age range (5.5-49.5 years). Volumes of the thalamus, hippocampus, amygdala, and anatomical subregions were estimated using FreeSurfer, and the linear effects of 22q11.2 gene dosage and non-linear effects of age were characterized with generalized additive mixed models (GAMMs). Positive gene dosage effects (volume increasing with copy number) were observed for total intracranial and whole hippocampus volumes, but not whole thalamus or amygdala volumes. Several amygdala subregions exhibited similar positive effects, with bi-directional effects found across thalamic nuclei. Distinct age-related trajectories were observed across the three groups. Notably, both 22qDel and 22qDup carriers exhibited flattened development of hippocampal CA2/3 subfields relative to TD controls. This study provides novel insights into the impact of 22q11.2 CNVs on subcortical brain structures and their developmental trajectories.

Higher evolutionary dynamics of gene copy number for Drosophila glue genes located near short repeat sequences.

BACKGROUND: During evolution, genes can experience duplications, losses, inversions and gene conversions. Why certain genes are more dynamic than others is poorly understood. Here we examine how several Sgs genes encoding glue proteins, which make up a bioadhesive that sticks the animal during metamorphosis, have evolved in Drosophila species. RESULTS: We examined high-quality genome assemblies of 24 Drosophila species to study the evolutionary dynamics of four glue genes that are present in D. melanogaster and are part of the same gene family - Sgs1, Sgs3, Sgs7 and Sgs8 - across approximately 30 millions of years. We annotated a total of 102 Sgs genes and grouped them into 4 subfamilies. We present here a new nomenclature for these Sgs genes based on protein sequence conservation, genomic location and presence/absence of internal repeats. Two types of glue genes were uncovered. The first category (Sgs1, Sgs3x, Sgs3e) showed a few gene losses but no duplication, no local inversion and no gene conversion. The second group (Sgs3b, Sgs7, Sgs8) exhibited multiple events of gene losses, gene duplications, local inversions and gene conversions. Our data suggest that the presence of short "new glue" genes near the genes of the latter group may have accelerated their dynamics. CONCLUSIONS: Our comparative analysis suggests that the evolutionary dynamics of glue genes is influenced by genomic context. Our molecular, phylogenetic and comparative analysis of the four glue genes Sgs1, Sgs3, Sgs7 and Sgs8 provides the foundation for investigating the role of the various glue genes during Drosophila life.