Search

Multiparametric magnetic resonance imaging (mpMRI) has transformed the management of localized prostate cancer by improving identification of clinically significant disease at diagnosis. Approximately 20% of primary prostate tumors are invisible to mpMRI, and we hypothesize that this invisibility reflects fundamental molecular properties of the tumor. We therefore profiled the genomes and transcriptomes of 40 International Society of Urological Pathology grade 2 tumors: 20 mpMRI-invisible (Prostate Imaging-Reporting and Data System [PI-RADS] v2

Large structural variants (SVs) in the human genome are difficult to detect and study by conventional sequencing technologies. With long-range genome analysis platforms, such as optical mapping, one can identify large SVs (>2 kb) across the genome in one experiment. Analyzing optical genome maps of 154 individuals from the 26 populations sequenced in the 1000 Genomes Project, we find that phylogenetic population patterns of large SVs are similar to those of single nucleotide variations in 86% of the human genome, while ~2% of the genome has high structural complexity. We are able to characterize SVs in many intractable regions of the genome, including segmental duplications and subtelomeric, pericentromeric, and acrocentric areas. In addition, we discover ~60 Mb of non-redundant genome content missing in the reference genome sequence assembly. Our results highlight the need for a comprehensive set of alternate haplotypes from different populations to represent SV patterns in the genome.

Mutations in the PARK2 gene are associated with early onset Parkinsonism. The Park2 -/- mouse, however, does not exhibit neurodegeneration or other Parkinsons disease (PD) phenotypes. Previously, we discovered that translation of Mcl-1, a pro-survival factor, is upregulated in the Park2 -/- mouse, suggesting a compensatory mechanism during development. Here we generated the Park2 -/- Mcl-1 +/- mouse and show that by reducing Mcl-1 gene dosage by 50%, the Park2 -/- genotype is sensitized, conferring both dopaminergic neuron loss and motor impairments. We propose that this murine model could be a useful tool for dissecting PD etiology and developing treatment strategies against this neurodegenerative disease.

Natural killer (NK) cells have increasingly become a target of interest for immunotherapies. NK cells express killer immunoglobulin-like receptors (KIRs), which play a vital role in immune response to tumors by detecting cellular abnormalities. The genomic region encoding the 16 KIR genes displays high polymorphic variability in human populations, making it difficult to resolve individual genotypes based on next generation sequencing data. As a result, the impact of polymorphic KIR variation on cancer phenotypes has been understudied. Currently, labor-intensive, experimental techniques are used to determine an individual's KIR gene copy number profile. Here, we develop an algorithm to determine the germline copy number of KIR genes from whole exome sequencing data and apply it to a cohort of nearly 5000 cancer patients. We use a k-mer based approach to capture sequences unique to specific genes, count their occurrences in the set of reads derived from an individual and compare the individual's k-mer distribution to that of the population. Copy number results demonstrate high concordance with population copy number expectations. Our method reveals that the burden of inhibitory KIR genes is associated with survival in two tumor types, highlighting the potential importance of KIR variation in understanding tumor development and response to immunotherapy.

Methyl-CpG-binding-Protein 2 (MeCP2) is an abundant nuclear protein highly enriched in neurons. Here we report live-cell single-molecule imaging studies of the kinetic features of mouse MeCP2 at high spatial-temporal resolution. MeCP2 displays dynamic features that are distinct from both highly mobile transcription factors and immobile histones. Stable binding of MeCP2 in living neurons requires its methyl-binding domain and is sensitive to DNA modification levels. Diffusion of unbound MeCP2 is strongly constrained by weak, transient interactions mediated primarily by its AT-hook domains, and varies with the level of chromatin compaction and cell type. These findings extend previous studies of the role of the MeCP2 MBD in high affinity DNA binding to living neurons, and identify a new role for its AT-hooks domains as critical determinants of its kinetic behavior. They suggest that limited nuclear diffusion of MeCP2 in live neurons contributes to its local impact on chromatin structure and gene expression.

Steroidogenic factor-1 (SF-1, also termed Ad4BP; NR5A1 in the official nomenclature) is a nuclear receptor transcription factor that plays a crucial role in the regulation of adrenal and gonadal development, function and maintenance. In addition to its classical role in regulating the expression of P450 steroid hydroxylases and other steroidogenic genes, involvement in other key processes such as cell survival/proliferation and cytoskeleton dynamics have also been highlighted for SF-1. SF-1 has a restricted pattern of expression, being expressed along the hypothalamic-pituitary axis and in steroidogenic organs since the time of their establishment. Reduced SF-1 expression affects proper gonadal and adrenal organogenesis and function. On the other hand, SF-1 overexpression is found in adrenocortical carcinoma and represents a prognostic marker for patients' survival. This review is focused on the current knowledge about SF-1 and the crucial importance of its dosage for adrenal gland development and function, from its involvement in adrenal cortex formation to tumorigenesis. Overall, data converge towards SF-1 being a key player in the complex network of transcriptional regulation within the adrenal gland in a dosage-dependent manner. [ABSTRACT FROM AUTHOR]

The safe production of β-alanine (BA) has attracted significant attention by its multifaceted applications in pharmaceutical, polymer, and nutrition. The extant high-yielding chemical and enzymatic methods of BA synthesis are handicapped by raw materials derived from petroleum resources, harsh reaction conditions, and catalyst instability. Consequently, this study explored a safe and alternative route via microbial fermentation, utilizing metabolic engineering of Bacillus megaterium to produce BA. The Bacillus subtilis panD gene (encoding L-aspartate-α-decarboxylase) was codon-optimized and overexpressed, which yielded 0.13 ± 0.05 g/L BA. Aspartate ammonia-lyase (AspA) and aspartate aminotransferase (AspB) based pathways were examined for BA production from glucose. NADH-dependent glutamate dehydrogenase (gdh) was used to regenerate the cofactor NAD+ in the pathway with AspB. Dosing of the rate liming panD showed a positive effect on BA production. The BA titer was further increased to 1.4 ± 0.06 g/L by over-expression of phosphoenolpyruvate carboxylase (PPC). Optimizing (NH4)2SO4 Pyridoxine, and NaHCO3 allowed the production of 2.41 ± 0.15 g/L BA. Fed-batch fermentation of the final strain allowed 17.60 ± 0.13 g/L BA production in 22 h. The present study has effectively unlocked the potential of engineering the B. megaterium for the sustainable production of the other ASP (L-aspartic acid) and BA-derived products at a large scale. [ABSTRACT FROM AUTHOR]

An understanding of how heterozygous loss-of-function mutations in autism spectrum disorder (ASD) risk genes, such as TBR1, contribute to ASD remains elusive. Conditional Tbr1 deletion during late mouse gestation in cortical layer 6 neurons (Tbr1layer6 mutants) provides novel insights into its function, including dendritic patterning, synaptogenesis, and cell-intrinsic physiology. These phenotypes occur in heterozygotes, providing insights into mechanisms that may underlie ASD pathophysiology. Restoring expression of Wnt7b largely rescues the synaptic deficit in Tbr1layer6 mutant neurons. Furthermore, Tbr1layer6 heterozygotes have increased anxiety-like behavior, a phenotype seen ASD. Integrating TBR1 chromatin immunoprecipitation sequencing (ChIP-seq) and RNA sequencing (RNA-seq) data from layer 6 neurons and activity of TBR1-bound candidate enhancers provides evidence for how TBR1 regulates layer 6 properties. Moreover, several putative TBR1 targets are ASD risk genes, placing TBR1 in a central position both for ASD risk and for regulating transcriptional circuits that control multiple steps in layer 6 development essential for the assembly of neural circuits.

Mitochondrial-derived peptides (MDPs) are rapidly emerging therapeutic targets to combat development of neurodegenerative diseases. SHLP2 (small humanin-like peptide 2) is a newly discovered MDP that is coded from the MT-RNR2 (Mitochondrially encoded 16S rRNA) gene in mitochondrial DNA (mtDNA). In the current study, we examined the biological consequences of treatment with exogenously-added SHLP2 in an in vitro human transmitochondrial age-related macular degeneration (AMD) ARPE-19 cell model. In AMD cells, we observed significant down-regulation of the MDP-coding MT-RNR2 gene, and remarkably reduced levels of all five oxidative phosphorylation (OXPHOS) complex I-V protein subunits that are involved in the electron transport chain; these results suggested mitochondrial toxicity and abnormal OXPHOS complex protein subunits' levels in AMD cells. However, treatment of AMD cells with SHLP2: (1) restored the normal levels of OXPHOS complex protein subunits, (2) prevented loss of viable cells and mitochondria, (3) increased the number of mtDNA copies, (4) induced anti-apoptotic effects, and (5) attenuated amyloid-β-induced cellular and mitochondrial toxicity. Cumulatively, our findings established the protective role of SHLP2 in AMD cells in vitro. In conclusion, this novel study supports the merit of SHLP2 in the treatment of AMD, a primary retinal disease that is a leading cause of blindness among the elderly population in the United States as well as worldwide.

The emergence of robust single-cell 'omics techniques enables studies of uncultivable species, allowing for the (re)discovery of diverse genomic features. In this study, we combine single-cell genomics and transcriptomics to explore genome evolution in ciliates (a > 1 Gy old clade). Analysis of the data resulting from these single-cell 'omics approaches show: 1) the description of the ciliates in the class Karyorelictea as "primitive" is inaccurate because their somatic macronuclei contain loci of varying copy number (i.e., they have been processed by genome rearrangements from the zygotic nucleus); 2) gene-sized somatic chromosomes exist in the class Litostomatea, consistent with Balbiani's (1890) observation of giant chromosomes in this lineage; and 3) gene scrambling exists in the underexplored Postciliodesmatophora (the classes Heterotrichea and Karyorelictea, abbreviated here as the Po-clade), one of two major clades of ciliates. Together these data highlight the complex evolutionary patterns underlying germline genome architectures in ciliates and provide a basis for further exploration of principles of genome evolution in diverse microbial lineages.

A fundamental question in the biology of sex differences has eluded direct study in humans: How does sex-chromosome dosage (SCD) shape genome function? To address this, we developed a systematic map of SCD effects on gene function by analyzing genome-wide expression data in humans with diverse sex-chromosome aneuploidies (XO, XXX, XXY, XYY, and XXYY). For sex chromosomes, we demonstrate a pattern of obligate dosage sensitivity among evolutionarily preserved X-Y homologs and update prevailing theoretical models for SCD compensation by detecting X-linked genes that increase expression with decreasing X- and/or Y-chromosome dosage. We further show that SCD-sensitive sex-chromosome genes regulate specific coexpression networks of SCD-sensitive autosomal genes with critical cellular functions and a demonstrable potential to mediate previously documented SCD effects on disease. These gene coexpression results converge with analysis of transcription factor binding site enrichment and measures of gene expression in murine knockout models to spotlight the dosage-sensitive X-linked transcription factor ZFX as a key mediator of SCD effects on wider genome expression. Our findings characterize the effects of SCD broadly across the genome, with potential implications for human phenotypic variation.

Terpenes are economically and ecologically important phytochemicals. Their synthesis is controlled by the terpene synthase (TPS) gene family, which is highly diversified throughout the plant kingdom. The plant family Myrtaceae are characterised by especially high terpene concentrations, and considerable variation in terpene profiles. Many Myrtaceae are grown commercially for terpene products including the eucalypts Corymbia and Eucalyptus. Eucalyptus grandis has the largest TPS gene family of plants currently sequenced, which is largely conserved in the closely related E. globulus. However, the TPS gene family has been well studied only in these two eucalypt species. The recent assembly of two Corymbia citriodora subsp. variegata genomes presents an opportunity to examine the conservation of this important gene family across more divergent eucalypt lineages. Manual annotation of the TPS gene family in C. citriodora subsp. variegata revealed a similar overall number, and relative subfamily representation, to that previously reported in E. grandis and E. globulus. Many of the TPS genes were in physical clusters that varied considerably between Eucalyptus and Corymbia, with several instances of translocation, expansion/contraction and loss. Notably, there was greater conservation in the subfamilies involved in primary metabolism than those involved in secondary metabolism, likely reflecting different selective constraints. The variation in cluster size within subfamilies and the broad conservation between the eucalypts in the face of this variation are discussed, highlighting the potential contribution of selection, concerted evolution and stochastic processes. These findings provide the foundation to better understand terpene evolution within the ecologically and economically important Myrtaceae.

We elucidated genomic and transcriptomic changes that accompany the evolution of melanoma from pre-malignant lesions by sequencing DNA and RNA from primary melanomas and their adjacent precursors, as well as matched primary tumors and regional metastases. In total, we analyzed 230 histopathologically distinct areas of melanocytic neoplasia from 82 patients. Somatic alterations sequentially induced mitogen-activated protein kinase (MAPK) pathway activation, upregulation of telomerase, modulation of the chromatin landscape, G1/S checkpoint override, ramp-up of MAPK signaling, disruption of the p53 pathway, and activation of the PI3K pathway; no mutations were specifically associated with metastatic progression, as these pathways were perturbed during the evolution of primary melanomas. UV radiation-induced point mutations steadily increased until melanoma invasion, at which point copy-number alterations also became prevalent.

In canines, transposon dynamics have been associated with a hyper-social behavioral syndrome, although the functional mechanism has yet to be described. We investigate the epigenetic and transcriptional consequences of these behavior-associated mobile element insertions (MEIs) in dogs and Yellowstone gray wolves. We posit that the transposons themselves may not be the causative feature; rather, their transcriptional regulation may exert the functional impact. We survey four outlier transposons associated with hyper-sociability, with the expectation that they are targeted for epigenetic silencing. We predict hyper-methylation of MEIs, suggestive that the epigenetic silencing of and not the MEIs themselves may be driving dysregulation of nearby genes. We found that transposon-derived sequences are significantly hyper-methylated, regardless of their copy number or species. Further, we have assessed transcriptome sequence data and found evidence that MEIs impact the expression levels of six genes (WBSCR17, LIMK1, GTF2I, WBSCR27, BAZ1B, and BCL7B), all of which have known roles in human Williams-Beuren syndrome due to changes in copy number, typically hemizygosity. Although further evidence is needed, our results suggest that a few insertions alter local expression at multiple genes, likely through a cis-regulatory mechanism that excludes proximal methylation.

Stably transmitted transgenes are indispensable for labeling cellular components and manipulating cellular functions. In Caenorhabditis elegans, transgenes are generally generated as inheritable multi-copy extrachromosomal arrays, which can be stabilized in the genome through a mutagenesis-mediated integration process. Standard methods to integrate extrachromosomal arrays primarily use protocols involving ultraviolet light plus trimethylpsoralen or gamma- or X-ray irradiation, which are laborious and time-consuming. Here, we describe a one-step integration method, following germline-mutagenesis induced by mini Singlet Oxygen Generator (miniSOG). Upon blue light treatment, miniSOG tagged to histone (Histone-miniSOG) generates reactive oxygen species (ROS) and induces heritable mutations, including DNA double-stranded breaks. We demonstrate that we can bypass the need to first establish extrachromosomal transgenic lines by coupling microinjection of desired plasmids with blue light illumination on Histone-miniSOG worms to obtain integrants in the F3 progeny. We consistently obtained more than one integrant from 12 injected animals in two weeks. This optogenetic approach significantly reduces the amount of time and labor for transgene integration. Moreover, it enables to generate stably expressed transgenes that cause toxicity in animal growth.

In canines, transposon dynamics have been associated with a hyper-social behavioral syndrome, although the functional mechanism has yet to be described. We investigate the epigenetic and transcriptional consequences of these behavior-associated mobile element insertions (MEIs) in dogs and Yellowstone gray wolves. We posit that the transposons themselves may not be the causative feature; rather, their transcriptional regulation may exert the functional impact. We survey four outlier transposons associated with hyper-sociability, with the expectation that they are targeted for epigenetic silencing. We predict hyper-methylation of MEIs, suggestive that the epigenetic silencing of and not the MEIs themselves may be driving dysregulation of nearby genes. We found that transposon-derived sequences are significantly hyper-methylated, regardless of their copy number or species. Further, we have assessed transcriptome sequence data and found evidence that MEIs impact the expression levels of six genes (WBSCR17, LIMK1, GTF2I, WBSCR27, BAZ1B, and BCL7B), all of which have known roles in human Williams-Beuren syndrome due to changes in copy number, typically hemizygosity. Although further evidence is needed, our results suggest that a few insertions alter local expression at multiple genes, likely through a cis-regulatory mechanism that excludes proximal methylation.

The end-Cretaceous extinction led to a massive faunal turnover, with placental mammals radiating in the wake of nonavian dinosaurs. Fossils indicate that Cretaceous stem placentals were generally insectivorous, whereas their earliest Cenozoic descendants occupied a variety of dietary niches. It is hypothesized that this dietary radiation resulted from the opening of niche space, following the extinction of dinosaurian carnivores and herbivores. We provide the first genomic evidence for the occurrence and timing of this dietary radiation in placental mammals. By comparing the genomes of 107 placental mammals, we robustly infer that chitinase genes (CHIAs), encoding enzymes capable of digesting insect exoskeletal chitin, were present as five functional copies in the ancestor of all placental mammals, and the number of functional CHIAs in the genomes of extant species positively correlates with the percentage of invertebrates in their diets. The diverse repertoire of CHIAs in early placental mammals corroborates fossil evidence of insectivory in Cretaceous eutherians, with descendant lineages repeatedly losing CHIAs beginning at the Cretaceous/Paleogene (K/Pg) boundary as they radiated into noninsectivorous niches. Furthermore, the timing of gene loss suggests that interordinal diversification of placental mammals in the Cretaceous predates the dietary radiation in the early Cenozoic, helping to reconcile a long-standing debate between molecular timetrees and the fossil record. Our results demonstrate that placental mammal genomes, including humans, retain a molecular record of the post-K/Pg placental adaptive radiation in the form of numerous chitinase pseudogenes.

Genome-wide association studies (GWAS) have identified over 100 risk loci for schizophrenia, but the causal mechanisms remain largely unknown. We performed a transcriptome-wide association study (TWAS) integrating a schizophrenia GWAS of 79,845 individuals from the Psychiatric Genomics Consortium with expression data from brain, blood, and adipose tissues across 3,693 primarily control individuals. We identified 157 TWAS-significant genes, of which 35 did not overlap a known GWAS locus. Of these 157 genes, 42 were associated with specific chromatin features measured in independent samples, thus highlighting potential regulatory targets for follow-up. Suppression of one identified susceptibility gene, mapk3, in zebrafish showed a significant effect on neurodevelopmental phenotypes. Expression and splicing from the brain captured most of the TWAS effect across all genes. This large-scale connection of associations to target genes, tissues, and regulatory features is an essential step in moving toward a mechanistic understanding of GWAS.

Variants of TREM2 are associated with Alzheimer's disease (AD). To study whether increasing TREM2 gene dosage could modify the disease pathogenesis, we developed BAC transgenic mice expressing human TREM2 (BAC-TREM2) in microglia. We found that elevated TREM2 expression reduced amyloid burden in the 5xFAD mouse model. Transcriptomic profiling demonstrated that increasing TREM2 levels conferred a rescuing effect, which includes dampening the expression of multiple disease-associated microglial genes and augmenting downregulated neuronal genes. Interestingly, 5xFAD/BAC-TREM2 mice showed further upregulation of several reactive microglial genes linked to phagocytosis and negative regulation of immune cell activation. Moreover, these mice showed enhanced process ramification and phagocytic marker expression in plaque-associated microglia and reduced neuritic dystrophy. Finally, elevated TREM2 gene dosage led to improved memory performance in AD models. In summary, our study shows that a genomic transgene-driven increase in TREM2 expression reprograms microglia responsivity and ameliorates neuropathological and behavioral deficits in AD mouse models.

rDNA loci, composed of hundreds of tandemly duplicated arrays of rRNA genes, are known to be among the most unstable genetic elements due to their repetitive nature. rDNA instability underlies aging (replicative senescence) in yeast cells, however, its contribution to the aging of multicellular organisms is poorly understood. In this study, we investigate the dynamics of rDNA loci during aging in the Drosophila male germline stem cell (GSC) lineage, and show that rDNA copy number decreases during aging. Our study further reveals that this age-dependent decrease in rDNA copy number is heritable from generation to generation, yet GSCs in young animals that inherited reduced rDNA copy number are capable of recovering normal rDNA copy number. Based on these findings, we propose that rDNA loci are dynamic genetic elements, where rDNA copy number changes dynamically yet is maintained through a recovery mechanism in the germline.

Dynamic evolutionary processes and complex structure make the Y chromosome among the most diverse and least understood regions in mammalian genomes. Here, we present an annotated assembly of the male specific region of the horse Y chromosome (eMSY), representing the first comprehensive Y assembly in odd-toed ungulates. The eMSY comprises single-copy, equine specific multi-copy, PAR transposed, and novel ampliconic sequence classes. The eMSY gene density approaches that of autosomes with the highest number of retained X-Y gametologs recorded in eutherians, in addition to novel Y-born and transposed genes. Horse, donkey and mule testis RNAseq reveals several candidate genes for stallion fertility. A novel testis-expressed XY ampliconic sequence class, ETSTY7, is shared with the parasite Parascaris genome, providing evidence for eukaryotic horizontal transfer and inter-chromosomal mobility. Our study highlights the dynamic nature of the Y and provides a reference sequence for improved understanding of equine male development and fertility.

Reducing tillage and growing cover crops, widely recommended practices for boosting soil health, have major impacts on soil communities. Surprisingly little is known about their impacts on soil microbial functional diversity, and especially so in irrigated Mediterranean ecosystems. In long-term experimental plots at the West Side Research and Extension Center in California's Central Valley, we characterized soil microbial communities in the presence or absence of physical disturbance due to tillage, in the presence or absence of cover crops, and at three depths: 0-5, 5-15 and 15-30 cm. This characterization included qPCR for bacterial and archaeal abundances, DNA sequencing of the 16S rRNA gene, and phylogenetic estimation of two ecologically important microbial traits (rRNA gene copy number and genome size). Total (bacterial + archaeal) diversity was higher in no-till than standard till; diversity increased with depth in no-till but decreased with depth in standard till. Total bacterial numbers were higher in cover cropped plots at all depths, while no-till treatments showed higher numbers in 0-5 cm but lower numbers at lower depths compared to standard tillage. Trait estimates suggested that different farming practices and depths favored distinctly different microbial life strategies. Tillage in the absence of cover crops shifted microbial communities towards fast growing competitors, while no-till shifted them toward slow growing stress tolerators. Across all treatment combinations, increasing depth resulted in a shift towards stress tolerators. Cover crops shifted the communities towards ruderals-organisms with wider metabolic capacities and moderate rates of growth. Overall, our results are consistent with decreasing nutrient availability with soil depth and under no-till treatments, bursts of nutrient availability and niche homogenization under standard tillage, and increases in C supply and variety provided by cover crops. Understanding how agricultural practices shift microbial abundance, diversity and life strategies, such as presented here, can assist with designing farming systems that can support high yields, while enhancing C sequestration and increasing resilience to climate change.

Gene dosage imbalance caused by copy number variations (CNVs) is a prominent contributor to brain disorders. In particular, 15q11.2 CNV duplications and deletions have been associated with autism spectrum disorder and schizophrenia, respectively. The mechanism underlying these diametric contributions remains unclear. We established both loss-of-function and gain-of-function mouse models of Cyfip1 , one of four genes within 15q11.2 CNVs. To assess the functional consequences of altered CYFIP1 levels, we performed systematic investigations on behavioral, electrophysiological, and biochemical phenotypes in both mouse models. In addition, we utilized RNA immunoprecipitation sequencing (RIP-seq) analysis to reveal molecular targets of CYFIP1 in vivo. Cyfip1 loss-of-function and gain-of function mouse models exhibited distinct and shared behavioral abnormalities related to autism spectrum disorder and schizophrenia. RIP-seq analysis identified messenger RNA targets of CYFIP1 in vivo, including postsynaptic NMDA receptor (NMDAR) complex components. In addition, these mouse models showed diametric changes in levels of postsynaptic NMDAR complex components at synapses because of dysregulated protein translation, resulting in bidirectional alteration of NMDAR-mediated signaling. Importantly, pharmacological balancing of NMDAR signaling in these mouse models with diametric Cyfip1 dosages rescues behavioral abnormalities. CYFIP1 regulates protein translation of NMDAR and associated complex components at synapses to maintain normal synaptic functions and behaviors. Our integrated analyses provide insight into how gene dosage imbalance caused by CNVs may contribute to divergent neuropsychiatric disorders. [ABSTRACT FROM AUTHOR]

Robustness to perturbation is a fundamental feature of complex organisms. Mutations are the raw material for evolution, yet robustness to their effects is required for species survival. The mechanisms that produce robustness are poorly understood. Nonlinearities are a ubiquitous feature of development that may link variation in development to phenotypic robustness. Here, we manipulate the gene dosage of a signaling molecule, Fgf8, a critical regulator of vertebrate development. We demonstrate that variation in Fgf8 expression has a nonlinear relationship to phenotypic variation, predicting levels of robustness among genotypes. Differences in robustness are not due to gene expression variance or dysregulation, but emerge from the nonlinearity of the genotype-phenotype curve. In this instance, embedded features of development explain robustness differences. How such features vary in natural populations and relate to genetic variation are key questions for unraveling the origin and evolvability of this feature of organismal development.

Transposons can impact the host genome by altering gene expression and participating in chromosome rearrangements. Therefore, organisms evolved different ways to minimize the level of transposition. In Saccharomyces cerevisiae and its close relative S. paradoxus, Ty1 copy number control (CNC) is mediated by the self-encoded restriction factor p22, which is derived from the GAG capsid gene and inhibits virus-like particle (VLP) assembly and function. Based on secondary screens of Ty1 cofactors, we identified LOC1, a RNA localization/ribosome biogenesis gene that affects Ty1 mobility predominantly in strains harboring Ty1 elements. Ribosomal protein mutants rps0bΔ and rpl7aΔ displayed similar CNC-specific phenotypes as loc1Δ, suggesting that ribosome biogenesis is critical for CNC. The level of Ty1 mRNA and Ty1 internal (Ty1i) transcripts encoding p22 was altered in these mutants, and displayed a trend where the level of Ty1i RNA increased relative to full-length Ty1 mRNA. The level of p22 increased in these mutants, and the half-life of p22 also increased in a loc1Δ mutant. Transcriptomic analyses revealed small changes in the level of Ty1 transcripts or efficiency of translation initiation in a loc1Δ mutant. Importantly, a loc1Δ mutant had defects in assembly of Gag complexes and packaging Ty1 RNA. Our results indicate that defective ribosome biogenesis enhances CNC by increasing the level of p22, and raise the possibility for versatile links between VLP assembly, its cytoplasmic environment, and a novel stress response.

The ecological concept of the r-K life history strategy is widely applied in macro-ecology to characterize functional traits of taxa. However, its adoption in microbial communities is limited, owing to the lack of a measureable, convenient functional trait for classification. In this study, we performed an experiment of stepwise organic amendments in triplicate anaerobic digesters. We found that high resource availability significantly favored microbial r-strategists such as Bacillus spp. Incremental resource availability heightened average rRNA operon copy number of microbial community, resulting in a strong, positive correlation (r>0.74, P

Shank2 is an abundant excitatory postsynaptic scaffolding protein that has been implicated in various neurodevelopmental and psychiatric disorders, including autism spectrum disorder (ASD), intellectual disability, attentiondeficit/hyperactivity disorder, and schizophrenia. Shank2-mutant mice show ASD-like behavioral deficits and altered synaptic and neuronal functions, but little is known about how different brain regions and gene dosages affect the transcriptomic phenotypes of these mice. Here, we performed RNA-Seqbased transcriptomic analyses of the prefrontal cortex, hippocampus, and striatum in adult Shank2 heterozygous (HT)- and homozygous (HM)-mutant mice lacking exons 6-7. The prefrontal cortical, hippocampal, and striatal regions showed distinct transcriptomic patterns associated with synapse, ribosome, mitochondria, spliceosome, and extracellular matrix (ECM). The three brain regions were also distinct in the expression of ASD-related and ASD-risk genes. These differential patterns were stronger in the prefrontal cortex where the HT transcriptome displayed increased synaptic gene expression and reverse-ASD patterns whereas the HM transcriptome showed decreased synaptic gene expression and ASD-like patterns. These results suggest brain region- and gene dosage-differential transcriptomic changes in Shank2-mutant mice. [ABSTRACT FROM AUTHOR]

Shank3 is an abundant excitatory postsynaptic scaffolding protein implicated in various neurodevelopmental disorders, including autism spectrum disorder (ASD), Phelan-McDermid syndrome, intellectual disability, and schizophrenia. Shank3-mutant mice show various molecular, synaptic, and behavioral deficits, but little is known about how transcriptomic phenotypes vary across different ages, brain regions, and gene dosages. Here, we report transcriptomic patterns in the forebrains of juvenile and adult homozygous Shank3-mutant mice that lack exons 14-16 and also the prefrontal, hippocampal, and striatal transcriptomes in adult heterozygous and homozygous Shank3-mutant mice. The juvenile and adult mutant transcriptomes show patterns opposite from and similar to those observed in ASD (termed reverse-ASD and ASD-like patterns), respectively. The juvenile transcriptomic changes accompany synaptic upregulations and ribosomal and mitochondrial downregulations, whereas the adult transcriptome show opposite changes. The prefrontal, hippocampal, and striatal transcriptomes show differential changes in ASD-related gene expressions and biological functions associated with synapse, ribosome, mitochondria, and spliceosome. These patterns also differ across heterozygous and homozygous Shank3- mutant mice. These results suggest age, brain region, and gene dosagedifferential transcriptomic changes in Shank3-mutant mice. [ABSTRACT FROM AUTHOR]

Latrophilin‐3 (LPHN3) is a brain specific G‐protein coupled receptor associated with increased risk of attention deficit hyperactivity disorder (ADHD) and cognitive deficits. CRISPR/Cas9 was used to generate a constitutive knockout (KO) rat of Lphn3 by deleting exon 3, based on human data that LPHN3 variants are associated with some cases of ADHD. Lphn3 KO rats are hyperactive with an attenuated response to ADHD medication and have cognitive deficits. Here, we tested KO, heterozygous (HET), and wildtype (WT) rats to determine if there was a gene‐dosage effect. We tested the rats in home‐cage activity starting at postnatal day (P)35 and P50, followed by tests of egocentric learning (Cincinnati water maze [CWM]), spatial learning (Morris water maze [MWM]), working memory (radial water maze [RWM]), incidental learning (novel object recognition [NOR]), acoustic startle response (ASR) habituation, tactile startle response (TSR) habituation, prepulse modification of acoustic startle, shuttle‐box passive avoidance, conditioned freezing, and a mirror image version of the CWM. KO and HET rats were hyperactive. KO and HET rats had egocentric (CWM) and spatial deficits (MWM), increased startle response, and KO rats showed less conditioned freezing on contextual and cued memory; there were no effects on working memory (RWM) or passive avoidance. The selective gene‐dosage effect in Lphn3 HET rats indicates that Lphn3 exhibits dominate expression on functions where it is most abundantly expressed (striatum, hippocampus) but not on behaviors mediated by regions of low expression. The data add further evidence to the impact of this synaptic protein on brain function and behavior. [ABSTRACT FROM AUTHOR]

Gene duplication events happen prevalently during evolution, and the mechanisms governing the loss or retention of duplicated genes are mostly elusive. Our genome scanning analysis revealed that trigger factor (TF), the one and only bacterial ribosome-associated molecular chaperone, is singly copied in virtually every bacterium except for a very few that possess two or more copies. However, even in these exceptions, only one complete TF copy exists, while other homologs lack the N-terminal domain that contains the conserved ribosome binding site (RBS) motif. Consistently, we demonstrated that the overproduction of the N-terminal complete TF proteins is detrimental to the cell, which can be rescued by removing the N-terminal domain. Our findings also indicated that TF overproduction leads to a decrease in protein productivity and profile changes in proteome due to its characteristic ribosome binding and holdase activities. Additionally, these N-terminal deficient TF homologs in bacteria with multiple TF homologs partition the function of TF via subfunctionalization. Our results revealed that TF is subjected to a dosage constraint that originates from its own intrinsic functions, which may drive the evolution and fates of duplicated TFs in bacteria., Importance: Gene duplication events presumably occur in tig , which encodes the ribosome-associated molecular chaperone trigger factor (TF). However, TF is singly copied in virtually every bacterium, and these exceptions with multiple TF homologs always retain only one complete copy while other homologs lack the N-terminal domain. Here, we reveal the manner and mechanism underlying the evolution and fates of TF duplicates in bacteria. We discovered that the mutation-to-loss or retention-to-sub/neofunctionalization of TF duplicates is associated with the dosage constraint of N-terminal complete TF. The dosage constraint of TF is attributed to its characteristic ribosome binding and substrate-holding activities, causing a decrease in protein productivity and profile changes in cellular proteome., Competing Interests: The authors declare no conflict of interest.

Background: The copy number status (CNS) of the survival motor neuron (SMN) gene may influence the risk and prognosis of amyotrophic lateral sclerosis (ALS) and lower motor neuron diseases (LMND) other than spinal muscular atrophy (SMA). However, previous studies of this association, mainly from Europe, have yielded controversial results, suggesting possible regional differences. Here, we investigated the effect of the SMN gene in Japanese patients with ALS and LMND., Methods: We examined the SMN copy numbers and clinical histories of 487 Japanese patients with sporadic ALS (281 men; mean age at onset 61.5 years), 50 with adult LMND (50 men; mean age at onset 58.4 years) and 399 Japanese controls (171 men; mean age 62.2 years). Patients with pathogenic mutations in ALS-causing genes were excluded. SMN1 and SMN2 copy numbers were determined using the droplet digital polymerase chain reaction., Results: The frequency of a copy number of one for the SMN2 gene was higher in patients with ALS (38.0%) than in healthy controls (30.8%) (odds ratio (OR) = 1.37, 95% confidence interval (CI) = 1.04-1.82, p < 0.05). The SMN2 copy number affected the survival time of patients with ALS (median time: 0 copies, 34 months; 1 copy, 39 months; 2 copies, 44 months; 3 copies, 54 months; log-rank test, p < 0.05). Cox regression analysis revealed that the SMN2 copy number was associated with increased mortality (hazard ratio = 0.84, 95% CI = 0.72-0.98, p < 0.05). Also, null SMN2 cases were significantly more frequent in the LMND group (12.0%) than in the control group (4.8%) (OR = 2.73, 95% CI = 1.06-6.98, p < 0.05)., Conclusions: Our findings suggest that SMN2 copy number reduction may adversely affect the onset and prognosis of MND, including ALS and LMND, in Japanese., (© 2024. The Author(s).)

Background: Wip1, is a p53-dependent Ser/Thr phosphatase involved in the timely termination of DDR. The PPM1D gene encoding Wip1 is deregulated and thus gained an oncogene character in common human solid tumors and cell lines. This study assessed the oncogenic potential of the PPM1D gene in human non- Hodgkin's lymphomas (NHL), the most common hematological malignancy worldwide., Methods and Results: FFPE human lymphoid hyperplasia (LH) (n = 17) and NHL tumor lymph node samples (n = 65) and human NHL cell lines were used to assess the oncogenic potential of the PPM1D gene in the present study. Copy number gain and mRNA expression analysis of the PPM1D/Wip1 gene were assessed by qRT-PCR analysis. Mutational analysis of Exon 6 of the PPM1D gene was performed by PCR amplification and Sanger sequencing. Expressions of Wip1 and p53 proteins were assessed by immunohistochemistry and Western blot analysis., Conclusions: We found that PPM1D gained gene copy number in NHL tumors by 0.7-8 times compared to the control (p < 0.01). Increased PPM1D/Wip1 gene copy number was associated with higher mRNA and protein expression in human NHL samples (p < 0.01). Overexpression of Wip1 in NHL tumors and NHL cell lines was associated with amplification level and was unaffected by p53 status. Furthermore, a heterozygous type insertion mutation was detected in exon 6 (c.1553&#95;1554insA) of the PPM1D gene particularly in DLBCL samples. Wip1 may have oncogenic potential, perhaps playing a role in the onset and progression of human NHL. The possible significance of Wip1 overexpression to chemotherapy response in NHL remains an intriguing question that requires more exploration., (© 2024. The Author(s), under exclusive licence to Springer Nature B.V.)

Protein copy numbers constrain systems-level properties of regulatory networks, but proportional proteomic data remain scarce compared to RNA-seq. We related mRNA to protein statistically using best-available data from quantitative proteomics and transcriptomics for 4366 genes in 369 cell lines. The approach starts with a protein's median copy number and hierarchically appends mRNA-protein and mRNA-mRNA dependencies to define an optimal gene-specific model linking mRNAs to protein. For dozens of cell lines and primary samples, these protein inferences from mRNA outmatch stringent null models, a count-based protein-abundance repository, empirical mRNA-to-protein ratios, and a proteogenomic DREAM challenge winner. The optimal mRNA-to-protein relationships capture biological processes along with hundreds of known protein-protein complexes, suggesting mechanistic relationships. We use the method to identify a viral-receptor abundance threshold for coxsackievirus B3 susceptibility from 1489 systems-biology infection models parameterized by protein inference. When applied to 796 RNA-seq profiles of breast cancer, inferred copy-number estimates collectively re-classify 26-29% of luminal tumors. By adopting a gene-centered perspective of mRNA-protein covariation across different biological contexts, we achieve accuracies comparable to the technical reproducibility of contemporary proteomics., (© 2024. The Author(s).)

SNX29 is a potential functional gene associated with meat production traits. Previous studies have shown that SNX29 copy number variation (CNV) could be implicated with phenotype in goats. However, in Diannan small-ear (DSE) pigs, the genetic impact of SNX29 CNV on growth traits remains unclear. Therefore, this study investigated the associations between SNX29 CNVs (CNV10810 and CNV10811) and growth traits in 415 DSE pigs. The results revealed that the CNV10810 mutation was significantly associated with backfat thickness in DSE pigs at 12 and 15 months old ( P  < 0.05), while the CNV10811 mutation had significant effects on various growth traits at 6 and 12 months old, particularly for body weight, body height, back height and backfat thickness ( P  < 0.05 or P  < 0.001). In conclusion, our results confirm that SNX29 CNV plays a role in regulating growth and development in pigs, thus suggesting its potential application for pig breeding programmes.

Sweet potato ( Ipomoea batatas ) is an important food crop that plays a pivotal role in preserving worldwide food security. Due to its polyploid genome, high heterogeneity, and phenotypic plasticity, sweet potato genetic characterization and breeding is challenging. Genome-wide association studies (GWASs) can provide important resources for breeders to improve breeding efficiency and effectiveness. GWASpoly was used to identify 28 single nucleotide polymorphisms (SNPs), comprising 21 unique genetic loci, associated with sweet potato storage root traits including dry matter (4 loci), subjective flesh color (5 loci), flesh hue angle (3 loci), and subjective skin color and skin hue angle (9 loci), in 384 accessions from the USDA sweet potato germplasm collection. The I. batatas 'Beauregard' and I. trifida reference genomes were utilized to identify candidate genes located within 100 kb from the SNPs that may affect the storage traits of dry matter, flesh color, and skin color. These candidate genes include transcription factors (especially Myb, bHLH, and WRKY family members), metabolite transporters, and metabolic enzymes and associated proteins involved in starch, carotenoid, and anthocyanin synthesis. A greater understanding of the genetic loci underlying sweet potato storage root traits will enable marker-assisted breeding of new varieties with desired traits. This study not only reinforces previous research findings on genes associated with dry matter and β-carotene content but also introduces novel genetic loci linked to these traits as well as other root characteristics.

The loss of a single chromosome in a diploid organism halves the dosage of many genes and is usually accompanied by a substantial decrease in fitness. We asked whether this decrease simply reflects the joint damage caused by individual gene dosage deficiencies. We measured the fitness effects of single heterozygous gene deletions in yeast and combined them for each chromosome. This predicted a negative growth rate, that is, lethality, for multiple monosomies. However, monosomic strains remained alive and grew as if much (often most) of the damage caused by single mutations had disappeared, revealing an exceptionally large and positive epistatic component of fitness. We looked for functional explanations by analyzing the transcriptomes. There was no evidence of increased (compensatory) gene expression on the monosomic chromosomes. Nor were there signs of the cellular stress response that would be expected if monosomy led to protein destabilization and thus cytotoxicity. Instead, all monosomic strains showed extensive upregulation of genes encoding ribosomal proteins, but in an indiscriminate manner that did not correspond to their altered dosage. This response did not restore the stoichiometry required for efficient biosynthesis, which probably became growth limiting, making all other mutation-induced metabolic defects much less important. In general, the modular structure of the cell leads to an effective fragmentation of the total mutational load. Defects outside the module(s) currently defining fitness lose at least some of their relevance, producing the epiphenomenon of positive interactions between individually negative effects., Competing Interests: HT, KT, AP, MM, RK No competing interests declared, (© 2023, Tutaj et al.)

A female protective effect has long been postulated as the primary explanation for the four-fold increase of autism spectrum disorder (ASD) diagnoses in males versus females. However, genetic and epidemiological investigations of this hypothesis have so far failed to explain the large difference in ASD prevalence between the sexes. To address this knowledge gap, we examined sex chromosome aneuploidy in a large ASD case-control cohort to evaluate the relationship between X and Y chromosome dosage and ASD risk. From these data, we modeled three relationships between sex chromosome dosage and ASD risk: the extra Y effect, the extra X effect, and sex chromosome haploinsufficiency. We found that the extra Y effect increased ASD risk significantly more than the extra X effect. Among females, we observed a large association between 45, X and ASD, confirming sex chromosome haploinsufficiency as a strong ASD risk factor. These results provide a framework for understanding the relationship between X and Y chromosome dosage on ASD, which may inform future research investigating genomic contributors to the observed sex difference., (© 2024. The Author(s).)

The Four Core Genotypes (FCG) is a mouse model system used to disentangle the function of sex chromosomes and hormones. We report that a copy of a 3.2 MB region of the X chromosome has translocated to the Y Sry- chromosome and thus increased the expression of X-linked genes including the single-stranded RNA sensor and autoimmune disease mediator Tlr7. This previously-unreported X-Y translocation complicates the interpretation of studies reliant on C57BL/6J FCG mice., (© 2024. The Author(s).)

Gaucher disease (GD) is caused by biallelic GBA1/Gba1 mutations that encode defective glucocerebrosidase (GCase). Progranulin (PGRN, encoded by GRN/Grn) is a modifier of GCase, but the interplay between PGRN and GCase, specifically GBA1/Gba1 mutations, contributing to GD severity is unclear. Mouse models were developed with various dosages of Gba1 D409V mutation against the PGRN deficiency (Grn-/-) [Grn-/-;Gba1D409V/WT (PG9Vwt), Grn-/-;Gba1D409V/D409V (PG9V), Grn-/-;Gba1D409V/Null (PG9VN)]. Disease progression in those mouse models was characterized by biochemical, pathological, transcriptomic, and neurobehavioral analyses. Compared to PG9Vwt, Grn-/-;Gba1WT/Null and Grn-/- mice that had a higher level of GCase activity and undetectable pathologies, homozygous or hemizygous D409V in PG9V or PG9VN, respectively, resulted in profound inflammation and neurodegeneration. PG9VN mice exhibited much earlier onset, shorter life span, tissue fibrosis, and more severe phenotypes than PG9V mice. Glycosphingolipid accumulation, inflammatory responses, lysosomal-autophagy dysfunction, microgliosis, retinal gliosis, as well as α-Synuclein increases were much more pronounced in PG9VN mice. Neurodegeneration in PG9VN was characterized by activated microglial phagocytosis of impaired neurons and programmed cell death due to necrosis and, possibly, pyroptosis. Brain transcriptomic analyses revealed the intrinsic relationship between D409V dosage, and the degree of altered gene expression related to lysosome dysfunction, microgliosis, and neurodegeneration in GD, suggesting the disease severity is dependent on a GCase activity threshold related to Gba1 D409V dosage and loss of PGRN. These findings contribute to a deeper understanding of GD pathogenesis by elucidating additional underlying mechanisms of interplay between PGRN and Gba1 mutation dosage in modulating GCase function and disease severity in GD and GBA1-associated neurodegenerative diseases., (© The Author(s) 2024. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

Tumors originating from thyroid follicular cells are the most common endocrine tumors, with rising incidence. Despite a generally good prognosis, up to 20% of patients experience recurrence and persistence, highlighting the need to identify the underlying molecular mechanisms. Dicer1 has been found to be altered in papillary thyroid cancer (PTC). Studies suggest that Dicer1 functions as a haploinsufficient tumor suppressor gene: partial loss promotes tumorigenesis, while complete loss prevents it. To investigate the effects of partial or total Dicer1 loss in PTC in vitro, we generated stable Dicer1 (+/-) cell lines from TPC1 using CRISPR-Cas9, though no Dicer1 (-/-) lines could be produced. Therefore, siRNA against Dicer1 was transfected into Dicer1 (+/-) cell lines to further decrease its expression. Transcriptomic analysis revealed changes in proliferation and cell locomotion. BrdU staining indicated a slow-down of the cell cycle, with fewer cells in S phase and more in G0-G1-phase. Additionally, transwell assays showed decreased invasion and migration after Dicer1 knockdown by siRNA. Moreover, Dicer1 overexpression led to decreased proliferation, invasion, and increased apoptosis. Our findings deepen the understanding of Dicer1 's role in thyroid cancer, demonstrating that both complete elimination and overexpression of Dicer1 inhibit thyroid oncogenesis, highlighting Dicer1 as a promising target for novel therapeutic strategies.

Following whole-genome duplication (WGD), duplicate gene pairs (homoeologs) can evolve varying degrees of expression divergence. However, the determinants influencing these relative expression level differences (RFPKM) between homoeologs remain elusive. In this study, we analyzed the RFPKM between homoeologs in 3 angiosperms, Nymphaea colorata, Nelumbo nucifera, and Acorus tatarinowii, all having undergone a single WGD since the origin of angiosperms. Our results show significant positive correlations in RFPKM of homoeologs among tissues within the same species, and among orthologs across these 3 species, indicating convergent expression balance/bias between homoeologous gene copies following independent WGDs. We linked RFPKM between homoeologs to gene attributes associated with dosage-balance constraints, such as protein-protein interactions, lethal-phenotype scores in Arabidopsis (Arabidopsis thaliana) orthologs, domain numbers, and expression breadth. Notably, homoeologs with lower RFPKM often had more interactions and higher lethal-phenotype scores, indicating selective pressures favoring balanced expression. Also, homoeologs with lower RFPKM were more likely to be retained after WGDs in angiosperms. Within Nelumbo, greater RFPKM between homoeologs correlated with increased cis- and trans-regulatory differentiation between species, highlighting the ongoing escalation of gene expression divergence. We further found that expression degeneration in 1 copy of homoeologs is inclined toward nonfunctionalization. Our research highlights the importance of balanced expression, shaped by dosage-balance constraints, in the evolutionary retention of homoeologs in plants., Competing Interests: Conflict of interest statement. None declared., (© The Author(s) 2024. Published by Oxford University Press on behalf of American Society of Plant Biologists. All rights reserved. For commercial re-use, please contact reprints@oup.com for reprints and translation rights for reprints. All other permissions can be obtained through our RightsLink service via the Permissions link on the article page on our site—for further information please contact journals.permissions@oup.com.)

Human salivary amylase secretion increases in response to stress; the activity has been reported to rise significantly with high-intensity exercise. The human salivary amylase gene (AMY1) has copy number variation, with the copy number correlating with salivary amylase activity. However, the relationship between individual AMY1 copy number and salivary amylase activity in response to exercise remains unclear. In this study, we investigated AMY1 copy number and fluctuations in amylase activity in 42 healthy university students (25 males and 17 females). Participants engaged in intermittent round-trip interval training on a basketball court. Saliva samples were collected pre- and post-exercise to measure amylase activity. DNA was extracted from the oral mucosa, and AMY1 copy number was quantified using RT-PCR. Results showed a significant increase in amylase activity postexercise. Additionally, amylase activity pre- and post-exercise was positively correlated with AMY1 copy number. The generalize linear model showed that the exercise-induced increase in amylase activity per AMY1 gene was negatively related to the AMY1 copy number and aerobic fitness. Gender has no effect on amylase activity. These results suggest a different mechanism for the constitutive and exercise-induced amylase secretion, while aerobic fitness may be independently involved in the secretion., (© 2024 The Author(s). Physiological Reports published by Wiley Periodicals LLC on behalf of The Physiological Society and the American Physiological Society.)

Polyploidization in plants often leads to increased cell size and grain size, which may be affected by the increased genome dosage and transcription abundance. The synthesized Triticum durum (AABB)-Haynaldia villosa (VV) amphiploid (AABBVV) has significantly increased grain size, especially grain length, than the tetraploid and diploid parents. To investigate how polyploidization affects grain development at the transcriptional level, we perform transcriptome analysis using the immature seeds of T. durum, H. villosa, and the amphiploid. The dosage effect genes are contributed more by differentially expressed genes from genome V of H. villosa. The dosage effect genes overrepresent grain development-related genes. Interestingly, the vernalization gene TaVRN1 is among the positive dosage effect genes in the T. durum‒H. villosa and T. turgidum‒Ae. tauschii amphiploids. The expression levels of TaVRN1 homologs are positively correlated with the grain size and weight. The TaVRN1-B1 or TaVRN1-D1 mutation shows delayed florescence, decreased cell size, grain size, and grain yield. These data indicate that dosage effect genes could be one of the important explanations for increased grain size by regulating grain development. The identification and functional validation of dosage effect genes may facilitate the finding of valuable genes for improving wheat yield., Competing Interests: Conflict of interest The authors declare no competing interests., (Copyright © 2024 Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, and Genetics Society of China. Published by Elsevier Ltd. All rights reserved.)

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., (© 2024. The Author(s).)

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., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

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., (© 2024 The Author(s). Molecular Plant Pathology published by British Society for Plant Pathology and John Wiley & Sons Ltd. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.)

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., (© 2024 The Author(s). Genes, Chromosomes and Cancer published by Wiley Periodicals LLC.)

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.