Structural variation and its effect on expression

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.

TADeus2: a web server facilitating the clinical diagnosis by pathogenicity assessment of structural variations disarranging 3D chromatin structure

In recent years great progress has been made in identification of structural variants (SV) in the human genome. However, the interpretation of SVs, especially located in non-coding DNA, remains challenging. One of the reasons stems in the lack of tools exclusively designed for clinical SVs evaluation acknowledging the 3D chromatin architecture. Therefore, we present TADeus2 a web server dedicated for a quick investigation of chromatin conformation changes, providing a visual framework for the interpretation of SVs affecting topologically associating domains (TADs). This tool provides a convenient visual inspection of SVs, both in a continuous genome view as well as from a rearrangement’s breakpoint perspective. Additionally, TADeus2 allows the user to assess the influence of analyzed SVs within flaking coding/non-coding regions based on the Hi-C matrix. Importantly, the SVs pathogenicity is quantified and ranked using TADA, ClassifyCNV tools and sampling-based P-value. TADeus2 is publicly available at https://tadeus2.mimuw.edu.pl.

Understanding the Early Evolutionary Stages of a Tandem Drosophilamelanogaster-Specific Gene Family: A Structural and Functional Population Study

Gene families underlie genetic innovation and phenotypic diversification. However, our understanding of the early genomic and functional evolution of tandemly arranged gene families remains incomplete as paralog sequence similarity hinders their accurate characterization. The Drosophila melanogaster-specific gene family Sdic is tandemly repeated and impacts sperm competition. We scrutinized Sdic in 20 geographically diverse populations using reference-quality genome assemblies, read-depth methodologies, and qPCR, finding that ∼90% of the individuals harbor 3-7 copies as well as evidence of population differentiation. In strains with reliable gene annotations, copy number variation (CNV) and differential transposable element insertions distinguish one structurally distinct version of the Sdic region per strain. All 31 annotated copies featured protein-coding potential and, based on the protein variant encoded, were categorized into 13 paratypes differing in their 3' ends, with 3-5 paratypes coexisting in any strain examined. Despite widespread gene conversion, the only copy present in all strains has functionally diverged at both coding and regulatory levels under positive selection. Contrary to artificial tandem duplications of the Sdic region that resulted in increased male expression, CNV in cosmopolitan strains did not correlate with expression levels, likely as a result of differential genome modifier composition. Duplicating the region did not enhance sperm competitiveness, suggesting a fitness cost at high expression levels or a plateau effect. Beyond facilitating a minimally optimal expression level, Sdic CNV acts as a catalyst of protein and regulatory diversity, showcasing a possible evolutionary path recently formed tandem multigene families can follow toward long-term consolidation in eukaryotic genomes.

Variation in the Evolution and Sequences of Proglucagon and the Receptors for Proglucagon-Derived Peptides in Mammals

The mammalian proglucagon gene (Gcg) encodes three glucagon like sequences, glucagon, glucagon-like peptide-1 (GLP-1), and glucagon-like peptide-2 that are of similar length and share sequence similarity, with these hormones having cell surface receptors, glucagon receptor (Gcgr), GLP-1 receptor (Glp1r), and GLP-2 receptor (Glp2r), respectively. Gcgr, Glp1r, and Glp2r are all class B1 G protein-coupled receptors (GPCRs). Despite their sequence and structural similarity, analyses of sequences from rodents have found differences in patterns of sequence conservation and evolution. To determine whether these were rodent-specific traits or general features of these genes in mammals I analyzed coding and protein sequences for proglucagon and the receptors for proglucagon-derived peptides from the genomes of 168 mammalian species. Single copy genes for each gene were found in almost all genomes. In addition to glucagon sequences within Hystricognath rodents (e.g., guinea pig), glucagon sequences from a few other groups (e.g., pangolins and some bats) as well as changes in the proteolytic processing of GLP-1 in some bats are suggested to have functional effects. GLP-2 sequences display increased variability but accepted few substitutions that are predicted to have functional consequences. In parallel, Glp2r sequences display the most rapid protein sequence evolution, and show greater variability in amino acids at sites involved in ligand interaction, however most were not predicted to have a functional consequence. These observations suggest that a greater diversity in biological functions for proglucagon-derived peptides might exist in mammals.

The Future of Livestock Management: A Review of Real-Time Portable Sequencing Applied to Livestock

Oxford Nanopore Technologies' MinION has proven to be a valuable tool within human and microbial genetics. Its capacity to produce long reads in real time has opened up unique applications for portable sequencing. Examples include tracking the recent African swine fever outbreak in China and providing a diagnostic tool for disease in the cassava plant in Eastern Africa. Here we review the current applications of Oxford Nanopore sequencing in livestock, then focus on proposed applications in livestock agriculture for rapid diagnostics, base modification detection, reference genome assembly and genomic prediction. In particular, we propose a future application: 'crush-side genotyping' for real-time on-farm genotyping for extensive industries such as northern Australian beef production. An initial in silico experiment to assess the feasibility of crush-side genotyping demonstrated promising results. SNPs were called from simulated Nanopore data, that included the relatively high base call error rate that is characteristic of the data, and calling parameters were varied to understand the feasibility of SNP calling at low coverages in a heterozygous population. With optimised genotype calling parameters, over 85% of the 10,000 simulated SNPs were able to be correctly called with coverages as low as 6×. These results provide preliminary evidence that Oxford Nanopore sequencing has potential to be used for real-time SNP genotyping in extensive livestock operations.

Vertebrate Chromosome Evolution

The study of chromosome evolution is undergoing a resurgence of interest owing to advances in DNA sequencing technology that facilitate the production of chromosome-scale whole-genome assemblies de novo. This review focuses on the history, methods, discoveries, and current challenges facing the field, with an emphasis on vertebrate genomes. A detailed examination of the literature on the biology of chromosome rearrangements is presented, specifically the relationship between chromosome rearrangements and phenotypic evolution, adaptation, and speciation. A critical review of the methods for identifying, characterizing, and visualizing chromosome rearrangements and computationally reconstructing ancestral karyotypes is presented. We conclude by looking to the future, identifying the enormous technical and scientific challenges presented by the accumulation of hundreds and eventually thousands of chromosome-scale assemblies.

Analysis of structural variants in four African cichlids highlights an association with developmental and immune related genes

BACKGROUND

East African lake cichlids are one of the most impressive examples of an adaptive radiation. Independently in Lake Victoria, Tanganyika, and Malawi, several hundreds of species arose within the last 10 million to 100,000 years. Whereas most analyses in cichlids focused on nucleotide substitutions across species to investigate the genetic bases of this explosive radiation, to date, no study has investigated the contribution of structural variants (SVs) in the evolution of adaptive traits across the three Great Lakes of East Africa.

RESULTS

Here, we annotate and characterize the repertoires and evolutionary potential of different SV classes (deletion, duplication, inversion, insertions and translocations) in four cichlid species: Haplochromis burtoni, Metriaclima zebra, Neolamprologus brichardi and Pundamilia nyererei. We investigate the patterns of gain and loss evolution for each SV type, enabling the identification of lineage specific events. Both deletions and inversions show a significant overlap with SINE elements, while inversions additionally show a limited, but significant association with DNA transposons. Inverted regions are enriched for genes regulating behaviour, or involved in skeletal and visual system development. We also find that duplicated regions show enrichment for genes associated with "antigen processing and presentation" and other immune related categories. Our pipeline and results were further tested by PCR validation of selected deletions and inversions, which confirmed respectively 7 out of 10 and 6 out of 9 events.

CONCLUSIONS

Altogether, we provide the first comprehensive overview of rearrangement evolution in East African cichlids, and some important insights into their likely contribution to adaptation.

Breakpoint Mapping of Symptomatic Balanced Translocations Links the EPHA6, KLF13 and UBR3 Genes to Novel Disease Phenotype

De novo balanced chromosomal aberrations (BCAs), such as reciprocal translocations and inversions, are genomic aberrations that, in approximately 25% of cases, affect the human phenotype. Delineation of the exact structure of BCAs may provide a precise diagnosis and/or point to new disease loci. We report on six patients with de novo balanced chromosomal translocations (BCTs) and one patient with a de novo inversion, in whom we mapped breakpoints to a resolution of 1 bp, using shallow whole-genome mate pair sequencing. In all seven cases, a disruption of at least one gene was found. In two patients, the phenotypic impact of the disrupted genes is well known (NFIA, ATP7A). In five patients, the aberration damaged genes: PARD3, EPHA6, KLF13, STK24, UBR3, MLLT10 and TLE3, whose influence on the human phenotype is poorly understood. In particular, our results suggest novel candidate genes for retinal degeneration with anophthalmia (EPHA6), developmental delay with speech impairment (KLF13), and developmental delay with brain dysembryoplastic neuroepithelial tumor (UBR3). In conclusion, identification of the exact structure of symptomatic BCTs using next generation sequencing is a viable method for both diagnosis and finding novel disease candidate genes in humans.

Characterization and evolutionary dynamics of complex regions in eukaryotic genomes

Complex regions in eukaryotic genomes are typically characterized by duplications of chromosomal stretches that often include one or more genes repeated in a tandem array or in relatively close proximity. Nevertheless, the repetitive nature of these regions, together with the often high sequence identity among repeats, have made complex regions particularly recalcitrant to proper molecular characterization, often being misassembled or completely absent in genome assemblies. This limitation has prevented accurate functional and evolutionary analyses of these regions. This is becoming increasingly relevant as evidence continues to support a central role for complex genomic regions in explaining human disease, developmental innovations, and ecological adaptations across phyla. With the advent of long-read sequencing technologies and suitable assemblers, the development of algorithms that can accommodate sample heterozygosity, and the adoption of a pangenomic-like view of these regions, accurate reconstructions of complex regions are now within reach. These reconstructions will finally allow for accurate functional and evolutionary studies of complex genomic regions, underlying the generation of genotype-phenotype maps of unprecedented resolution.

Local and global chromatin interactions are altered by large genomic deletions associated with human brain development

Large copy number variants (CNVs) in the human genome are strongly associated with common neurodevelopmental, neuropsychiatric disorders such as schizophrenia and autism. Here we report on the epigenomic effects of the prominent large deletion CNVs on chromosome 22q11.2 and on chromosome 1q21.1. We use Hi-C analysis of long-range chromosome interactions, including haplotype-specific Hi-C analysis, ChIP-Seq analysis of regulatory histone marks, and RNA-Seq analysis of gene expression patterns. We observe changes on all the levels of analysis, within the deletion boundaries, in the deletion flanking regions, along chromosome 22q, and genome wide. We detect gene expression changes as well as pronounced and multilayered effects on chromatin states, chromosome folding and on the topological domains of the chromatin, that emanate from the large CNV locus. These findings suggest basic principles of how such large genomic deletions can alter nuclear organization and affect genomic molecular activity.

Copy number variants in the human genome (CNVs) are associated with neurodevelopmental and psychiatric disorders such as schizophrenia and autism. Here the authors investigate how  the large deletion CNV on chromosome 22q11.2 alters chromatin organization.

Mapping of breakpoints in balanced chromosomal translocations by shallow whole-genome sequencing points to EFNA5, BAHD1 and PPP2R5E as novel candidates for genes causing human Mendelian disorders

BACKGROUND

Mapping the breakpoints in de novo balanced chromosomal translocations (BCT) in symptomatic individuals provides a unique opportunity to identify in an unbiased way the likely causative genetic defect and thus find novel human disease candidate genes. Our aim was to fine-map breakpoints of de novo BCTs in a case series of nine patients.

METHODS

Shallow whole-genome mate pair sequencing (SGMPS) together with long-range PCR and Sanger sequencing. In one case (BCT disrupting BAHD1 and RET) cDNA analysis was used to verify expression of a fusion transcript in cultured fibroblasts.

RESULTS

In all nine probands 11 disrupted genes were found, that is, EFNA5, EBF3, LARGE, PPP2R5E, TXNDC5, ZNF423, NIPBL, BAHD1, RET, TRPS1 and SLC4A10. Five subjects had translocations that disrupted genes with so far unknown (EFNA5, BAHD1, PPP2R5E, TXNDC5) or poorly delineated impact on the phenotype (SLC4A10, two previous reports of BCT disrupting the gene). The four genes with no previous disease associations (EFNA5, BAHD1, PPP2R5E, TXNDC5), when compared with all human genes by a bootstrap test, had significantly higher pLI (p<0.017) and DOMINO (p<0.02) scores indicating enrichment in genes likely to be intolerant to single copy damage. Inspection of individual pLI and DOMINO scores, and local topologically associating domain structure suggested that EFNA5, BAHD1 and PPP2R5E were particularly good candidates for novel disease loci. The pathomechanism for BAHD1 may involve deregulation of expression due to fusion with RET promoter.

CONCLUSION

SGMPS in symptomatic carriers of BCTs is a powerful approach to delineate novel human gene-disease associations.

Genotype/phenotype correlation in a female patient with 21q22.3 and 12p13.33 duplications

Many chromosomal rearrangements that lead to copy-number gains or losses have been shown to cause distinctive and recognizable clinical phenotypes. Conventional cytogenetic analysis can detect many, but not all, rearrangements depending on its power of resolution. The wide use of whole-genome array-based comparative genomic hybridization (array-CGH) techniques has allowed the detection of novel syndromes and to establish genotype-phenotype correlations by delineating at high resolution the regions involved in specific chromosomal aberrations. We report on a two and half-year-old female patient with intellectual disability and distinctive phenotypic features resulting from a de novo duplication of about 0.3 Mb in 21q22.3 associated with duplication of about 0.3 Mb in 12p13.33. The patient's chromosomal abnormalities were identified at the cytogenetic molecular level, using SNP array analysis, while GTG banding technique revealed a normal karyotype. Clinical findings of the patient were compared with Down syndrome and 12p duplication syndrome. This study suggests that an area of contiguous genes on the distal part of chromosome 21 (21q22.3) contribute to the Down syndrome phenotype and indicates that genes in the distal region of 12p (12p13.33) account for many facial characteristics and hypotonia of trisomy 12p syndrome.

Identification of structural variation in mouse genomes

Structural variation is variation in structure of DNA regions affecting DNA sequence length and/or orientation. It generally includes deletions, insertions, copy-number gains, inversions, and transposable elements. Traditionally, the identification of structural variation in genomes has been challenging. However, with the recent advances in high-throughput DNA sequencing and paired-end mapping (PEM) methods, the ability to identify structural variation and their respective association to human diseases has improved considerably. In this review, we describe our current knowledge of structural variation in the mouse, one of the prime model systems for studying human diseases and mammalian biology. We further present the evolutionary implications of structural variation on transposable elements. We conclude with future directions on the study of structural variation in mouse genomes that will increase our understanding of molecular architecture and functional consequences of structural variation.

The impact of chromosomal rearrangements on regulation of gene expression

The effects that coding region single-nucleotide polymorphisms or mutations have on gene expression have been well documented, predominantly owing to their association with disease. The effects of structural chromosomal rearrangements are also receiving increasing attention with the development of new techniques that allow accurate, high-resolution data, whether genomic interaction or transcriptome data, to be generated right down to the single-cell level. Over the past 18 months, these advances in experimental techniques have been used to further confirm and delineate the substantial effects that chromosome rearrangements can have on the regulation of gene expression and provide evidence of direct links between the two.

Evolution of genomic structural variation and genomic architecture in the adaptive radiations of African cichlid fishes

African cichlid fishes are an ideal system for studying explosive rates of speciation and the origin of diversity in adaptive radiation. Within the last few million years, more than 2000 species have evolved in the Great Lakes of East Africa, the largest adaptive radiation in vertebrates. These young species show spectacular diversity in their coloration, morphology and behavior. However, little is known about the genomic basis of this astonishing diversity. Recently, five African cichlid genomes were sequenced, including that of the Nile Tilapia (Oreochromis niloticus), a basal and only relatively moderately diversified lineage, and the genomes of four representative endemic species of the adaptive radiations, Neolamprologus brichardi, Astatotilapia burtoni, Metriaclima zebra, and Pundamila nyererei. Using the Tilapia genome as a reference genome, we generated a high-resolution genomic variation map, consisting of single nucleotide polymorphisms (SNPs), short insertions and deletions (indels), inversions and deletions. In total, around 18.8, 17.7, 17.0, and 17.0 million SNPs, 2.3, 2.2, 1.4, and 1.9 million indels, 262, 306, 162, and 154 inversions, and 3509, 2705, 2710, and 2634 deletions were inferred to have evolved in N. brichardi, A. burtoni, P. nyererei, and M. zebra, respectively. Many of these variations affected the annotated gene regions in the genome. Different patterns of genetic variation were detected during the adaptive radiation of African cichlid fishes. For SNPs, the highest rate of evolution was detected in the common ancestor of N. brichardi, A. burtoni, P. nyererei, and M. zebra. However, for the evolution of inversions and deletions, we found that the rates at the terminal taxa are substantially higher than the rates at the ancestral lineages. The high-resolution map provides an ideal opportunity to understand the genomic bases of the adaptive radiation of African cichlid fishes.

Impacts of variation in the human genome on gene regulation

Recent advances in fast and inexpensive DNA sequencing have enabled the extensive study of genomic and transcriptomic variation in humans. Human genomic variation is composed of sequence and structural changes including single-nucleotide and multinucleotide variants, short insertions or deletions (indels), larger copy number variants, and similarly sized copy neutral inversions and translocations. It is now well established that any two genomes differ extensively and that structural changes constitute a prominent source of this variation. There have also been major technological advances in RNA sequencing to globally quantify and describe diversity in transcripts. Large consortia such as the 1000 Genomes Project and the ENCODE (ENCyclopedia Of DNA Elements) Project are producing increasingly comprehensive maps outlining the regions of the human genome containing variants and functional elements, respectively. Integration of genetic variation data and extensive annotation of functional genomic elements, along with the ability to measure global transcription, allow the impacts of genetic variants on gene expression to be resolved. There are several well-established models by which genetic variants affect gene regulation depending on the type, nature, and position of the variant with respect to the affected genes. These effects can be manifested in two ways: changes to transcript sequences and isoforms by coding variants, and changes to transcript abundance by dosage or regulatory variants. Here, we review the current state of how genetic variations impact gene regulation locally and globally in the human genome.

Duplication of 20p12.3 associated with familial Wolff-Parkinson-White syndrome

Wolff-Parkinson-White (WPW) syndrome is caused by preexcitation of the ventricular myocardium via an accessory pathway which increases the risk for paroxysmal supraventricular tachycardia. The condition is often sporadic and of unknown etiology in the majority of cases. Autosomal dominant inheritance and association with congenital heart defects or ventricular hypertrophy were described. Microdeletions of 20p12.3 have been associated with WPW syndrome with either cognitive dysfunction or Alagille syndrome. Here, we describe the association of 20p12.3 duplication with WPW syndrome in a patient who presented with non-immune hydrops. Her paternal uncle carries the duplication and has attention-deficit hyperactivity disorder and electrocardiographic findings consistent with WPW. The 769 kb duplication was detected by the Affymetrix Whole Genome-Human SNP Array 6.0 and encompasses two genes and the first two exons of a third gene. We discuss the potential role of the genes in the duplicated region in the pathogenesis of WPW and possible neurobehavioral abnormalities. Our data provide additional support for a significant role of 20p12.3 chromosomal rearrangements in the etiology of WPW syndrome.

Genome-wide patterns of standing genetic variation in a marine population of three-spined sticklebacks

Since the end of the Pleistocene, the three-spined stickleback (Gasterosteus aculeatus) has repeatedly colonized and adapted to various freshwater habitats probably originating from ancestral marine populations. Standing genetic variation and the underlying genomic architecture both have been speculated to contribute to recent adaptive radiations of sticklebacks. Here, we expand on the current genomic resources of this fish by providing extensive genome-wide variation data from six individuals from a marine (North Sea) stickleback population. Using next-generation sequencing and a combination of paired-end and mate-pair libraries, we detected a wide size range of genetic variation. Among the six individuals, we found more than 7% of the genome is polymorphic, consisting of 2599111 SNPs, 233464 indels and structural variation (SV) (>50 bp) such as 1054 copy-number variable regions (deletions and duplications) and 48 inversions. Many of these polymorphisms affect gene and coding sequences. Based on SNP diversity, we determined outlier regions concordant with signatures expected under adaptive evolution. As some of these outliers overlap with pronounced regions of copy-number variation, we propose the consideration of such SV when analysing SNP data from re-sequencing approaches. We further discuss the value of this resource on genome-wide variation for further investigation upon the relative contribution of standing variation on the parallel evolution of sticklebacks and the importance of the genomic architecture in adaptive radiation.