Copy number variation in human health, disease, and evolution

Copy number variations within fibroblast growth factor 13 gene influence growth traits and alternative splicing in cattle

Previous researches revealed a copy number variation (CNV) region in the bovine fibroblast growth factor 13 (FGF13) gene. However, its effects remain unknown. This study detected the various copy number types in seven Chinese cattle breeds and analysed their population genetic characteristics and effects on growth traits and transcription levels. Copy number Loss was more frequent in Caoyuan Red cattle and Xianan cattle than in the other breeds. Association analysis between CNV and growth traits of Qinchuan indicated that the CNV was significantly related to chest depth, hip width and hucklebone width (P < 0.05). Additionally, the growth traits of individuals with copy number Loss were significantly inferior to those with copy number Gain or Median (P < 0.05). Besides, we found two splicing isoforms, AS1 and AS2, in FGF13 gene, which resulted from alternative 5' splicing sites of intron 1. These isoforms showed varied expression levels in various tissues. Moreover, CNV was significantly and negatively associated with the mRNA expression of AS1 (r = -0.525, P < 0.05). The CNVs in bovine FGF13 gene negatively regulated growth traits and gene transcription. These observations provide new insights into bovine FGF13 gene, delivering potentially useful information for future Chinese cattle breeding programs.

Advancing human genotyping: The Infinium HTS iSelect Custom microarray panel (Rita) development study

Single Nucleotide Polymorphisms (SNPs), as the most prevalent type of variation in the human genome, play a pivotal role in influencing human traits. They are extensively utilized in diverse fields such as population genetics, forensic science, and genetic medicine. This study focuses on the 'Rita' BeadChip, a custom SNP microarray panel developed using Illumina Infinium HTS technology. Designed for high-throughput genotyping, the panel facilitates the analysis of over 4000 markers efficiently and cost-effectively. After careful clustering performed on a set of 1000 samples, an evaluation of the Rita panel was undertaken, assessing its sensitivity, repeatability, reproducibility, precision, accuracy, and resistance to contamination. The panel's performance was evaluated in various scenarios, including sex estimation and parental relationship assessment, using GenomeStudio data analysis software. Findings show that over 95 % of the custom BeadChip assay markers were successful, with better performance of transitions over other mutations, and a considerably lower success rate for Y chromosome loci. An exceptional call rate exceeding 99 % was demonstrated for control samples, even with DNA input as low as 0.781 ng. Call rates above 80 % were still obtained with DNA quantities under 0.1 ng, indicating high sensitivity and suitability for forensic applications where DNA quantity is often limited. Repeatability, reproducibility, and precision studies revealed consistency of the panel's performance across different batches and operators, with no significant deviations in call rates or genotyping results. Accuracy assessments, involving comparison with multiple available genetic databases, including the 1000 Genome Project and HapMap, denoted over 99 % concordance, establishing the Rita panel's reliability in genotyping. The contamination study revealed insights into background noise and allowed the definition of thresholds for sample quality evaluation. Multiple metrics for differentiating between negative controls and true samples were highlighted, increasing the reliability of the obtained results. The sex estimation tool in GenomeStudio proved highly effective, correctly assigning sex in all samples with autosomal loci call rates above 97 %. The parental relationship assessment of family trios highlighted the utility of GenomeStudio in identifying genotyping errors or potential Mendelian inconsistencies, promoting the application of arrays such as Rita in kinship testing. Overall, this evaluation confirms the Rita microarray as a robust, high-throughput genotyping tool, underscoring its potential in genetic research and forensic applications. With its custom content and adaptable design, it not only meets current genotyping demands but also opens avenues for further research and application expansion in the field of genetic analysis.

Amplification editing enables efficient and precise duplication of DNA from short sequence to megabase and chromosomal scale

Duplication is a foundation of molecular evolution and a driver of genomic and complex diseases. Here, we develop a genome editing tool named Amplification Editing (AE) that enables programmable DNA duplication with precision at chromosomal scale. AE can duplicate human genomes ranging from 20 bp to 100 Mb, a size comparable to human chromosomes. AE exhibits activity across various cell types, encompassing diploid, haploid, and primary cells. AE exhibited up to 73.0% efficiency for 1 Mb and 3.4% for 100 Mb duplications, respectively. Whole-genome sequencing and deep sequencing of the junctions of edited sequences confirm the precision of duplication. AE can create chromosomal microduplications within disease-relevant regions in embryonic stem cells, indicating its potential for generating cellular and animal models. AE is a precise and efficient tool for chromosomal engineering and DNA duplication, broadening the landscape of precision genome editing from an individual genetic locus to the chromosomal scale.

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.

Development of DNA Insertion-specific Markers Based on the Intergenic Region of Oplegnathus punctatus Cdkn1/srsf3 for Sex Identification

Spotted knifejaw (Oplegnathus punctatus) is a marine economic fish with high food and ecological value, and its growth process has obvious male and female sexual dimorphism, with males growing significantly faster than females. However, the current sex identification technology is not yet mature, which will limit the growth rate of O. punctatus aquaculture and the efficiency of separate sex breeding, so the development of efficient sex molecular markers is imperative. This study identified a 926 bp DNA insertion fragment in the cdkn1/srsf3 intergenic region of O. punctatus males through whole-genome scanning, comparative genomics, and structural variant analysis. A pair of primers was designed based on the insertion information of the Y chromosome intergenic region in male individuals. Agarose gel electrophoresis revealed the amplification of two DNA fragments, 1118 bp and 192 bp, in male O. punctatus individuals. The 926 bp fragment was identified as the insertion in the intergenic region of cdkn1/srsf3 in males, while only a single 192 bp DNA fragment was amplified in females. The biological sex of the individuals identified in this manner was consistent with their known phenotypic sex. In this study, we developed a method to detect DNA insertion variants in the intergenic region of O. punctatus. Additionally, we introduced a new DNA marker for the rapid identification of the sex of O. punctatus, which enhances detection efficiency. The text has important reference significance and application value in sex identification, all-male breeding, and lineage selection. It provides new insights into the regulation of variation in the intergenic region of cdkn1/srsf3 genes and the study of RNA shearing.

Clinical Relevance of the Systematic Analysis of Copy Number Variants in the Genetic Study of Cardiomyopathies

Cardiomyopathies (CMs), one of the main causes of sudden death among the young population, are a heterogeneous group of myocardial diseases, usually with a genetic cause. Next-Generation Sequencing (NGS) has expanded the genes studied for CMs; however, the yield is still around 50%. The systematic study of Copy Number Variants (CNVs) could contribute to improving our diagnostic capacity. These alterations have already been described as responsible for cardiomyopathies in some cases; however, their impact has been rarely assessed. We analyzed the clinical significance of CNVs in cardiomyopathies by studying 11,647 affected patients, many more than those considered in previously published studies. We evaluated the yield of the systematic study of CNVs in a production context using NGS and a novel CNV detection software tool v2.0 that has demonstrated great efficacy, maximizing sensitivity and avoiding false positives. We obtained a CNV analysis yield of 0.8% that fluctuated depending on the type of cardiomyopathy studied (0.29% HCM, 1.41% DCM, 1.88% ARVC, 1.8% LVNC, 1.45% RCM), and we present the frequency of occurrence for 18 genes that agglutinate the 95 pathogenic/likely pathogenic CNVs detected. We conclude the importance of including in diagnostic tests a systematic study of these genetic alterations for the different cardiomyopathies.

Genetic analysis of 106 sporadic cases with hearing loss in the UAE population

BACKGROUND

Hereditary hearing loss is a rare hereditary condition that has a significant presence in consanguineous populations. Despite its prevalence, hearing loss is marked by substantial genetic diversity, which poses challenges for diagnosis and screening, particularly in cases with no clear family history or when the impact of the genetic variant requires functional analysis, such as in the case of missense mutations and UTR variants. The advent of next-generation sequencing (NGS) has transformed the identification of genes and variants linked to various conditions, including hearing loss. However, there remains a high proportion of undiagnosed patients, attributable to various factors, including limitations in sequencing coverage and gaps in our knowledge of the entire genome, among other factors. In this study, our objective was to comprehensively identify the spectrum of genes and variants associated with hearing loss in a cohort of 106 affected individuals from the UAE.

RESULTS

In this study, we investigated 106 sporadic cases of hearing impairment and performed genetic analyses to identify causative mutations. Screening of the GJB2 gene in these cases revealed its involvement in 24 affected individuals, with specific mutations identified. For individuals without GJB2 mutations, whole exome sequencing (WES) was conducted. WES revealed 33 genetic variants, including 6 homozygous and 27 heterozygous DNA changes, two of which were previously implicated in hearing loss, while 25 variants were novel. We also observed multiple potential pathogenic heterozygous variants across different genes in some cases. Notably, a significant proportion of cases remained without potential pathogenic variants.

CONCLUSIONS

Our findings confirm the complex genetic landscape of hearing loss and the limitations of WES in achieving a 100% diagnostic rate, especially in conditions characterized by genetic heterogeneity. These results contribute to our understanding of the genetic basis of hearing loss and emphasize the need for further research and comprehensive genetic analyses to elucidate the underlying causes of this condition.

Identification of novel driver risk genes in CNV loci associated with neurodevelopmental disorders

Copy-number variants (CNVs) are genome-wide structural variations involving the duplication or deletion of large nucleotide sequences. While these types of variations can be commonly found in humans, large and rare CNVs are known to contribute to the development of various neurodevelopmental disorders (NDDs), including autism spectrum disorder (ASD). Nevertheless, given that these NDD-risk CNVs cover broad regions of the genome, it is particularly challenging to pinpoint the critical gene(s) responsible for the manifestation of the phenotype. In this study, we performed a meta-analysis of CNV data from 11,614 affected individuals with NDDs and 4,031 control individuals from SFARI database to identify 41 NDD-risk CNV loci, including 24 novel regions. We also found evidence for dosage-sensitive genes within these regions being significantly enriched for known NDD-risk genes and pathways. In addition, a significant proportion of these genes was found to (1) converge in protein-protein interaction networks, (2) be among most expressed genes in the brain across all developmental stages, and (3) be hit by deletions that are significantly over-transmitted to individuals with ASD within multiplex ASD families from the iHART cohort. Finally, we conducted a burden analysis using 4,281 NDD cases from Decipher and iHART cohorts, and 2,504 neurotypical control individuals from 1000 Genomes and iHART, which resulted in the validation of the association of 162 dosage-sensitive genes driving risk for NDDs, including 22 novel NDD-risk genes. Importantly, most NDD-risk CNV loci entail multiple NDD-risk genes in agreement with a polygenic model associated with the majority of NDD cases.

Precision medicine in colorectal cancer: Leveraging multi-omics, spatial omics, and artificial intelligence

Colorectal cancer (CRC) is a leading cause of cancer-related deaths. Recent advancements in genomic technologies and analytical approaches have revolutionized CRC research, enabling precision medicine. This review highlights the integration of multi-omics, spatial omics, and artificial intelligence (AI) in advancing precision medicine for CRC. Multi-omics approaches have uncovered molecular mechanisms driving CRC progression, while spatial omics have provided insights into the spatial heterogeneity of gene expression in CRC tissues. AI techniques have been utilized to analyze complex datasets, identify new treatment targets, and enhance diagnosis and prognosis. Despite the tumor's heterogeneity and genetic and epigenetic complexity, the fusion of multi-omics, spatial omics, and AI shows the potential to overcome these challenges and advance precision medicine in CRC. The future lies in integrating these technologies to provide deeper insights and enable personalized therapies for CRC patients.

Effects of urban-induced mutations on ecology, evolution and health

Increasing evidence suggests that urbanization is associated with higher mutation rates, which can affect the health and evolution of organisms that inhabit cities. Elevated pollution levels in urban areas can induce DNA damage, leading to de novo mutations. Studies on mutations induced by urban pollution are most prevalent in humans and microorganisms, whereas studies of non-human eukaryotes are rare, even though increased mutation rates have the potential to affect organisms and their populations in contemporary time. Our Perspective explores how higher mutation rates in urban environments could impact the fitness, ecology and evolution of populations. Most mutations will be neutral or deleterious, and higher mutation rates associated with elevated pollution in urban populations can increase the risk of cancer in humans and potentially other species. We highlight the potential for urban-driven increased deleterious mutational loads in some organisms, which could lead to a decline in population growth of a wide diversity of organisms. Although beneficial mutations are expected to be rare, we argue that higher mutation rates in urban areas could influence adaptive evolution, especially in organisms with short generation times. Finally, we explore avenues for future research to better understand the effects of urban-induced mutations on the fitness, ecology and evolution of city-dwelling organisms.

Omics-based molecular classifications empowering in precision oncology

BACKGROUND

In the past decades, cancer enigmatical heterogeneity at distinct expression levels could interpret disparities in therapeutic response and prognosis. It built hindrances to precision medicine, a tactic to tailor customized treatment informed by the tumors' molecular profile. Single-omics analysis dissected the biological features associated with carcinogenesis to some extent but still failed to revolutionize cancer treatment as expected. Integrated omics analysis incorporated tumor biological networks from diverse layers and deciphered a holistic overview of cancer behaviors, yielding precise molecular classification to facilitate the evolution and refinement of precision medicine.

CONCLUSION

This review outlined the biomarkers at multiple expression layers to tutor molecular classification and pinpoint tumor diagnosis, and explored the paradigm shift in precision therapy: from single- to multi-omics-based subtyping to optimize therapeutic regimens. Ultimately, we firmly believe that by parsing molecular characteristics, omics-based typing will be a powerful assistant for precision oncology.

Novel and Recurrent Copy Number Variants in ABCA4-Associated Retinopathy

ABCA4 is the most frequently mutated gene leading to inherited retinal disease (IRD) with over 2200 pathogenic variants reported to date. Of these, ~1% are copy number variants (CNVs) involving the deletion or duplication of genomic regions, typically >50 nucleotides in length. An in-depth assessment of the current literature based on the public database LOVD, regarding the presence of known CNVs and structural variants in ABCA4, and additional sequencing analysis of ABCA4 using single-molecule Molecular Inversion Probes (smMIPs) for 148 probands highlighted recurrent and novel CNVs associated with ABCA4-associated retinopathies. An analysis of the coverage depth in the sequencing data led to the identification of eleven deletions (six novel and five recurrent), three duplications (one novel and two recurrent) and one complex CNV. Of particular interest was the identification of a complex defect, i.e., a 15.3 kb duplicated segment encompassing exon 31 through intron 41 that was inserted at the junction of a downstream 2.7 kb deletion encompassing intron 44 through intron 47. In addition, we identified a 7.0 kb tandem duplication of intron 1 in three cases. The identification of CNVs in ABCA4 can provide patients and their families with a genetic diagnosis whilst expanding our understanding of the complexity of diseases caused by ABCA4 variants.

Genome-wide association analyses reveals copy number variant regions associated with reproduction and disease traits in Canadian Holstein cattle

This study aimed to evaluate the impact of copy number variants (CNVs) on 13 reproduction and 12 disease traits in Holstein cattle. Intensity signal files containing Log R ratio and B allele frequency information from 13,730 Holstein animals genotyped with a 95K SNP panel, and 8,467 Holstein animals genotyped with a 50K SNP panel were used to identify the CNVs. Subsequently, the identified CNVs were validated using whole genome sequence data from 126 animals, resulting in 870 high-confidence CNV regions (CNVRs) on 12,131 animals. Out of these, 54 CNVRs had frequencies higher than or equal to 1% in the population and were used in the genome-wide association analysis (one CNVR at a time, including the G matrix). Results revealed that 4 CNVRs were significantly (p-value < 3.7 × 10-5) associated with at least one of the traits analyzed in this study. Specifically, 2 CNVRs were associated with 3 reproduction traits (i.e., calf survival, first service to conception, and non-return rate), and 2 CNVRs were associated with 2 disease traits (i.e., metritis and retained placenta). These CNVRs harbored genes implicated in immune response, cellular signaling, and neuronal development, supporting their potential involvement in these traits. Further investigations to unravel the mechanistic and functional implications of these CNVRs on the mentioned traits are warranted.

Structural deviations of the posterior fossa and the cerebellum and their cognitive links in a neurodevelopmental deletion syndrome

High-impact genetic variants associated with neurodevelopmental disorders provide biologically-defined entry points for mechanistic investigation. The 3q29 deletion (3q29Del) is one such variant, conferring a 40-100-fold increased risk for schizophrenia, as well as high risk for autism and intellectual disability. However, the mechanisms leading to neurodevelopmental disability remain largely unknown. Here, we report the first in vivo quantitative neuroimaging study in individuals with 3q29Del (N = 24) and neurotypical controls (N = 1608) using structural MRI. Given prior radiology reports of posterior fossa abnormalities in 3q29Del, we focused our investigation on the cerebellum and its tissue-types and lobules. Additionally, we compared the prevalence of cystic/cyst-like malformations of the posterior fossa between 3q29Del and controls and examined the association between neuroanatomical findings and quantitative traits to probe gene-brain-behavior relationships. 3q29Del participants had smaller cerebellar cortex volumes than controls, before and after correction for intracranial volume (ICV). An anterior-posterior gradient emerged in finer grained lobule-based and voxel-wise analyses. 3q29Del participants also had larger cerebellar white matter volumes than controls following ICV-correction and displayed elevated rates of posterior fossa arachnoid cysts and mega cisterna magna findings independent of cerebellar volume. Cerebellar white matter and subregional gray matter volumes were associated with visual-perception and visual-motor integration skills as well as IQ, while cystic/cyst-like malformations yielded no behavioral link. In summary, we find that abnormal development of cerebellar structures may represent neuroimaging-based biomarkers of cognitive and sensorimotor function in 3q29Del, adding to the growing evidence identifying cerebellar pathology as an intersection point between syndromic and idiopathic forms of neurodevelopmental disabilities.

Exome copy number variant detection, analysis, and classification in a large cohort of families with undiagnosed rare genetic disease

Copy number variants (CNVs) are significant contributors to the pathogenicity of rare genetic diseases and, with new innovative methods, can now reliably be identified from exome sequencing. Challenges still remain in accurate classification of CNV pathogenicity. CNV calling using GATK-gCNV was performed on exomes from a cohort of 6,633 families (15,759 individuals) with heterogeneous phenotypes and variable prior genetic testing collected at the Broad Institute Center for Mendelian Genomics of the Genomics Research to Elucidate the Genetics of Rare Diseases consortium and analyzed using the seqr platform. The addition of CNV detection to exome analysis identified causal CNVs for 171 families (2.6%). The estimated sizes of CNVs ranged from 293 bp to 80 Mb. The causal CNVs consisted of 140 deletions, 15 duplications, 3 suspected complex structural variants (SVs), 3 insertions, and 10 complex SVs, the latter two groups being identified by orthogonal confirmation methods. To classify CNV variant pathogenicity, we used the 2020 American College of Medical Genetics and Genomics/ClinGen CNV interpretation standards and developed additional criteria to evaluate allelic and functional data as well as variants on the X chromosome to further advance the framework. We interpreted 151 CNVs as likely pathogenic/pathogenic and 20 CNVs as high-interest variants of uncertain significance. Calling CNVs from existing exome data increases the diagnostic yield for individuals undiagnosed after standard testing approaches, providing a higher-resolution alternative to arrays at a fraction of the cost of genome sequencing. Our improvements to the classification approach advances the systematic framework to assess the pathogenicity of CNVs.

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.

Genetics of ischemic stroke functional outcome

Ischemic stroke, which accounts for 87% of cerebrovascular accidents, is responsible for massive global burden both in terms of economic cost and personal hardship. Many stroke survivors face long-term disability-a phenotype associated with an increasing number of genetic variants. While clinical variables such as stroke severity greatly impact recovery, genetic polymorphisms linked to functional outcome may offer physicians a unique opportunity to deliver personalized care based on their patient's genetic makeup, leading to improved outcomes. A comprehensive catalogue of the variants at play is required for such an approach. In this review, we compile and describe the polymorphisms associated with outcome scores such as modified Rankin Scale and Barthel Index. Our search identified 74 known genetic polymorphisms spread across 48 features associated with various poststroke disability metrics. The known variants span diverse biological systems and are related to inflammation, vascular homeostasis, growth factors, metabolism, the p53 regulatory pathway, and mitochondrial variation. Understanding how these variants influence functional outcome may be helpful in maximizing poststroke recovery.

Copy number variation analysis identifies MIR9-3 and MIR1299 as novel miRNA candidate genes for CAKUT

BACKGROUND

Congenital anomalies of the kidney and urinary tract (CAKUT) represent a frequent cause of pediatric kidney failure. CNVs, as a major class of genomic variations, can also affect miRNA regions. Common CNV corresponding miRNAs (cCNV-miRNAs) are functional variants regulating crucial processes which could affect urinary system development. Thus, we hypothesize that cCNV-miRNAs are associated with CAKUT occurrence and its expressivity.

METHODS

The extraction and filtering of common CNVs, identified in control samples deposited in publicly available databases gnomAD v2.1 and dbVar, were coupled with mapping of miRNA sequences using UCSC Genome Browser. After verification of the mapped miRNAs using referent miRBase V22.1, prioritization of cCNV-miRNA candidates has been performed using bioinformatic annotation and literature research. Genotyping of miRNA gene copy numbers for MIR9-3, MIR511, and MIR1299, was conducted on 221 CAKUT patients and 192 controls using TaqMan™ technology.

RESULTS

We observed significantly different MIR9-3 and MIR1299 gene copy number distribution between CAKUT patients and controls (Chi-square, P = 0.006 and P = 0.0002, respectively), while difference of MIR511 copy number distribution showed nominal significance (Chi-square, P = 0.027). The counts of less and more than two of MIR1299 copy numbers were more frequent within CAKUT patients compared to controls (P = 0.01 and P = 0.008, respectively) and also in cohort of patients with anomalies of the urinary tract compared to controls (P = 0.016 and P = 0.003, respectively).

CONCLUSIONS

Copy number variations of miRNA genes represent a novel avenue in clarification of the inheritance complexity in CAKUT and provide potential evidence about the association of common genetic variation with CAKUT phenotypes.

Detection and characterization of copy number variation in three differentially-selected Nellore cattle populations

Introduction: Nellore cattle (Bos taurus indicus) is the main beef cattle breed raised in Brazil. This breed is well adapted to tropical conditions and, more recently, has experienced intensive genetic selection for multiple performance traits. Over the past 43 years, an experimental breeding program has been developed in the Institute of Animal Science (IZ, Sertaozinho, SP, Brazil), which resulted in three differentially-selected lines known as Nellore Control (NeC), Nellore Selection (NeS), and Nellore Traditional (NeT). The primary goal of this selection experiment was to determine the response to selection for yearling weight (YW) and residual feed intake (RFI) on Nellore cattle. The main objectives of this study were to: 1) identify copy number variation (CNVs) in Nellore cattle from three selection lines; 2) identify and characterize CNV regions (CNVR) on these three lines; and 3) perform functional enrichment analyses of the CNVR identified. Results: A total of 14,914 unique CNVs and 1,884 CNVRs were identified when considering all lines as a single population. The CNVRs were non-uniformly distributed across the chromosomes of the three selection lines included in the study. The NeT line had the highest number of CNVRs (n = 1,493), followed by the NeS (n = 823) and NeC (n = 482) lines. The CNVRs covered 23,449,890 bp (0.94%), 40,175,556 bp (1.61%), and 63,212,273 bp (2.54%) of the genome of the NeC, NeS, and NeT lines, respectively. Two CNVRs were commonly identified between the three lines, and six, two, and four exclusive regions were identified for NeC, NeS, and NeT, respectively. All the exclusive regions overlap with important genes, such as SMARCD3, SLC15A1, and MAPK1. Key biological processes associated with the candidate genes were identified, including pathways related to growth and metabolism. Conclusion: This study revealed large variability in CNVs and CNVRs across three Nellore lines differentially selected for YW and RFI. Gene annotation and gene ontology analyses of the exclusive CNVRs to each line revealed specific genes and biological processes involved in the expression of growth and feed efficiency traits. These findings contribute to the understanding of the genetic mechanisms underlying the phenotypic differences among the three Nellore selection lines.

A novel method for simultaneous detection of hematological tumors and infectious pathogens by metagenomic next generation sequencing of plasma

BACKGROUND

Metagenomic next-generation sequencing (mNGS) is valuable for pathogen identification; however, distinguishing between infectious diseases and conditions with potentially similar clinical manifestations, including malignant tumors, is challenging. Therefore, we developed a method for simultaneous detection of infectious pathogens and cancer in blood samples.

METHODS

Plasma samples (n = 244) were collected from 150 and 94 patients with infections and hematological malignancies, respectively, and analyzed by mNGS for pathogen detection, alongside human tumor chromosomal copy number variation (CNV) analysis (≥5Mbp or 10Mbp CNV region). Further, an evaluation set, comprising 87 plasma samples, was analyzed by mNGS and human CNV analysis, to validate the feasibility of the method.

RESULTS

Among 94 patients with hematological malignancy, sensitivity values of CNV detection for tumor diagnosis were 69.15 % and 32.98 % for CNV region 5Mbp and 10Mbp, respectively, with corresponding specificities of 92.62 % and 100 % in the infection group. Area under the ROC curve (AUC) values for 5Mbp and 10Mbp region were 0.825 and 0.665, respectively, which was a significant difference of 0.160 (95 % CI: 0.110-0.210; p < 0.001), highlighting the superiority of 5Mbp output region data. Six patients with high-risk CNV results were identified in the validation study: three with history of tumor treatment, two eventually newly-diagnosed with hematological malignancies, and one with indeterminate final diagnosis.

CONCLUSIONS

Concurrent CNV analysis alongside mNGS for infection diagnosis is promising for detecting malignant tumors. We recommend adopting a CNV region of 10Mbp over 5Mbp for our model, because of the lower false-positive rate (FPR).

Human embryonic genetic mosaicism and its effects on development and disease

Nearly every mammalian cell division is accompanied by a mutational event that becomes fixed in a daughter cell. When carried forward to additional cell progeny, a clone of variant cells can emerge. As a result, mammals are complex mosaics of clones that are genetically distinct from one another. Recent high-throughput sequencing studies have revealed that mosaicism is common, clone sizes often increase with age and specific variants can affect tissue function and disease development. Variants that are acquired during early embryogenesis are shared by multiple cell types and can affect numerous tissues. Within tissues, variant clones compete, which can result in their expansion or elimination. Embryonic mosaicism has clinical implications for genetic disease severity and transmission but is likely an under-recognized phenomenon. To better understand its implications for mosaic individuals, it is essential to leverage research tools that can elucidate the mechanisms by which expanded embryonic variants influence development and disease.

(Zebra)fishing for nephrogenesis genes

Kidney disease is a devastating condition affecting millions of people worldwide, where over 100,000 patients in the United States alone remain waiting for a lifesaving organ transplant. Concomitant with a surge in personalized medicine, single-gene mutations, and polygenic risk alleles have been brought to the forefront as core causes of a spectrum of renal disorders. With the increasing prevalence of kidney disease, it is imperative to make substantial strides in the field of kidney genetics. Nephrons, the core functional units of the kidney, are epithelial tubules that act as gatekeepers of body homeostasis by absorbing and secreting ions, water, and small molecules to filter the blood. Each nephron contains a series of proximal and distal segments with explicit metabolic functions. The embryonic zebrafish provides an ideal platform to systematically dissect the genetic cues governing kidney development. Here, we review the use of zebrafish to discover nephrogenesis genes.

Integrated genomic and transcriptomic analysis reveals the activation of PI3K signaling pathway in HPV-independent cervical cancers

BACKGROUND

HPV-independent cervical cancers (HPV-ind CCs) are uncommon with worse prognosis and poorly understood. This study investigated the molecular characteristics of HPV-ind CCs, aiming to explore new strategies for HPV-ind CCs.

METHODS

HPV status of 1010 cervical cancer patients were detected by RT-PCR, PCR and RNA-sequencing (RNA-seq). Whole exome sequencing (WES) and RNA-seq were performed in identified HPV-ind CCs. The efficacy of PI3Kα inhibitor BYL719 in HPV-ind CCs was evaluated in cell lines, patient-derived organoids (PDOs) and patient-derived xenografts (PDXs).

RESULTS

Twenty-five CCs were identified as HPV-ind, which were more common seen in older, adenocarcinoma patients and exhibited poorer prognosis as well as higher tumor mutation burden compared to HPV-associated CCs. HPV-ind CCs were featured with highly activated PI3K/AKT signaling pathway, particularly, PIK3CA being the most predominant genomic alteration (36%). BYL719 demonstrated superior tumor suppression in vitro and in vivo. Furthermore, HPV-ind CCs were classified into two subtypes according to distinct prognosis by gene expression profiles, the metabolism subtype and immune subtype.

CONCLUSIONS

This study reveals the prevalence, clinicopathology, and molecular features of HPV-ind CCs and emphasizes the importance of PIK3CA mutations and PI3K pathway activation in tumorigenesis, which suggests the potential significance of PI3Kα inhibitors in HPV-ind CC patients.

Validation of low-pass genome sequencing for prenatal diagnosis

OBJECTIVE

Chromosomal microarray (CMA), while considered the gold standard for detecting copy number variants (CNVs) in prenatal diagnostics, has its limitations, including the necessity to replace aging microarray equipment, low throughput, a static design, and an inefficient multi-day workflow. This study evaluates the feasibility of low-pass genome sequencing (LP-GS) as a potential replacement for CMA in prenatal diagnostics.

METHODS

We comprehensively compared LP-GS at 10x and 5x average depths with CMA in a prenatal laboratory. We examined parameters, including concordance, sensitivity, specificity, workflow efficiency, and cost-effectiveness.

RESULTS

We found a high degree of agreement between LP-GS and CMA for detecting CNVs and absence of heterozygosity. Furthermore, compared to CMA, LP-GS increased workflow efficiency and proved to be cost-neutral at 10x and cost-effective at 5x.

CONCLUSION

Our study suggests that LP-GS is a promising alternative to CMA in prenatal diagnostics, offering advantages, including a more efficient workflow and scalability for larger testing volumes. Importantly, for clinical laboratories that have adopted next-generation sequencing in a separate capacity, LP-GS facilitates a unified NGS-centric approach, enabling workflow consolidation. By offering a single, streamlined platform for detecting a broad range of genetic variants, LP-GS may represent a critical step toward enhancing the diagnostic capabilities of prenatal laboratories.

Dynamic interplay between human alpha-satellite DNA structure and centromere functions

Maintenance of genome stability relies on functional centromeres for correct chromosome segregation and faithful inheritance of the genetic information. The human centromere is the primary constriction within mitotic chromosomes made up of repetitive alpha-satellite DNA hierarchically organized in megabase-long arrays of near-identical higher order repeats (HORs). Centromeres are epigenetically specified by the presence of the centromere-specific histone H3 variant, CENP-A, which enables the assembly of the kinetochore for microtubule attachment. Notably, centromeric DNA is faithfully inherited as intact haplotypes from the parents to the offspring without intervening recombination, yet, outside of meiosis, centromeres are akin to common fragile sites (CFSs), manifesting crossing-overs and ongoing sequence instability. Consequences of DNA changes within the centromere are just starting to emerge, with unclear effects on intra- and inter-generational inheritance driven by centromere's essential role in kinetochore assembly. Here, we review evidence of meiotic selection operating to mitigate centromere drive, as well as recent reports on centromere damage, recombination and repair during the mitotic cell division. We propose an antagonistic pleiotropy interpretation to reconcile centromere DNA instability as both driver of aneuploidy that underlies degenerative diseases, while also potentially necessary for the maintenance of homogenized HORs for centromere function. We attempt to provide a framework for this conceptual leap taking into consideration the structural interface of centromere-kinetochore interaction and present case scenarios for its malfunctioning. Finally, we offer an integrated working model to connect DNA instability, chromatin, and structural changes with functional consequences on chromosome integrity.

Optical Genome Mapping Reveals Genomic Alterations upon Gene Editing in hiPSCs: Implications for Neural Tissue Differentiation and Brain Organoid Research

Genome editing, notably CRISPR (cluster regularly interspaced short palindromic repeats)/Cas9 (CRISPR-associated protein 9), has revolutionized genetic engineering allowing for precise targeted modifications. This technique's combination with human induced pluripotent stem cells (hiPSCs) is a particularly valuable tool in cerebral organoid (CO) research. In this study, CRISPR/Cas9-generated fluorescently labeled hiPSCs exhibited no significant morphological or growth rate differences compared with unedited controls. However, genomic aberrations during gene editing necessitate efficient genome integrity assessment methods. Optical genome mapping, a high-resolution genome-wide technique, revealed genomic alterations, including chromosomal copy number gain and losses affecting numerous genes. Despite these genomic alterations, hiPSCs retain their pluripotency and capacity to generate COs without major phenotypic changes but one edited cell line showed potential neuroectodermal differentiation impairment. Thus, this study highlights optical genome mapping in assessing genome integrity in CRISPR/Cas9-edited hiPSCs emphasizing the need for comprehensive integration of genomic and morphological analysis to ensure the robustness of hiPSC-based models in cerebral organoid research.

A Pilot Study of Multiplex Ligation-Dependent Probe Amplification Evaluation of Copy Number Variations in Romanian Children with Congenital Heart Defects

Congenital heart defects (CHDs) have had an increasing prevalence over the last decades, being one of the most common congenital defects. Their etiopathogenesis is multifactorial in origin. About 10-15% of all CHD can be attributed to copy number variations (CNVs), a type of submicroscopic structural genetic alterations. The aim of this study was to evaluate the involvement of CNVs in the development of congenital heart defects. We performed a cohort study investigating the presence of CNVs in the 22q11.2 region and GATA4, TBX5, NKX2-5, BMP4, and CRELD1 genes in patients with syndromic and isolated CHDs. A total of 56 patients were included in the study, half of them (28 subjects) being classified as syndromic. The most common heart defect in our study population was ventricular septal defect (VSD) at 39.28%. There were no statistically significant differences between the two groups in terms of CHD-type distribution, demographical, and clinical features, with the exceptions of birth length, weight, and length at the time of blood sampling, that were significantly lower in the syndromic group. Through multiplex ligation-dependent probe amplification (MLPA) analysis, we found two heterozygous deletions in the 22q11.2 region, both in patients from the syndromic group. No CNVs involving GATA4, NKX2-5, TBX5, BMP4, and CRELD1 genes were identified in our study. We conclude that the MLPA assay may be used as a first genetic test in patients with syndromic CHD and that the 22q11.2 region may be included in the panels used for screening these patients.

Progress in single-cell multimodal sequencing and multi-omics data integration

With the rapid advance of single-cell sequencing technology, cell heterogeneity in various biological processes was dissected at different omics levels. However, single-cell mono-omics results in fragmentation of information and could not provide complete cell states. In the past several years, a variety of single-cell multimodal omics technologies have been developed to jointly profile multiple molecular modalities, including genome, transcriptome, epigenome, and proteome, from the same single cell. With the availability of single-cell multimodal omics data, we can simultaneously investigate the effects of genomic mutation or epigenetic modification on transcription and translation, and reveal the potential mechanisms underlying disease pathogenesis. Driven by the massive single-cell omics data, the integration method of single-cell multi-omics data has rapidly developed. Integration of the massive multi-omics single-cell data in public databases in the future will make it possible to construct a cell atlas of multi-omics, enabling us to comprehensively understand cell state and gene regulation at single-cell resolution. In this review, we summarized the experimental methods for single-cell multimodal omics data and computational methods for multi-omics data integration. We also discussed the future development of this field.

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.

On the core segmentation algorithms of copy number variation detection tools

Shotgun sequencing is a high-throughput method used to detect copy number variants (CNVs). Although there are numerous CNV detection tools based on shotgun sequencing, their quality varies significantly, leading to performance discrepancies. Therefore, we conducted a comprehensive analysis of next-generation sequencing-based CNV detection tools over the past decade. Our findings revealed that the majority of mainstream tools employ similar detection rationale: calculates the so-called read depth signal from aligned sequencing reads and then segments the signal by utilizing either circular binary segmentation (CBS) or hidden Markov model (HMM). Hence, we compared the performance of those two core segmentation algorithms in CNV detection, considering varying sequencing depths, segment lengths and complex types of CNVs. To ensure a fair comparison, we designed a parametrical model using mainstream statistical distributions, which allows for pre-excluding bias correction such as guanine-cytosine (GC) content during the preprocessing step. The results indicate the following key points: (1) Under ideal conditions, CBS demonstrates high precision, while HMM exhibits a high recall rate. (2) For practical conditions, HMM is advantageous at lower sequencing depths, while CBS is more competitive in detecting small variant segments compared to HMM. (3) In case involving complex CNVs resembling real sequencing, HMM demonstrates more robustness compared with CBS. (4) When facing large-scale sequencing data, HMM costs less time compared with the CBS, while their memory usage is approximately equal. This can provide an important guidance and reference for researchers to develop new tools for CNV detection.

Advancements in copy number variation screening in herbivorous livestock genomes and their association with phenotypic traits

Copy number variations (CNVs) have garnered increasing attention within the realm of genetics due to their prevalence in human, animal, and plant genomes. These structural genetic variations have demonstrated associations with a broad spectrum of phenotypic diversity, economic traits, environmental adaptations, epidemics, and other essential aspects of both plants and animals. Furthermore, CNVs exhibit extensive sequence variability and encompass a wide array of genomes. The advancement and maturity of microarray and sequencing technologies have catalyzed a surge in research endeavors pertaining to CNVs. This is particularly prominent in the context of livestock breeding, where molecular markers have gained prominence as a valuable tool in comparison to traditional breeding methods. In light of these developments, a contemporary and comprehensive review of existing studies on CNVs becomes imperative. This review serves the purpose of providing a brief elucidation of the fundamental concepts underlying CNVs, their mutational mechanisms, and the diverse array of detection methods employed to identify these structural variations within genomes. Furthermore, it seeks to systematically analyze the recent advancements and findings within the field of CNV research, specifically within the genomes of herbivorous livestock species, including cattle, sheep, horses, and donkeys. The review also highlighted the role of CNVs in shaping various phenotypic traits including growth traits, reproductive traits, pigmentation and disease resistance etc., in herbivorous livestock. The main goal of this review is to furnish readers with an up-to-date compilation of knowledge regarding CNVs in herbivorous livestock genomes. By integrating the latest research findings and insights, it is anticipated that this review will not only offer pertinent information but also stimulate future investigations into the realm of CNVs in livestock. In doing so, it endeavors to contribute to the enhancement of breeding strategies, genomic selection, and the overall improvement of herbivorous livestock production and resistance to diseases.

Rare copy-number variants as modulators of common disease susceptibility

BACKGROUND

Copy-number variations (CNVs) have been associated with rare and debilitating genomic disorders (GDs) but their impact on health later in life in the general population remains poorly described.

METHODS

Assessing four modes of CNV action, we performed genome-wide association scans (GWASs) between the copy-number of CNV-proxy probes and 60 curated ICD-10 based clinical diagnoses in 331,522 unrelated white British UK Biobank (UKBB) participants with replication in the Estonian Biobank.

RESULTS

We identified 73 signals involving 40 diseases, all of which indicating that CNVs increased disease risk and caused earlier onset. We estimated that 16% of these associations are indirect, acting by increasing body mass index (BMI). Signals mapped to 45 unique, non-overlapping regions, nine of which being linked to known GDs. Number and identity of genes affected by CNVs modulated their pathogenicity, with many associations being supported by colocalization with both common and rare single-nucleotide variant association signals. Dissection of association signals provided insights into the epidemiology of known gene-disease pairs (e.g., deletions in BRCA1 and LDLR increased risk for ovarian cancer and ischemic heart disease, respectively), clarified dosage mechanisms of action (e.g., both increased and decreased dosage of 17q12 impacted renal health), and identified putative causal genes (e.g., ABCC6 for kidney stones). Characterization of the pleiotropic pathological consequences of recurrent CNVs at 15q13, 16p13.11, 16p12.2, and 22q11.2 in adulthood indicated variable expressivity of these regions and the involvement of multiple genes. Finally, we show that while the total burden of rare CNVs-and especially deletions-strongly associated with disease risk, it only accounted for ~ 0.02% of the UKBB disease burden. These associations are mainly driven by CNVs at known GD CNV regions, whose pleiotropic effect on common diseases was broader than anticipated by our CNV-GWAS.

CONCLUSIONS

Our results shed light on the prominent role of rare CNVs in determining common disease susceptibility within the general population and provide actionable insights for anticipating later-onset comorbidities in carriers of recurrent CNVs.

Template switching between the leading and lagging strands at replication forks generates inverted copy number variants through hairpin-capped extrachromosomal DNA

Inherited and germ-line de novo copy number variants (CNVs) are increasingly found to be correlated with human developmental and cancerous phenotypes. Several models for template switching during replication have been proposed to explain the generation of these gross chromosomal rearrangements. We proposed a model of template switching (ODIRA-origin dependent inverted repeat amplification) in which simultaneous ligation of the leading and lagging strands at diverging replication forks could generate segmental inverted triplications through an extrachromosomal inverted circular intermediate. Here, we created a genetic assay using split-ura3 cassettes to trap the proposed inverted intermediate. However, instead of recovering circular inverted intermediates, we found inverted linear chromosomal fragments ending in native telomeres-suggesting that a template switch had occurred at the centromere-proximal fork of a replication bubble. As telomeric inverted hairpin fragments can also be created through double strand breaks we tested whether replication errors or repair of double stranded DNA breaks were the most likely initiating event. The results from CRISPR/Cas9 cleavage experiments and growth in the replication inhibitor hydroxyurea indicate that it is a replication error, not a double stranded break that creates the inverted junctions. Since inverted amplicons of the SUL1 gene occur during long-term growth in sulfate-limited chemostats, we sequenced evolved populations to look for evidence of linear intermediates formed by an error in replication. All of the data are compatible with a two-step version of the ODIRA model in which sequential template switching at short inverted repeats between the leading and lagging strands at a replication fork, followed by integration via homologous recombination, generates inverted interstitial triplications.

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.

Small Indels in the Androgen Receptor Gene: Phenotype Implications and Mechanisms of Mutagenesis

CONTEXT

Despite high abundance of small indels in human genomes, their precise roles and underlying mechanisms of mutagenesis in Mendelian disorders require further investigation.

OBJECTIVE

To profile the distribution, functional implications, and mechanisms of small indels in the androgen receptor (AR) gene in individuals with androgen insensitivity syndrome (AIS).

METHODS

We conducted a systematic review of previously reported indels within the coding region of the AR gene, including 3 novel indels. Distribution throughout the AR coding region was examined and compared with genomic population data. Additionally, we assessed their impact on the AIS phenotype and investigated potential mechanisms driving their occurrence.

RESULTS

A total of 82 indels in AIS were included. Notably, all frameshift indels exhibited complete AIS. The distribution of indels across the AR gene showed a predominance in the N-terminal domain, most leading to frameshift mutations. Small deletions accounted for 59.7%. Most indels occurred in nonrepetitive sequences, with 15.8% situated within triplet regions. Gene burden analysis demonstrated significant enrichment of frameshift indels in AIS compared with controls (P < .00001), and deletions were overrepresented in AIS (P < .00001).

CONCLUSION

Our findings underscore a robust genotype-phenotype relationship regarding small indels in the AR gene in AIS, with a vast majority presenting complete AIS. Triplet regions and homopolymeric runs emerged as prone loci for small indels within the AR. Most were frameshift indels, with polymerase slippage potentially explaining half of AR indel occurrences. Complex frameshift indels exhibited association with palindromic runs. These discoveries advance understanding of the genetic basis of AIS and shed light on potential mechanisms underlying pathogenic small indel events.

Comparative genomic diversity analysis of copy number variations (CNV) in indicine and taurine cattle thriving in Europe and Indian subcontinent

Copy number variations (CNVs) include deletions, duplications, and insertions that are larger than 50 bp in size causing structural variation responsible for diversity, adaptation, and breed development. Indian cattle breeds are highly diverse from the taurine breeds. The pattern of CNVRs in 191 animals belonging to 39 cattle breeds (four Indicine and 35 Taurine) was studied based on Illumina 777K BovineHD chip data. The Indicine breeds revealed 2590 CNVs and 335 copy number variation regions (CNVRs) in autosomes. Out of the identified CNVs, 50 were found to be novel. Structure analysis revealed admixed nature of Siri. Neighbor joining tree from CNVR data showed that hot (Kankrej and Hallikar) and cold (Ladakhi and Siri) adapted cattle breeds clustered separately. CNVR of Indian and European breeds revealed that Balkan and Italian breeds of Podolian group are admixed with Indian cattle breeds corroborating indicine introgression (6.1-13.5%). CNVRs spanning the regions of olfactory receptors and immune system genes were identified. AMOVA revealed 9% variation among populations which is 2% greater than SNP based studies showing higher inclusion of variation by CNVR. Detailed analysis of CNVs/CNVRs in Indian cattle adapted to hot and cold climate, and their diversity among worldwide cattle is presented in this study.

High population frequencies of MICA copy number variations originate from independent recombination events

MICA is a stress-induced ligand of the NKG2D receptor that stimulates NK and T cell responses and was identified as a key determinant of anti-tumor immunity. The MICA gene is located inside the MHC complex and is in strong linkage disequilibrium with HLA-B. While an HLA-B*48-linked MICA deletion-haplotype was previously described in Asian populations, little is known about other MICA copy number variations. Here, we report the genotyping of more than two million individuals revealing high frequencies of MICA duplications (1%) and MICA deletions (0.4%). Their prevalence differs between ethnic groups and can rise to 2.8% (Croatia) and 9.2% (Mexico), respectively. Targeted sequencing of more than 70 samples indicates that these copy number variations originate from independent nonallelic homologous recombination events between segmental duplications upstream of MICA and MICB. Overall, our data warrant further investigation of disease associations and consideration of MICA copy number data in oncological study protocols.

Rare Genomic Copy Number Variants Implicate New Candidate Genes for Bicuspid Aortic Valve

Bicuspid aortic valve (BAV), the most common congenital heart defect, is a major cause of aortic valve disease requiring valve interventions and thoracic aortic aneurysms predisposing to acute aortic dissections. The spectrum of BAV ranges from early onset valve and aortic complications (EBAV) to sporadic late onset disease. Rare genomic copy number variants (CNVs) have previously been implicated in the development of BAV and thoracic aortic aneurysms. We determined the frequency and gene content of rare CNVs in EBAV probands (n = 272) using genome-wide SNP microarray analysis and three complementary CNV detection algorithms (cnvPartition, PennCNV, and QuantiSNP). Unselected control genotypes from the Database of Genotypes and Phenotypes were analyzed using identical methods. We filtered the data to select large genic CNVs that were detected by multiple algorithms. Findings were replicated in cohorts with late onset sporadic disease (n = 5040). We identified 34 large and rare (< 1:1000 in controls) CNVs in EBAV probands. The burden of CNVs intersecting with genes known to cause BAV when mutated was increased in case-control analysis. CNVs intersecting with GATA4 and DSCAM were enriched in cases, recurrent in other datasets, and segregated with disease in families. In total, we identified potentially pathogenic CNVs in 8% of EBAV cases, implicating alterations of candidate genes at these loci in the pathogenesis of BAV.

Genome-wide identification of candidate copy number polymorphism genes associated with complex traits of Tibetan-sheep

Copy number variation (CNV) is a genetic structural polymorphism important for phenotypic diversity and important economic traits of livestock breeds, and it plays an important role in the desired genetic variation. This study used whole genome sequencing to detect the CNV variation in the genome of 6 local Tibetan sheep groups. We detected 69,166 CNV events and 7230 copy number variable regions (CNVRs) after merging the overlapping CNVs, accounting for 2.72% of the reference genome. The CNVR length detected ranged from 1.1 to 1693.5 Kb, with a total length of 118.69 Mb and an average length of 16.42 Kb per CNVR. Functional GO cluster analysis showed that the CNVR genes were mainly involved in sensory perception systems, response to stimulus, and signal transduction. Through CNVR-based Vst analysis, we found that the CACNA2D3 and CTBP1 genes related to hypoxia adaptation, the HTR1A gene related to coat color, and the TRNAS-GGA and PIK3C3 genes related to body weight were all strongly selected. The findings of our study will contribute novel insights into the genetic structural variation underlying hypoxia adaptation and economically important traits in Tibetan sheep.

Exome copy number variant detection, analysis and classification in a large cohort of families with undiagnosed rare genetic disease

Copy number variants (CNVs) are significant contributors to the pathogenicity of rare genetic diseases and with new innovative methods can now reliably be identified from exome sequencing. Challenges still remain in accurate classification of CNV pathogenicity. CNV calling using GATK-gCNV was performed on exomes from a cohort of 6,633 families (15,759 individuals) with heterogeneous phenotypes and variable prior genetic testing collected at the Broad Institute Center for Mendelian Genomics of the GREGoR consortium. Each family's CNV data was analyzed using the seqr platform and candidate CNVs classified using the 2020 ACMG/ClinGen CNV interpretation standards. We developed additional evidence criteria to address situations not covered by the current standards. The addition of CNV calling to exome analysis identified causal CNVs for 173 families (2.6%). The estimated sizes of CNVs ranged from 293 bp to 80 Mb with estimates that 44% would not have been detected by standard chromosomal microarrays. The causal CNVs consisted of 141 deletions, 15 duplications, 4 suspected complex structural variants (SVs), 3 insertions and 10 complex SVs, the latter two groups being identified by orthogonal validation methods. We interpreted 153 CNVs as likely pathogenic/pathogenic and 20 CNVs as high interest variants of uncertain significance. Calling CNVs from existing exome data increases the diagnostic yield for individuals undiagnosed after standard testing approaches, providing a higher resolution alternative to arrays at a fraction of the cost of genome sequencing. Our improvements to the classification approach advances the systematic framework to assess the pathogenicity of CNVs.

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.

Integrating whole-genome sequencing and transcriptomic findings in the diagnosis and management of Coffin-Siris syndrome

INTRODUCTION

Although the whole-exome sequencing (WES) approach has been widely used in clinic, many rare diseases with syndromic and nonsyndromic neurological manifestations remain undiagnosed. Coffin-Siris syndrome (CSS) is a rare autosomal dominant genetic disease characterized by neurodevelopmental delay. A suspected diagnosis can be made based on the typical CSS clinical features; however, molecular genetic testing is necessary for a confirmed diagnosis.

OBJECTIVES

Three CSS-like patients with negative results in the WES and chromosomal microarray analysis (CMA) were recruited in this study.

METHODS

We used whole-genome sequencing (WGS) technology to sequence the peripheral blood of the three families. To further explore the possible pathogenesis of CSS, we performed RNA-sequencing (RNA-seq).

RESULTS

WGS identified the three CSS patients were carrying de novo copy number variants of the ARID1B gene, which have not been reported before. RNA-seq identified 184 differentially expressed genes (DEGs), with 116 up-regulated and 68 down-regulated. Functional annotation of DEGs showed that two biological processes (immune response, chemokine activity) and two signaling pathways (cytokine-cytokine receptor interaction, chemokine activity) were highlighted. We speculated that ARID1B deficiency might trigger abnormal immune responses, which may be involved in the pathophysiologic mechanisms of CSS.

CONCLUSION

Our research provided further support for WGS application in CSS diagnosis and made an investigational approach for the underlying mechanisms of CSS.

Regulatory processes that maintain or alter ribosomal DNA stability during the repair of programmed DNA double-strand breaks

Organisms have evolved elaborate mechanisms that maintain genome stability. Deficiencies in these mechanisms result in changes to the nucleotide sequence as well as copy number and structural variations in the genome. Genome instability has been implicated in numerous human diseases. However, genomic alterations can also be beneficial as they are an essential part of the evolutionary process. Organisms sometimes program genomic changes that drive genetic and phenotypic diversity. Therefore, genome alterations can have both positive and negative impacts on cellular growth and functions, which underscores the need to control the processes that restrict or induce such changes to the genome. The ribosomal RNA gene (rDNA) is highly abundant in eukaryotic genomes, forming a cluster where numerous rDNA copies are tandemly arrayed. Budding yeast can alter the stability of its rDNA cluster by changing the rDNA copy number within the cluster or by producing extrachromosomal rDNA circles. Here, we review the mechanisms that regulate the stability of the budding yeast rDNA cluster during repair of DNA double-strand breaks that are formed in response to programmed DNA replication fork arrest.

Simulating Space Conditions Evokes Different DNA Damage Responses in Immature and Mature Cells of the Human Hematopoietic System

The impact of space radiation and microgravity on DNA damage responses has been discussed controversially, largely due to the variety of model systems engaged. Here, we performed side-by-side analyses of human hematopoietic stem/progenitor cells (HSPC) and peripheral blood lymphocytes (PBL) cultivated in a 2D clinostat to simulate microgravity before, during and after photon and particle irradiation. We demonstrate that simulated microgravity (SMG) accelerates the early phase of non-homologous end joining (NHEJ)-mediated repair of simple, X-ray-induced DNA double-strand breaks (DSBs) in PBL, while repair kinetics in HSPC remained unaltered. Repair acceleration was lost with increasing LET of ion exposures, which increases the complexity of DSBs, precluding NHEJ and requiring end resection for successful repair. Such cell-type specific effect of SMG on DSB repair was dependent on the NF-кB pathway pre-activated in PBL but not HSPC. Already under unperturbed growth conditions HSPC and PBL suffered from SMG-induced replication stress associated with accumulation of single-stranded DNA and DSBs, respectively. We conclude that in PBL, SMG-induced DSBs promote repair of radiation-induced damage in an adaptive-like response. HSPC feature SMG-induced single-stranded DNA and FANCD2 foci, i.e., markers of persistent replication stress and senescence that may contribute to a premature decline of the immune system in space.

Copy Number Variations in Neuropsychiatric Disorders

Neuropsychiatric disorders are complex conditions that represent a significant global health burden with complex and multifactorial etiologies. Technological advances in recent years have improved our understanding of the genetic architecture of the major neuropsychiatric disorders and the genetic loci involved. Previous studies mainly investigated genome-wide significant SNPs to elucidate the cross-disorder and disorder-specific genetic basis of neuropsychiatric disorders. Although copy number variations represent a major source of genetic variations, they are known risk factors in developing a variety of human disorders, including certain neuropsychiatric diseases. In this review, we demonstrate the current understanding of CNVs contributing to liability for schizophrenia, bipolar disorder, and major depressive disorder.

Copy number variation as a tool for implementing pregnancy as an aging model

Copy number variations (CNV) are a major contributor to genome variability and have been linked to aging and other degradable phenotypes such as pregnancy physiology. To demonstrate how pregnancy can be used as a model of aging, we used CNVs from pregnant mice. Candidate CNVs were selected by applying case-control analysis in human centenarians compared with control groups. These CNVs were aligned with the mouse genome and their copy variation was assessed using qRT-PCR in liver and blood tissue samples from pregnant mice throughout pregnancy (baseline; first, second, and third trimester; post-partum). Eight of the ten selected CNVs demonstrated a significant decline/increase trend throughout the pregnancy followed by opposite direction soon after delivery in the liver and blood of the mouse tissues. Furthermore, significant differential expression was detected among the candidate CNVs' close vicinity genes (APA2A, LSS, RBDHF1, PLAAT1, and SCL17A2), but not in the WSCD2 gene. Establishing a genetic link between longevity and pregnancy is a significant step toward implementing the pregnancy process as a model for aging. These results in pregnant mice highlight the mechanism and similarities between pregnancy and aging. Investigating the mechanisms that cause such rejuvenation after labor could change our aging treatment paradigm.

3' UTR Deletion of FBXO28 in a Patient with Brain Abnormalities and Developmental Delay

Constitutional deletions of chromosome 1q42 region are rare. The phenotype spectrum associated with this copy number change is variable, including developmental delay, intellectual disability, seizures, and dysmorphology. This study describes a patient with developmental delays and brain abnormalities. G-banded karyotype, FISH, SNP oligonucleotide microarray analysis (SOMA), and whole exome sequencing analysis were performed. Postnatal reanalysis of prenatal SOMA and follow-up parental testing revealed a paternally inherited 63 kb deletion at 1q42.11 in the patient. We characterized the clinical features of this patient, providing insight into the clinical phenotype associated with deletions of the 1q42.11 sub-band. Our study provides new evidence supporting the potential functional importance of the FBXO28 3' UTR region and the hypothesis that FBXO28 is a critical gene in the pathogenesis of chromosome 1q41q42 microdeletion syndrome. It also highlights the different goals and reporting criteria between prenatal and postnatal microarray tests.

Genome sequencing of 2000 canids by the Dog10K consortium advances the understanding of demography, genome function and architecture

BACKGROUND

The international Dog10K project aims to sequence and analyze several thousand canine genomes. Incorporating 20 × data from 1987 individuals, including 1611 dogs (321 breeds), 309 village dogs, 63 wolves, and four coyotes, we identify genomic variation across the canid family, setting the stage for detailed studies of domestication, behavior, morphology, disease susceptibility, and genome architecture and function.

RESULTS

We report the analysis of > 48 M single-nucleotide, indel, and structural variants spanning the autosomes, X chromosome, and mitochondria. We discover more than 75% of variation for 239 sampled breeds. Allele sharing analysis indicates that 94.9% of breeds form monophyletic clusters and 25 major clades. German Shepherd Dogs and related breeds show the highest allele sharing with independent breeds from multiple clades. On average, each breed dog differs from the UU_Cfam_GSD_1.0 reference at 26,960 deletions and 14,034 insertions greater than 50 bp, with wolves having 14% more variants. Discovered variants include retrogene insertions from 926 parent genes. To aid functional prioritization, single-nucleotide variants were annotated with SnpEff and Zoonomia phyloP constraint scores. Constrained positions were negatively correlated with allele frequency. Finally, the utility of the Dog10K data as an imputation reference panel is assessed, generating high-confidence calls across varied genotyping platform densities including for breeds not included in the Dog10K collection.

CONCLUSIONS

We have developed a dense dataset of 1987 sequenced canids that reveals patterns of allele sharing, identifies likely functional variants, informs breed structure, and enables accurate imputation. Dog10K data are publicly available.

Exome sequencing identified rare recurrent copy number variants and hereditary breast cancer susceptibility

Copy number variants (CNVs) are a major source of genetic variation and can disrupt genes or affect gene dosage. They are known to be causal or underlie predisposition to various diseases. However, the role of CNVs in inherited breast cancer susceptibility has not been thoroughly investigated. To address this, we performed whole-exome sequencing based analysis of rare CNVs in 98 high-risk Northern Finnish breast cancer cases. After filtering, selected candidate alleles were validated and characterized with a combination of orthogonal methods, including PCR-based approaches, optical genome mapping and long-read sequencing. This revealed three recurrent alterations: a 31 kb deletion co-occurring with a retrotransposon insertion (delins) in RAD52, a 13.4 kb deletion in HSD17B14 and a 64 kb partial duplication of RAD51C. Notably, all these genes encode proteins involved in pathways previously identified as essential for breast cancer development. Variants were genotyped in geographically matched cases and controls (altogether 278 hereditary and 1983 unselected breast cancer cases, and 1229 controls). The RAD52 delins and HSD17B14 deletion both showed significant enrichment among cases with indications of hereditary disease susceptibility. RAD52 delins was identified in 7/278 cases (2.5%, P = 0.034, OR = 2.86, 95% CI = 1.10-7.45) and HSD17B14 deletion in 8/278 cases (2.9%, P = 0.014, OR = 3.28, 95% CI = 1.31-8.23), the frequency of both variants in the controls being 11/1229 (0.9%). This suggests a role for RAD52 and HSD17B14 in hereditary breast cancer susceptibility. The RAD51C duplication was very rare, identified only in 2/278 of hereditary cases and 2/1229 controls (P = 0.157, OR = 4.45, 95% CI = 0.62-31.70). The identification of recurrent CNVs in these genes, and especially the relatively high frequency of RAD52 and HSD17B14 alterations in the Finnish population, highlights the importance of studying CNVs alongside single nucleotide variants when searching for genetic factors underlying hereditary disease predisposition.

Applications of advanced technologies for detecting genomic structural variation

Chromosomal structural variation (SV) encompasses a heterogenous class of genetic variants that exerts strong influences on human health and disease. Despite their importance, many structural variants (SVs) have remained poorly characterized at even a basic level, a discrepancy predicated upon the technical limitations of prior genomic assays. However, recent advances in genomic technology can identify and localize SVs accurately, opening new questions regarding SV risk factors and their impacts in humans. Here, we first define and classify human SVs and their generative mechanisms, highlighting characteristics leveraged by various SV assays. We next examine the first-ever gapless assembly of the human genome and the technical process of assembling it, which required third-generation sequencing technologies to resolve structurally complex loci. The new portions of that "telomere-to-telomere" and subsequent pangenome assemblies highlight aspects of SV biology likely to develop in the near-term. We consider the strengths and limitations of the most promising new SV technologies and when they or longstanding approaches are best suited to meeting salient goals in the study of human SV in population-scale genomics research, clinical, and public health contexts. It is a watershed time in our understanding of human SV when new approaches are expected to fundamentally change genomic applications.