Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Preprint”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 451 records · Page 25Linked to original sources

Evidence of parental care as a newly identified reproductive isolating barrier.

Variation in behavior can contribute to reproductive isolation by preventing gene flow among populations. Here, we tested the novel hypothesis that parental care, when dysregulated, can function as a reproductive isolating mechanism in three-spined stickleback fish (Gasterosteus aculeatus). In the typical "common" stickleback ecotype, males provide care to their offspring by fanning with their pectoral fins and defending their nest. In contrast, a divergent "white" stickleback ecotype has evolutionarily lost care and disperses embryos into the surrounding environment. We examined how paternal care from common, white, and F1 hybrid fathers influenced offspring survival. We detected no intrinsic incompatibilities in embryos, but F1 hybrid fathers exhibited dysregulated parental care that coincided with decreased survival of parented offspring. The increased offspring mortality may be explained by further dysregulation of feeding and parenting circuitry, as F2 hybrids exhibited significantly higher rates of post-fertilization filial cannibalism than common or white fathers. Additionally, despite strong divergence in nesting-building and courtship behavior, F1 hybrids achieved mating success at a similar rate as male commons and whites, suggesting that prezygotic barriers against hybrids may be weak. The observed postzygotic isolation, and potentially weak prezygotic isolation in lab-based studies, suggest that hybridization is likely occurring at low rates in the wild. Population genetic analysis supported this, as low proportions of putative hybrids were detected in sympatric sites. Together, these results may help explain why genetic divergence between these phenotypically distinct ecotypes is low and provide evidence that dysregulated parental behavior can act as a newly discovered postzygotic reproductive isolating barrier.

Reproductive isolation↗

The landscape of structural variation in pediatric cancer.

Structural variants (SVs) account for over 60% of the driver variants in pediatric cancer, and in many cases act as the cancer initiating event. To study SVs from a pan-cancer perspective, we analyzed 1,616 pediatric cancer genomes in 16 major cancer types of hematological malignancies (n = 908), brain tumors (n = 183), and solid tumors (n = 525) and compared their profiles to those of 2,203 adult cancers. The SV burden varied ~100-fold across pediatric cancer types and demonstrated an 8- to 16-fold reduction compared to adult brain and solid tumors but was comparable in pediatric versus adult hematological malignancies. Recurrent SV hotspots occurred uniquely in pediatric acute lymphoblastic leukemias (ALLs) in proximity to RAG-mediated recombination signal sequences (RSS) and disrupted multiple immune-related loci as well as 69 genes, which often involved cryptic RSS sites. By contrast, such hotspots affected only immune-related loci but not driver genes in adult lymphoid cancers. Eight SV signatures extracted from the cohort had varying distributions across cancer types, with clustered translocations reflecting templated insertions in osteosarcoma, and medium-sized deletions (10 kb to 1 Mb) enriched in cancers with RAG-mediated deletions. Intra-patient evolutionary analysis in 13 patients with multiple spatiotemporally distinct samples revealed that RAG-mediated recombination in leukemia and complex rearrangements in solid tumors occurred both early in disease initiation and continuously during later diversification, contributing to clonal heterogeneity. Finally, we found that both driver genes and fragile sites were the two genomic regions most frequently disrupted by SVs. The unique and diverse SV landscapes that emerged from this comprehensive analysis expand the scope of RSS-mediated mutagenesis in pediatric ALL and will be a valuable resource for guiding future functional studies and the design of clinical genomic testing in pediatric cancer.

Journal Article↗

Tulp3 quantitative alleles titrate requirements for viability, brain development, and kidney homeostasis but do not suppress Zfp423 mutations in mice.

Tubby-like protein 3 (TULP3) regulates receptor trafficking in primary cilia and antagonizes SHH signaling. Tulp3 knockout mice are embryonic lethal with developmental abnormalities in multiple organs, while tissue-specific knockouts and viable missense alleles cause polycystic kidney disease. Human patients with TULP3 mutations present with variable, but often multi-organ fibrotic disease. We previously showed that mouse and human Tulp3 expression is negatively regulated by ZNF423, which is required for SHH sensitivity in some progenitor cell models. The level of TULP3 function required to prevent mutant phenotypes has not been known. Here we report a Tulp3 quantitative allelic series, designed by targeting the polypyrimidine tract 5' to the splice acceptor of a critical exon, that shows distinct dose-response effects on viability, brain overgrowth, weight gain, and cystic kidney disease. We find limited evidence for genetic interaction with Zfp423 null or hypomorphic mutations. Together, these results establish an approach to developing quantitative allelic series by exon exclusion, rank-order dose-sensitivity of Tulp3 phenotypes, and model thresholds for TULP3 function to prevent severe outcomes.

Journal Article↗

Separating direct, indirect and parent-of-origin genetic effects in the human population.

Here, we present a novel approach to estimate the degree to which the phenotypic effect of a DNA locus is attributable to four components: alleles in the child (direct genetic effects), alleles in the mother and the father (indirect genetic effects), or is dependent upon the parent from which it is inherited (parent-of-origin, PofO effects). Applying our model, JODIE, to 30,000 child-mother-father trios with phased DNA information from the Estonian Biobank (EstBB) and the Norwegian Mother, Father, Child Cohort (MoBa), we jointly estimate the phenotypic variance attributable to these four effects unbiased of assortative mating (AM) for height, body mass index (BMI) and childhood educational test score (EA). For all three traits, direct effects make the largest contribution to the genetic effect variance. But we find that parental indirect genetic effects make an equivalent combined contribution, and that there is a non-zero PofO effect variance for all traits. We calculate the heritability that would be obtained at the population-level in the absence of AM for common DNA loci, and show that the proportional contribution of direct effects to these heritability values can be calculated as 64.0% for EA in MoBa, 77.1% and 63.4% for height in MoBa and EstBB, and 81.2% and 88.0% for BMI in MoBa and EstBB. Additionally, using within-family genome-wide association testing, we identify 276 independently associated DNA regions that replicate across two additional biobanks, which all show a genotype-phenotype relationship that reflects an interplay of direct, indirect and PofO effects. Determining how direct, parental and PofO genetic effects combine across loci genome-wide to influence human phenotypic variation requires joint modeling of parental and child genotypes alongside the parental origin of loci and here, we make the first attempt to do this in the human population.

EstBB↗

Maternal Immune Activation Disrupts Epigenomic and Functional Maturation of Cortical Excitatory Neurons.

Elevated levels of maternal pro-inflammatory cytokines during gestation can disrupt offspring neural development, increasing the risk of neurodevelopmental disorders. We studied the effects of Poly(I:C)-induced maternal immune activation (PIC-MIA) during mid-gestation on developing cortical excitatory neurons' DNA methylation and transcriptome. PIC-MIA disrupted the developmental regulation of synapse-related genes and of genes implicated in autism spectrum disorders. Genomic regions that gain or lose DNA methylation during normal development were altered following PIC-MIA, including neurodevelopmental transcription factor binding sites. The DNA methylation and transcriptional changes were consistent with a delay in excitatory neuron maturation. Whole-cell recordings showed that PIC-MIA preferentially altered the physiological development of layer 5 excitatory neurons. Taken together, present results suggest that alterations in the epigenome, through the disruption of circuit formation, may drive the long-term consequences of maternal infection during gestation.

DNA methylation↗

Lamin A/C Deficiency Drives Genomic Instability and Poor Survival in Small-Cell Lung Cancer through Increased R-loop Accumulation.

Lamin A/C (LMNA), a key component of the nuclear envelope, is essential for maintaining nuclear integrity and genome organization [1]. While LMNA dysregulation has been implicated in genomic instability across cancer and aging, the underlying mechanisms remain poorly understood [2]. Here, we investigate LMNA's role in small-cell lung cancer (SCLC), a highly aggressive malignancy characterized by extreme genomic instability [3, 4]. We demonstrate that LMNA depletion promotes R-loop accumulation, transcription-replication conflicts, replication stress, DNA breaks, and micronuclei formation. Mechanistically, LMNA loss disrupts nuclear pore complex distribution, reducing phenylalanine-glycine (FG)-nucleoporin incorporation and impairing RNA export efficiency. Furthermore, we show that LMNA expression is epigenetically repressed by EZH2 during SCLC differentiation from neuroendocrine (NE) to non-NE states. Clinically, low LMNA levels correlate with significantly worse survival in SCLC patients. These findings uncover a novel role for LMNA in safeguarding genome integrity and shaping tumor heterogeneity, with broad implications for cancer and aging.

Biological Sciences↗

Analysis of duplication and possible sub-functionalization of wing gene network components in pea aphids.

A fundamental focus of evolutionary-developmental biology is uncovering the genetic mechanisms responsible for the gain and loss of characters. One approach to this question is to investigate changes in the coordinated expression of a group of genes important for the development of a character of interest (a gene regulatory network). Here we consider the possibility that modifications to the wing gene regulatory network (wGRN), as defined by work primarily done in Drosophila melanogaster, were involved in the evolution of wing dimorphisms of the pea aphid (Acyrthosiphon pisum). We hypothesize that this may have occurred via changes in expression levels or duplication followed by sub-functionalization of wGRN components. To test this, we annotated members of the wGRN in the pea aphid genome and assessed their expression levels in first and third nymphal instars of winged and wingless morphs of males and asexual females. We find that only two of the 32 assessed genes exhibit morph-biased expression. We also find that three wing genes (apterous (ap), warts (wts), and decapentaplegic (dpp)) have undergone gene duplication. In each case, the resulting paralogs show signs of functional divergence, exhibiting either sex-, morph-, or stage-specific expression. Two gene duplicates, wts2 and dpp3, are of particular interest with respect to wing dimorphism, as they exhibit a wingless male-specific isoform and wingless male-biased expression, respectively. These results supplement our understanding of trends in developmental gene network evolution, such as side-stepping pleiotropic constraint via duplication and sub-functionalization, underlying the emergence of novel phenotypes.

dimorphism↗

Improving RNA Secondary Structure Prediction Through Expanded Training Data.

In recent years, deep learning has revolutionized protein structure prediction, achieving remarkable speed and accuracy. RNA structure prediction, however, has lagged behind. Although several methods have shown some success in predicting RNA secondary and tertiary structures, none have reached the accuracy observed with contemporary protein models. The lack of success of these RNA structure prediction models has been proposed to be due to limited high-quality structural information that can be used as training data. To probe this proposed limitation, we developed a large and diverse dataset comprising paired RNA sequences and their corresponding secondary structures. We assess the utility of this enhanced dataset by retraining on a deep learning model, SincFold. We find that SincFold exhibited improved generalization to some previously unseen RNA families, enhancing its capability to predict accurate de novo RNA secondary structures. The RNASSTR dataset provides a substantial advance for RNA structure modeling, laying a strong foundation for the development of future RNA secondary structure prediction algorithms.

Journal Article↗

Androgens mediate sexual dimorphism in Pilarowski-Bjornsson Syndrome.

Sex-specific penetrance in autosomal dominant Mendelian conditions is largely understudied. The neurodevelopmental disorder Pilarowski-Bjornsson syndrome (PILBOS) was initially described in females. Here, we describe the clinical and genetic characteristics of the largest PILBOS cohort to date, showing that both sexes can exhibit PILBOS features, although males are overrepresented. A mouse model carrying a human-derived Chd1 missense variant (Chd1 R616Q/+) displays female-restricted phenotypes, including growth deficiency, anxiety and hypotonia. Orchiectomy unmasks a growth deficiency phenotype in male Chd1 R616Q/+ mice, while testosterone rescues the phenotype in females, implicating androgens in phenotype modulation. In the gnomAD and UK Biobank databases, rare missense variants in CHD1 are overrepresented in males, supporting a male protective effect. We identify 33 additional highly constrained autosomal genes with missense variant overrepresentation in males. Our results support androgen-regulated sexual dimorphism in PILBOS and open novel avenues to understand the mechanistic basis of sexual dimorphism in other autosomal Mendelian disorders.

CHD1↗

Disruption of Microhomology-mediated End-joining in Ewing Sarcoma.

Ewing sarcoma (EwS) is a group of bone and soft tissue cancers in children and young adults. Since EwS cells have pronounced sensitivity to radiation and chemotherapy-induced DNA damage, the role of the oncoprotein, EWS-FLI1, in DNA repair is likely. Here, we demonstrate that EWS-FLI1 causes a defect in microhomology-mediated end-joining (MMEJ) repair. EWSR1 is a splicing factor that promotes the faithful splicing of the POLQ pre-mRNA, required for the expression of POLΘ, a critical protein in the MMEJ pathway. Expression of EWS-FLI1, or loss of EWSR1, causes exon 25 skipping of the POLQ transcript, decreased POLΘ expression, impaired MMEJ, and cellular sensitivity to inhibitors of the Fanconi Anemia (FA), NHEJ, or HR pathways, through the mechanism of synthetic lethality. Knockdown of EWS-FLI1 expression restores POLΘ mitotic foci and increases MMEJ activity. Inhibitors of the FA, NHEJ, or HR therefore may provide a targeted therapy for patients with EwS.

Alternative end-joining↗

RNA Pol I activity is required for meiotic chromatin organization and the H3K4me3 gradient essential for oogenesis, independent of ribosome synthesis.

Oogenesis requires extensive and dynamic chromatin remodeling that primes gene promoters for later transcriptional activation during embryonic development. Here, we uncover a pivotal, non-canonical role for RNA Polymerase I (Pol I) in driving these chromatin state transitions during Caenorhabditis elegans oogenesis. Using the auxin-inducible degron system to selectively deplete either Pol I catalytic subunits or ribosome assembly factors, we disentangle the consequences of impaired nucleolar integrity from reductions in ribosome biogenesis. Strikingly, although disrupting ribosome assembly caused minimal effects on oocyte production, loss of Pol I activity led to widespread changes in chromatin accessibility, a dampening of the distal-proximal H3K4me3 gradient required for oogenesis, reduced synapsis, and elevated ATM/ATR phosphorylation, resulting in fewer but significantly larger oocytes. Despite their promoters becoming more accessible, oogenesis genes did not show large changes in steady-state mRNA, consistent with transcriptional repression prior to fertilization. Instead, Pol I depletion prematurely remodeled oogenic chromatin, through a misdirection of H3K4me3 deposition towards promoters normally primed for zygotic genome activation. These findings reveal an epigenetic gating function for nucleolar integrity in oocyte maturation: Pol I preserves three-dimensional chromatin organization and maintains proper spatiotemporal regulation of histone modifications, independent of ribosome production. Given the evolutionary conservation of nucleolar dynamics and histone modifications during gametogenesis, our work suggests that nucleolar stress, whether from environmental factors, aging, or genetic disorders, could broadly compromise fertility by disrupting oogenic chromatin priming.

Journal Article↗

Evaluating selection at intermediate scales within genes provides robust identification of genes under positive selection in M. tuberculosis clinical isolates.

Multiple studies have reported genes in the M. tuberculosis (Mtb) genome that are under diversifying selection, based on genetic variants among Mtb clinical isolates. These might reflect adaptions to selection pressures associated with modern clinical treatment of TB. Many, but not all, of these genes under selection are related to drug resistance. Most of these studies have evaluated selection at the gene-level. However, positive selection can be evaluated on different scales, including individual sites (codons) and local regions within an ORF. In this paper, we use GenomegaMap, a Bayesian method for calculating selection, to evaluate selection of genes in the Mtb genome at all three levels. We present evidence that the intermediate analysis (windows of codons) yields the most credible list of candidate genes under selection (excluding PPE and PE_PGRS genes, which are predicted less reliably due to frequent sequencing errors). A further advantage of this approach is that it identifies specific regions within proteins that are under selective pressure, which is useful for structural and functional interpretation. In an analysis of two separate collections of Mtb clinical isolates (from Moldova; and a globally-representative set), we observed 53 and 173 significant genes under selection, with 36% overlap. The lists of genes under selection include many drug-resistance genes, as well as other genes that have previously been reported to be under selection (resR, phoR). The specific regions under selection identified within drug-resistance genes are shown to correspond to protein structural features known to be involved in resistance, supporting accuracy of the method. Positive selection in several ESX-1-related genes was also observed, suggesting adaptation to immune pressure.

adaptation↗

Efficient evidence-based genome annotation with EviAnn.

For many years, machine learning-based ab initio gene finding approaches have been central components of eukaryotic genome annotation pipelines, and they remain so today. The reliance on these approaches was originally sustained by the high cost and low availability of gene expression data, a primary source of evidence for gene annotation along with protein homology. However, innovations in modern sequencing technologies have revolutionized the acquisition of gene expression data, allowing scientists to rely more heavily on this class of evidence. In addition, proteins found in a multitude of well-annotated genomes represent another invaluable resource for gene annotation. Existing annotation packages often underutilize these data sources, which prompted us to develop EviAnn (Evidence-based Annotator), a novel evidence-based eukaryotic gene annotation system. EviAnn takes a strongly data-driven approach, building the exon-intron structure of genes from transcript alignments or protein-sequence homology rather than from purely ab initio gene finding techniques. We show that when provided with the same input data, EviAnn consistently outperforms current state-of-the-art packages including BRAKER3, MAKER2, and FINDER, while utilizing considerably less computer time. Annotation of a mammalian genome can be completed in less than an hour on a single multi-core server. EviAnn is freely available under an open-source license from https://github.com/alekseyzimin/EviAnn_release and from Bioconda as "eviann".

Journal Article↗

Genome-wide mapping of cyclic AMP receptor protein binding in Enteroaggregative Escherichia coli reveals targeting of virulence-associated genes.

Bacterial pathogens use a wide array of virulence factors to colonise and subsequently elicit disease in their host. These factors are often subject to extensive regulation at the transcriptional level, to ensure that their expression is timely. Although many pathogens use bespoke transcription factors that primarily target virulence genes, global transcription factors also sometimes play a role in controlling these genes. Enteroaggregative Escherichia coli (EAEC) is a significant cause of watery and mucoid diarrhoea globally. The organism colonises the small intestine before producing toxins that elicit disease, using a multitude of virulence factors that are encoded both chromosomally and on virulence plasmids. In this work, we have studied the cAMP Receptor Protein (CRP), a well-characterised bacterial global transcription factor, focusing on its role in pathogenicity of the prototype EAEC strain 042. We show that, although most functional CRP binding sites on the chromosome are conserved between E. coli K-12 and 042, CRP has been co-opted to couple the expression of some virulence genes to the nutritional state of the cell. We report novel mechanisms for CRP-dependent regulation of genes, whose products contribute to adhesion, production of a bacterial antibiotic, and export of a polysaccharide capsule.

CRP↗

Alzheimer's subtypes A supervised, unsupervised, multimodal, multilayered embedded recursive (SUMMER) AI study.

Since Alzheimer's disease (AD) is a heterogeneous disease, different subtypes may have distinct biological, genetic, and clinical characteristics, requiring tailored interventions. While several proposed subtypes of AD exist, there is still no clear consensus on a definitive classification. By leveraging complementary AI approaches, including supervised and unsupervised learning, within a recursive pipeline (SUMMER) that integrates multimodal datasets encompassing MRI measurements, phenotypes, and genetic data, our goal was to generate robust scientific evidence for identifying AD subtypes. Data was downloaded from the Alzheimer's Disease Neuroimaging Initiative (ADNI) database and included neuroimaging data (MRI), genetics (SNPs), clinical diagnosis, and demographics. 1133 European American participants' images, aged 55-95, were included in this study. The analysis was multi-fold, where the first step involved applying an unsupervised application to a subset of the MRI sample (AD + cognitively normal (CN) aged matched groups, 100 men aged 68-85 years, and 76 women aged 68-85 years). The MRI brain gray matter was segmented into 44 regions of interest (ROIs) according to a standard atlas, and 618 features were extracted, including ROI voxel intensity measurements such as minimum, maximum, and histogram variables. Results identified a cluster of subtype AD men and a cluster of subtype AD women that were distinct from the rest of their respective samples. In the next step, the integrity of the identified subtype AD clusters was investigated using the XGBoost supervised machine learning application with genetic features (SNPs, N=36,724) and labels: the identified subtype AD cluster vs. the rest of the sample, stratified by sex. A significant AD subtype men model (accuracy=0.85, F1=0.72, AUC=0.83) and a significant women AD subtype model (accuracy=0.81, F1=0.81, AUC=0.81) were built, confirming the homogeneity of the isolated AD subtype clusters. Discriminative biomarkers were extracted from the significant models, including selected ROIs and SNPs. Finally, the subtype models were tested on an unseen subset of ADNI data. The genetic-based models identified clusters of AD subtype participants consisting of 34% of the men AD group and 47% of the women AD group. Phenotypic analysis indicates that lower body weight was associated with the women's AD subtype. Complex diseases like AD demand a sophisticated, multimodal approach for precise diagnosis. Effectively identifying disease subtypes enhances the potential for personalized treatment, ultimately improving patient outcomes.

Journal Article↗

Mapping the transcriptional regulatory network of a fungal pathogen by exploiting transcription factor perturbation.

Cryptococcus neoformans is a deadly fungal pathogen. Upon entering a mammalian host, it deploys a voluminous polysaccharide capsule that is necessary for it to survive host defenses and maintain an infection. Capsule expansion is regulated transcriptionally, as deletion of many transcription factors (TFs) alters capsule. Thus, we set out to map the transcriptional regulatory network of C. neoformans - that is, to identify the TFs that directly regulate each gene in the genome. First, we carried out RNA-seq of 120 single-TF-deletion strains, together with wild-type controls. We then applied NetProphet3, a TF network mapping algorithm, to predict the direct functional targets of each TF. Unexpectedly, analysis of this network indicated that there are no TFs that primarily regulate genes involved in capsule formation. Rather, the TFs that play a role in deploying capsule also regulate many other genes and processes. Comparison to a TF network map we built for Saccharomyces cerevisiae, a distantly related model yeast, identified pairs of TFs that are functionally orthologous - that is, their targets are enriched for orthologous genes. In many cases, these pairs are different from the ones identified by sequence homology alone. We suggest that network analyses should be used to complement sequence comparison when searching for functionally orthologous transcription factors. Our network map can be searched and visualized at http://cryptococcus.net.

Journal Article↗

SLC35G3 is a UDP-N-acetylglucosamine transporter for sperm glycoprotein formation and underpins male fertility in mice.

Despite the recognized importance of glycans in biological phenomena, their complex roles in spermatogenesis and sperm function remain unclear. SLC35G3, a 10-transmembrane protein specifically found in early round spermatids, belongs to the sugar-nucleotide transporter family, indicating its involvement in glycan formation. In this study, we found that Slc35g3 knockout male mice were sterile due to impaired sperm functions in uterotubal junction passage, zona pellucida binding, and oocyte fusion. Mouse SLC35G3 has UDP-GlcNAc transporter activity, and its ablation caused abnormal processing of the sperm plasma membrane and acrosome membrane proteins. Reported human SLC35G3 mutations (F267L and T179HfsTer27) diminished the UDP-GlcNAc transporter activity of SLC35G3, implying infertility risks in males carrying these mutations. Our findings unveil the vital roles of SLC35G3 in the glycan formation of sperm membrane proteins critical for sperm fertilizing ability.

Biological Sciences↗

The distribution of fitness effects varies phylogenetically across animals.

The distribution of fitness effects (DFE) describes the selection coefficients () of newly arising mutations and fundamentally influences population genetic processes. However, the extent and mechanisms of DFE variation have not been systematically investigated across species with divergent phylogenetic histories and ecological functions. Here, we inferred the DFE in natural populations of eleven animal (sub)species, including humans, mice, fin whales, vaquitas, wolves, collared flycatchers, pied flycatchers, halictid bees, Drosophila, and mosquitoes. We find that the DFE co-varies with phylogeny, where the expected mutation effects are more similar in closely related species (). Additionally, mammals have a higher proportion of strongly deleterious mutations (22% to 47% in mammals; 0.0% to 5.4% in insects and birds) and a lower proportion of weakly deleterious mutations than insects and birds. Population size is significantly negatively correlated with the expected impact of new deleterious mutations (), and the proportion of new beneficial mutations (). These findings align with Fisher's Geometric Model (FGM), which defines organismal complexity as the number of phenotypes under selection. Consistent with the FGM's predictions, we observe that mutations are more deleterious in complex organisms, while beneficial mutations occur more frequently in smaller populations to compensate for the drift load. Our study demonstrates strong phylogenetic constraints in the evolution of a fundamental population genetics parameter, and proposes that, through mechanisms of global epistasis, long-term population size and organismal complexity drive variation in the DFE across animals.

Fisher’s geometric model↗