Segmental copy number variation shapes tissue transcriptomes

Association analysis of the sorting nexin 29 (SNX29) gene copy number variations with growth traits in Diannan small-ear (DSE) pigs

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

MIR937 amplification potentiates ovarian cancer progression by attenuating FBXO16 inhibition on ULK1-mediated autophagy

High-grade serous ovarian carcinoma (HGSOC) is one of the most lethal gynecological cancer. Genetic studies have revealed gene copy number alterations (CNAs) frequently occurred in HGSOC pathogenesis, however the function and mechanism of CNAs for microRNAs are still not fully understood. Here, we show the dependence on gene copy number amplification of MIR937 that enhances cell autophagy and dictates HGSOC proliferative activity. Data mining of TCGA database revealed MIR937 amplification is correlated with increased MIR937 expression and cell proliferation of HGSOC. Deletion of MIR937 in HGSOC cells led to impaired autophagy and retarded cell proliferation, and the extent for its inhibitory effects scaled with the degree of MIR937 copy loss. Rescue assay confirmed miR-937-5p, a mature product of MIR937, was sufficient to restore its oncogenic function. Mechanistically, MIR937 amplification raised the expression of miR-937-5p, enhanced its binding to 3' UTR of FBXO16 transcript, and thereby restricting FBXO16 degradative effects on ULK1. Our results demonstrate that MIR937 amplification augments cell autophagy and proliferation, and suggest an alternative strategy of MIR937/FBXO16/ULK1 targeting for HGSOC treatment.

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.

Allele-specific transcriptional effects of subclonal copy number alterations enable genotype-phenotype mapping in cancer cells

Subclonal copy number alterations are a prevalent feature in tumors with high chromosomal instability and result in heterogeneous cancer cell populations with distinct phenotypes. However, the extent to which subclonal copy number alterations contribute to clone-specific phenotypes remains poorly understood. We develop TreeAlign, which computationally integrates independently sampled single-cell DNA and RNA sequencing data from the same cell population. TreeAlign accurately encodes dosage effects from subclonal copy number alterations, the impact of allelic imbalance on allele-specific transcription, and obviates the need to define genotypic clones from a phylogeny a priori, leading to highly granular definitions of clones with distinct expression programs. These improvements enable clone-clone gene expression comparisons with higher resolution and identification of expression programs that are genomically independent. Our approach sets the stage for dissecting the relative contribution of fixed genomic alterations and dynamic epigenetic processes on gene expression programs in cancer.

Circulating tumor DNA-based copy-number profiles enable monitoring treatment effects during therapy in high-grade serous carcinoma

Circulating tumor DNA (ctDNA) analysis has emerged as a promising tool for detecting and profiling longitudinal genomics changes in cancer. While copy-number alterations (CNAs) play a major role in cancers, treatment effect monitoring using copy-number profiles has received limited attention as compared to mutations. A major reason for this is the insensitivity of CNA analysis for the real-life tumor-fraction ctDNA samples. We performed copy-number analysis on 152 plasma samples obtained from 29 patients with high-grade serous ovarian cancer (HGSC) using a sequencing panel targeting over 500 genes. Twenty-one patients had temporally matched tissue and plasma sample pairs, which enabled assessing concordance with tissues sequenced with the same panel or whole-genome sequencing and to evaluate sensitivity. Our approach could detect concordant CNA profiles in most plasma samples with as low as 5% tumor content and highly amplified regions in samples with ∼1% of tumor content. Longitudinal profiles showed changes in the CNA profiles in seven out of 11 patients with high tumor-content plasma samples at relapse. These changes included focal acquired or lost copy-numbers, even though most of the genome remained stable. Two patients displayed major copy-number profile changes during therapy. Our analysis revealed ctDNA-detectable subclonal selection resulting from both surgical operations and chemotherapy. Overall, longitudinal ctDNA data showed acquired and diminished CNAs at relapse when compared to pre-treatment samples. These results highlight the importance of genomic profiling during treatment as well as underline the usability of ctDNA.

Cell-free multi-omics analysis reveals potential biomarkers in gastrointestinal cancer patients' blood

During cancer progression, tumorigenic and immune signals are spread through circulating molecules, such as cell-free DNA (cfDNA) and cell-free RNA (cfRNA) in the blood. So far, they have not been comprehensively investigated in gastrointestinal cancers. Here, we profile 4 categories of cell-free omics data from patients with colorectal cancer and patients with stomach adenocarcinoma and then assay 15 types of genomic, epigenomic, and transcriptomic variations. We find that multi-omics data are more appropriate for detection of cancer genes compared with single-omics data. In particular, cfRNAs are more sensitive and informative than cfDNAs in terms of detection rate, enriched functional pathways, etc. Moreover, we identify several peripheral immune signatures that are suppressed in patients with cancer. Specifically, we establish a γδ-T cell score and a cancer-associated-fibroblast (CAF) score, providing insights into clinical statuses like cancer stage and survival. Overall, we reveal a cell-free multi-molecular landscape that is useful for blood monitoring in personalized cancer treatment.

Exploiting allele-specific transcriptional effects of subclonal copy number alterations for genotype-phenotype mapping in cancer cell populations

Somatic copy number alterations drive aberrant gene expression in cancer cells. In tumors with high levels of chromosomal instability, subclonal copy number alterations (CNAs) are a prevalent feature which often result in heterogeneous cancer cell populations with distinct phenotypes1. However, the extent to which subclonal CNAs contribute to clone-specific phenotypes remains poorly understood, in part due to the lack of methods to quantify how CNAs influence gene expression at a subclone level. We developed TreeAlign, which computationally integrates independently sampled single-cell DNA and RNA sequencing data from the same cell population and explicitly models gene dosage effects from subclonal alterations. We show through quantitative benchmarking data and application to human cancer data with single cell DNA and RNA libraries that TreeAlign accurately encodes clone-specific transcriptional effects of subclonal CNAs, the impact of allelic imbalance on allele-specific transcription, and obviates the need to arbitrarily define genotypic clones from a phylogenetic tree a priori. Combined, these advances lead to highly granular definitions of clones with distinct copy-number driven expression programs with increased resolution and accuracy over competing methods. The resulting improvement in assignment of transcriptional phenotypes to genomic clones enables clone-clone gene expression comparisons and explicit inference of genes that are mechanistically altered through CNAs, and identification of expression programs that are genomically independent. Our approach sets the stage for dissecting the relative contribution of fixed genomic alterations and dynamic epigenetic processes on gene expression programs in cancer.

Recent advances and current challenges in population genomics of structural variation in animals and plants

The field of population genomics has seen a surge of studies on genomic structural variation over the past two decades. These studies witnessed that structural variation is taxonomically ubiquitous and represent a dominant form of genetic variation within species. Recent advances in technology, especially the development of long-read sequencing platforms, have enabled the discovery of structural variants (SVs) in previously inaccessible genomic regions which unlocked additional structural variation for population studies and revealed that more SVs contribute to evolution than previously perceived. An increasing number of studies suggest that SVs of all types and sizes may have a large effect on phenotype and consequently major impact on rapid adaptation, population divergence, and speciation. However, the functional effect of the vast majority of SVs is unknown and the field generally lacks evidence on the phenotypic consequences of most SVs that are suggested to have adaptive potential. Non-human genomes are heavily under-represented in population-scale studies of SVs. We argue that more research on other species is needed to objectively estimate the contribution of SVs to evolution. We discuss technical challenges associated with SV detection and outline the most recent advances towards more representative reference genomes, which opens a new era in population-scale studies of structural variation.

Genome-wide association analysis reveals 6 copy number variations associated with the number of cervical vertebrae in Pekin ducks

As a critical developmental stage in vertebrates, the vertebral column formation process is under strict control; however, we observed variations in the number of cervical vertebrae in duck populations in our previous study. Here, we further explored the variations in the number of vertebrae in two duck populations: 421 Pekin duck × mallard F2 ducks and 850 Pekin ducks. Using resequencing data of 125 Pekin ducks with different numbers of cervical vertebrae and 352 Pekin duck × mallard F2 ducks with different numbers of thoracic vertebrae, we detected whole-genome copy number variations (CNVs) and implemented a genome-wide association study (GWAS) to identify the genetic variants related to the traits. The findings verified the existence of variations in the number of cervical vertebrae in duck populations. The number of cervical vertebrae in most ducks was 15, while that in a small number of the ducks was 14 or 16. The number of cervical vertebrae had a positive influence on the neck production, and one cervical vertebra addition could increase 11 g or 2 cm of duck neck. Genome-wide CNV association analysis identified six CNVs associated with the number of cervical vertebrae, and the associated CNV regions covered 15 genes which included WNT10A and WNT6. These findings improve our understanding of the variations in the number of vertebrae in ducks and lay a foundation for future duck breeding.

Copy number variant-based genome wide association study reveals immune-related genes associated with parasite resistance in a heritage sheep breed from the United States

Florida Native is a heritage sheep breed in the United States and expresses superior ability to regulate gastrointestinal nematodes. The objective of the present study was to investigate the importance of copy number variants (CNVs) on resistance to natural Haemonchus contortus infections. A total of 300 Florida Native sheep were evaluated. Phenotypic records included fecal egg count (FEC, eggs/gram), FAMACHA© score, percentage cell volume (PCV, %), body condition score (BCS) and average daily gain (ADG, kg). Sheep were genotyped using the GGP Ovine 50K single nucleotide polymorphism (SNP) chip. Log ratios from 45.2 k SNP markers spanning the entire genome were utilized for CNV detection. After quality control, 261 animals with CNVs and phenotypic records were used for the association testing. Association tests were carried out using correlation-trend test and principal component analysis correction to identify CNVs associated with FEC, FAMACHA©, PCV, BCS and ADG. Significant CNVs were detected when their adjusted p-value was <.05 after FDR correction. A total of 8124 CNVs were identified, which gave 246 non-overlapping CNVs. Fourteen CNVs were significantly associated with FEC and PCV. CNVs associated with FEC overlapped 14 Quantitative Trait Locus previously associated with H. contortus resistance. Our study demonstrated for the first time that CNVs could be potentially involved with parasite resistance in Florida Native sheep. Immune-related genes such as CCL1, CCL2, CCL8, CCL11, NOS2, TNF, CSF3 and STAT3 genes could play an important role for controlling H. contortus resistance. These genes could be potentially utilized as candidate markers for selection of parasite resistance in this breed.

Clonal somatic copy number altered driver events inform drug sensitivity in high-grade serous ovarian cancer

Chromosomal instability is a major challenge to patient stratification and targeted drug development for high-grade serous ovarian carcinoma (HGSOC). Here we show that somatic copy number alterations (SCNAs) in frequently amplified HGSOC cancer genes significantly correlate with gene expression and methylation status. We identify five prevalent clonal driver SCNAs (chromosomal amplifications encompassing MYC, PIK3CA, CCNE1, KRAS and TERT) from multi-regional HGSOC data and reason that their strong selection should prioritise them as key biomarkers for targeted therapies. We use primary HGSOC spheroid models to test interactions between in vitro targeted therapy and SCNAs. MYC chromosomal copy number is associated with in-vitro and clinical response to paclitaxel and in-vitro response to mTORC1/2 inhibition. Activation of the mTOR survival pathway in the context of MYC-amplified HGSOC is statistically associated with increased prevalence of SCNAs in genes from the PI3K pathway. Co-occurrence of amplifications in MYC and genes from the PI3K pathway is independently observed in squamous lung cancer and triple negative breast cancer. In this work, we show that identifying co-occurrence of clonal driver SCNA genes could be used to tailor therapeutics for precision medicine.

Robustness of the Autophagy Pathway to Somatic Copy Number Losses

Autophagy allows cells to temporarily tolerate energy stress by replenishing critical metabolites through self-digestion, thereby attenuating the cytotoxic effects of anticancer drugs that target tumor metabolism. Autophagy defects could therefore mark a metabolically vulnerable cancer state and open a therapeutic window. While mutations of autophagy genes (ATGs) are notably rare in cancer, haploinsufficiency network analyses across many cancers have shown that the autophagy pathway is frequently hit by somatic copy number losses of ATGs such as MAP1LC3B/ATG8F (LC3), BECN1/ATG6 (Beclin-1), and ATG10. Here, we used CRISPR/Cas9 technology to delete increasing numbers of copies of one or more of these ATGs in non-small cell lung cancer cells and examined the effects on sensitivity to compounds targeting aerobic glycolysis, a hallmark of cancer metabolism. Whereas the complete knockout of one ATG blocked autophagy and led to profound metabolic vulnerability, this was not the case for combinations of different nonhomozygous deletions. In cancer patients, the effect of ATG copy number loss was blunted at the protein level and did not lead to the accumulation of p62 as a sign of reduced autophagic flux. Thus, the autophagy pathway is shown to be markedly robust and resilient, even with the concomitant copy number loss of key autophagy genes.

Differential Regulation of NK Cell Receptors in Acute Lymphoblastic Leukemia

Cancer immunotherapies are preferred over conventional treatments which are highly cytotoxic to normal cells. Focus has been on T cells but natural killer (NK) cells have equal potential. Concepts in cancer control and influence of sex require further investigation to improve successful mobilization of immune cells in cancer patients. Acute lymphoblastic leukemia (ALL) is a hematological malignancy mainly of B cell (B-ALL) and T cell (T-ALL) subtypes. Influence of ALL on NK cell is still unclear. Targeted next-generation sequencing was conducted on 62 activating/inhibitory receptors, ligands, effector, and exhaustion molecules on T-ALL (6 males) and normal controls (NC) (4 males and 4 females). Quantitative PCR (q-PCR) further investigated copy number variation (CNV), methylation index (MI), and mRNA expression of significant genes in T-ALL (14 males), NC (12 males and 12 females), and B-ALL samples (N = 12 males and 12 females). Bioinformatics revealed unique variants particularly rs2253849 (T>C) in KLRC1 and rs1141715 (A>G) in KLRC2 only among T-ALL (allele frequency 0.8-1.0). Gene amplification was highest in female B-ALL compared to male B-ALL (KLRC2, KLRC4, and NCR3, p < 0.05) and lowest in male T-ALL cumulating in deletion of KLRD1 and CD69. MI was higher in male ALL of both subtypes compared to normal (KIR2DL1-2 and 4 and KIR2DS2 and 4, p < 0.05) as well as to female B-ALL (KIR3DL2 and KIR2DS2, p < 0.05). mRNA expressions were low. Thus, ALL subtypes potentially regulated NK cell suppression by different mechanisms which should be considered in future immunotherapies for ALL.

Copy-number intratumor heterogeneity increases the risk of relapse in chemotherapy-naive stage II colon cancer

Optimal selection of high-risk patients with stage II colon cancer is crucial to ensure clinical benefit of adjuvant chemotherapy. Here, we investigated the prognostic value of genomic intratumor heterogeneity and aneuploidy for disease recurrence. We combined targeted sequencing, SNP arrays, fluorescence in situ hybridization, and immunohistochemistry on a retrospective cohort of 84 untreated stage II colon cancer patients. We assessed the clonality of copy-number alterations (CNAs) and mutations, CD8⁺ lymphocyte infiltration, and their association with time to recurrence. Prognostic factors were included in machine learning analysis to evaluate their ability to predict individual relapse risk. Tumors from recurrent patients displayed a greater proportion of CNAs compared with non-recurrent (mean 31.3% versus 23%, respectively; p = 0.014). Furthermore, patients with elevated tumor CNA load exhibited a higher risk of recurrence compared with those with low levels [p = 0.038; hazard ratio (HR) 2.46], which was confirmed in an independent cohort (p = 0.004; HR 3.82). Candidate chromosome-specific aberrations frequently observed in recurrent cases included gain of the chromosome arm 13q (p = 0.02; HR 2.67) and loss of heterozygosity at 17q22-q24.3 (p = 0.05; HR 2.69). CNA load positively correlated with intratumor heterogeneity (R = 0.52; p < 0.0001). Consistently, incremental subclonal CNAs were associated with an elevated risk of relapse (p = 0.028; HR 2.20), which we did not observe for subclonal single-nucleotide variants and small insertions and deletions. The clinico-genomic model rated an area under the curve of 0.83, achieving a 10% incremental gain compared with clinicopathological markers (p = 0.047). In conclusion, tumor aneuploidy and copy-number intratumor heterogeneity were predictive of poor outcome and improved discriminative performance in early-stage colon cancer. © 2022 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

G-Quadruplex Matters in Tissue-Specific Tumorigenesis by BRCA1 Deficiency

How and why distinct genetic alterations, such as BRCA1 mutation, promote tumorigenesis in certain tissues, but not others, remain an important issue in cancer research. The underlying mechanisms may reveal tissue-specific therapeutic vulnerabilities. Although the roles of BRCA1, such as DNA damage repair and stalled fork stabilization, obviously contribute to tumor suppression, these ubiquitously important functions cannot explain tissue-specific tumorigenesis by BRCA1 mutations. Recent advances in our understanding of the cancer genome and fundamental cellular processes on DNA, such as transcription and DNA replication, have provided new insights regarding BRCA1-associated tumorigenesis, suggesting that G-quadruplex (G4) plays a critical role. In this review, we summarize the importance of G4 structures in mutagenesis of the cancer genome and cell type-specific gene regulation, and discuss a recently revealed molecular mechanism of G4/base excision repair (BER)-mediated transcriptional activation. The latter adequately explains the correlation between the accumulation of unresolved transcriptional regulatory G4s and multi-level genomic alterations observed in BRCA1-associated tumors. In summary, tissue-specific tumorigenesis by BRCA1 deficiency can be explained by cell type-specific levels of transcriptional regulatory G4s and the role of BRCA1 in resolving it. This mechanism would provide an integrated understanding of the initiation and development of BRCA1-associated tumors.

Pan-Cancer Study on Protein Kinase C Family as a Potential Biomarker for the Tumors Immune Landscape and the Response to Immunotherapy

The protein kinase C (PKC) family has been described with its role in some cancers, either as a promoter or suppressor. PKC signaling also regulates a molecular switch between transactivation and transrepression activity of the peroxisome proliferator-activated receptor alpha (PPARalpha). However, the role of different PKC enzymes in tumor immunity remains poorly defined. This study aims to investigate the correlation between PKC genes and tumor immunity, in addition to studying the probability of their use as predictive biomarkers for tumor immunity and immunotherapeutic response. The ssGSEA and the ESTIMATE methods were used to assess 28 tumor-infiltrating lymphocytes (TILs) and the immune component of each cancer, then correlated with PKC levels. Prediction of PKC levels-dependent immunotherapeutic response was based on human leukocytic antigen (HLA) gene enrichment scores and programmed cell death 1 ligand (PD-L1) expression. Univariate and multivariate Cox analysis was performed to evaluate the prognostic role of PKC genes in cancers. Methylation level and CNAs could drive the expression levels of some PKC members, especially PRKCI, whose CNGs are predicted to elevate their level in many cancer types. The most crucial finding in this study was that PKC isoenzymes are robust biomarkers for the tumor immune status, PRKCB, PRKCH, and PRKCQ as stimulators, while PRKCI and PRKCZ as inhibitors in most cancers. Also, PKC family gene levels can be used as predictors for the response to immunotherapies, especially HLAs dependent and PD-L1 blockade-dependent ones. In addition to its prognostic function, all PKC family enzymes are promising tumor immunity biomarkers and can help select suitable immune therapy in different cancers.

Copy number variations of the KAT6A gene are associated with body measurements of Chinese sheep breeds

Copy number variation (CNV) is one kind of genomic structure variations and presents as gains and losses of genomic fragments. More recently, we have made an atlas of CNV maps for livestock. In the future, it is a primary focus to determine the phenotypic effects of candidate CNVs. Lysine Acetyltransferase 6 A (KAT6A) is a protein coding gene and plays a critical role in many cellular processes. However, the effects of KAT6A CNVs on sheep body measurements remains unknown. In this study, we performed quantitative polymerase chain reaction (qPCR) to detect the presences and distributions of three CNV regions within KAT6A gene in 672 sheep from four Chinese breeds. Association analysis indicated that the three CNVs of KAT6A gene were significantly associated with body measurement(s) in Small-tailed Han sheep (STH) and Hu sheep (HU) (p < 0.05), while no effects on Large-tailed Han sheep (LTH) were observed (p > 0.05) were observed. Additionally, only one CNV was significantly associated with body measurement (body length) in Chaka sheep (CK) (p < 0.05). Our study provided evidence that the CNV(s) of KAT6A gene could be used as candidate marker(s) for molecular breedings of STH, HU, and CK breeds.

Gv1, a Zinc Finger Gene Controlling Endogenous MLV Expression

The genomes of inbred mice harbor around 50 endogenous murine leukemia virus (MLV) loci, although the specific complement varies greatly between strains. The Gv1 locus is known to control the transcription of endogenous MLVs and to be the dominant determinant of cell-surface presentation of MLV envelope, the G_IX antigen. Here, we identify a single Krüppel-associated box zinc finger protein (ZFP) gene, Zfp998, as Gv1 and show it to be necessary and sufficient to determine the GIX+ phenotype. By long-read sequencing of bacterial artificial chromosome clones from 129 mice, the prototypic GIX+ strain, we reveal the source of sufficiency and deficiency as splice-acceptor variations and highlight the varying origins of the chromosomal region encompassing Gv1. Zfp998 becomes the second identified ZFP gene responsible for epigenetic suppression of endogenous MLVs in mice and further highlights the prominent role of this gene family in control of endogenous retroviruses.

The Breeding of Winter-Hardy Malting Barley

In breeding winter malting barley, one recurring strategy is to cross a current preferred spring malting barley to a winter barley. This is because spring malting barleys have the greatest amalgamation of trait qualities desirable for malting and brewing. Spring barley breeding programs can also cycle their material through numerous generations each year-some managing even six-which greatly accelerates combining desirable alleles to generate new lines. In a winter barley breeding program, a single generation per year is the limit when the field environment is used and about two generations per year if vernalization and greenhouse facilities are used. However, crossing the current favored spring malting barley to a winter barley may have its downsides, as winter-hardiness too may be an amalgamation of desirable alleles assembled together that confers the capacity for prolonged cold temperature conditions. In this review I touch on some general criteria that give a variety the distinction of being a malting barley and some of the general trends made in the breeding of spring malting barleys. But the main objective of this review is to pull together different aspects of what we know about winter-hardiness from the seemingly most essential aspect, which is survival in the field, to molecular genetics and gene regulation, and then finish with ideas that might help further our insight for predictability purposes.

Importance of determining variations in the number of copies in newborns with autosomal aneuploidies

Introducción.

Las aneuploidías son trastornos genéticos frecuentes en la práctica clínica; sin embargo, se conoce poco sobre las otras variantes genéticas que modifican el fenotipo final.

Objetivo.

Determinar las variantes en el número de copias y las regiones con pérdida de heterocigosidad autosómica mayor de 0,5 % o de regiones mayores de 10 Mb en neonatos con aneuploidías autosómicas.

Materiales y métodos.

Se hizo el análisis cromosómico por micromatrices a los neonatos con aneuploidías autosómicas (n=7), trisomía 21 (n=5) y trisomía 18 (n=2) evaluados en los hospitales Antonio Lorena y Regional de Cusco, Perú, en el 2018.

Resultados.

En dos neonatos se encontraron variantes en el número de copias, patogénicas o probablemente patogénicas, en regiones diferentes al cromosoma 21 o al 18. Además, se observaron dos variantes del número de copias con más de 500 kpb de patogenia desconocida.

Conclusiones.

Si bien el número de pacientes era muy reducido, es importante resaltar que se encontraron otras variantes en el número de copias que se han descrito asociadas con trastornos del neurodesarrollo, varias anomalías congénitas, hipoacusia y talla baja o alta, entre otras, lo que probablemente influye negativamente en el fenotipo de este grupo de pacientes.

Copy number variation: Characteristics, evolutionary and pathological aspects

Copy number variants (CNVs) were the subject of extensive research in the past years. They are common features of the human genome that play an important role in evolution, contribute to population diversity, development of certain diseases, and influence host–microbiome interactions. CNVs have found application in the molecular diagnosis of many diseases and in non-invasive prenatal care, but their full potential is only emerging. CNVs are expected to have a tremendous impact on screening, diagnosis, prognosis, and monitoring of several disorders, including cancer and cardiovascular disease. Here, we comprehensively review basic definitions of the term CNV, outline mechanisms and factors involved in CNV formation, and discuss their evolutionary and pathological aspects. We suggest a need for better defined distinguishing criteria and boundaries between known types of CNVs.

A Mechanism Leading to Changes in Copy Number Variations Affected by Transcriptional Level Might Be Involved in Evolution, Embryonic Development, Senescence, and Oncogenesis Mediated by Retrotransposons

In recent years, more and more evidence has emerged showing that changes in copy number variations (CNVs) correlated with the transcriptional level can be found during evolution, embryonic development, and oncogenesis. However, the underlying mechanisms remain largely unknown. The success of the induced pluripotent stem cell suggests that genome changes could bring about transformations in protein expression and cell status; conversely, genome alterations generated during embryonic development and senescence might also be the result of genome changes. With rapid developments in science and technology, evidence of changes in the genome affected by transcriptional level has gradually been revealed, and a rational and concrete explanation is needed. Given the preference of the HIV-1 genome to insert into transposons of genes with high transcriptional levels, we propose a mechanism based on retrotransposons facilitated by specific pre-mRNA splicing style and homologous recombination (HR) to explain changes in CNVs in the genome. This mechanism is similar to that of the group II intron that originated much earlier. Under this proposed mechanism, CNVs on genome are dynamically and spontaneously extended in a manner that is positively correlated with transcriptional level or contract as the cell divides during evolution, embryonic development, senescence, and oncogenesis, propelling alterations in them. Besides, this mechanism explains several critical puzzles in these processes. From evidence collected to date, it can be deduced that the message contained in genome is not just three-dimensional but will become four-dimensional, carrying more genetic information.

Core transcriptional networks in Williams syndrome: IGF1-PI3K-AKT-mTOR, MAPK and actin signaling at the synapse echo autism

Gene networks for disorders of social behavior provide the mechanisms critical for identifying therapeutic targets and biomarkers. Large behavioral phenotypic effects of small human deletions make the positive sociality of Williams syndrome (WS) ideal for determining transcriptional networks for social dysfunction currently based on DNA variations for disorders such as autistic spectrum disorder (ASD) and schizophrenia (SCHZ). Consensus on WS networks has been elusive due to the need for larger cohort size, sensitive genome-wide detection and analytic tools. We report a core set of WS network perturbations in a cohort of 58 individuals (34 with typical, 6 atypical deletions and 18 controls). Genome-wide exon-level expression arrays robustly detected changes in differentially expressed gene (DEG) transcripts from WS deleted genes that ranked in the top 11 of 12 122 transcripts, validated by quantitative reverse transcription PCR, RNASeq and western blots. WS DEG's were strictly dosed in the full but not the atypical deletions that revealed a breakpoint position effect on non-deleted CLIP2, a caveat for current phenotypic mapping based on copy number variants. Network analyses tested the top WS DEG's role in the dendritic spine, employing GeneMANIA to harmonize WS DEGs with comparable query gene-sets. The results indicate perturbed actin cytoskeletal signaling analogous to the excitatory dendritic spines. Independent protein-protein interaction analyses of top WS DEGs generated a 100-node graph annotated topologically revealing three interacting pathways, MAPK, IGF1-PI3K-AKT-mTOR/insulin and actin signaling at the synapse. The results indicate striking similarity of WS transcriptional networks to genome-wide association study-based ASD and SCHZ risk suggesting common network dysfunction for these disorders of divergent sociality.

Characterization of Copy-Number Variations and Possible Candidate Genes in Recurrent Pregnancy Losses

It is well established that embryonic chromosomal abnormalities (both in the number of chromosomes and the structure) account for 50% of early pregnancy losses. However, little is known regarding the potential differences in the incidence and distribution of chromosomal abnormalities between patients with sporadic abortion (SA) and recurrent pregnancy loss (RPL), let alone the role of submicroscopic copy-number variations (CNVs) in these cases. The aim of the present study was to systematically evaluate the role of embryonic chromosomal abnormalities and CNVs in the etiology of RPL compared with SA. Over a 3-year period, 1556 fresh products of conception (POCs) from miscarriage specimens were investigated using single nucleotide polymorphism array (SNP-array) and CNV sequencing (CNV-seq) in this study, along with further functional enrichment analysis. Chromosomal abnormalities were identified in 57.52% (895/1556) of all cases. Comparisons of the incidence and distributions of chromosomal abnormalities within the SA group and RPL group and within the different age groups were performed. Moreover, 346 CNVs in 173 cases were identified, including 272 duplications, 2 deletions and 72 duplications along with deletions. Duplications in 16q24.3 and 16p13.3 were significantly more frequent in RPL cases, and thereby considered to be associated with RPL. There were 213 genes and 131 signaling pathways identified as potential RPL candidate genes and signaling pathways, respectively, which were centered primarily on six functional categories. The results of the present study may improve our understanding of the etiologies of RPL and assist in the establishment of a population-based diagnostic panel of genetic markers for screening RPL amongst Chinese women.

Metabarcoding From Microbes to Mammals: Comprehensive Bioassessment on a Global Scale

Global biodiversity loss is unprecedented, and threats to existing biodiversity are growing. Given pervasive global change, a major challenge facing resource managers is a lack of scalable tools to rapidly and consistently measure Earth's biodiversity. Environmental genomic tools provide some hope in the face of this crisis, and DNA metabarcoding, in particular, is a powerful approach for biodiversity assessment at large spatial scales. However, metabarcoding studies are variable in their taxonomic, temporal, or spatial scope, investigating individual species, specific taxonomic groups, or targeted communities at local or regional scales. With the advent of modern, ultra-high throughput sequencing platforms, conducting deep sequencing metabarcoding surveys with multiple DNA markers will enhance the breadth of biodiversity coverage, enabling comprehensive, rapid bioassessment of all the organisms in a sample. Here, we report on a systematic literature review of 1,563 articles published about DNA metabarcoding and summarize how this approach is rapidly revolutionizing global bioassessment efforts. Specifically, we quantify the stakeholders using DNA metabarcoding, the dominant applications of this technology, and the taxonomic groups assessed in these studies. We show that while DNA metabarcoding has reached global coverage, few studies deliver on its promise of near-comprehensive biodiversity assessment. We then outline how DNA metabarcoding can help us move toward real-time, global bioassessment, illustrating how different stakeholders could benefit from DNA metabarcoding. Next, we address barriers to widespread adoption of DNA metabarcoding, highlighting the need for standardized sampling protocols, experts and computational resources to handle the deluge of genomic data, and standardized, open-source bioinformatic pipelines. Finally, we explore how technological and scientific advances will realize the promise of total biodiversity assessment in a sample—from microbes to mammals—and unlock the rich information genomics exposes, opening new possibilities for merging whole-system DNA metabarcoding with (1) abundance and biomass quantification, (2) advanced modeling, such as species occupancy models, to improve species detection, (3) population genetics, (4) phylogenetics, and (5) food web and functional gene analysis. While many challenges need to be addressed to facilitate widespread adoption of environmental genomic approaches, concurrent scientific and technological advances will usher in methods to supplement existing bioassessment tools reliant on morphological and abiotic data. This expanded toolbox will help ensure that the best tool is used for the job and enable exciting integrative techniques that capitalize on multiple tools. Collectively, these new approaches will aid in addressing the global biodiversity crisis we now face.

Neurodevelopmental trajectory and modifiers of 16p11.2 microdeletion: A follow-up study of four Chinese children carriers

Background

Neurodevelopmental disorders (NDDs) are a group of disorders with high genetic and phenotypic heterogeneities. The 16p11.2 microdeletion has been implicated as an important genetic risk factor for NDDs.

Methods

Multiple genetic tests were used to detect the 16p11.2 microdeletion from 918 Chinese children with NDDs. Targeted sequencing of genes in the 16p11.2 interval was performed in all carriers of the 16p11.2 microdeletion, and whole‐genome expression profiling analysis was performed for the patient carriers and normal carriers in their intra‐family.

Results

Three patients carrying the 16p11.2 microdeletion were screened out, indicating a frequency of 0.33% for the 16p11.2 microdeletion in this cohort. We reviewed the neurodevelopmental trajectories of the 16p11.2 microdeletion carriers from childhood to puberty and confirmed that this microdeletion was associated with abnormal neurodevelopment, with varied neurodevelopmental phenotypes. A differential PRRT2 genotype (rs10204, T>C) was identified between patients and normal carriers of the 16p11.2 microdeletion. Moreover, the determination of differential whole‐genome expression profiling demonstrated the destruction of the top‐ranked network in neurogenesis and accounted for observation of abnormal neurodevelopmental phenotypes in the 16p11.2 microdeletion carriers.

Conclusions

We have provided the frequency of the 16p11.2 microdeletion in a Chinese pediatric NDD cohort with a variable NDD phenotype from childhood to puberty, which is useful for Chinese geneticists/pediatricians to conduct the 16p11.2 microdeletion testing in children with NDDs.

A de novo 13q31.3 microduplication encompassing the miR-17 ~ 92 cluster results in features mirroring those associated with Feingold syndrome 2

Hemizygosity of the MIR17HG gene encoding the miR-17 ~ 92 cluster is associated with Feingold syndrome 2 characterized by intellectual disability, skeletal abnormalities, short stature, and microcephaly. Here, we report on a female with a de novo 13q31.3 microduplication encompassing MIR17HG but excluding GPC5. She presented developmental delay, skeletal and digital abnormalities, and features such as tall stature and macrocephaly mirroring those of Feingold syndrome 2 patients. The limited extent of the proband's rearrangement to the miR cluster and the corresponding normal expression level of the neighboring GPC5 in her cells, together with previously described data on affected individuals of two families carrying overlapping duplications of the miR-17 ~ 92 cluster that comprise part of GPC5, who likewise presented macrocephaly, developmental delay, as well as skeletal, digital and stature abnormalities, allow to define a new syndrome due to independent microduplication of the miR-17 ~ 92 cluster.

A genome-wide survey of copy number variations reveals an asymmetric evolution of duplicated genes in rice

Background

Copy number variations (CNVs) are an important type of structural variations in the genome that usually affect gene expression levels by gene dosage effect. Understanding CNVs as part of genome evolution may provide insights into the genetic basis of important agricultural traits and contribute to the crop breeding in the future. While available methods to detect CNVs utilizing next-generation sequencing technology have helped shed light on prevalence and effects of CNVs, the complexity of crop genomes poses a major challenge and requires development of additional tools.

Results

Here, we generated genomic and transcriptomic data of 93 rice (Oryza sativa L.) accessions and developed a comprehensive pipeline to call CNVs in this large-scale dataset. We analyzed the correlation between CNVs and gene expression levels and found that approximately 13% of the identified genes showed a significant correlation between their expression levels and copy numbers. Further analysis showed that about 36% of duplicate pairs were involved in pseudogenetic events while only 5% of them showed functional differentiation. Moreover, the offspring copy mainly contributed to the expression levels and seemed more likely to become a pseudogene, whereas the parent copy tended to maintain the function of ancestral gene.

Conclusion

We provide a high-accuracy CNV dataset that will contribute to functional genomics studies and molecular breeding in rice. We also showed that gene dosage effect of CNVs in rice is not exponential or linear. Our work demonstrates that the evolution of duplicated genes is asymmetric in both expression levels and gene fates, shedding a new insight into the evolution of duplicated genes.

Sensitivity to gene dosage and gene expression affects genes with copy number variants observed among neuropsychiatric diseases

BACKGROUND

Copy number variants (CNVs) have been reported to be associated with diseases, traits, and evolution. However, it is hard to determine which gene should have priority as a target for further functional experiments if a CNV is rare or a singleton. In this study, we attempted to overcome this issue by using two approaches: by assessing the influences of gene dosage sensitivity and gene expression sensitivity. Dosage sensitive genes derived from two-round whole-genome duplication in previous studies. In addition, we proposed a cross-sectional omics approach that utilizes open data from GTEx to assess the effect of whole-genome CNVs on gene expression.

METHODS

Affymetrix Genome-Wide SNP Array 6.0 was used to detect CNVs by PennCNV and CNV Workshop. After quality controls for population stratification, family relationship and CNV detection, 287 patients with narcolepsy, 133 patients with essential hypersomnia, 380 patients with panic disorders, 164 patients with autism, 784 patients with Alzheimer disease and 1280 healthy individuals remained for the enrichment analysis.

RESULTS

Overall, significant enrichment of dosage sensitive genes was found across patients with narcolepsy, panic disorders and autism. Particularly, significant enrichment of dosage-sensitive genes in duplications was observed across all diseases except for Alzheimer disease. For deletions, less or no enrichment of dosage-sensitive genes with deletions was seen in the patients when compared to the healthy individuals. Interestingly, significant enrichments of genes with expression sensitivity in brain were observed in patients with panic disorder and autism. While duplications presented a higher burden, deletions did not cause significant differences when compared to the healthy individuals. When we assess the effect of sensitivity to genome dosage and gene expression at the same time, the highest ratio of enrichment was observed in the group including dosage-sensitive genes and genes with expression sensitivity only in brain. In addition, shared CNV regions among the five neuropsychiatric diseases were also investigated.

CONCLUSIONS

This study contributed the evidence that dosage-sensitive genes are associated with CNVs among neuropsychiatric diseases. In addition, we utilized open data from GTEx to assess the effect of whole-genome CNVs on gene expression. We also investigated shared CNV region among neuropsychiatric diseases.

Transcriptional effects of copy number alterations in a large set of human cancers

Copy number alterations (CNAs) can promote tumor progression by altering gene expression levels. Due to transcriptional adaptive mechanisms, however, CNAs do not always translate proportionally into altered expression levels. By reanalyzing >34,000 gene expression profiles, we reveal the degree of transcriptional adaptation to CNAs in a genome-wide fashion, which strongly associate with distinct biological processes. We then develop a platform-independent method-transcriptional adaptation to CNA profiling (TACNA profiling)-that extracts the transcriptional effects of CNAs from gene expression profiles without requiring paired CNA profiles. By applying TACNA profiling to >28,000 patient-derived tumor samples we define the landscape of transcriptional effects of CNAs. The utility of this landscape is demonstrated by the identification of four genes that are predicted to be involved in tumor immune evasion when transcriptionally affected by CNAs. In conclusion, we provide a novel tool to gain insight into how CNAs drive tumor behavior via altered expression levels.

Copy number variation profiling in pharmacogenes using panel-based exome resequencing and correlation to human liver expression

Structural variants including copy number variations (CNV) have gained widespread attention, especially in pharmacogenomics but for several genes functional relevance and clinical evidence are still lacking. Detection of CNVs in next-generation sequencing data is challenging but offers widespread applications. We developed a cohort-based CNV detection workflow to extract CNVs from read counts of targeted NGS of 340 genes involved in absorption, distribution, metabolism and excretion (ADME) of drugs. We applied our method to 150 human liver tissue samples and correlated identified CNVs to mRNA expression levels. In total, we identified 445 deletions (73%) and 167 duplications (27%) in 36 pharmacogenes including all well-known CNVs of CYPs, GSTs, SULTs, UGTs, numerous described rare CNVs of CYP2E1, SLC16A3 or UGT2B15 as well as novel observations, e.g., for SLC22A12, SLC22A17 and GPS2 (G Protein Pathway Suppressor 2). We were able to fine-map complex CNVs of CYP2A6 and CYP2D6 with exon resolution. Correlation analysis confirmed known expression patterns for common CNVs and suggested an influence on expression variability for some rare CNVs. Our straightforward CNV detection workflow can be easily applied to any NGS coverage data and helped to analyze CNVs in an ADME-NGS panel of 340 pharmacogenes to improve genotype-phenotype correlations.

Copy Number Variation of the CADM2 Gene and Its Association with Growth Traits in Yak

Simple Summary

Cell adhesion molecule 2 (CADM2), also known as synaptic cell adhesion molecule 2 (SYNCAM2), is the mediator of synaptic signals enriched in the brain. Overlaps between copy number variation (CNV) regions in CADM2 and quantitative trait loci (QTL) related to body weight have been clarified in a previous study. In this study, two loci were amplified in the CADM2 gene (CNV1: 235,915 bp, exon 1 and partial intron 1; CNV2: 60,430 bp, intron 9) to explore the relationship between CNV types in the CADM2 gene and growth traits in 350 Ashidan yaks. Association analysis illustrated that no significant effect was found on growth traits in CNV1. However, the CNV2 mutation had a significant effect on body weight at the sixth month (p < 0.05). Individuals with the gain-type copy number variation CNV2 were significantly superior to those with loss- or normal-type in terms of body weight (p < 0.05). In summary, this study confirmed that CADM2-CNVs affect growth traits in yaks, and may be candidate genes for successful yak breeding and genetics projects.

Abstract

Copy number variation (CNV) is currently accepted as a common source of genetic variation. It is reported that CNVs may influence the resistance to disease and complex economic traits, such as residual feed intake, muscle formation, and fat deposition in livestock. Cell adhesion molecule 2 (CADM2) is expressed widely in the brain and adipose tissue and can regulate body weight through the central nervous system. Growth traits are important economic traits for animal selection. In this study, we aimed to explore the effect of CADM2 gene copy number variants on yak growth traits. Here, two CNVs in the CADM2 gene were investigated using the quantitative polymerase chain reaction (qPCR), and the association of the CNVs with growth traits in yak was analyzed using statistical methods by SPSS software. Differences were considered significant if the p value was < 0.05. Statistical analysis indicated significant association of CADM2-CNV2 with the body weight of the Chinese Ashidan yak. A significant effect of CNV2 (p < 0.05) was found on body weight at 6 months. In CNV2, the gain-type copy number variation exhibited greater performance than the other variants, with greater body weight observed at 6 months (p < 0.05). To the best of our knowledge, this is the first attempt to investigate the function of CADM2-CNVs and their association with growth traits in animals. This may be a useful candidate marker in marker-assisted selection of yaks.

Bovine milk transcriptome analysis reveals microRNAs and RNU2 involved in mastitis

Mastitis is a common inflammatory infectious disease in dairy cows. To understand the microRNA (miRNA) expression profile changes during bovine mastitis, we undertook a genome-wide miRNA study of normal milk and milk that tested positive on the California mastitis test for bovine mastitis (CMT+). Twenty-five miRNAs were differentially expressed (23 miRNAs upregulated and two downregulated) during bovine mastitis relative to their expression in normal milk. Upregulated mature miR-1246 probably derived from a U2 small nuclear RNA rather than an miR-1246 precursor. The significantly upregulated miRNA precursors and RNU2 were significantly enriched on bovine chromosome 19, which is homologous to human chromosome 17. A gene ontology analysis of the putative mRNA targets of the significantly upregulated miRNAs showed that these miRNAs were involved in binding target mRNA transcripts and regulating target gene expression, and a Kyoto Encyclopedia of Genes and Genomes pathway analysis showed that the upregulated miRNAs were predominantly related to cancer and immune system pathways. Three novel miRNAs were associated with bovine mastitis and were relatively highly expressed in milk. We confirmed that one of the novel mastitis-related miRNAs was significantly upregulated using a digital PCR system. The differentially expressed miRNAs were involved in human cancers, infections, and immune-related diseases. The genome-wide analysis of miRNA profiles in this study provides insight into bovine mastitis and inflammatory diseases. DATABASES: The miRNAseq generated for this study can be found in the Sequence Read Archive (SRA) under BioProject Number PRJNA421075 and SRA Study Number SRP126134 (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA421075).

Greater genetic and regulatory plasticity of retained duplicates in Epichloë endophytic fungi

Gene duplicates can act as a source of genetic material from which new functions arise. Most duplicated genes revert to single copy genes and only a small proportion are retained. However, it remains unclear why some duplicate genes persist in the genome for an extended time. We investigate this question by analysing retained gene duplicates in the fungal genus Epichloë, ascomycete fungi that form close endophytic symbioses with their host grasses. Retained duplicates within this genus have two independent origins, but both long pre-date the origin and diversification of the genus Epichloë. We find that loss of retained duplicates within the genus is frequent and often associated with speciation. Retained duplicates have faster evolutionary rates (Ka) and show relaxed selection (Ka/Ks) compared to single copy genes. Both features are time-dependent. Through comparison of conspecific strains, we find greater evolutionary rates in coding regions and sequence divergence in regulatory regions of retained duplicates than single copy genes, with this pattern more pronounced for strains adapted to different grass host species. Consistent with this sequence divergence in regulatory regions, transcriptome analyses show greater expression variation of retained duplicates than single copy genes. This suggest that cis-regulatory changes make important contributions to the expression patterns of retained duplicates. Coupled with supporting observations from the model yeast Saccharomyces cerevisiae, these data suggest that genetic robustness and regulatory plasticity are common drivers behind the retention of duplicated genes in fungi.

C57BL/6 and 129 inbred mouse strains differ in Gbp2 and Gbp2b expression in response to inflammatory stimuli in vivo

Background: Infections cause the production of inflammatory cytokines such as Interferon gamma (IFNγ). IFNγ in turn prompts the upregulation of a range of host defence proteins including members of the family of guanylate binding proteins (Gbps). In humans and mice alike, GBPs restrict the intracellular replication of invasive microbes and promote inflammation. To study the physiological functions of Gbp family members, the most commonly chosen in vivo models are mice harbouring loss-of-function mutations in either individual Gbp genes or the entire Gbp gene cluster on mouse chromosome 3. Individual Gbp deletion strains differ in their design, as some strains exist on a pure C57BL/6 genetic background, while other strains contain a 129-derived genetic interval encompassing the Gbp gene cluster on an otherwise C57BL/6 genetic background. Methods: To determine whether the presence of 129 alleles of paralogous Gbps could influence the phenotypes of 129-congenic Gbp-deficient strains, we studied the expression of Gbps in both C57BL/6J and 129/Sv mice following in vivo stimulation with adjuvants and after infection with either Toxoplasma gondii or Shigella flexneri. Results: We show that C57BL/6J relative to 129/Sv mice display moderately elevated expression of Gbp2, but more prominently, are also defective for Gbp2b (formerly Gbp1) mRNA induction upon immune priming. Notably, Toxoplasma infections induce robust Gbp2b protein expression in both strains of mice, suggestive of a Toxoplasma-activated mechanism driving Gbp2b protein translation. We further find that the higher expression of Gbp2b mRNA in 129/Sv mice correlates with a gene duplication event at the Gbp2b locus resulting in two copies of the Gbp2b gene on the haploid genome of the 129/Sv strain. Conclusions: Our findings demonstrate functional differences between 129 and C57BL/6 Gbp alleles which need to be considered in the design and interpretation of studies utilizing mouse models, particularly for phenotypes influenced by Gbp2 or Gbp2b expression.

A Survey of Copy Number Variation in the Porcine Genome Detected From Whole-Genome Sequence

Copy number variations (CNVs) are gains and losses of large regions of genomic sequence between individuals of a species. Although CNVs have been associated with various phenotypic traits in humans and other species, the extent to which CNVs impact phenotypic variation remains unclear. In swine, as well as many other species, relatively little is understood about the frequency of CNV in the genome, sizes, locations, and other chromosomal properties. In this work, we identified and characterized CNV by utilizing whole-genome sequence from 240 members of an intensely phenotyped experimental swine herd at the U.S. Meat Animal Research Center (USMARC). These animals included all 24 of the purebred founding boars (12 Duroc and 12 Landrace), 48 of the founding Yorkshire-Landrace composite sows, 109 composite animals from generations 4 through 9, 29 composite animals from generation 15, and 30 purebred industry boars (15 Landrace and 15 Yorkshire) used as sires in generations 10 through 15. Using a combination of split reads, paired-end mapping, and read depth approaches, we identified a total of 3,538 copy number variable regions (CNVRs), including 1,820 novel CNVRs not reported in previous studies. The CNVRs covered 0.94% of the porcine genome and overlapped 1,401 genes. Gene ontology analysis identified that CNV-overlapped genes were enriched for functions related to organism development. Additionally, CNVRs overlapped with many known quantitative trait loci (QTL). In particular, analysis of QTL previously identified in the USMARC herd showed that CNVRs were most overlapped with reproductive traits, such as age of puberty and ovulation rate, and CNVRs were significantly enriched for reproductive QTL.

Copy Number Variation of the SHE Gene in Sheep and Its Association with Economic Traits

Simple Summary

Src Homology 2 Domain Containing E (SHE) is a protein coding gene, and member of the SH2 domain-containing family. Sequencing revealed a 2000 bp copy number variation in the SHE gene. There is overlap between SHE copy number variation (CNV) and quantitative trait loci related to milk fat percentage and bone density. A total of 750 sheep, including Chaka sheep (CKS), Hu sheep (HS), Small Tail Han sheep (STHS) and Large Tail Han sheep (LTHS) were available to detect the CNV of the SHE gene and correlate these gene variations with economic traits. The results showed that there were more individuals with SHE copy number loss in CKS and HS than in STHS and LTHS. Association analysis showed that gain and normal copy number types performed better in body length (p < 0.05), circumference of cannon bone (p < 0.05), heart girth (p < 0.05), chest width (p < 0.05) and high at the cross (p < 0.05) in CKS, HS and STHS. Chi-square analyses found significant variation in the CNV of the SHE gene, so it varies greatly between varieties. These findings clarified the relationship between the CNV of the SHE gene and the economic traits in these four kinds of sheep, and provide a reference for sheep breeding.

Abstract

Copy number variation (CNV) caused by gene rearrangement is an important part of genomic structural variation. We found that the copy number variation region of the Src Homology 2 Domain Containing E (SHE) gene correlates with a quantitative trait locus of sheep related to milk fat percentage and bone density. The aim of our study was to detect the copy number variation of the SHE gene in four sheep breeds and to conduct a correlation analysis with economic traits, hoping to provide some reference for sheep breeding. In this study, we examined 750 sheep from four Chinese breeds: Chaka sheep (CKS), Hu sheep (HS), Large Tail Han sheep (LTHS) and Small Tail Han sheep (STHS). We used qPCR to evaluate the copy number of the SHE gene, and then used general linear models to analyze the associations between CNV and economic traits. The results showed that there were more individuals with SHE copy number loss in CKS and HS than in STHS and LTHS individuals. Association analyses showed that gain and normal copy number types were correlated to body length, circumference of cannon bone, heart girth, chest width and high at the cross in CKS, HS and STHS (p < 0.05), but this association was not observed for LTHS. Chi-square values (χ²) found prominent differences in CNV distribution among the studied breeds. Overall, the CNV of the SHE gene may be an important consideration for sheep molecular breeding.

An evolving view of copy number variants

Copy number variants (CNVs) are regions of the genome that vary in integer copy number. CNVs, which comprise both amplifications and deletions of DNA sequence, have been identified across all domains of life, from bacteria and archaea to plants and animals. CNVs are an important source of genetic diversity, and can drive rapid adaptive evolution and progression of heritable and somatic human diseases, such as cancer. However, despite their evolutionary importance and clinical relevance, CNVs remain understudied compared to single-nucleotide variants (SNVs). This is a consequence of the inherent difficulties in detecting CNVs at low-to-intermediate frequencies in heterogeneous populations of cells. Here, we discuss molecular methods used to detect CNVs, the limitations associated with using these techniques, and the application of new and emerging technologies that present solutions to these challenges. The goal of this short review and perspective is to highlight aspects of CNV biology that are understudied and define avenues for further research that address specific gaps in our knowledge of these complex alleles. We describe our recently developed method for CNV detection in which a fluorescent gene functions as a single-cell CNV reporter and present key findings from our evolution experiments in Saccharomyces cerevisiae. Using a CNV reporter, we found that CNVs are generated at a high rate and undergo selection with predictable dynamics across independently evolving replicate populations. Many CNVs appear to be generated through DNA replication-based processes that are mediated by the presence of short, interrupted, inverted-repeat sequences. Our results have important implications for the role of CNVs in evolutionary processes and the molecular mechanisms that underlie CNV formation. We discuss the possible extension of our method to other applications, including tracking the dynamics of CNVs in models of human tumors.

Computational detection and experimental validation of segmental duplications and associated copy number variations in water buffalo ( Bubalus bubalis )

Duplicated sequences are an important source of gene evolution and structural variation within mammalian genomes. Using a read depth approach based on next-generation sequencing, we performed a genome-wide analysis of segmental duplications (SDs) and associated copy number variations (CNVs) in the water buffalo (Bubalus bubalis). By aligning short reads of Olimpia (the reference water buffalo) to the UMD3.1 cattle genome, we identified 1,038 segmental duplications comprising 44.6 Mb (equivalent to ~1.73% of the cattle genome) of the autosomal and X chromosomal sequence in the buffalo genome. We experimentally validated 70.3% (71/101) of these duplications using fluorescent in situ hybridization. We also detected a total of 1,344 CNV regions across 14 additional water buffaloes, amounting to 59.8 Mb of variable sequence or the equivalent of 2.2% of the cattle genome. The CNV regions overlap 1,245 genes that are significantly enriched for specific biological functions including immune response, oxygen transport, sensory system and signal transduction. Additionally, we performed array Comparative Genomic Hybridization (aCGH) experiments using the 14 water buffaloes as test samples and Olimpia as the reference. Using a linear regression model, a high Pearson correlation (r = 0.781) was observed between the log₂ ratios between copy number estimates and the log₂ ratios of aCGH probes. We further designed Quantitative PCR assays to confirm CNV regions within or near annotated genes and found 74.2% agreement with our CNV predictions. These results confirm sub-chromosome-scale structural rearrangements present in the cattle and water buffalo. The information on genome variation that will be of value for evolutionary and phenotypic studies, and may be useful for selective breeding of both species.

Probe-based association analysis identifies several deletions associated with average daily gain in beef cattle

BACKGROUND

Average daily gain (ADG) is an important trait that contributes to the production efficiency and economic benefits in the beef cattle industry. The molecular mechanisms of ADG have not yet been fully explored because most recent association studies for ADG are based on SNPs or haplotypes. We reported a systematic CNV discovery and association analysis for ADG in Chinese Simmental beef cattle.

RESULTS

Our study identified 4912 nonredundant CNVRs with a total length of ~ 248.7 Mb, corresponding to ~ 8.9% of the cattle genome. Using probe-based CNV association, we identified 24 and 12 significant SNP probes within five deletions and two duplications for ADG, respectively. Among them, we found one common deletion with 89 kb imbedded in LHFPL Tetraspan Subfamily Member 6 (LHFPL6) at 22.9 Mb on BTA12, which has high frequency (12.9%) dispersing across population. CNV selection test using V_ST statistic suggested this common deletion may be under positive selection in Chinese Simmental cattle. Moreover, this deletion was not overlapped with any candidate SNP for ADG compared with previous SNPs-based association studies, suggesting its important role for ADG. In addition, we identified one rare deletion near gene Growth Factor Receptor-bound Protein 10 (GRB10) at 5.1 Mb on BTA4 for ADG using both probe-based association and region-based approaches.

CONCLUSIONS

Our results provided some valuable insights to elucidate the genetic basis of ADG in beef cattle, and these findings offer an alternative perspective to understand the genetic mechanism of complex traits in terms of copy number variations in farm animals.

A copy number variant scan in the autochthonous Valdostana Red Pied cattle breed and comparison with specialized dairy populations

Copy number variants (CNVs) are an important source of genomic structural variation, recognized to influence phenotypic variation in many species. Many studies have focused on identifying CNVs within and between human and livestock populations alike, but only few have explored population-genetic properties in cattle based on CNVs derived from a high-density SNP array. We report a high-resolution CNV scan using Illumina’s 777k BovineHD Beadchip for Valdostana Red Pied (VRP), an autochthonous Italian dual-purpose cattle population reared in the Alps that did not undergo strong selection for production traits. After stringent quality control and filtering, CNVs were called across 108 bulls using the PennCNV software. A total of 6,784 CNVs were identified, summarized to 1,723 CNV regions (CNVRs) on 29 autosomes covering a total of ~59 Mb of the UMD3.1 assembly. Among the mapped CNVRs, there were 812 losses, 832 gains and 79 complexes. We subsequently performed a comparison of CNVs detected in the VRP and those available from published studies in the Italian Brown Swiss (IBS) and Mexican Holstein (HOL). A total of 171 CNVRs were common to all three breeds. Between VRP and IBS, 474 regions overlapped, while only 313 overlapped between VRP and HOL, indicating a more similar genetic background among populations with common origins, i.e. the Alps. The principal component, clustering and admixture analyses showed a clear separation of the three breeds into three distinct clusters. In order to describe the distribution of CNVs within and among breeds we used the pair V_ST statistic, considering only the CNVRs shared to more than 5 individuals (within breed). We identified unique and highly differentiated CNVs (n = 33), some of which could be due to specific breed selection and adaptation. Genes and QTL within these regions were characterized.

Structural Genome Variations Related to Craniosynostosis

Craniosynostosis refers to a condition during early development in which one or more of the fibrous sutures of the skull prematurely fuse by turning into bone, which produces recognizable patterns of cranial shape malformations depending on which suture(s) are affected. In addition to cases with isolated cranial dysmorphologies, craniosynostosis appears in syndromes that include skeletal features of the eyes, nose, palate, hands, and feet as well as impairment of vision, hearing, and intellectual development. Approximately 85% of the cases are nonsyndromic sporadic and emerge after de novo structural genome rearrangements or single nucleotide variation, while the remainders consist of syndromic cases following mendelian inheritance. By karyotyping, genome wide linkage, and CNV analyses as well as by whole exome and whole genome sequencing, numerous candidate genes for craniosynostosis belonging to the FGF, Wnt, BMP, Ras/ERK, ephrin, hedgehog, STAT, and retinoic acid signaling pathways have been identified. Many of the craniosynostosis-related candidate genes form a functional network based upon protein-protein or protein-DNA interactions. Depending on which node of this craniosynostosis-related network is affected by a gene mutation or a change in gene expression pattern, a distinct craniosynostosis syndrome or set of phenotypes ensues. Structural variations may alter the dosage of one or several genes or disrupt the genomic architecture of genes and their regulatory elements within topologically associated chromatin domains. These may exert dominant effects by either haploinsufficiency, dominant negative partial loss of function, gain of function, epistatic interaction, or alteration of levels and patterns of gene expression during development. Molecular mechanisms of dominant modes of action of these mutations may include loss of one or several binding sites for cognate protein partners or transcription factor binding sequences. Such losses affect interactions within functional networks governing development and consequently result in phenotypes such as craniosynostosis. Many of the novel variants identified by genome wide CNV analyses, whole exome and whole genome sequencing are incorporated in recently developed diagnostic algorithms for craniosynostosis.

Mutational and transcriptional landscape of spontaneous gene duplications and deletions in Caenorhabditis elegans

Gene duplication and deletion are pivotal processes shaping the structural and functional repertoire of genomes, with implications for disease, adaptation, and evolution. We employed a mutation accumulation (MA) framework partnered with high-throughput genomics to assess the molecular and transcriptional characteristics of newly arisen gene copy-number variants (CNVs) in Caenorhabditis elegans populations subjected to varying intensity of selection. Here, we report a direct spontaneous genome-wide rate of gene duplication of 2.9 × 10^-5/gene per generation in C. elegans, the highest for any species to date. The rate of gene deletion is sixfold lower (5 × 10^-6/gene per generation). Deletions of highly expressed genes are particularly deleterious, given their paucity in even the N = 1 lines with minimal efficacy of selection. The increase in average transcript abundance of new duplicates arising under minimal selection is significantly greater than twofold compared with single copies of the same gene, suggesting that genes in segmental duplications are frequently overactive at inception. The average increase in transcriptional activity of gene duplicates is greater in the N = 1 MA lines than in MA lines with larger population bottlenecks. There is an inverse relationship between the ancestral transcription levels of new gene duplicates and population size, with duplicate copies of highly expressed genes less likely to accumulate in larger populations. Our results demonstrate a fitness cost of increased transcription following duplication, which results in purifying selection against new gene duplicates. However, on average, duplications also provide a significant increase in gene expression that can facilitate adaptation to novel environmental challenges.

Divergent patterns of genic copy number variation in KCNIP1 gene reveal risk locus of type 2 diabetes in Chinese population

Copy number variation (CNV) has emerged as another important genetic marker in addition to SNP for understanding etiology of complex disease. Kv channel interacting protein 1 (KCNIP1) is a Ca²⁺-dependent transcriptional modulator that contributes to the regulation of insulin secretion. Previous genome-wide CNV assay identified the KCNIP1 gene encompassing a CNV region, however, its further effect and risk rate on type 2 diabetes (T2D) have rarely been addressed, especially in Chinese population. The current study aims to detect and excavate genetic distribution profile of KCNIP1 CNV in Chinese T2D and control populations, and further to investigate the associations with clinical characteristics. Divergent patterns of the KCNIP1 CNV were identified (p < 0.01), in which the copy number gain was predominant in T2D, while the copy number normal accounted for the most in control group. Consistently, the individuals with copy number gain showed significant risk on T2D (OR = 4.550, p < 0.01). The KCNIP1 copy numbers presented significantly positive correlations with fasting plasma glucose and glycated hemoglobin in T2D. For OGTT test, the T2D patients with copy number gain had remarkably elevated glucose contents (60, 120, 180-min, p < 0.05 or p < 0.01) and diminished insulin levels (60, 120-min, p < 0.05) than those with copy number loss and normal, which suggested that the KCNIP1 CNV was correlated with the glucose and insulin action. This is the first CNV association study of the KCNIP1 gene in Chinese population, and these data indicated that KCNIP1 might function as a T2D-susceptibility gene whose dysregulation alters insulin production.

Whole-Body Single-Cell Sequencing Reveals Transcriptional Domains in the Annelid Larval Body

Animal bodies comprise diverse arrays of cells. To characterize cellular identities across an entire body, we have compared the transcriptomes of single cells randomly picked from dissociated whole larvae of the marine annelid Platynereis dumerilii. We identify five transcriptionally distinct groups of differentiated cells, each expressing a unique set of transcription factors and effector genes that implement cellular phenotypes. Spatial mapping of cells into a cellular expression atlas, and wholemount in situ hybridization of group-specific genes reveals spatially coherent transcriptional domains in the larval body, comprising, for example, apical sensory-neurosecretory cells versus neural/epidermal surface cells. These domains represent new, basic subdivisions of the annelid body based entirely on differential gene expression, and are composed of multiple, transcriptionally similar cell types. They do not represent clonal domains, as revealed by developmental lineage analysis. We propose that the transcriptional domains that subdivide the annelid larval body represent families of related cell types that have arisen by evolutionary diversification. Their possible evolutionary conservation makes them a promising tool for evo-devo research.

Comparative sequence alignment reveals River Buffalo genomic structural differences compared with cattle

This study sought to characterize differences in gene content, regulation and structure between taurine cattle and river buffalo (one subspecies of domestic water buffalo) using the extensively annotated UMD3.1 cattle reference genome as a basis for comparisons. We identified 127 deletion CNV regions in river buffalo representing 5 annotated cattle genes. We also characterized 583 merged mobile element insertion (MEI) events within the upstream regions of annotated cattle genes. Transcriptome analysis in various tissue types on river buffalo confirmed the absence of four cattle genes. Four genes which may be related to phenotypic differences in meat quality and color, had upstream MEI predictions and were found to have significantly elevated expression in river buffalo compared with cattle. Our comparative alignment approach and gene expression analyses suggested a functional role for many genomic structural variations, which may contribute to the unique phenotypes of river buffalo.