Mechanisms of change in gene copy number

Copy number variations affecting growth curve parameters in a crossbred chicken population

Copy number variations (CNVs) are key structural variations in the genome and may contribute to phenotypic differences. In this study, we used a F2 chicken population created from reciprocal crossing between fast-growing Arian broiler line and Urmia native chickens. The chickens were genotyped by 60 K SNP BeadChip, and PennCNV algorithm was used to detect genome-wide CNVs. The growth curve parameters of W0, k, L, Wf, Wi, ti and average GR were used as phenotypic data. The association between CNV and growth curve parameters was carried out using the CNVRanger R/Bioconductor package. Five CNV regions (CNVRs) were chosen for the validation experiment using qPCR. Gene enrichment analysis was done using WebGestalt. The STRING database was used to search for significant pathways. The results identified 966 CNVs and 600 CNVRs including 468 gains, 67 losses, and 65 both events on autosomal chromosomes. Validation of the CNVRs obtained from the qPCR assay were 79 % consistent with the prediction by PennCNV. A total of 43 significant CNVs were obtained for the seven growth curve parameters. The 416 genes annotated for significant CNVs. Six genes out of 416 genes were most related to growth curve parameters. These genes were LCP2, Dock2, CD80, CYFIP1, NIPA1 and NIPA2. Some of these genes in their biological process were associated with the growth, reproduction and development of cells or organs that ultimately lead to the growth of the body. The results of the study could pave the way for better understanding the molecular process of CNVs and growth curve parameters in birds.

Identification of novel compound heterozygous variants of the ALMS1 gene in a child with Alström syndrome by whole genome sequencing

BACKGROUND

Alström syndrome (ALMS), a rare recessively inherited ciliopathy caused by mutations in ALMS1, is characterized by retinal dystrophy, childhood obesity, sensorineural hearing loss, and type 2 diabetes mellitus. The majority of pathogenic variants in ALMS1 are nonsense and frameshift mutations, which would lead to premature protein truncation, whereas copy number variants are seldom reported.

METHODS

Herein, we present a 10-year-old Chinese girl with ALMS. The potential causative genetic variant was confirmed through whole genome sequencing, quantitative real-time PCR analysis, and Sanger sequencing. Additionally, breakpoint analysis was performed to determine the exact breakpoint site of the large deletion and elucidate its probable formation mechanism.

RESULTS

The patient had a cor triatriatum sinister (CTS) structure. Genetic analysis identified novel compound heterozygous variants in the patient, consisting of a frameshift variant c.4414_4415delGT (p.V1472Nfs*26) in ALMS1 and a novel large deletion at chr2:73,612,355-73,626,339, which encompasses exon 1 of the ALMS1 gene. Moreover, breakpoint analysis revealed that the large deletion probably formed through the microhomology-mediated end joining (MMEJ) mechanism due to the 6-bp microhomologies (TCCTTC) observed at both ends of the breakpoints.

CONCLUSIONS

In this study, novel compound heterozygous variants in the ALMS1 gene were identified in an ALMS patient with a CTS structure. The molecular confirmation of these variants expands the mutational spectrum of ALMS1, while the manifestation of ALMS in the patient provides additional clinical insights into this syndrome.

Modelling DNA damage-repair and beyond

The paper presents a review of mechanistic modelling studies of DNA damage and DNA repair, and consequences to follow in mammalian cell nucleus. We hypothesize DNA deletions are consequences of repair of double strand breaks leading to the modifications of genome that play crucial role in long term development of genetic inheritance and diseases. The aim of the paper is to review formation mechanisms underlying naturally occurring DNA deletions in the human genome and their potential relevance for bridging the gap between induced DNA double strand breaks and deletions in damaged human genome from endogenous and exogenous events. The model of the cell nucleus presented enables simulation of DNA damage at molecular level identifying the spectrum of damage induced in all chromosomal territories and loops. Our mechanistic modelling of DNA repair for double stand breaks (DSB), single strand breaks (SSB) and base damage (BD), shows the complexity of DNA damage is responsible for the longer repair times and the reason for the biphasic feature of mammalian cells repair curves. In the absence of experimentally determined data, the mechanistic model of repair predicts the in vivo rate constants for the proteins involved in the repair of DSB, SSB, and of BD.

Genome-wide association study of copy number variations with shank traits in a F2 crossbred chicken population

Copy number variations (CNVs) are large-scale changes in the DNA sequence that can affect the genetic structure and phenotype of an organism. The purpose of this study was to investigate the existing CNVs and their associations with the shank diameter (ShD) and shank length (ShL) traits using data from an F2 crossbred chicken population. To carry out the study, 312 chickens were genotyped using the Illumina 60k SNP Beadchip. The shank traits of the birds were measured from day 1 to 12 weeks of age. penncnv and cnvruler tools were used to find copy numbers and regions with copy number changes (CNVR), respectively. The CNVRanger package was used to perform a genome-wide association study between shank traits and CNVs. Gene ontology research in CNVRs was carried out using the david database. In this investigation, 966 CNVs and 606 regions with copy number changes were discovered. The copy number states and variations were randomly distributed along the length of the autosomal chromosomes. Weeks 1-4, 9 and 12 of growth revealed a significant association of copy number variations with shank traits, false discovery rate (FDR-corrected p-value < 0.01), and the majority of CNVs that were statistically significant were found on chromosomes 1-3. These CNV segments are nearby genes such as KCNJ12, FGF6 and MYF5, which are fundamental to growth and development. In addition, gene set analyses revealed terms related to muscle physiology, regulation of cellular processes and potassium channels.

afila, the origin and nature of a major innovation in the history of pea breeding

The afila (af) mutation causes the replacement of leaflets by a branched mass of tendrils in the compound leaves of pea - Pisum sativum L. This mutation was first described in 1953, and several reports of spontaneous af mutations and induced mutants with a similar phenotype exist. Despite widespread introgression into breeding material, the nature of af and the origin of the alleles used remain unknown. Here, we combine comparative genomics with reverse genetic approaches to elucidate the genetic determinants of af. We also investigate haplotype diversity using a set of AfAf and afaf cultivars and breeding lines and molecular markers linked to seven consecutive genes. Our results show that deletion of two tandemly arranged genes encoding Q-type Cys(2)His(2) zinc finger transcription factors, PsPALM1a and PsPALM1b, is responsible for the af phenotype in pea. Eight haplotypes were identified in the af-harbouring genomic region on chromosome 2. These haplotypes differ in the size of the deletion, covering more or less genes. Diversity at the af locus is valuable for crop improvement and sheds light on the history of pea breeding for improved standing ability. The results will be used to understand the function of PsPALM1a/b and to transfer the knowledge for innovation in related crops.

Prevalent fast evolution of genes involved in heterochromatin functions

Heterochromatin is a gene-poor and repeat-rich genomic compartment universally found in eukaryotes. Despite its low transcriptional activity, heterochromatin plays important roles in maintaining genome stability, organizing chromosomes, and suppressing transposable elements (TEs). Given the importance of these functions, it is expected that the 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 (IDRs), purifying selection may have maintained the proportions of IDRs of these proteins. Together with the observed negative associations between evolutionary rates of these genes and genomic TE abundance, 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 TEs 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.

A comprehensive review of deep learning-based variant calling methods

Genome sequencing data have become increasingly important in the field of personalized medicine and diagnosis. However, accurately detecting genomic variations remains a challenging task. Traditional variation detection methods rely on manual inspection or predefined rules, which can be time-consuming and prone to errors. Consequently, deep learning-based approaches for variation detection have gained attention due to their ability to automatically learn genomic features that distinguish between variants. In our review, we discuss the recent advancements in deep learning-based algorithms for detecting small variations and structural variations in genomic data, as well as their advantages and limitations.

Large TRAPPC11 gene deletions as a cause of muscular dystrophy and their estimated genesis

BACKGROUND

Transport protein particle (TRAPP) is a multiprotein complex that functions in localising proteins to the Golgi compartment. The TRAPPC11 subunit has been implicated in diseases affecting muscle, brain, eye and to some extent liver. We present three patients who are compound heterozygotes for a missense variant and a structural variant in the TRAPPC11 gene. TRAPPC11 structural variants have not yet been described in association with a disease. In order to reveal the estimated genesis of identified structural variants, we performed sequencing of individual breakpoint junctions and analysed the extent of homology and the presence of repetitive elements in and around the breakpoints.

METHODS

Biochemical methods including isoelectric focusing on serum transferrin and apolipoprotein C-III, as well as mitochondrial respiratory chain complex activity measurements, were used. Muscle biopsy samples underwent histochemical analysis. Next-generation sequencing was employed for identifying sequence variants associated with neuromuscular disorders, and Sanger sequencing was used to confirm findings.

RESULTS

We suppose that non-homologous end joining is a possible mechanism of deletion origin in two patients and non-allelic homologous recombination in one patient. Analyses of mitochondrial function performed in patients' skeletal muscles revealed an imbalance of mitochondrial metabolism, which worsens with age and disease progression.

CONCLUSION

Our results contribute to further knowledge in the field of neuromuscular diseases and mutational mechanisms. This knowledge is important for understanding the molecular nature of human diseases and allows us to improve strategies for identifying disease-causing mutations.

Copy number variation introduced by a massive mobile element facilitates global thermal adaptation in a fungal wheat pathogen

Copy number variation (CNV) can drive rapid evolution in changing environments. In microbial pathogens, such adaptation is a key factor underpinning epidemics and colonization of new niches. However, the genomic determinants of such adaptation remain poorly understood. Here, we systematically investigate CNVs in a large genome sequencing dataset spanning a worldwide collection of 1104 genomes from the major wheat pathogen Zymoseptoria tritici. We found overall strong purifying selection acting on most CNVs. Genomic defense mechanisms likely accelerated gene loss over episodes of continental colonization. Local adaptation along climatic gradients was likely facilitated by CNVs affecting secondary metabolite production and gene loss in general. One of the strongest loci for climatic adaptation is a highly conserved gene of the NAD-dependent Sirtuin family. The Sirtuin CNV locus localizes to an ~68-kb Starship mobile element unique to the species carrying genes highly expressed during plant infection. The element has likely lost the ability to transpose, demonstrating how the ongoing domestication of cargo-carrying selfish elements can contribute to selectable variation within populations. Our work highlights how standing variation in gene copy numbers at the global scale can be a major factor driving climatic and metabolic adaptation in microbial species.

Low-input PacBio sequencing generates high-quality individual fly genomes and characterizes mutational processes

Long-read sequencing, exemplified by PacBio, revolutionizes genomics, overcoming challenges like repetitive sequences. However, the high DNA requirement ( > 1 µg) is prohibitive for small organisms. We develop a low-input (100 ng), low-cost, and amplification-free library-generation method for PacBio sequencing (LILAP) using Tn5-based tagmentation and DNA circularization within one tube. We test LILAP with two Drosophila melanogaster individuals, and generate near-complete genomes, surpassing preexisting single-fly genomes. By analyzing variations in these two genomes, we characterize mutational processes: complex transpositions (transposon insertions together with extra duplications and/or deletions) prefer regions characterized by non-B DNA structures, and gene conversion of transposons occurs on both DNA and RNA levels. Concurrently, we generate two complete assemblies for the endosymbiotic bacterium Wolbachia in these flies and similarly detect transposon conversion. Thus, LILAP promises a broad PacBio sequencing adoption for not only mutational studies of flies and their symbionts but also explorations of other small organisms or precious samples.

Regulators of rDNA array morphology in fission yeast

Nucleolar morphology is a well-established indicator of ribosome biogenesis activity that has served as the foundation of many screens investigating ribosome production. Missing from this field of study is a broad-scale investigation of the regulation of ribosomal DNA morphology, despite the essential role of rRNA gene transcription in modulating ribosome output. We hypothesized that the morphology of rDNA arrays reflects ribosome biogenesis activity. We established GapR-GFP, a prokaryotic DNA-binding protein that recognizes transcriptionally-induced overtwisted DNA, as a live visual fluorescent marker for quantitative analysis of rDNA organization in Schizosaccharomyces pombe. We found that the morphology-which we refer to as spatial organization-of the rDNA arrays is dynamic throughout the cell cycle, under glucose starvation, RNA pol I inhibition, and TOR activation. Screening the haploid S. pombe Bioneer deletion collection for spatial organization phenotypes revealed large ribosomal protein (RPL) gene deletions that alter rDNA organization. Further work revealed RPL gene deletion mutants with altered rDNA organization also demonstrate resistance to the TOR inhibitor Torin1. A genetic analysis of signaling pathways essential for this resistance phenotype implicated many factors including a conserved MAPK, Pmk1, previously linked to extracellular stress responses. We propose RPL gene deletion triggers altered rDNA morphology due to compensatory changes in ribosome biogenesis via multiple signaling pathways, and we further suggest compensatory responses may contribute to human diseases such as ribosomopathies. Altogether, GapR-GFP is a powerful tool for live visual reporting on rDNA morphology under myriad conditions.

The black honey bee genome: insights on specific structural elements and a first step towards pangenomes

BACKGROUND

The honey bee reference genome, HAv3.1, was produced from a commercial line sample that was thought to have a largely dominant Apis mellifera ligustica genetic background. Apis mellifera mellifera, often referred to as the black bee, has a separate evolutionary history and is the original type in western and northern Europe. Growing interest in this subspecies for conservation and non-professional apicultural practices, together with the necessity of deciphering genome backgrounds in hybrids, triggered the necessity for a specific genome assembly. Moreover, having several high-quality genomes is becoming key for taking structural variations into account in pangenome analyses.

RESULTS

Pacific Bioscience technology long reads were produced from a single haploid black bee drone. Scaffolding contigs into chromosomes was done using a high-density genetic map. This allowed for re-estimation of the recombination rate, which was over-estimated in some previous studies due to mis-assemblies, which resulted in spurious inversions in the older reference genomes. The sequence continuity obtained was very high and the only limit towards continuous chromosome-wide sequences seemed to be due to tandem repeat arrays that were usually longer than 10 kb and that belonged to two main families, the 371 and 91 bp repeats, causing problems in the assembly process due to high internal sequence similarity. Our assembly was used together with the reference genome to genotype two structural variants by a pangenome graph approach with Graphtyper2. Genotypes obtained were either correct or missing, when compared to an approach based on sequencing depth analysis, and genotyping rates were 89 and 76% for the two variants.

CONCLUSIONS

Our new assembly for the Apis mellifera mellifera honey bee subspecies demonstrates the utility of multiple high-quality genomes for the genotyping of structural variants, with a test case on two insertions and deletions. It will therefore be an invaluable resource for future studies, for instance by including structural variants in GWAS. Having used a single haploid drone for sequencing allowed a refined analysis of very large tandem repeat arrays, raising the question of their function in the genome. High quality genome assemblies for multiple subspecies such as presented here, are crucial for emerging projects using pangenomes.

HBV integrations reshaping genomic structures promote hepatocellular carcinoma

OBJECTIVE

Hepatitis B virus (HBV)-related hepatocellular carcinoma (HCC), mostly characterised by HBV integrations, is prevalent worldwide. Previous HBV studies mainly focused on a few hotspot integrations. However, the oncogenic role of the other HBV integrations remains unclear. This study aimed to elucidate HBV integration-induced tumourigenesis further.

DESIGN

Here, we illuminated the genomic structures encompassing HBV integrations in 124 HCCs across ages using whole genome sequencing and Nanopore long reads. We classified a repertoire of integration patterns featured by complex genomic rearrangement. We also conducted a clustered regularly interspaced short palindromic repeat (CRISPR)-based gain-of-function genetic screen in mouse hepatocytes. We individually activated each candidate gene in the mouse model to uncover HBV integration-mediated oncogenic aberration that elicits tumourigenesis in mice.

RESULTS

These HBV-mediated rearrangements are significantly enriched in a bridge-fusion-bridge pattern and interchromosomal translocations, and frequently led to a wide range of aberrations including driver copy number variations in chr 4q, 5p (TERT), 6q, 8p, 16q, 9p (CDKN2A/B), 17p (TP53) and 13q (RB1), and particularly, ultra-early amplifications in chr8q. Integrated HBV frequently contains complex structures correlated with the translocation distance. Paired breakpoints within each integration event usually exhibit different microhomology, likely mediated by different DNA repair mechanisms. HBV-mediated rearrangements significantly correlated with young age, higher HBV DNA level and TP53 mutations but were less prevalent in the patients subjected to prior antiviral therapies. Finally, we recapitulated the TONSL and TMEM65 amplification in chr8q led by HBV integration using CRISPR/Cas9 editing and demonstrated their tumourigenic potentials.

CONCLUSION

HBV integrations extensively reshape genomic structures and promote hepatocarcinogenesis (graphical abstract), which may occur early in a patient's life.

Analysis of copy number variants detected by sequencing in spontaneous abortion

BACKGROUND

The incidence of spontaneous abortion (SA), which affects approximately 15-20% of pregnancies, is the most common complication of early pregnancy. Pathogenic copy number variations (CNVs) are recognized as potential genetic causes of SA. However, CNVs of variants of uncertain significance (VOUS) have been identified in products of conceptions (POCs), and their correlation with SA remains uncertain.

RESULTS

Of 189 spontaneous abortion cases, trisomy 16 was the most common numerical chromosome abnormality, followed by monosomy X. CNVs most often occurred on chromosomes 4 and 8. Gene Ontology and signaling pathway analysis revealed significant enrichment of genes related to nervous system development, transmembrane transport, cell adhesion, and structural components of chromatin. Furthermore, genes within the VOUS CNVs were screened by integrating human placental expression profiles, PhyloP scores, and Residual Variance Intolerance Score (RVIS) percentiles to identify potential candidate genes associated with spontaneous abortion. Fourteen potential candidate genes (LZTR1, TSHZ1, AMIGO2, H1-4, H2BC4, H2AC7, H3C8, H4C3, H3C6, PHKG2, PRR14, RNF40, SRCAP, ZNF629) were identified. Variations in LZTR1, TSHZ1, and H4C3 may contribute to embryonic lethality.

CONCLUSIONS

CNV sequencing (CNV-seq) analysis is an effective technique for detecting chromosomal abnormalities in POCs and identifying potential candidate genes for SA.

Genome integrity as a potential index of longevity in Ashkenazi Centenarian's families

The aging process, or senescence, is characterized by age-specific decline in physical and physiological function, and increased frailty and genomic changes, including mutation accumulation. However, the mechanisms through which changes in genomic architecture influence human longevity have remained obscure. Copy number variants (CNVs), an abundant class of genomic variants, offer unique opportunities for understanding age-related genomic changes. Here we report the spectrum of CNVs in a cohort of 670 Ashkenazi Jewish centenarians, their progeny, and unrelated controls. The average ages of these groups were 97.4 ± 2.8, 69.2 ± 9.2, and 66.5 ± 7.0 respectively. For the first time, we compared different size classes of CNVs, from 1 kB to 100 MB in size. Using a high-resolution custom Affymetrix array, targeting 44,639 genomic regions, we identified a total of 12,166, 22,188, and 10,285 CNVs in centenarians, their progeny, and control groups, respectively. Interestingly, the offspring group showed the highest number of unique CNVs, followed by control and centenarians. While both gains and losses were found in all three groups, centenarians showed a significantly higher average number of both total gains and losses relative to their controls (p < 0.0327, 0.0182, respectively). Moreover, centenarians showed a lower total length of genomic material lost, suggesting that they may maintain superior genomic integrity over time. We also observe a significance fold increase of CNVs among the offspring, implying greater genomic integrity and a putative mechanism for longevity preservation. Genomic regions that experienced loss or gains appear to be distributed across many sites in the genome and contain genes involved in DNA transcription, cellular transport, developmental pathways, and metabolic functions. Our findings suggest that the exceptional longevity observed in centenarians may be attributed to the prolonged maintenance of functionally important genes. These genes are intrinsic to specific genomic regions as well as to the overall integrity of the genomic architecture. Additionally, a strong association between longer CNVs and differential gene expression observed in this study supports the notion that genomic integrity could positively influence longevity.

Evaluation of 100 Dutch cases with 16p11.2 deletion and duplication syndromes; from clinical manifestations towards personalized treatment options

The 16p11.2 deletion syndrome is a clinically heterogeneous disorder, characterized by developmental delay, intellectual disability, hyperphagia, obesity, macrocephaly and psychiatric problems. Cases with 16p11.2 duplication syndrome have similar neurodevelopmental problems, but typically show a partial 'mirror phenotype' with underweight and microcephaly. Various copy number variants (CNVs) of the chromosomal 16p11.2 region have been described. Most is known about the 'typical' 16p11.2 BP4-BP5 (29.6-30.2 Mb; ~600 kb) deletions and duplications, but there are also several published cohorts with more distal 16p11.2 BP2-BP3 CNVs (28.8-29.0 Mb; ~220 kb), who exhibit clinical overlap. We assessed 100 cases with various pathogenic 16p11.2 CNVs and compared their clinical characteristics to provide more clear genotype-phenotype correlations and raise awareness of the different 16p11.2 CNVs. Neurodevelopmental and weight issues were reported in the majority of cases. Cases with distal 16p11.2 BP2-BP3 deletion showed the most severe obesity phenotype (73.7% obesity, mean BMI SDS 3.2). In addition to the more well defined typical 16p11.2 BP4-BP5 and distal 16p11.2 BP2-BP3 CNVs, we describe the clinical features of five cases with other, overlapping, 16p11.2 CNVs in more detail. Interestingly, four cases had a second genetic diagnosis and 18 cases an additional gene variant of uncertain significance, that could potentially help explain the cases' phenotypes. In conclusion, we provide an overview of our Dutch cohort of cases with various pathogenic 16p11.2 CNVs and relevant second genetic findings, that can aid in adequately recognizing, diagnosing and counseling of individuals with 16p11.2 CNVs, and describe the personalized medicine for cases with these conditions.

White horses - non-coding sequences drive premature hair greying and predisposition to melanoma

The Grey allele in horses is causing premature hair greying and susceptibility to melanoma. The causal mutation is a 4.6 kb tandem duplication in intron 6 of the Syntaxin 17 gene. A recent study demonstrated that the most common allele at the Grey locus (G3) involves three tandem copies of this sequence, whilst a more rare allele (G2) has two tandem copies and the wild-type allele (G1) only one copy. The G3 allele is causing fast greying and high incidence of skin melanoma, whereas the G2 allele is causing slow greying and no obvious increase in melanoma incidence. Further somatic copy number expansion has been documented in melanoma tissue from Grey horses. Functional studies showed that this intronic sequence acts as a weak melanocyte-specific enhancer that becomes substantially stronger by the copy number expansion. The Grey mutation is associated with upregulated expression of both Syntaxin 17 and the neighbouring NR4A3 gene in Grey horse melanomas. It is still an open question which of these genes is most important for the phenotypic effects or if causality is due to the combined effect of upregulation of both genes. Interestingly, RNAseq data in the Human Protein Atlas give support for a possible role of NR4A3 because it is particularly upregulated in human skin cancer, and it belongs to a cluster of genes associated with skin cancer and melanin biosynthesis. The Grey mutation and its association with melanoma provide a possibility to study the path to tumour development in numerous Grey horses carrying exactly the same predisposing mutation.

TMBstable: a variant caller controls performance variation across heterogeneous sequencing samples

In cancer genomics, variant calling has advanced, but traditional mean accuracy evaluations are inadequate for biomarkers like tumor mutation burden, which vary significantly across samples, affecting immunotherapy patient selection and threshold settings. In this study, we introduce TMBstable, an innovative method that dynamically selects optimal variant calling strategies for specific genomic regions using a meta-learning framework, distinguishing it from traditional callers with uniform sample-wide strategies. The process begins with segmenting the sample into windows and extracting meta-features for clustering, followed by using a pre-trained meta-model to select suitable algorithms for each cluster, thereby addressing strategy-sample mismatches, reducing performance fluctuations and ensuring consistent performance across various samples. We evaluated TMBstable using both simulated and real non-small cell lung cancer and nasopharyngeal carcinoma samples, comparing it with advanced callers. The assessment, focusing on stability measures, such as the variance and coefficient of variation in false positive rate, false negative rate, precision and recall, involved 300 simulated and 106 real tumor samples. Benchmark results showed TMBstable's superior stability with the lowest variance and coefficient of variation across performance metrics, highlighting its effectiveness in analyzing the counting-based biomarker. The TMBstable algorithm can be accessed at https://github.com/hello-json/TMBstable for academic usage only.

Genomic and metabolic instability during long-term fermentation of an industrial Saccharomyces cerevisiae strain engineered for C5 sugar utilization

The genetic stability and metabolic robustness of production strains is one of the key criteria for the production of bio-based products by microbial fermentation on an industrial scale. These criteria were here explored in an industrial ethanol-producer strain of Saccharomyces cerevisiae able to co-ferment D-xylose and L-arabinose with glucose through the chromosomal integration of several copies of pivotal genes for the use of these pentose (C5) sugars. Using batch sequential cultures in a controlled bioreactor that mimics long-term fermentation in an industrial setting, this strain was found to exhibit significant fluctuations in D-xylose and L-arabinose consumption as early as the 50th generation and beyond. These fluctuations seem not related to the few low-consumption C5 sugar clones that appeared throughout the sequential batch cultures at a frequency lower than 1.5% and that were due to the reduction in the number of copies of transgenes coding for C5 sugar assimilation enzymes. Also, subpopulations enriched with low or high RAD52 expression, whose expression level was reported to be proportional to homologous recombination rate did not exhibit defect in C5-sugar assimilation, arguing that other mechanisms may be responsible for copy number variation of transgenes. Overall, this work highlighted the existence of genetic and metabolic instabilities in an industrial yeast which, although modest in our conditions, could be more deleterious in harsher industrial conditions, leading to reduced production performance.

Large-scale duplication events underpin population-level flexibility in tRNA gene copy number in Pseudomonas fluorescens SBW25

The complement of tRNA genes within a genome is typically considered to be a (relatively) stable characteristic of an organism. Here, we demonstrate that bacterial tRNA gene set composition can be more flexible than previously appreciated, particularly regarding tRNA gene copy number. We report the high-rate occurrence of spontaneous, large-scale, tandem duplication events in laboratory populations of the bacterium Pseudomonas fluorescens SBW25. The identified duplications are up to ∼1 Mb in size (∼15% of the wildtype genome) and are predicted to change the copy number of up to 917 genes, including several tRNA genes. The observed duplications are inherently unstable: they occur, and are subsequently lost, at extremely high rates. We propose that this unusually plastic type of mutation provides a mechanism by which tRNA gene set diversity can be rapidly generated, while simultaneously preserving the underlying tRNA gene set in the absence of continued selection. That is, if a tRNA set variant provides no fitness advantage, then high-rate segregation of the duplication ensures the maintenance of the original tRNA gene set. However, if a tRNA gene set variant is beneficial, the underlying duplication fragment(s) may persist for longer and provide raw material for further, more stable, evolutionary change.

Evolving copy number gains promote tumor expansion and bolster mutational diversification

The timing and fitness effect of somatic copy number alterations (SCNA) in cancer evolution remains poorly understood. Here we present a framework to determine the timing of a clonal SCNA that encompasses multiple gains. This involves calculating the proportion of time from its last gain to the onset of population expansion (lead time) as well as the proportion of time prior to its first gain (initiation time). Our method capitalizes on the observation that a genomic segment, while in a specific copy number (CN) state, accumulates point mutations proportionally to its CN. Analyzing 184 whole genome sequenced samples from 75 patients across five tumor types, we commonly observe late gains following early initiating events, occurring just before the clonal expansion relevant to the sampling. These include gains acquired after genome doubling in more than 60% of cases. Notably, mathematical modeling suggests that late clonal gains may contain final-expansion drivers. Lastly, SCNAs bolster mutational diversification between subpopulations, exacerbating the circle of proliferation and increasing heterogeneity.

The Complex Landscape of Structural Divergence Between the Drosophila pseudoobscura and D. persimilis Genomes

Structural genomic variants are key drivers of phenotypic evolution. They can span hundreds to millions of base pairs and can thus affect large numbers of genetic elements. Although structural variation is quite common within and between species, its characterization depends upon the quality of genome assemblies and the proportion of repetitive elements. Using new high-quality genome assemblies, we report a complex and previously hidden landscape of structural divergence between the genomes of Drosophila persimilis and D. pseudoobscura, two classic species in speciation research, and study the relationships among structural variants, transposable elements, and gene expression divergence. The new assemblies confirm the already known fixed inversion differences between these species. Consistent with previous studies showing higher levels of nucleotide divergence between fixed inversions relative to collinear regions of the genome, we also find a significant overrepresentation of INDELs inside the inversions. We find that transposable elements accumulate in regions with low levels of recombination, and spatial correlation analyses reveal a strong association between transposable elements and structural variants. We also report a strong association between differentially expressed (DE) genes and structural variants and an overrepresentation of DE genes inside the fixed chromosomal inversions that separate this species pair. Interestingly, species-specific structural variants are overrepresented in DE genes involved in neural development, spermatogenesis, and oocyte-to-embryo transition. Overall, our results highlight the association of transposable elements with structural variants and their importance in driving evolutionary divergence.

A view of the pan-genome of domesticated Cowpea (Vigna unguiculata [L.] Walp.)

Cowpea, Vigna unguiculata L. Walp., is a diploid warm-season legume of critical importance as both food and fodder in sub-Saharan Africa. This species is also grown in Northern Africa, Europe, Latin America, North America, and East to Southeast Asia. To capture the genomic diversity of domesticates of this important legume, de novo genome assemblies were produced for representatives of six subpopulations of cultivated cowpea identified previously from genotyping of several hundred diverse accessions. In the most complete assembly (IT97K-499-35), 26,026 core and 4963 noncore genes were identified, with 35,436 pan genes when considering all seven accessions. GO terms associated with response to stress and defense response were highly enriched among the noncore genes, while core genes were enriched in terms related to transcription factor activity, and transport and metabolic processes. Over 5 million single nucleotide polymorphisms (SNPs) relative to each assembly and over 40 structural variants >1 Mb in size were identified by comparing genomes. Vu10 was the chromosome with the highest frequency of SNPs, and Vu04 had the most structural variants. Noncore genes harbor a larger proportion of potentially disruptive variants than core genes, including missense, stop gain, and frameshift mutations; this suggests that noncore genes substantially contribute to diversity within domesticated cowpea.

Inherent genome instability underlies trisomy 21-associated myeloid malignancies

Constitutional trisomy 21 (T21) is a state of aneuploidy associated with high incidence of childhood acute myeloid leukemia (AML). T21-associated AML is preceded by transient abnormal myelopoiesis (TAM), which is triggered by truncating mutations in GATA1 generating a short GATA1 isoform (GATA1s). T21-associated AML emerges due to secondary mutations in hematopoietic clones bearing GATA1s. Since aneuploidy generally impairs cellular fitness, the paradoxically elevated risk of myeloid malignancy in T21 is not fully understood. We hypothesized that individuals with T21 bear inherent genome instability in hematopoietic lineages that promotes leukemogenic mutations driving the genesis of TAM and AML. We found that individuals with T21 show increased chromosomal copy number variations (CNVs) compared to euploid individuals, suggesting that genome instability could be underlying predisposition to TAM and AML. Acquisition of GATA1s enforces myeloid skewing and maintenance of the hematopoietic progenitor state independently of T21; however, GATA1s in T21 hematopoietic progenitor cells (HPCs) further augments genome instability. Increased dosage of the chromosome 21 (chr21) gene DYRK1A impairs homology-directed DNA repair as a mechanism of elevated mutagenesis. These results posit a model wherein inherent genome instability in T21 drives myeloid malignancy in concert with GATA1s mutations.

Scrambling the genome in cancer: causes and consequences of complex chromosome rearrangements

Complex chromosome rearrangements, known as chromoanagenesis, are widespread in cancer. Based on large-scale DNA sequencing of human tumours, the most frequent type of complex chromosome rearrangement is chromothripsis, a massive, localized and clustered rearrangement of one (or a few) chromosomes seemingly acquired in a single event. Chromothripsis can be initiated by mitotic errors that produce a micronucleus encapsulating a single chromosome or chromosomal fragment. Rupture of the unstable micronuclear envelope exposes its chromatin to cytosolic nucleases and induces chromothriptic shattering. Found in up to half of tumours included in pan-cancer genomic analyses, chromothriptic rearrangements can contribute to tumorigenesis through inactivation of tumour suppressor genes, activation of proto-oncogenes, or gene amplification through the production of self-propagating extrachromosomal circular DNAs encoding oncogenes or genes conferring anticancer drug resistance. Here, we discuss what has been learned about the mechanisms that enable these complex genomic rearrangements and their consequences in cancer.

Characterizing PALB2 intragenic duplication breakpoints in a triple-negative breast cancer case using long-read sequencing

Introduction

Accurate identification and characterization of Large Genomic Rearrangements (LGR), especially duplications, are crucial for precise diagnosis and risk assessment. In this report, we characterized an intragenic duplication breakpoint of PALB2 to determine its pathogenicity significance.

Methods

A 52-year-old female with triple-negative breast cancer was diagnosed with a novel PALB2 LGR. An efficient and accurate methodology was applied, combining long-read sequencing and transcript analysis for the rapid characterization of the duplication.

Results

Duplication of exons 5 and 6 of PALB2 was validated by transcript analysis. Long-read sequencing enabled the localization of breakpoints within Alu elements, providing insights into the mechanism of duplication via non-allelic homologous recombination.

Conclusion

Using our combined methodology, we reclassified the PALB2 duplication as a pathogenic variant. This reclassification suggests a possible causative link between this specific genetic alteration and the aggressive phenotype of the patient.

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

BACKGROUND

During evolution, genes can experience duplications, losses, inversions and gene conversions. Why certain genes are more dynamic than others is poorly understood. Here we examine how several Sgs genes encoding glue proteins, which make up a bioadhesive that sticks the animal during metamorphosis, have evolved in Drosophila species.

RESULTS

We examined high-quality genome assemblies of 24 Drosophila species to study the evolutionary dynamics of four glue genes that are present in D. melanogaster and are part of the same gene family - Sgs1, Sgs3, Sgs7 and Sgs8 - across approximately 30 millions of years. We annotated a total of 102 Sgs genes and grouped them into 4 subfamilies. We present here a new nomenclature for these Sgs genes based on protein sequence conservation, genomic location and presence/absence of internal repeats. Two types of glue genes were uncovered. The first category (Sgs1, Sgs3x, Sgs3e) showed a few gene losses but no duplication, no local inversion and no gene conversion. The second group (Sgs3b, Sgs7, Sgs8) exhibited multiple events of gene losses, gene duplications, local inversions and gene conversions. Our data suggest that the presence of short "new glue" genes near the genes of the latter group may have accelerated their dynamics.

CONCLUSIONS

Our comparative analysis suggests that the evolutionary dynamics of glue genes is influenced by genomic context. Our molecular, phylogenetic and comparative analysis of the four glue genes Sgs1, Sgs3, Sgs7 and Sgs8 provides the foundation for investigating the role of the various glue genes during Drosophila life.

Acquired copy number variation in prostate tumours: a review of common somatic copy number alterations, how they are formed and their clinical utility

Prostate cancer is one of the most commonly diagnosed cancers in men and unfortunately, disease will progress in up to a third of patients despite primary treatment. Currently, there is a significant lack of prognostic tests that accurately predict disease course; however, the acquisition of somatic chromosomal variation in the form of DNA copy number variants may help understand disease progression. Notably, studies have found that a higher burden of somatic copy number alterations (SCNA) correlates with more aggressive disease, recurrence after surgery and metastasis. Here we will review the literature surrounding SCNA formation, including the roles of key tumour suppressors and oncogenes (PTEN, BRCA2, NKX3.1, ERG and AR), and their potential to inform diagnostic and prognostic clinical testing to improve predictive value. Ultimately, SCNAs, or inherited germline alterations that predispose to SCNAs, could have significant clinical utility in diagnostic and prognostic tests, in addition to guiding therapeutic selection.

Conumee 2.0: enhanced copy-number variation analysis from DNA methylation arrays for humans and mice

MOTIVATION

Copy-number variations (CNVs) are common genetic alterations in cancer and their detection may impact tumor classification and therapeutic decisions. However, detection of clinically relevant large and focal CNVs remains challenging when sample material or resources are limited. This has motivated us to create a software tool to infer CNVs from DNA methylation arrays which are often generated as part of clinical routines and in research settings.

RESULTS

We present our R package, conumee 2.0, that combines tangent normalization, an adjustable genomic binning heuristic, and weighted circular binary segmentation to utilize DNA methylation arrays for CNV analysis and mitigate technical biases and batch effects. Segmentation results were validated in a lung squamous cell carcinoma dataset from TCGA (n = 367 samples) by comparison to segmentations derived from genotyping arrays (Pearson's correlation coefficient of 0.91). We further introduce a segmented block bootstrapping approach to detect focal alternations that achieved 60.9% sensitivity and 98.6% specificity for deletions affecting CDKN2A/B (60.0% and 96.9% for RB1, respectively) in a low-grade glioma cohort from TCGA (n = 239 samples). Finally, our tool provides functionality to detect and summarize CNVs across large sample cohorts.

AVAILABILITY AND IMPLEMENTATION

Conumee 2.0 is available under open-source license at: https://github.com/hovestadtlab/conumee2.

Atypical cortical networks in children at high-genetic risk of psychiatric and neurodevelopmental disorders

Although many genetic risk factors for psychiatric and neurodevelopmental disorders have been identified, the neurobiological route from genetic risk to neuropsychiatric outcome remains unclear. 22q11.2 deletion syndrome (22q11.2DS) is a copy number variant (CNV) syndrome associated with high rates of neurodevelopmental and psychiatric disorders including autism spectrum disorder (ASD), attention deficit hyperactivity disorder (ADHD) and schizophrenia. Alterations in neural integration and cortical connectivity have been linked to the spectrum of neuropsychiatric disorders seen in 22q11.2DS and may be a mechanism by which the CNV acts to increase risk. In this study, magnetoencephalography (MEG) was used to investigate electrophysiological markers of local and global network function in 34 children with 22q11.2DS and 25 controls aged 10-17 years old. Resting-state oscillatory activity and functional connectivity across six frequency bands were compared between groups. Regression analyses were used to explore the relationships between these measures, neurodevelopmental symptoms and IQ. Children with 22q11.2DS had altered network activity and connectivity in high and low frequency bands, reflecting modified local and long-range cortical circuitry. Alpha and theta band connectivity were negatively associated with ASD symptoms while frontal high frequency (gamma band) activity was positively associated with ASD symptoms. Alpha band activity was positively associated with cognitive ability. These findings suggest that haploinsufficiency at the 22q11.2 locus impacts short and long-range cortical circuits, which could be a mechanism underlying neurodevelopmental and psychiatric vulnerability in this high-risk group.

Copy Number Variations in Pancreatic Cancer: From Biological Significance to Clinical Utility

Pancreatic ductal adenocarcinoma (PDAC) is the most common type of pancreatic cancer, characterized by high tumor heterogeneity and a poor prognosis. Inter- and intra-tumoral heterogeneity in PDAC is a major obstacle to effective PDAC treatment; therefore, it is highly desirable to explore the tumor heterogeneity and underlying mechanisms for the improvement of PDAC prognosis. Gene copy number variations (CNVs) are increasingly recognized as a common and heritable source of inter-individual variation in genomic sequence. In this review, we outline the origin, main characteristics, and pathological aspects of CNVs. We then describe the occurrence of CNVs in PDAC, including those that have been clearly shown to have a pathogenic role, and further highlight some key examples of their involvement in tumor development and progression. The ability to efficiently identify and analyze CNVs in tumor samples is important to support translational research and foster precision oncology, as copy number variants can be utilized to guide clinical decisions. We provide insights into understanding the CNV landscapes and the role of both somatic and germline CNVs in PDAC, which could lead to significant advances in diagnosis, prognosis, and treatment. Although there has been significant progress in this field, understanding the full contribution of CNVs to the genetic basis of PDAC will require further research, with more accurate CNV assays such as single-cell techniques and larger cohorts than have been performed to date.

Tandem gene duplications contributed to high-level azole resistance in a rapidly expanding Candida tropicalis population

Invasive diseases caused by the globally distributed commensal yeast Candida tropicalis are associated with mortality rates of greater than 50%. Notable increases of azole resistance have been observed in this species, particularly within Asia-Pacific regions. Here, we carried out a genetic population study on 1571 global C. tropicalis isolates using multilocus sequence typing (MLST). In addition, whole-genome sequencing (WGS) analysis was conducted on 629 of these strains, comprising 448 clinical invasive strains obtained in this study and 181 genomes sourced from public databases. We found that MLST clade 4 is the predominant azole-resistant clone. WGS analyses demonstrated that dramatically increasing rates of azole resistance are associated with a rapid expansion of cluster AZR, a sublineage of clade 4. Cluster AZR isolates exhibited a distinct high-level azole resistance, which was induced by tandem duplications of the ERG11A395T gene allele. Ty3/gypsy-like retrotransposons were found to be highly enriched in this population. The alarming expansion of C. tropicalis cluster AZR population underscores the urgent need for strategies against growing threats of antifungal resistance.

Rare X-linked variants carry predominantly male risk in autism, Tourette syndrome, and ADHD

Autism spectrum disorder (ASD), Tourette syndrome (TS), and attention-deficit/hyperactivity disorder (ADHD) display strong male sex bias, due to a combination of genetic and biological factors, as well as selective ascertainment. While the hemizygous nature of chromosome X (Chr X) in males has long been postulated as a key point of "male vulnerability", rare genetic variation on this chromosome has not been systematically characterized in large-scale whole exome sequencing studies of "idiopathic" ASD, TS, and ADHD. Here, we take advantage of informative recombinations in simplex ASD families to pinpoint risk-enriched regions on Chr X, within which rare maternally-inherited damaging variants carry substantial risk in males with ASD. We then apply a modified transmission disequilibrium test to 13,052 ASD probands and identify a novel high confidence ASD risk gene at exome-wide significance (MAGEC3). Finally, we observe that rare damaging variants within these risk regions carry similar effect sizes in males with TS or ADHD, further clarifying genetic mechanisms underlying male vulnerability in multiple neurodevelopmental disorders that can be exploited for systematic gene discovery.

Mock community experiments can inform on the reliability of eDNA metabarcoding data: a case study on marine phytoplankton

Environmental DNA metabarcoding is increasingly implemented in biodiversity monitoring, including phytoplankton studies. Using 21 mock communities composed of seven unicellular diatom and dinoflagellate algae, assembled with different composition and abundance by controlling the number of cells, we tested the accuracy of an eDNA metabarcoding protocol in reconstructing patterns of alpha and beta diversity. This approach allowed us to directly evaluate both qualitative and quantitative metabarcoding estimates. Our results showed non-negligible rates (17-25%) of false negatives (i.e., failure to detect a taxon in a community where it was included), for three taxa. This led to a statistically significant underestimation of metabarcoding-derived alpha diversity (Wilcoxon p = 0.02), with the detected species richness being lower than expected (based on cell numbers) in 8/21 mock communities. Considering beta diversity, the correlation between metabarcoding-derived and expected community dissimilarities was significant but not strong (R2 = 0.41), indicating suboptimal accuracy of metabarcoding results. Average biovolume and rDNA gene copy number were estimated for the seven taxa, highlighting a potential, though not exhaustive, role of the latter in explaining the recorded biases. Our findings highlight the importance of mock communities for assessing the reliability of phytoplankton eDNA metabarcoding studies and identifying their limitations.

Association between the copy number variation of CCSER1 gene and growth traits in Chinese Capra hircus (goat) populations

Recently, Coiled-coil serine-rich protein 1 (CCSER1) gene is reported to be related to economic traits in livestock, and become a hotspot. In our study, we detected CCSER1 gene CNV in 693 goats from six breeds (GZB, GZW, AN, BH, HG, TH) by quantitative real-time PCR (qPCR) and the association analysis between the types of CNV and growth traits. Then, CCSER1 gene expression pattern was discovered in seven tissues from NB goats. Our results showed that the CCSER1 gene copy numbers were distributed differently in the aforementioned six breeds. The type of CCSER1 gene CNV was significantly associated with body weight and heart girth traits in GZW goat, in which individuals with deletion type were dominant in body weight trait (P < 0.05), while the normal type individuals were more advantageous in heart girth trait (P < 0.01); and there was a significant association with heart girth in TH goat (P < 0.05), which normal type was the dominant one. The expression profile revealed that CCSER1 gene has the highest level in the lung, followed by the small intestine and heart. In conclusion, our result is dedicated to an in-depth study of the novel CCSER1 gene CNV site and to provide essential information for Chinese goats molecular selective breeding in the future.

A pan-cancer landscape of pathogenic somatic copy number variations

OBJECTIVE

Copy number variations (CNVs) play crucial roles in physiological and pathological processes, including cancer. However, the functional implications of somatic CNVs in tumor progression and evolution remain unclear. This study focuses on identifying CNV alterations with high pathogenic potential that drive and sustain tumorigenesis, distinguishing them from passenger alterations that accumulate during tumor growth. Our goal is to explore the variability of CNVs across different tumor types and infer their impact on tumor cell functions.

METHODS

Starting from 7352 copy number profiles across 33 different cancer types, we infer the pathogenicity of each CNV and perform both intra- and inter-tumor analyses to predict the functional impact of different genomic patterns. We evaluate the actionability of genes belonging to altered regions and we correlate the presence of pathogenic regions with genome instability patterns and patients' survival.

RESULTS

Our analysis uncovered large heterogeneity among different tumors suggesting in many cases distinct genetic drivers of tumorigenesis. Recurrent genomic alterations frequently coincide with dysfunctional homologous recombination pathways and negative regulation of the immune system. In certain tumors, the number of pathogenic CNVs emerged as a prognostic biomarker, highlighting their significance in cancer progression.

CONCLUSION

This study contributes to elucidate the functional impact of pathogenic CNVs in tumor progression and sheds light on their potential as prognostic markers in specific cancer types.

DNA cytosine deamination is associated with recurrent Somatic Copy Number Alterations in stomach adenocarcinoma

Stomach Adenocarcinoma (STAD) is a leading cause of death worldwide. Somatic Copy Number Alterations (SCNAs), which result in Homologous recombination (HR) deficiency in double-strand break repair, are associated with the progression of STAD. However, the landscape of frequent breakpoints of SCNAs (hotspots) and their functional impacts remain poorly understood. In this study, we aimed to explore the frequency and impact of these hotspots in 332 STAD patients and 1,043 cancer cells using data from the Cancer Genome Atlas (TCGA) and Cancer Cell Line Encyclopedia (CCLE). We studied the rates of DSB (Double-Strand Breaks) loci in STAD patients by employing the Non-Homogeneous Poisson Distribution (λ), based on which we identified 145 DSB-hotspots with genes affected. We further verified DNA cytosine deamination as a critical process underlying the burden of DSB in STAD. Finally, we illustrated the clinical impact of the significant biological processes. Our findings highlighted the relationship between DNA cytosine deamination and SCNA in cancer was associated with recurrent Somatic Copy Number Alterations in STAD.

Deep multiple-instance learning accurately predicts gene haploinsufficiency and deletion pathogenicity

Copy number losses (deletions) are a major contributor to the etiology of severe genetic disorders. Although haploinsufficient genes play a critical role in deletion pathogenicity, current methods for deletion pathogenicity prediction fail to integrate multiple lines of evidence for haploinsufficiency at the gene level, limiting their power to pinpoint deleterious deletions associated with genetic disorders. Here we introduce DosaCNV, a deep multiple-instance learning framework that, for the first time, models deletion pathogenicity jointly with gene haploinsufficiency. By integrating over 30 gene-level features potentially predictive of haploinsufficiency, DosaCNV shows unmatched performance in prioritizing pathogenic deletions associated with a broad spectrum of genetic disorders. Furthermore, DosaCNV outperforms existing methods in predicting gene haploinsufficiency even though it is not trained on known haploinsufficient genes. Finally, DosaCNV leverages a state-of-the-art technique to quantify the contributions of individual gene-level features to haploinsufficiency, allowing for human-understandable explanations of model predictions. Altogether, DosaCNV is a powerful computational tool for both fundamental and translational research.

Disruption of regulatory domains and novel transcripts as disease-causing mechanisms

Deletions, duplications, insertions, inversions, and translocations, collectively called structural variations (SVs), affect more base pairs of the genome than any other sequence variant. The recent technological advancements in genome sequencing have enabled the discovery of tens of thousands of SVs per human genome. These SVs primarily affect non-coding DNA sequences, but the difficulties in interpreting their impact limit our understanding of human disease etiology. The functional annotation of non-coding DNA sequences and methodologies to characterize their three-dimensional (3D) organization in the nucleus have greatly expanded our understanding of the basic mechanisms underlying gene regulation, thereby improving the interpretation of SVs for their pathogenic impact. Here, we discuss the various mechanisms by which SVs can result in altered gene regulation and how these mechanisms can result in rare genetic disorders. Beyond changing gene expression, SVs can produce novel gene-intergenic fusion transcripts at the SV breakpoints.

Cracking the pattern of tumor evolution based on single-cell copy number alterations

Copy number alterations (CNAs) are a key characteristic of tumor development and progression. The accumulation of various CNAs during tumor development plays a critical role in driving tumor evolution. Heterogeneous clones driven by distinct CNAs have different selective advantages, leading to differential patterns of tumor evolution that are essential for developing effective cancer therapies. Recent advances in single-cell sequencing technology have enabled genome-wide copy number profiling of tumor cell populations at single-cell resolution. This has made it possible to explore the evolutionary patterns of CNAs and accurately discover the mechanisms of intra-tumor heterogeneity. Here, we propose a two-step statistical approach that distinguishes neutral, linear, branching and punctuated evolutionary patterns for a tumor cell population based on single-cell copy number profiles. We assessed our approach using a variety of simulated and real single-cell genomic and transcriptomic datasets, demonstrating its high accuracy and robustness in predicting tumor evolutionary patterns. We applied our approach to single-cell DNA sequencing data from 20 breast cancer patients and observed that punctuated evolution is the dominant evolutionary pattern in breast cancer. Similar conclusions were drawn when applying the approach to single-cell RNA sequencing data obtained from 132 various cancer patients. Moreover, we found that differential immune cell infiltration is associated with specific evolutionary patterns. The source code of our study is available at https://github.com/FangWang-SYSU/PTEM.

The Dynamic Fungal Genome: Polyploidy, Aneuploidy and Copy Number Variation in Response to Stress

Fungal species have dynamic genomes and often exhibit genomic plasticity in response to stress. This genome plasticity often comes with phenotypic consequences that affect fitness and resistance to stress. Fungal pathogens exhibit genome plasticity in both clinical and agricultural settings and often during adaptation to antifungal drugs, posing significant challenges to human health. Therefore, it is important to understand the rates, mechanisms, and impact of large genomic changes. This review addresses the prevalence of polyploidy, aneuploidy, and copy number variation across diverse fungal species, with special attention to prominent fungal pathogens and model species. We also explore the relationship between environmental stress and rates of genomic changes and highlight the mechanisms underlying genotypic and phenotypic changes. A comprehensive understanding of these dynamic fungal genomes is needed to identify novel solutions for the increase in antifungal drug resistance.

The rate of de novo structural variation is increased in in vitro-produced offspring and preferentially affects the paternal genome

Assisted reproductive technologies (ARTs), including in vitro maturation and fertilization (IVF), are increasingly used in human and animal reproduction. Whether these technologies directly affect the rate of de novo mutation (DNM), and to what extent, has been a matter of debate. Here we take advantage of domestic cattle, characterized by complex pedigrees that are ideally suited to detect DNMs and by the systematic use of ART, to study the rate of de novo structural variation (dnSV) in this species and how it is impacted by IVF. By exploiting features of associated de novo point mutations (dnPMs) and dnSVs in clustered DNMs, we provide strong evidence that (1) IVF increases the rate of dnSV approximately fivefold, and (2) the corresponding mutations occur during the very early stages of embryonic development (one- and two-cell stage), yet primarily affect the paternal genome.

Integrative genomic and transcriptomic analysis reveals genetic alterations associated with the early progression of follicular lymphoma

Follicular lymphoma (FL), the most common indolent lymphoma, is a clinically and genetically heterogeneous disease. However, the prognostic value of driver gene mutations and copy number alterations has not been systematically assessed. Here, we analysed the clinical-biological features of 415 FL patients to identify variables associated with disease progression within 24 months of first-line therapy (POD24). Patients with B symptoms, elevated lactate dehydrogenase and β2-microglobulin levels, unfavourable baseline haemoglobin levels, advanced stage, and high-risk FL International Prognostic Index (FLIPI) scores had an increased risk of POD24, with FLIPI being the most important factor in logistic regression. HIST1H1D, identified as a driver mutation, was correlated with POD24. Gains of 6p22.2 (HIST1H1D) and 18q21.33 (BCL2) and loss of 1p36.13 (NBPF1) predicted POD24 independent of FLIPI. Gene expression profiling of FL samples showed that the POD24 cohort was significantly enriched in the inflammatory response (mediated by interferon and tumour necrosis factor), cell cycle regulation (transcription, replication and proliferation) sets and PI3K-AKT-mTOR signalling. This result was further validated with transcriptome-wide information provided by RNA-seq at single-cell resolution. Our study, performed on a large cohort of FL patients, highlights the importance of distinctive genetic alterations and gene expression relevant to disease diagnosis and early progression.

A chromosome-scale genome assembly of the grape powdery mildew pathogen Erysiphe necator reveals its genomic architecture and previously unknown features of its biology

Erysiphe necator is an obligate fungal pathogen that causes grape powdery mildew, globally the most important disease on grapevines. Previous attempts to obtain a quality genome assembly for this pathogen were hindered by its high repetitive DNA content. Here, chromatin conformation capture (Hi-C) with long-read PacBio sequencing was combined to obtain a chromosome-scale assembly and a high-quality annotation for E. necator isolate EnFRAME01. The resulting 81.1 Mb genome assembly is 98% complete and consists of 34 scaffolds, 11 of which represent complete chromosomes. All chromosomes contain large centromeric-like regions and lack synteny to the 11 chromosomes of the cereal PM pathogen Blumeria graminis. Further analysis of their composition showed that repeats and transposable elements (TEs) occupy 62.7% of their content. TEs were almost evenly interspersed outside centromeric and telomeric regions and massively overlapped with regions of annotated genes, suggesting that they could have a significant functional impact. Abundant gene duplicates were observed as well, particularly in genes encoding candidate secreted effector proteins. Moreover, younger in age gene duplicates exhibited more relaxed selection pressure and were more likely to be located physically close in the genome than older duplicates. A total of 122 genes with copy number variations among six isolates of E. necator were also identified and were enriched in genes that were duplicated in EnFRAME01, indicating they may reflect an adaptive variation. Taken together, our study illuminates higher-order genomic architectural features of E. necator and provides a valuable resource for studying genomic structural variations in this pathogen. IMPORTANCE Grape powdery mildew caused by the ascomycete fungus Erysiphe necator is economically the most important and recurrent disease in vineyards across the world. The obligate biotrophic nature of E. necator hinders the use of typical genetic methods to elucidate its pathogenicity and adaptation to adverse conditions, and thus comparative genomics has been a major method to study its genome biology. However, the current reference genome of E. necator isolate C-strain is highly fragmented with many non-coding regions left unassembled. This incompleteness prohibits in-depth comparative genomic analyses and the study of genomic structural variations (SVs) that are known to affect several aspects of microbial life, including fitness, virulence, and host adaptation. By obtaining a chromosome-scale genome assembly and a high-quality gene annotation for E. necator, we reveal the organization of its chromosomal content, unearth previously unknown features of its biology, and provide a reference for studying genomic SVs in this pathogen.

Deregulated miRNA Expression in Triple-Negative Breast Cancer of Ancestral Genomic-Characterized Latina Patients

Among patients with triple-negative breast cancer (TNBC), several studies have suggested that deregulated microRNA (miRNA) expression may be associated with a more aggressive phenotype. Although tumor molecular signatures may be race- and/or ethnicity-specific, there is limited information on the molecular profiles in women with TNBC of Hispanic and Latin American ancestry. We simultaneously profiled TNBC biopsies for the genome-wide copy number and miRNA global expression from 28 Latina women and identified a panel of 28 miRNAs associated with copy number alterations (CNAs). Four selected miRNAs (miR-141-3p, miR-150-5p, miR-182-5p, and miR-661) were validated in a subset of tumor and adjacent non-tumor tissue samples, with miR-182-5p being the most discriminatory among tissue groups (AUC value > 0.8). MiR-141-3p up-regulation was associated with increased cancer recurrence; miR-661 down-regulation with larger tumor size; and down-regulation of miR-150-5p with larger tumor size, high p53 expression, increased cancer recurrence, presence of distant metastasis, and deceased status. This study reinforces the importance of integration analysis of CNAs and miRNAs in TNBC, allowing for the identification of interactions among molecular mechanisms. Additionally, this study emphasizes the significance of considering the patients ancestral background when examining TNBC, as it can influence the relationship between intrinsic tumor molecular characteristics and clinical manifestations of the disease.

Copy Number Variation Discovery in South African Nguni-Sired and Bonsmara-Sired Crossbred Cattle

Crossbreeding forms part of Climate-Smart beef production and is one of the strategies to mitigate the effects of climate change. Two Nguni-sired and three Bonsmara-sired crossbred animals underwent whole genome sequencing. Following quality control and file preparation, the sequence data were investigated for genome-wide copy number variation (CNV) using the panelcn.MOPS tool. A total of 355 CNVs were identified in the crossbreds, of which 274 were unique in Bonsmara-sired crossbreds and 81 unique in the Nguni-sired crossbreds. Genes that differed in copy number in both crossbreds included genes related to growth (SCRN2, LOC109572916) and fertility-related factors (RPS28, LOC1098562432, LOC109570037). Genes that were present only in the Bonsmara-sired crossbreds included genes relating to lipid metabolism (MAF1), olfaction (LOC109569114), body size (HES7), immunity (LOC10957335, LOC109877039) and disease (DMBT1). Genes that were present only in the Nguni-sired crossbreds included genes relating to ketosis (HMBOX1) and amino acid transport (LOC109572916). Results of this study indicate that Nguni and Bonsmara cattle can be utilized in crossbreeding programs as they may enhance the presence of economically important traits associated with both breeds. This will produce crossbred animals that are good meat producers, grow faster, have high fertility, strong immunity and a better chance of producing in South Africa's harsh climate conditions. Ultimately, this study provides new genetic insights into the adaptability of Nguni and Bonsmara crossbred cattle.

Kinetoplast Genome of Leishmania spp. Is under Strong Purifying Selection

Instability is an intriguing characteristic of many protist genomes, and trypanosomatids are not an exception in this respect. Some regions of trypanosomatid genomes evolve fast. For instance, the trypanosomatid mitochondrial (kinetoplast) genome consists of fairly conserved maxicircle and minicircle molecules that can, nevertheless, possess high nucleotide substitution rates between closely related strains. Recent experiments have demonstrated that rapid laboratory evolution can result in the non-functionality of multiple genes of kinetoplast genomes due to the accumulation of mutations or loss of critical genomic components. An example of a loss of critical components is the reported loss of entire minicircle classes in Leishmania tarentolae during laboratory cultivation, which results in an inability to generate some correctly encoded genes. In the current work, we estimated the evolutionary rates of mitochondrial and nuclear genome regions of multiple natural Leishmania spp. We analyzed synonymous and non-synonymous substitutions and, rather unexpectedly, found that the coding regions of kinetoplast maxicircles are among the most variable regions of both genomes. In addition, we demonstrate that synonymous substitutions greatly predominate among maxicircle coding regions and that most maxicircle genes show signs of purifying selection. These results imply that maxicircles in natural Leishmania populations remain functional despite their high mutation rate.

Mycoplasma DnaK increases DNA copy number variants in vivo

The human microbiota affects critical cellular functions, although the responsible mechanism(s) is still poorly understood. In this regard, we previously showed that Mycoplasma fermentans DnaK, an HSP70 chaperone protein, hampers the activity of important cellular proteins responsible for DNA integrity. Here, we describe a novel DnaK knock-in mouse model generated in our laboratory to study the effect of M. fermentans DnaK expression in vivo. By using an array-based comparative genomic hybridization assay, we demonstrate that exposure to DnaK was associated with a higher number of DNA copy number variants (CNVs) indicative of unbalanced chromosomal alterations, together with reduced fertility and a high rate of fetal abnormalities. Consistent with their implication in genetic disorders, one of these CNVs caused a homozygous Grid2 deletion, resulting in an aberrant ataxic phenotype that recapitulates the extensive biallelic deletion in the Grid2 gene classified in humans as autosomal recessive spinocerebellar ataxia 18. Our data highlight a connection between components of the human urogenital tract microbiota, namely Mycoplasmas, and genetic abnormalities in the form of DNA CNVs, with obvious relevant medical, diagnostic, and therapeutic implications.

Mechanisms of insertions at a DNA double-strand break

Insertions and deletions (indels) are common sources of structural variation, and insertions originating from spontaneous DNA lesions are frequent in cancer. We developed a highly sensitive assay called insertion and deletion sequencing (Indel-seq) to monitor rearrangements in human cells at the TRIM37 acceptor locus that reports indels stemming from experimentally induced and spontaneous genome instability. Templated insertions, which derive from sequences genome wide, require contact between donor and acceptor loci, require homologous recombination, and are stimulated by DNA end-processing. Insertions are facilitated by transcription and involve a DNA/RNA hybrid intermediate. Indel-seq reveals that insertions are generated via multiple pathways. The broken acceptor site anneals with a resected DNA break or invades the displaced strand of a transcription bubble or R-loop, followed by DNA synthesis, displacement, and then ligation by non-homologous end joining. Our studies identify transcription-coupled insertions as a critical source of spontaneous genome instability that is distinct from cut-and-paste events.