CNVpytor: a tool for copy number variation detection and analysis from read depth and allele imbalance in whole-genome sequencing

Genomic analyses reveal high diversity and rapid evolution of Pichia kudriavzevii within a neonatal intensive care unit in Delhi, India

Pichia kudriavzevii causes life-threatening infections in immunocompromised hosts, including hospitalized neonates. This pathogen is intrinsically resistant to fluconazole, while uncommon P. kudriavzevii strains resistant to multiple antifungal drugs, including voriconazole, amphotericin B, and echinocandins, have also been reported from healthcare environments. Thus, understanding how P. kudriavzevii spread, persist, and adapt to healthcare settings could help us develop better infection management strategies. In this study, whole genome sequencing identifies multiple outbreaks of bloodstream infections in a single neonatal intensive care unit (NICU) over 5 years caused by genetically diverse strains of P. kudriavzevii. Interestingly, two genetically distinct clusters of P. kudriavzevii strains showed frequent loss of heterozygosity (LOH) events between two temporal samples. The first outbreak cluster (2015-2016) showed LOH at chromosomes 1, 4, and 5, and the other outbreak cluster (2020) exhibited LOH at chromosome 2. The circulation of two separate strain clusters of P. kudriavzevii suggests nosocomial transmission in the NICU in different time periods. Furthermore, we compared the transcriptomic profiles of three isolates of clusters I and II that exhibited distinct fluconazole and itraconazole MICs. While no significant difference in gene expression was found at the azole-target gene ERG11 or the ATP-binding cassette (ABC) transporter genes, such differences were found in genes involved in cell division and filamentation, such as SIR2 (sirtuin deacetylase) and RFA1 (replication factor A). Interestingly, increased filamentation was observed in clade I isolate exhibiting high fluconazole MICs. Together, our study indicates significant diversity, persistence, and rapid evolution of P. kudriavzevii within a single NICU.

Unveiling the fitness of Saccharomyces cerevisiae strains for lignocellulosic bioethanol: a genomic exploration through fermentation stress tests

Lignocellulosic biomass holds significant promise as a substrate for bioethanol production, yet the financial viability of lignocellulosic fermentation poses challenges. The pre-treatment step needed for lignocellulosic substrates generates inhibitors that impede Saccharomyces cerevisiae growth, affecting the fermentation process and overall yield. In modern sugarcane-to-ethanol plants, a rapid succession of yeast strains occurs, with dominant strains prevailing. Therefore, yeast strains with both dominance potential and inhibitor tolerance are crucial towards the development of superior strains with industrial fitness. This study adopted a hybrid approach combining biotechnology and bioinformatics to explore a cluster of 20 S. cerevisiae strains, including industrial and oenological strains exhibiting diverse phenotypic features. In-depth genomic analyses focusing on gene copy number variations (CNVs) and single nucleotide polymorphisms (SNPs) were conducted and compared with results from fermentation tests once inoculated in multiple strains kinetics under stressing conditions such as low nitrogen availability and high formic or acetic acid levels. Some strains showed high resistance to biotic stress and acetic acid. Moreover, four out of 20 strains - namely S. cerevisiae YI30, Fp89, Fp90 and CESPLG05 - displayed promising resistance also to formic acid, the most impactful weak acids in pre-treated lignocellulosic biomass. These strains have the potential to be used for the development of superior S. cerevisiae strains tailored for lignocellulosic bioethanol production.

Comprehensive analysis of somatic mutations and structural variations in domestic pig

Understanding somatic mutations and structural variations in domestic pigs (Sus scrofa domestica) is critical due to their increasing importance as model organisms in biomedical research. In this study, we conducted a comprehensive analysis through whole-genome sequencing of skin, organs, and blood samples. By examining two pig pedigrees, we investigated the inheritance and sharedness of structural variants among fathers, mothers, and offsprings. Utilizing single-cell clonal expansion techniques, we observed significant variations in the number of somatic mutations across different tissues. An in-house developed pipeline enabled precise filtering and analysis of these mutations, resulting in the construction of individual phylogenetic trees for two pigs. These trees explored the developmental relationships between different tissues, revealing insights into clonal expansions from various anatomical locations. This study enhances the understanding of pig genomes, affirming their increasing value in clinical and genomic research, and provides a foundation for future studies in other animals, paralleling previous studies in mice and humans. This approach not only deepens our understanding of mammalian genomic variations but also strengthens the role of pigs as a crucial model in human health and disease research.

Genome wide detection of CNV and their association with body size in Danzhou chickens

Copy number variation (CNV) is a crucial component of genetic diversity in the genome, serving as the foundation for the genetic architecture and phenotypic variability of complex traits. In this study, we examined CNVs in the Danzhou (DZ) chicken, an indigenous breed exclusive to Hainan Province, China. By employing whole-genome resequencing data from 200 DZ chickens, we conducted a comprehensive genome-wide analysis of CNVs using CNVpytor and performed CNV-based genome-wide association studies (GWAS) on 6 body size traits, including body slope length (BSL), keel length (KeL), tibial length (TiL), tibial circumference (TiC), chest width (ChW), and chest depth (ChD) utilizing linear mixed model methods considering a genomic relationship matrix. We identified a total of 144,265 autosomal CNVs among the 200 individuals, comprising 67,818 deletions and 76,447 duplications. After merging these variants together, we obtained 4,824 distinct copy number variant regions, which accounted for approximately 20% of the chicken autosomal genome. Furthermore, we discovered several significantly associated CNV segments with body size traits located proximal to genes such as IHH, WNT6, WNT10A, LPR4, FZD2, WNT7B, and GNAS that have been extensively implicated in skeletal development and growth processes. These findings enhance our understanding of CNVs in chickens and their potential impact on body size traits by revealing candidate genes involved in the regulation of these traits. This establishes a solid framework for future studies and may prove particularly beneficial for exploring genetic structural variation in chickens.

Leveraging the T2T assembly to resolve rare and pathogenic inversions in reference genome gaps

Chromosomal inversions (INVs) are particularly challenging to detect due to their copy-number neutral state and association with repetitive regions. Inversions represent about 1/20 of all balanced structural chromosome aberrations and can lead to disease by gene disruption or altering regulatory regions of dosage-sensitive genes in cis Short-read genome sequencing (srGS) can only resolve ∼70% of cytogenetically visible inversions referred to clinical diagnostic laboratories, likely due to breakpoints in repetitive regions. Here, we study 12 inversions by long-read genome sequencing (lrGS) (n = 9) or srGS (n = 3) and resolve nine of them. In four cases, the inversion breakpoint region was missing from at least one of the human reference genomes (GRCh37, GRCh38, T2T-CHM13) and a reference agnostic analysis was needed. One of these cases, an INV9 mappable only in de novo assembled lrGS data using T2T-CHM13 disrupts EHMT1 consistent with a Mendelian diagnosis (Kleefstra syndrome 1; MIM#610253). Next, by pairwise comparison between T2T-CHM13, GRCh37, and GRCh38, as well as the chimpanzee and bonobo, we show that hundreds of megabases of sequence are missing from at least one human reference, highlighting that primate genomes contribute to genomic diversity. Aligning population genomic data to these regions indicated that these regions are variable between individuals. Our analysis emphasizes that T2T-CHM13 is necessary to maximize the value of lrGS for optimal inversion detection in clinical diagnostics. These results highlight the importance of leveraging diverse and comprehensive reference genomes to resolve unsolved molecular cases in rare diseases.

Structural Variants in COL1A1 and COL1A2 in Osteogenesis Imperfecta

Osteogenesis Imperfecta (OI) is a heterogeneous skeletal dysplasia characterized by bone fragility, skeletal deformities, and short stature. Most commonly, it is caused by autosomal dominant variants in the type I collagen genes, COL1A1 or COL1A2. Type I collagen is the main protein of the extracellular matrix in the skeleton and changes in its structure or quantity may lead to OI. 85%-90% of OI cases occur due to sequence variants in type I collagen genes, while OI caused by structural abnormalities in type I collagen genes is less common. In most cases, haploinsufficiency of type I collagen is associated with a milder OI phenotype. Large genomic deletions often involve several genes within the same chromosomal region, leading to microdeletion syndromes with OI features. Here, we report eight Swedish patients from five unrelated families with OI due to structural variants in the COL1A1 and COL1A2 genes. One patient with OI type III had a complex rearrangement with a deletion and duplication event in COL1A2, leading to reduced COL1A2 expression. Three other patients from two different families with OI type I had whole gene deletions involving COL1A1. In one family, three affected individuals with OI type I had a small intragenic deletion of exons 11-12 in COL1A2. One patient had a 2.1 Mb de novo deletion encompassing COL1A1 and DLX3 genes and features of OI and tricho-dento-osseous syndrome. Overall, this study highlights the importance of investigating gene dosage abnormalities in patients with OI and further delineates clinical and genetic variability of OI caused by structural variants in type I collagen genes.

High level of aneuploidy and recurrent loss of chromosome 11 as relevant features of somatotroph pituitary tumors

BACKGROUND

Somatotroph neuroendocrine pituitary tumors (sPitNET) are a subtype of pituitary tumors that commonly cause acromegaly. Our study aimed to determine the spectrum of DNA copy number abnormalities (CNAs) in sPitNETs and their relevance.

METHODS

A landscape of CNAs in sPitNETs was determined using combined whole-genome approaches involving low-pass whole genome sequencing and SNP microarrays. Fluorescent in situ hybridization (FISH) was used for microscopic validation of CNAs. The tumors were also subjected to transcriptome and DNA methylation analyses with RNAseq and microarrays, respectively.

RESULTS

We observed a wide spectrum of cytogenetic changes ranging from multiple deletions, recurrent chromosome 11 loss, stable genomes, to duplication of the majority of the chromosomes. The identified CNAs were confirmed with FISH. sPitNETs with multiple duplications were characterized by intratumoral heterogeneity in chromosome number variation in individual tumor cells, as determined with FISH. These tumors were separate CNA-related sPitNET subtype in clustering analyses with CNA signature specific for whole genome doubling-related etiology. This subtype encompassed GNAS-wild type, mostly densely granulated tumors with favorable expression level of known prognosis-related genes, notably enriched with POUF1/NR5A1-double positive PitNETs. Chromosomal deletions in sPitNETs are functionally relevant. They occurred in gene-dense DNA regions and were related to genes downregulation and increased DNA methylation in the CpG island and promoter regions in the affected regions. Recurrent loss of chromosome 11 was reflected by lowered MEN1 and AIP. No such unequivocal relevance was found for chromosomal gains. Comparisons of transcriptomes of selected most cytogenetically stable sPitNETs with tumors with recurrent loss of chromosome 11 showed upregulation of processes related to gene dosage compensation mechanism in tumors with deletion. Comparison of stable tumors with those with multiple duplications showed upregulation of processes related to mitotic spindle, DNA repair, and chromatin organization. Both comparisons showed upregulation of the processes related to immune infiltration in cytogenetically stable tumors and deconvolution of DNA methylation data indicated a higher content of specified immune cells and lower tumor purity in these tumors.

CONCLUSIONS

sPitNETs fall into three relevant cytogenetic groups: highly aneuploid tumors characterized by known prognostically favorable features and low aneuploidy tumors including specific subtype with chromosome 11 loss.

A comprehensive map of copy number variations in dromedary camels based on whole genome sequence data

Copy number variants (CNVs) are structural variants within the eukaryotic genome that vary among individuals of a species. These variants have been associated with different phenotypic traits, making them a valuable consideration as markers for designing breeding programmes. In this study, whole genome sequence data of 60 dromedary camel samples originating from the Arabian Peninsula were analyzed to construct a comprehensive dromedary CNV map. Utilizing four CNV callers employing read-depth, split-read and paired-end mapping approaches, a total of 37,519 CNV events (17,847 deletions and 19,672 duplications) were called on the dromedary autosomes. These CNV events were merged into 2,557 regions, categorized as 1,322 losses, 122 gains, and 1,113 "mixed regions" comprising both types. The cumulative size of the CNV regions amounted to 22.5 Mb, covering roughly 1.16% of the dromedary autosomes. Approximately 32% of the defined CNV regions (comprising 60% losses, 18% gains, and 0.27% mixed regions) were found in ≥ 90% of the dromedary samples, classifying them as prevalent regions. Genes with biological functions related to the different adaptive physiologies of dromedary camels, such as fertility, heat stress, musculoskeletal development, and fat metabolism, were overlapping with or in close proximity to ~ 68% of the defined CNV regions, demonstrating their potential role in dromedaries' physiology. This study presents the first comprehensive CNV map of dromedary camels and builds on the present knowledge in understanding the genetic structure of this species.

Somatic mosaicism in schizophrenia brains reveals prenatal mutational processes

Germline mutations modulate the risk of developing schizophrenia (SCZ). Much less is known about the role of mosaic somatic mutations in the context of SCZ. Deep (239×) whole-genome sequencing (WGS) of brain neurons from 61 SCZ cases and 25 controls postmortem identified mutations occurring during prenatal neurogenesis. SCZ cases showed increased somatic variants in open chromatin, with increased mosaic CpG transversions (CpG>GpG) and T>G mutations at transcription factor binding sites (TFBSs) overlapping open chromatin, a result not seen in controls. Some of these variants alter gene expression, including SCZ risk genes and genes involved in neurodevelopment. Although these mutational processes can reflect a difference in factors indirectly involved in disease, increased somatic mutations at developmental TFBSs could also potentially contribute to SCZ.

SeqVerify: An accessible analysis tool for cell line genomic integrity, contamination, and gene editing outcomes

Over the last decade, advances in genome editing and pluripotent stem cell (PSC) culture have let researchers generate edited PSC lines to study a wide variety of biological questions. However, abnormalities in cell lines such as aneuploidy, mutations, on-target and off-target editing errors, and microbial contamination can arise during PSC culture or due to undesired editing outcomes. The ongoing decline of next-generation sequencing prices has made whole-genome sequencing (WGS) a promising option for detecting these abnormalities. However, this approach has been held back by a lack of easily usable data analysis software. Here, we present SeqVerify, a computational pipeline designed to take raw WGS data and a list of intended genome edits, and verify that the edits are present and that there are no abnormalities. We anticipate that SeqVerify will be a useful tool for researchers generating edited PSCs, and more broadly, for cell line quality control in general.

A novel copy number variant in the murine Cdh23 gene gives rise to profound deafness and vestibular dysfunction

Hearing loss is the most common congenital sensory deficit worldwide and exhibits high genetic heterogeneity, making molecular diagnoses elusive for most individuals. Detecting novel mutations that contribute to hearing loss is crucial to providing accurate personalized diagnoses, tailored interventions, and improving prognosis. Copy number variants (CNVs) are structural mutations that are understudied, potential contributors to hearing loss. Here, we present the Abnormal Wobbly Gait (AWG) mouse, the first documented mutant exhibiting waltzer-like locomotor dysfunction, hyperactivity, circling behaviour, and profound deafness caused by a spontaneous CNV deletion in cadherin 23 (Cdh23). We were unable to identify the causative mutation through a conventional whole-genome sequencing (WGS) and variant detection pipeline, but instead found a linked variant in hexokinase 1 (Hk1) that was insufficient to recapitulate the AWG phenotype when introduced into C57BL/6J mice using CRISPR-Cas9. Investigating nearby deafness-associated genes revealed a pronounced downregulation of Cdh23 mRNA and a complete absence of full-length CDH23 protein, which is critical for the development and maintenance of inner ear hair cells, in whole head extracts from AWG neonates. Manual inspection of WGS read depth plots of the Cdh23 locus revealed a putative 10.4 kb genomic deletion of exons 11 and 12 that was validated by PCR and Sanger sequencing. This study underscores the imperative to refine variant detection strategies to permit identification of pathogenic CNVs easily missed by conventional variant calling to enhance diagnostic precision and ultimately improve clinical outcomes for individuals with genetically heterogenous disorders such as hearing loss.

Genetic architecture of long-distance migration and population genomics of the endangered Japanese eel

The Japanese eel (Anguilla japonica), a flagship anguillid species for conservation, is known for its long-distance-oriented migration. However, our understanding of the genetic architecture underlying long-distance migration and population genomic characteristics of A. japonica is still limited. Here, we generated a high-quality chromosome-level genome assembly and conducted whole-genome resequencing of 218 individuals to explore these aspects. Strong signals of selection were found on genes involved in long-distance aerobic exercise and navigation, which might be associated with evolutionary adaptation to long-distance migrations. Low genetic diversity was detected, which might result from genetic drift associated with demographic declines. Both mitochondrial and nuclear genomic datasets supported the existence of a single panmictic population for Japanese eel, despite signals of single-generation selection. Candidate genes for local selection involved in functions like development and circadian rhythm. The findings can provide insights to adaptative evolution to long-distance migration and inform conservation efforts for A. japonica.

Genome-wide analysis and visualization of copy number with CNVpytor in igv.js

SUMMARY

Copy number variation (CNV) and alteration (CNA) analysis is a crucial component in many genomic studies and its applications span from basic research to clinic diagnostics and personalized medicine. CNVpytor is a tool featuring a read depth-based caller and combined read depth and B-allele frequency (BAF) based 2D caller to find CNVs and CNAs. The tool stores processed intermediate data and CNV/CNA calls in a compact HDF5 file-pytor file. Here, we describe a new track in igv.js that utilizes pytor and whole genome variant files as input for on-the-fly read depth and BAF visualization, CNV/CNA calling and analysis. Embedding into HTML pages and Jupiter Notebooks enables convenient remote data access and visualization simplifying interpretation and analysis of omics data.

AVAILABILITY AND IMPLEMENTATION

The CNVpytor track is integrated with igv.js and available at https://github.com/igvteam/igv.js. The documentation is available at https://github.com/igvteam/igv.js/wiki/cnvpytor. Usage can be tested in the IGV-Web app at https://igv.org/app and also on https://github.com/abyzovlab/CNVpytor.

A case of an Angelman-syndrome caused by an intragenic duplication of UBE3A uncovered by adaptive nanopore sequencing

Adaptive nanopore sequencing as a diagnostic method for imprinting disorders and episignature analysis revealed an intragenic duplication of Exon 6 and 7 in UBE3A (NM_000462.5) in a patient with relatively mild Angelman-like syndrome. In an all-in-one nanopore sequencing analysis DNA hypomethylation of the SNURF:TSS-DMR, known contributing deletions on the maternal allele and point mutations in UBE3A could be ruled out as disease drivers. In contrast, breakpoints and orientation of the tandem duplication could clearly be defined. Segregation analysis in the family showed that the duplication derived de novo in the maternal grandfather. Our study shows the benefits of an all-in-one nanopore sequencing approach for the diagnostics of Angelman syndrome and other imprinting disorders.

Investigating copy number variants in schizophrenia pedigrees using a new consensus pipeline called PECAN

Copy number variants (CNVs) have been implicated in many human diseases, including psychiatric disorders. Whole genome sequencing offers advantages in CNV calling compared to previous array-based methods. Here we present a robust and transparent CNV calling pipeline, PECAN (PEdigree Copy number vAriaNt calling), for short-read, whole genome sequencing data, comprised of a novel combination of four calling methods and structural variant genotyping. This method is scalable and can incorporate pedigree information to retain lower-confidence CNVs that would otherwise be discarded. We have robustly benchmarked PECAN using gold-standard CNV calls for two well-established evaluation samples, NA12878 and HG002, showing that PECAN performs with high precision and recall on both datasets, outperforming another pedigree-based CNV calling pipeline. As part of this work, we provide a list of high-confidence gold standard CNVs for the NA12878 reference sample, curated from multiple studies. We applied PECAN to a collection of pedigrees multiply affected with schizophrenia and identified a rare deletion that perfectly co-segregates with schizophrenia in one of the pedigrees. The CNV overlaps the gene PITRM1, which has been implicated in a complex phenotype including ataxia, developmental delay, and schizophrenia-like episodes in affected adults.

High-depth whole-genome sequencing identifies structure variants, copy number variants and short tandem repeats associated with Parkinson's disease

While numerous single nucleotide variants and small indels have been identified in Parkinson's disease (PD), the contribution of structural variants (SVs), copy number variants (CNVs), and short tandem repeats (STRs) remains poorly understood. Here we investigated the association using the high-depth whole-genome sequencing data from 466 Chinese PD patients and 513 controls. Totally, we identified 29,561 SVs, 32,153 CNVs, and 174,905 STRs, and found that CNV deletions were significantly enriched in the end-proportion of autosomal chromosomes in PD. After genome-wide association analysis and replication in an external cohort of 352 cases and 547 controls, we validated that the 1.6 kb-deletion neighboring MUC19, 12.4kb-deletion near RXFP1 and GGGAAA repeats in SLC2A13 were significantly associated with PD. Moreover, the MUC19 deletion and the SLC2A13 5-copy repeat reduced the penetrance of the LRRK2 G2385R variant. Moreover, genes with these variants were dosage-sensitive. These data provided novel insights into the genetic architecture of PD.

An assessment of the genomic structural variation landscape in Sub-Saharan African populations

Structural variants are responsible for a large part of genomic variation between individuals and play a role in both common and rare diseases. Databases cataloguing structural variants notably do not represent the full spectrum of global diversity, particularly missing information from most African populations. To address this representation gap, we analysed 1,091 high-coverage African genomes, 545 of which are public data sets, and 546 which have been analysed for structural variants for the first time. Variants were called using five different tools and datasets merged and jointly called using SURVIVOR. We identified 67,795 structural variants throughout the genome, with 10,421 genes having at least one variant. Using a conservative overlap in merged data, 6,414 of the structural variants (9.5%) are novel compared to the Database of Genomic Variants. This study contributes to knowledge of the landscape of structural variant diversity in Africa and presents a reliable dataset for potential applications in population genetics and health-related research.

Assembly and Genome Annotation of Different Strains of Apple Fruit Moth Virus (Cydia pomonella granulovirus)

Cydia pomonella granulovirus is a natural pathogen for Cydia pomonella that is used as a biocontrol agent of insect populations. The study of granulovirus virulence is of particular interest since the development of resistance in natural populations of C. pomonella has been observed during the long-term use of the Mexican isolate CpGV. In our study, we present the genomes of 18 CpGV strains endemic to southern Russia and from Kazakhstan, as well as a strain included in the commercial preparation "Madex Twin", which were sequenced and analyzed. We performed comparative genomic analysis using several tools. From comparisons at the level of genes and protein products that are involved in the infection process of virosis, synonymous and missense substitution variants have been identified. The average nucleotide identity has demonstrated a high similarity with other granulovirus genomes of different geographic origins. Whole-genome alignment of the 18 genomes relative to the reference revealed regions of low similarity. Analysis of gene repertoire variation has shown that BZR GV 4, BZR GV 6, and BZR GV L-7 strains have been the closest in gene content to the commercial "Madex Twin" strain. We have confirmed two deletions using read depth coverage data in regions lacking genes shown by homology analysis for granuloviruses BZR GV L-4 and BZR GV L-6; however, they are not related to the known genes causing viral pathogenicity. Thus, we have isolated novel CpGV strains and analyzed their potential as strains producing highly effective bioinsecticides against C. pomonella.

Limited Added Diagnostic Value of Whole Genome Sequencing in Genetic Testing of Inherited Retinal Diseases in a Swiss Patient Cohort

The purpose of this study was to assess the added diagnostic value of whole genome sequencing (WGS) for patients with inherited retinal diseases (IRDs) who remained undiagnosed after whole exome sequencing (WES). WGS was performed for index patients in 66 families. The datasets were analyzed according to GATK's guidelines. Additionally, DeepVariant was complemented by GATK's workflow, and a novel structural variant pipeline was developed. Overall, a molecular diagnosis was established in 19/66 (28.8%) index patients. Pathogenic deletions and one deep-intronic variant contributed to the diagnostic yield in 4/19 and 1/19 index patients, respectively. The remaining diagnoses (14/19) were attributed to exonic variants that were missed during WES analysis due to bioinformatic limitations, newly described loci, or unclear pathogenicity. The added diagnostic value of WGS equals 5/66 (9.6%) for our cohort, which is comparable to previous studies. This figure would decrease further to 1/66 (1.5%) with a standardized and reliable copy number variant workflow during WES analysis. Given the higher costs and limited added value, the implementation of WGS as a first-tier assay for inherited eye disorders in a diagnostic laboratory remains untimely. Instead, progress in bioinformatic tools and communication between diagnostic and clinical teams have the potential to ameliorate diagnostic yields.

The clinical utility and diagnostic implementation of human subject cell transdifferentiation followed by RNA sequencing

RNA sequencing (RNA-seq) has recently been used in translational research settings to facilitate diagnoses of Mendelian disorders. A significant obstacle for clinical laboratories in adopting RNA-seq is the low or absent expression of a significant number of disease-associated genes/transcripts in clinically accessible samples. As this is especially problematic in neurological diseases, we developed a clinical diagnostic approach that enhanced the detection and evaluation of tissue-specific genes/transcripts through fibroblast-to-neuron cell transdifferentiation. The approach is designed specifically to suit clinical implementation, emphasizing simplicity, cost effectiveness, turnaround time, and reproducibility. For clinical validation, we generated induced neurons (iNeurons) from 71 individuals with primary neurological phenotypes recruited to the Undiagnosed Diseases Network. The overall diagnostic yield was 25.4%. Over a quarter of the diagnostic findings benefited from transdifferentiation and could not be achieved by fibroblast RNA-seq alone. This iNeuron transcriptomic approach can be effectively integrated into diagnostic whole-transcriptome evaluation of individuals with genetic disorders.

The rise of baobab trees in Madagascar

The baobab trees (genus Adansonia) have attracted tremendous attention because of their striking shape and distinctive relationships with fauna1. These spectacular trees have also influenced human culture, inspiring innumerable arts, folklore and traditions. Here we sequenced genomes of all eight extant baobab species and argue that Madagascar should be considered the centre of origin for the extant lineages, a key issue in their evolutionary history2,3. Integrated genomic and ecological analyses revealed the reticulate evolution of baobabs, which eventually led to the species diversity seen today. Past population dynamics of Malagasy baobabs may have been influenced by both interspecific competition and the geological history of the island, especially changes in local sea levels. We propose that further attention should be paid to the conservation status of Malagasy baobabs, especially of Adansonia suarezensis and Adansonia grandidieri, and that intensive monitoring of populations of Adansonia za is required, given its propensity for negatively impacting the critically endangered Adansonia perrieri.

The prevalence of copy number increase at multiallelic copy number variants associated with cave colonization

Copy number variation is a common contributor to phenotypic diversity, yet its involvement in ecological adaptation is not easily discerned. Instances of parallelly evolving populations of the same species in a similar environment marked by strong selective pressures present opportunities to study the role of copy number variants (CNVs) in adaptation. By identifying CNVs that repeatedly occur in multiple populations of the derived ecotype and are not (or are rarely) present in the populations of the ancestral ecotype, the association of such CNVs with adaptation to the novel environment can be inferred. We used this paradigm to identify CNVs associated with recurrent adaptation of the Mexican tetra (Astyanax mexicanus) to cave environment. Using a read-depth approach, we detected CNVs from previously re-sequenced genomes of 44 individuals belonging to two ancestral surfaces and three derived cave populations. We identified 102 genes and 292 genomic regions that repeatedly diverge in copy number between the two ecotypes and occupy 0.8% of the reference genome. Functional analysis revealed their association with processes previously recognized to be relevant for adaptation, such as vision, immunity, oxygen consumption, metabolism, and neural function and we propose that these variants have been selected for in the cave or surface waters. The majority of the ecotype-divergent CNVs are multiallelic and display copy number increases in cavefish compared to surface fish. Our findings suggest that multiallelic CNVs - including gene duplications - and divergence in copy number provide a fast route to produce novel phenotypes associated with adaptation to subterranean life.

Mind the gap: the relevance of the genome reference to resolve rare and pathogenic inversions

Both long-read genome sequencing (lrGS) and the recently published Telomere to Telomere (T2T) reference genome provide increased coverage and resolution across repetitive regions promising heightened structural variant detection and improved mapping. Inversions (INV), intrachromosomal segments which are rotated 180° and inserted back into the same chromosome, are a class of structural variants particularly challenging to detect due to their copy-number neutral state and association with repetitive regions. Inversions represent about 1/20 of all balanced structural chromosome aberrations and can lead to disease by gene disruption or altering regulatory regions of dosage sensitive genes in cis . Here we remapped the genome data from six individuals carrying unsolved cytogenetically detected inversions. An INV6 and INV10 were resolved using GRCh38 and T2T-CHM13. Finally, an INV9 required optical genome mapping, de novo assembly of lrGS data and T2T-CHM13. This inversion disrupted intron 25 of EHMT1, confirming a diagnosis of Kleefstra syndrome 1 (MIM#610253). These three inversions, only mappable in specific references, prompted us to investigate the presence and population frequencies of differential reference regions (DRRs) between T2T-CHM13, GRCh37, GRCh38, the chimpanzee and bonobo, and hundreds of megabases of DRRs were identified. Our results emphasize the significance of the chosen reference genome and the added benefits of lrGS and optical genome mapping in solving rearrangements in challenging regions of the genome. This is particularly important for inversions and may impact clinical diagnostics.

Association of Structural Forms of 17q21.31 with the Risk of Progressive Supranuclear Palsy and MAPT Sub-haplotypes

Importance

The chromosome 17q21.31 region, containing a 900 Kb inversion that defines H1 and H2 haplotypes, represents the strongest genetic risk locus in progressive supranuclear palsy (PSP). In addition to H1 and H2, various structural forms of 17q21.31, characterized by the copy number of α, β, and γ duplications, have been identified. However, the specific effect of each structural form on the risk of PSP has never been evaluated in a large cohort study.

Objective

To assess the association of different structural forms of 17q.21.31, defined by the copy numbers of α, β, and γ duplications, with the risk of PSP and MAPT sub-haplotypes.

Design setting and participants

Utilizing whole genome sequencing data of 1,684 (1,386 autopsy confirmed) individuals with PSP and 2,392 control subjects, a case-control study was conducted to investigate the association of copy numbers of α, β, and γ duplications and structural forms of 17q21.31 with the risk of PSP. All study subjects were selected from the Alzheimer's Disease Sequencing Project (ADSP) Umbrella NG00067.v7. Data were analyzed between March 2022 and November 2023.

Main outcomes and measures

The main outcomes were the risk (odds ratios [ORs]) for PSP with 95% CIs. Risks for PSP were evaluated by logistic regression models.

Results

The copy numbers of α and β were associated with the risk of PSP only due to their correlation with H1 and H2, while the copy number of γ was independently associated with the increased risk of PSP. Each additional duplication of γ was associated with 1.10 (95% CI, 1.04-1.17; P = 0.0018) fold of increased risk of PSP when conditioning H1 and H2. For the H1 haplotype, addition γ duplications displayed a higher odds ratio for PSP: the odds ratio increases from 1.21 (95%CI 1.10-1.33, P = 5.47 × 10-5) for H1β1γ1 to 1.29 (95%CI 1.16-1.43, P = 1.35 × 10-6) for H1β1γ2, 1.45 (95%CI 1.27-1.65, P = 3.94 × 10-8) for H1β1γ3, and 1.57 (95%CI 1.10-2.26, P = 1.35 × 10-2) for H1β1γ4. Moreover, H1β1γ3 is in linkage disequilibrium with H1c (R2 = 0.31), a widely recognized MAPT sub-haplotype associated with increased risk of PSP. The proportion of MAPT sub-haplotypes associated with increased risk of PSP (i.e., H1c, H1d, H1g, H1o, and H1h) increased from 34% in H1β1γ1 to 77% in H1β1γ4.

Conclusions and relevance

This study revealed that the copy number of γ was associated with the risk of PSP independently from H1 and H2. The H1 haplotype with more γ duplications showed a higher odds ratio for PSP and were associated with MAPT sub-haplotypes with increased risk of PSP. These findings expand our understanding of how the complex structure at 17q21.31 affect the risk of PSP.

Increased privatization of a public resource leads to spread of cooperation in a microbial population

The phenomenon of cooperation is prevalent at all levels of life. In one such manifestation of cooperation in microbial communities, some cells produce costly extracellular resources that are freely available to others. These resources are referred to as public goods. Saccharomyces cerevisiae secretes invertase (public good) in the periplasm to hydrolyze sucrose into glucose and fructose, which are then imported by the cells. After hydrolysis of sucrose, a cooperator retains only 1% of the monosaccharides, while 99% of the monosaccharides diffuse into the environment and can be utilized by any cell. The non-producers of invertase (cheaters) exploit the invertase-producing cells (cooperators) by utilizing the monosaccharides and not paying the metabolic cost of producing the invertase. In this work, we investigate the evolutionary dynamics of this cheater-cooperator system. In a co-culture, if cheaters are selected for their higher fitness, the population will collapse. On the other hand, for cooperators to survive in the population, a strategy to increase fitness would likely be required. To understand the adaptation of cooperators in sucrose, we performed a coevolution experiment in sucrose. Our results show that cooperators increase in fitness as the experiment progresses. This phenomenon was not observed in environments which involved a non-public good system. Genome sequencing reveals duplication of several HXT transporters in the evolved cooperators. Based on these results, we hypothesize that increased privatization of the monosaccharides is the most likely explanation of spread of cooperators in the population.IMPORTANCEHow is cooperation, as a trait, maintained in a population? In order to answer this question, we perform a coevolution experiment between two strains of yeast-one which produces a public good to release glucose and fructose in the media, thus generating a public resource, and the other which does not produce public resource and merely benefits from the presence of the cooperator strain. What is the outcome of this coevolution experiment? We demonstrate that after ~200 generations of coevolution, cooperators increase in frequency in the co-culture. Remarkably, in all parallel lines of our experiment, this is obtained via duplication of regions which likely allow greater privatization of glucose and fructose. Thus, increased privatization, which is intuitively thought to be a strategy against cooperation, enables spread of cooperation.

Generation of a Commercial-Scale Founder Population of Porcine Reproductive and Respiratory Syndrome Virus Resistant Pigs Using CRISPR-Cas

Disease resistance genes in livestock provide health benefits to animals and opportunities for farmers to meet the growing demand for affordable, high-quality protein. Previously, researchers used gene editing to modify the porcine CD163 gene and demonstrated resistance to a harmful virus that causes porcine reproductive and respiratory syndrome (PRRS). To maximize potential benefits, this disease resistance trait needs to be present in commercially relevant breeding populations for multiplication and distribution of pigs. Toward this goal, a first-of-its-kind, scaled gene editing program was established to introduce a single modified CD163 allele into four genetically diverse, elite porcine lines. This effort produced healthy pigs that resisted PRRS virus infection as determined by macrophage and animal challenges. This founder population will be used for additional disease and trait testing, multiplication, and commercial distribution upon regulatory approval. Applying CRISPR-Cas to eliminate a viral disease represents a major step toward improving animal health.

Rediscovering tandem repeat variation in schizophrenia: challenges and opportunities

Tandem repeats (TRs) are prevalent throughout the genome, constituting at least 3% of the genome, and often highly polymorphic. The high mutation rate of TRs, which can be orders of magnitude higher than single-nucleotide polymorphisms and indels, indicates that they are likely to make significant contributions to phenotypic variation, yet their contribution to schizophrenia has been largely ignored by recent genome-wide association studies (GWAS). Tandem repeat expansions are already known causative factors for over 50 disorders, while common tandem repeat variation is increasingly being identified as significantly associated with complex disease and gene regulation. The current review summarizes key background concepts of tandem repeat variation as pertains to disease risk, elucidating their potential for schizophrenia association. An overview of next-generation sequencing-based methods that may be applied for TR genome-wide identification is provided, and some key methodological challenges in TR analyses are delineated.

The landscape of genomic structural variation in Indigenous Australians

Indigenous Australians harbour rich and unique genomic diversity. However, Aboriginal and Torres Strait Islander ancestries are historically under-represented in genomics research and almost completely missing from reference datasets1-3. Addressing this representation gap is critical, both to advance our understanding of global human genomic diversity and as a prerequisite for ensuring equitable outcomes in genomic medicine. Here we apply population-scale whole-genome long-read sequencing4 to profile genomic structural variation across four remote Indigenous communities. We uncover an abundance of large insertion-deletion variants (20-49 bp; n = 136,797), structural variants (50  b-50 kb; n = 159,912) and regions of variable copy number (>50 kb; n = 156). The majority of variants are composed of tandem repeat or interspersed mobile element sequences (up to 90%) and have not been previously annotated (up to 62%). A large fraction of structural variants appear to be exclusive to Indigenous Australians (12% lower-bound estimate) and most of these are found in only a single community, underscoring the need for broad and deep sampling to achieve a comprehensive catalogue of genomic structural variation across the Australian continent. Finally, we explore short tandem repeats throughout the genome to characterize allelic diversity at 50 known disease loci5, uncover hundreds of novel repeat expansion sites within protein-coding genes, and identify unique patterns of diversity and constraint among short tandem repeat sequences. Our study sheds new light on the dimensions and dynamics of genomic structural variation within and beyond Australia.

BRACNAC: A BRCA1 and BRCA2 Copy Number Alteration Caller from Next-Generation Sequencing Data

Detecting copy number variations (CNVs) and alterations (CNAs) in the BRCA1 and BRCA2 genes is essential for testing patients for targeted therapy applicability. However, the available bioinformatics tools were initially designed for identifying CNVs/CNAs in whole-genome or -exome (WES) NGS data or targeted NGS data without adaptation to the BRCA1/2 genes. Most of these tools were tested on sample cohorts of limited size, with their use restricted to specific library preparation kits or sequencing platforms. We developed BRACNAC, a new tool for detecting CNVs and CNAs in the BRCA1 and BRCA2 genes in NGS data of different origin. The underlying mechanism of this tool involves various coverage normalization steps complemented by CNV probability evaluation. We estimated the sensitivity and specificity of our tool to be 100% and 94%, respectively, with an area under the curve (AUC) of 94%. The estimation was performed using the NGS data obtained from 213 ovarian and prostate cancer samples tested with in-house and commercially available library preparation kits and additionally using multiplex ligation-dependent probe amplification (MLPA) (12 CNV-positive samples). Using freely available WES and targeted NGS data from other research groups, we demonstrated that BRACNAC could also be used for these two types of data, with an AUC of up to 99.9%. In addition, we determined the limitations of the tool in terms of the minimum number of samples per NGS run (≥20 samples) and the minimum expected percentage of CNV-negative samples (≥80%). We expect that our findings will improve the efficacy of BRCA1/2 diagnostics. BRACNAC is freely available at the GitHub server.

Genomic data resources of the Brain Somatic Mosaicism Network for neuropsychiatric diseases

Somatic mosaicism is defined as an occurrence of two or more populations of cells having genomic sequences differing at given loci in an individual who is derived from a single zygote. It is a characteristic of multicellular organisms that plays a crucial role in normal development and disease. To study the nature and extent of somatic mosaicism in autism spectrum disorder, bipolar disorder, focal cortical dysplasia, schizophrenia, and Tourette syndrome, a multi-institutional consortium called the Brain Somatic Mosaicism Network (BSMN) was formed through the National Institute of Mental Health (NIMH). In addition to genomic data of affected and neurotypical brains, the BSMN also developed and validated a best practices somatic single nucleotide variant calling workflow through the analysis of reference brain tissue. These resources, which include >400 terabytes of data from 1087 subjects, are now available to the research community via the NIMH Data Archive (NDA) and are described here.

SeqVerify: An accessible analysis tool for cell line genomic integrity, contamination, and gene editing outcomes

1Over the last decade, advances in genome editing and pluripotent stem cell (PSC) culture have let researchers generate edited PSC lines to study a wide variety of biological questions. However, abnormalities in cell lines such as aneuploidy, on-target and off-target editing errors, and microbial contamination can arise during PSC culture or due to undesired editing outcomes. Any of these abnormalities can invalidate experiments, so detecting them is crucial. The ongoing decline of next-generation sequencing prices has made whole genome sequencing (WGS) an effective quality control option, since WGS can detect any abnormality involving changes to DNA sequences or presence of unwanted sequences. However, this approach has suffered from a lack of easily usable data analysis software. Here, we present SeqVerify, a computational pipeline designed to take raw WGS data and a list of intended edits, and verify that the edits are present and that there are no abnormalities. We anticipate that SeqVerify will be a useful tool for researchers generating edited PSCs, and more broadly, for cell line quality control in general.

Read-depth based approach on whole genome resequencing data reveals important insights into the copy number variation (CNV) map of major global buffalo breeds

BACKGROUND

Elucidating genome-wide structural variants including copy number variations (CNVs) have gained increased significance in recent times owing to their contribution to genetic diversity and association with important pathophysiological states. The present study aimed to elucidate the high-resolution CNV map of six different global buffalo breeds using whole genome resequencing data at two coverages (10X and 30X). Post-quality control, the sequence reads were aligned to the latest draft release of the Bubaline genome. The genome-wide CNVs were elucidated using a read-depth approach in CNVnator with different bin sizes. Adjacent CNVs were concatenated into copy number variation regions (CNVRs) in different breeds and their genomic coverage was elucidated.

RESULTS

Overall, the average size of CNVR was lower at 30X coverage, providing finer details. Most of the CNVRs were either deletion or duplication type while the occurrence of mixed events was lesser in number on a comparative basis in all breeds. The average CNVR size was lower at 30X coverage (0.201 Mb) as compared to 10X (0.013 Mb) with the finest variants in Banni buffaloes. The maximum number of CNVs was observed in Murrah (2627) and Pandharpuri (25,688) at 10X and 30X coverages, respectively. Whereas the minimum number of CNVs were scored in Surti at both coverages (2092 and 17,373). On the other hand, the highest and lowest number of CNVRs were scored in Jaffarabadi (833 and 10,179 events) and Surti (783 and 7553 events) at both coverages. Deletion events overnumbered duplications in all breeds at both coverages. Gene profiling of common overlapped genes and longest CNVRs provided important insights into the evolutionary history of these breeds and indicate the genomic regions under selection in respective breeds.

CONCLUSION

The present study is the first of its kind to elucidate the high-resolution CNV map in major buffalo populations using a read-depth approach on whole genome resequencing data. The results revealed important insights into the divergence of major global buffalo breeds along the evolutionary timescale.

Phenotypic variation and genomic variation in insect virulence traits reveal patterns of intraspecific diversity in a locust-specific fungal pathogen

Intraspecific pathogen diversity is crucial for understanding the evolution and maintenance of adaptation in host-pathogen interactions. Traits associated with virulence are often a significant source of variation directly impacted by local selection pressures. The specialist fungal entomopathogen, Metarhizium acridum, has been widely implemented as a biological control agent of locust pests in tropical regions of the world. However, few studies have accounted for natural intraspecific phenotypic and genetic variation. Here, we examine the diversity of nine isolates of M. acridum spanning the known geographic distribution, in terms of (1) virulence towards two locust species, (2) growth rates on three diverse nutrient sources, and (3) comparative genomics to uncover genomic variability. Significant variability in patterns of virulence and growth was shown among the isolates, suggesting intraspecific ecological specialization. Different patterns of virulence were shown between the two locust species, indicative of potential host preference. Additionally, a high level of diversity among M. acridum isolates was observed, revealing increased variation in subtilisin-like proteases from the Pr1 family. These results culminate in the first in-depth analysis regarding multiple facets of natural variation in M. acridum, offering opportunities to understand critical evolutionary drivers of intraspecific diversity in pathogens.

Modeling idiopathic autism in forebrain organoids reveals an imbalance of excitatory cortical neuron subtypes during early neurogenesis

Idiopathic autism spectrum disorder (ASD) is highly heterogeneous, and it remains unclear how convergent biological processes in affected individuals may give rise to symptoms. Here, using cortical organoids and single-cell transcriptomics, we modeled alterations in the forebrain development between boys with idiopathic ASD and their unaffected fathers in 13 families. Transcriptomic changes suggest that ASD pathogenesis in macrocephalic and normocephalic probands involves an opposite disruption of the balance between excitatory neurons of the dorsal cortical plate and other lineages such as early-generated neurons from the putative preplate. The imbalance stemmed from divergent expression of transcription factors driving cell fate during early cortical development. While we did not find genomic variants in probands that explained the observed transcriptomic alterations, a significant overlap between altered transcripts and reported ASD risk genes affected by rare variants suggests a degree of gene convergence between rare forms of ASD and the developmental transcriptome in idiopathic ASD.

Establishment and characterization of patient-derived xenograft of a rare pediatric anaplastic pleomorphic xanthoastrocytoma (PXA) bearing a CDC42SE2-BRAF fusion

Pleomorphic xanthoastrocytoma (PXA) is a rare subset of primary pediatric glioma with 70% 5-year disease free survival. However, up to 20% of cases present with local recurrence and malignant transformation into more aggressive type anaplastic PXA (AXPA) or glioblastoma. The understanding of disease etiology and mechanisms driving PXA and APXA are limited, and there is no standard of care. Therefore, development of relevant preclinical models to investigate molecular underpinnings of disease and to guide novel therapeutic approaches are of interest. Here, for the first time we established, and characterized a patient-derived xenograft (PDX) from a leptomeningeal spread of a patient with recurrent APXA bearing a novel CDC42SE2-BRAF fusion. An integrated -omics analysis was conducted to assess model fidelity of the genomic, transcriptomic, and proteomic/phosphoproteomic landscapes. A stable xenoline was derived directly from the patient recurrent tumor and maintained in 2D and 3D culture systems. Conserved histology features between the PDX and matched APXA specimen were maintained through serial passages. Whole exome sequencing (WES) demonstrated a high degree of conservation in the genomic landscape between PDX and matched human tumor, including small variants (Pearson's r = 0.794-0.839) and tumor mutational burden (~ 3 mutations/MB). Large chromosomal variations including chromosomal gains and losses were preserved in PDX. Notably, chromosomal gain in chromosomes 4-9, 17 and 18 and loss in the short arm of chromosome 9 associated with homozygous 9p21.3 deletion involving CDKN2A/B locus were identified in both patient tumor and PDX sample. Moreover, chromosomal rearrangement involving 7q34 fusion; CDC42SE-BRAF t (5;7) (q31.1, q34) (5:130,721,239, 7:140,482,820) was identified in the PDX tumor, xenoline and matched human tumor. Transcriptomic profile of the patient's tumor was retained in PDX (Pearson r = 0.88) and in xenoline (Pearson r = 0.63) as well as preservation of enriched signaling pathways (FDR Adjusted P < 0.05) including MAPK, EGFR and PI3K/AKT pathways. The multi-omics data of (WES, transcriptome, and reverse phase protein array (RPPA) was integrated to deduce potential actionable pathways for treatment (FDR < 0.05) including KEGG01521, KEGG05202, and KEGG05200. Both xenoline and PDX were resistant to the MEK inhibitors trametinib or mirdametinib at clinically relevant doses, recapitulating the patient's resistance to such treatment in the clinic. This set of APXA models will serve as a preclinical resource for developing novel therapeutic regimens for rare anaplastic PXAs and pediatric high-grade gliomas bearing BRAF fusions.

WGS Revealed Novel BBS5 Pathogenic Variants, Missed by WES, Causing Ciliary Structure and Function Defects

Bardet-Biedl syndrome (BBS) is an autosomal recessive ciliopathy that affects multiple organs, leading to retinitis pigmentosa, polydactyly, obesity, renal anomalies, cognitive impairment, and hypogonadism. Until now, biallelic pathogenic variants have been identified in at least 24 genes delineating the genetic heterogeneity of BBS. Among those, BBS5 is a minor contributor to the mutation load and is one of the eight subunits forming the BBSome, a protein complex implied in protein trafficking within the cilia. This study reports on a European BBS5 patient with a severe BBS phenotype. Genetic analysis was performed using multiple next-generation sequencing (NGS) tests (targeted exome, TES and whole exome, WES), and biallelic pathogenic variants could only be identified using whole-genome sequencing (WGS), including a previously missed large deletion of the first exons. Despite the absence of family samples, the biallelic status of the variants was confirmed. The BBS5 protein's impact was confirmed on the patient's cells (presence/absence and size of the cilium) and ciliary function (Sonic Hedgehog pathway). This study highlights the importance of WGS and the challenge of reliable structural variant detection in patients' genetic explorations as well as functional tests to assess a variant's pathogenicity.

Clonally Selected Lines After CRISPR-Cas Editing Are Not Isogenic

The CRISPR-Cas9 system has enabled researchers to precisely modify/edit the sequence of a genome. A typical editing experiment consists of two steps: (1) editing cultured cells; (2) cell cloning and selection of clones with and without intended edit, presumed to be isogenic. The application of CRISPR-Cas9 system may result in off-target edits, whereas cloning will reveal culture-acquired mutations. We analyzed the extent of the former and the latter by whole genome sequencing in three experiments involving separate genomic loci and conducted by three independent laboratories. In all experiments we hardly found any off-target edits, whereas detecting hundreds to thousands of single nucleotide mutations unique to each clone after relatively short culture of 10-20 passages. Notably, clones also differed in copy number alterations (CNAs) that were several kb to several mb in size and represented the largest source of genomic divergence among clones. We suggest that screening of clones for mutations and CNAs acquired in culture is a necessary step to allow correct interpretation of DNA editing experiments. Furthermore, since culture associated mutations are inevitable, we propose that experiments involving derivation of clonal lines should compare a mix of multiple unedited lines and a mix of multiple edited lines.

Analysis of Genome Structure and Its Variations in Potato Cultivars Grown in Russia

Solanum tuberosum L. (common potato) is one of the most important crops produced almost all over the world. Genomic sequences of potato opens the way for studying the molecular variations related to diversification. We performed a reconstruction of genomic sequences for 15 tetraploid potato cultivars grown in Russia using short reads. Protein-coding genes were identified; conserved and variable parts of pan-genome and the repertoire of the NBS-LRR genes were characterized. For comparison, we used additional genomic sequences for twelve South American potato accessions, performed analysis of genetic diversity, and identified the copy number variations (CNVs) in two these groups of potato. Genomes of Russian potato cultivars were more homogeneous by CNV characteristics and have smaller maximum deletion size in comparison with South American ones. Genes with different CNV occurrences in two these groups of potato accessions were identified. We revealed genes of immune/abiotic stress response, transport and five genes related to tuberization and photoperiod control among them. Four genes related to tuberization and photoperiod were investigated in potatoes previously (phytochrome A among them). A novel gene, homologous to the poly(ADP-ribose) glycohydrolase (PARG) of Arabidopsis, was identified that may be involved in circadian rhythm control and contribute to the acclimatization processes of Russian potato cultivars.

ConanVarvar: a versatile tool for the detection of large syndromic copy number variation from whole-genome sequencing data

BACKGROUND

A wide range of tools are available for the detection of copy number variants (CNVs) from whole-genome sequencing (WGS) data. However, none of them focus on clinically-relevant CNVs, such as those that are associated with known genetic syndromes. Such variants are often large in size, typically 1-5 Mb, but currently available CNV callers have been developed and benchmarked for the discovery of smaller variants. Thus, the ability of these programs to detect tens of real syndromic CNVs remains largely unknown.

RESULTS

Here we present ConanVarvar, a tool which implements a complete workflow for the targeted analysis of large germline CNVs from WGS data. ConanVarvar comes with an intuitive R Shiny graphical user interface and annotates identified variants with information about 56 associated syndromic conditions. We benchmarked ConanVarvar and four other programs on a dataset containing real and simulated syndromic CNVs larger than 1 Mb. In comparison to other tools, ConanVarvar reports 10-30 times less false-positive variants without compromising sensitivity and is quicker to run, especially on large batches of samples.

CONCLUSIONS

ConanVarvar is a useful instrument for primary analysis in disease sequencing studies, where large CNVs could be the cause of disease.

Diagnosis of Prader-Willi syndrome and Angelman syndrome by targeted nanopore long-read sequencing

The CpG island flanking the promoter region of SNRPN on chromosome 15q11.2 contains CpG sites that are completely methylated in the maternally derived allele and unmethylated in the paternally derived allele. Both unmethylated and methylated alleles are observed in normal individuals. Only the methylated allele is observed in patients with Prader-Willi syndrome, whereas only the unmethylated allele is observed in those with Angelman syndrome. Hence, detection of aberrant methylation at the differentially methylated region is fundamental to the molecular diagnosis of Prader-Willi syndrome and Angelman syndromes. Traditionally, bisulfite treatment and methylation-sensitive restriction enzyme treatment or methylation-specific multiplex ligation-dependent probe amplification (MS-MLPA) have been used. We here developed a long-read sequencing assay that can distinguish methylated and unmethylated CpG sites at 15q11.2 by the difference in current intensity generated from nanopore reads. We successfully diagnosed 4 Prader-Willi syndrome patients and 3 Angelman syndrome patients by targeting differentially methylated regions. Concurrent copy number analysis, homozygosity analysis, and structural variant analysis also allowed us to precisely delineate the underlying pathogenic mechanisms, including gross deletion, uniparental heterodisomy, uniparental isodisomy, or imprinting defect. Furthermore, we showed allele-specific methylation in imprinting-related differentially methylated regions on chromosomes 6, 7, 11, 14, and 20 in a normal individual together with 4 Prader-Willi patients and 3 Angelman syndrome patients. Hence, presently reported method is likely to be applicable to the diagnosis of imprinting disorders other than Prader-Willi syndrome and Angelman syndrome as well.

Performance comparisons of methylation and structural variants from low-input whole-genome methylation sequencing

Aim: Whole-genome methylation sequencing carries both DNA methylation and structural variant information (single nucleotide variant [SNV]; copy number variant [CNV]); however, limited data is available on the reliability of obtaining this information simultaneously from low-input DNA using various library preparation and sequencing protocols. Methods: A HapMap NA12878 sample was sequenced with three protocols (EM-sequencing, QIA-sequencing and Swift-sequencing) and their performance was compared on CpG methylation measurement and SNV and CNV detection. Results: At low DNA input (10-25 ng), EM-sequencing was superior in almost all metrics except CNV detection where all protocols were similar. EM-sequencing captured the highest number of CpGs and true SNVs. Conclusion: EM-sequencing is suitable to detect methylation, SNVs and CNVs from single sequencing with low-input DNA.

Chromosome-level assembly of the Phytophthora agathidicida genome reveals adaptation in effector gene families

Phytophthora species are notorious plant pathogens, with some causing devastating tree diseases that threaten the survival of their host species. One such example is Phytophthora agathidicida, the causal agent of kauri dieback - a root and trunk rot disease that kills the ancient, iconic and culturally significant tree species, Agathis australis (New Zealand kauri). A deeper understanding of how Phytophthora pathogens infect their hosts and cause disease is critical for the development of effective treatments. Such an understanding can be gained by interrogating pathogen genomes for effector genes, which are involved in virulence or pathogenicity. Although genome sequencing has become more affordable, the complete assembly of Phytophthora genomes has been problematic, particularly for those with a high abundance of repetitive sequences. Therefore, effector genes located in repetitive regions could be truncated or missed in a fragmented genome assembly. Using a combination of long-read PacBio sequences, chromatin conformation capture (Hi-C) and Illumina short reads, we assembled the P. agathidicida genome into ten complete chromosomes, with a genome size of 57 Mb including 34% repeats. This is the first Phytophthora genome assembled to chromosome level and it reveals a high level of syntenic conservation with the complete genome of Peronospora effusa, the only other completely assembled genome sequence of an oomycete. All P. agathidicida chromosomes have clearly defined centromeres and contain candidate effector genes such as RXLRs and CRNs, but in different proportions, reflecting the presence of gene family clusters. Candidate effector genes are predominantly found in gene-poor, repeat-rich regions of the genome, and in some cases showed a high degree of duplication. Analysis of candidate RXLR effector genes that occur in multicopy gene families indicated half of them were not expressed in planta. Candidate CRN effector gene families showed evidence of transposon-mediated recombination leading to new combinations of protein domains, both within and between chromosomes. Further analysis of this complete genome assembly will help inform new methods of disease control against P. agathidicida and other Phytophthora species, ultimately helping decipher how Phytophthora pathogens have evolved to shape their effector repertoires and how they might adapt in the future.

Genome-wide detection of copy number variation in American mink using whole-genome sequencing

BACKGROUND

Copy number variations (CNVs) represent a major source of genetic diversity and contribute to the phenotypic variation of economically important traits in livestock species. In this study, we report the first genome-wide CNV analysis of American mink using whole-genome sequence data from 100 individuals. The analyses were performed by three complementary software programs including CNVpytor, DELLY and Manta.

RESULTS

A total of 164,733 CNVs (144,517 deletions and 20,216 duplications) were identified representing 5378 CNV regions (CNVR) after merging overlapping CNVs, covering 47.3 Mb (1.9%) of the mink autosomal genome. Gene Ontology and KEGG pathway enrichment analyses of 1391 genes that overlapped CNVR revealed potential role of CNVs in a wide range of biological, molecular and cellular functions, e.g., pathways related to growth (regulation of actin cytoskeleton, and cAMP signaling pathways), behavior (axon guidance, circadian entrainment, and glutamatergic synapse), lipid metabolism (phospholipid binding, sphingolipid metabolism and regulation of lipolysis in adipocytes), and immune response (Wnt signaling, Fc receptor signaling, and GTPase regulator activity pathways). Furthermore, several CNVR-harbored genes associated with fur characteristics and development (MYO5A, RAB27B, FGF12, SLC7A11, EXOC2), and immune system processes (SWAP70, FYN, ORAI1, TRPM2, and FOXO3).

CONCLUSIONS

This study presents the first genome-wide CNV map of American mink. We identified 5378 CNVR in the mink genome and investigated genes that overlapped with CNVR. The results suggest potential links with mink behaviour as well as their possible impact on fur quality and immune response. Overall, the results provide new resources for mink genome analysis, serving as a guideline for future investigations in which genomic structural variations are present.

Prognostic significance of copy number variation in B-cell acute lymphoblastic leukemia

Copy number variations (CNVs) are widespread in both pediatric and adult cases of B-cell acute lymphoblastic leukemia (B-ALL); however, their clinical significance remains unclear. This review primarily discusses the most prevalent CNVs in B-ALL to elucidate their clinical value and further personalized management of this population. The discovery of the molecular mechanism of gene deletion and the development of targeted drugs will further enhance the clinical prognosis of B-ALL.

Analysis of somatic mutations in 131 human brains reveals aging-associated hypermutability

We analyzed 131 human brains (44 neurotypical, 19 with Tourette syndrome, 9 with schizophrenia, and 59 with autism) for somatic mutations after whole genome sequencing to a depth of more than 200×. Typically, brains had 20 to 60 detectable single-nucleotide mutations, but ~6% of brains harbored hundreds of somatic mutations. Hypermutability was associated with age and damaging mutations in genes implicated in cancers and, in some brains, reflected in vivo clonal expansions. Somatic duplications, likely arising during development, were found in ~5% of normal and diseased brains, reflecting background mutagenesis. Brains with autism were associated with mutations creating putative transcription factor binding motifs in enhancer-like regions in the developing brain. The top-ranked affected motifs corresponded to MEIS (myeloid ectopic viral integration site) transcription factors, suggesting a potential link between their involvement in gene regulation and autism.

A wMel Wolbachia variant in Aedes aegypti from field-collected Drosophila melanogaster with increased phenotypic stability under heat stress

Mosquito‐borne diseases remain a major cause of morbidity and mortality. Population replacement strategies involving the wMel strain of Wolbachia are being used widely to control mosquito‐borne diseases. However, these strategies may be influenced by temperature because wMel is vulnerable to heat. wMel infections in Drosophila melanogaster are genetically diverse, but few transinfections of wMel variants have been generated in Aedes aegypti. Here, we successfully transferred a wMel variant (termed wMelM) originating from a field‐collected D. melanogaster into Ae. aegypti. The new wMelM variant (clade I) is genetically distinct from the original wMel transinfection (clade III), and there are no genomic differences between wMelM in its original and transinfected host. We compared wMelM with wMel in its effects on host fitness, temperature tolerance, Wolbachia density, vector competence, cytoplasmic incompatibility and maternal transmission under heat stress in a controlled background. wMelM showed a higher heat tolerance than wMel, likely due to higher overall densities within the mosquito. Both wMel variants had minimal host fitness costs, complete cytoplasmic incompatibility and maternal transmission, and dengue virus blocking under laboratory conditions. Our results highlight phenotypic differences between Wolbachia variants and wMelM shows potential as an alternative strain in areas with strong seasonal temperature fluctuations.

Identification of Candidate Genes Controlling Soybean Cyst Nematode Resistance in "Handou 10" Based on Genome and Transcriptome Analyzes

Soybean cyst nematode (SCN; Heterodera glycines Ichinohe) is a highly destructive pathogen for soybean production worldwide. The use of resistant varieties is the most effective way of preventing yield loss. Handou 10 is a commercial soybean variety with desirable agronomic traits and SCN resistance, however genes underlying the SCN resistance in the variety are unknown. An F_2:8 recombinant inbred line (RIL) population derived from a cross between Zheng 9525 (susceptible) and Handou 10 was developed and its resistance to SCN HG type 2.5.7 (race 1) and 1.2.5.7 (race 2) was identified. We identified seven quantitative trait loci (QTLs) with additive effects. Among these, three QTLs on Chromosomes 7, 8, and 18 were resistant to both races. These QTLs could explain 1.91-7.73% of the phenotypic variation of SCN's female index. The QTLs on chromosomes 8 and 18 have already been reported and were most likely overlapped with rhg1 and Rhg4 loci, respectively. However, the QTL on chromosome 7 was novel. Candidate genes for the three QTLs were predicted through genes functional analysis and transcriptome analysis of infected roots of Handou 10 vs. Zheng 9525. Transcriptome analysis performed also indicated that the plant-pathogen interaction played an important role in the SCN resistance for Handou 10. The information will facilitate SCN-resistant gene cloning, and the novel resistant gene will be a source for improving soybeans' resistance to SCN.

CNVpytor: a tool for copy number variation detection and analysis from read depth and allele imbalance in whole-genome sequencing

BACKGROUND

Detecting copy number variations (CNVs) and copy number alterations (CNAs) based on whole-genome sequencing data is important for personalized genomics and treatment. CNVnator is one of the most popular tools for CNV/CNA discovery and analysis based on read depth.

FINDINGS

Herein, we present an extension of CNVnator developed in Python-CNVpytor. CNVpytor inherits the reimplemented core engine of its predecessor and extends visualization, modularization, performance, and functionality. Additionally, CNVpytor uses B-allele frequency likelihood information from single-nucleotide polymorphisms and small indels data as additional evidence for CNVs/CNAs and as primary information for copy number-neutral losses of heterozygosity.

CONCLUSIONS

CNVpytor is significantly faster than CNVnator-particularly for parsing alignment files (2-20 times faster)-and has (20-50 times) smaller intermediate files. CNV calls can be filtered using several criteria, annotated, and merged over multiple samples. Modular architecture allows it to be used in shared and cloud environments such as Google Colab and Jupyter notebook. Data can be exported into JBrowse, while a lightweight plugin version of CNVpytor for JBrowse enables nearly instant and GUI-assisted analysis of CNVs by any user. CNVpytor release and the source code are available on GitHub at https://github.com/abyzovlab/CNVpytor under the MIT license.