Primate segmental duplications: crucibles of evolution, diversity and disease

EFL1 mutations impair eIF6 release to cause Shwachman-Diamond syndrome

Shwachman-Diamond syndrome (SDS) is a recessive disorder typified by bone marrow failure and predisposition to hematological malignancies. SDS is predominantly caused by deficiency of the allosteric regulator Shwachman-Bodian-Diamond syndrome that cooperates with elongation factor-like GTPase 1 (EFL1) to catalyze release of the ribosome antiassociation factor eIF6 and activate translation. Here, we report biallelic mutations in EFL1 in 3 unrelated individuals with clinical features of SDS. Cellular defects in these individuals include impaired ribosomal subunit joining and attenuated global protein translation as a consequence of defective eIF6 eviction. In mice, Efl1 deficiency recapitulates key aspects of the SDS phenotype. By identifying biallelic EFL1 mutations in SDS, we define this leukemia predisposition disorder as a ribosomopathy that is caused by corruption of a fundamental, conserved mechanism, which licenses entry of the large ribosomal subunit into translation.

Structural variants identified by Oxford Nanopore PromethION sequencing of the human genome

We sequenced the genome of the Yoruban reference individual NA19240 on the long-read sequencing platform Oxford Nanopore PromethION for evaluation and benchmarking of recently published aligners and germline structural variant calling tools, as well as a comparison with the performance of structural variant calling from short-read sequencing data. The structural variant caller Sniffles after NGMLR or minimap2 alignment provides the most accurate results, but additional confidence or sensitivity can be obtained by a combination of multiple variant callers. Sensitive and fast results can be obtained by minimap2 for alignment and a combination of Sniffles and SVIM for variant identification. We describe a scalable workflow for identification, annotation, and characterization of tens of thousands of structural variants from long-read genome sequencing of an individual or population. By discussing the results of this well-characterized reference individual, we provide an approximation of what can be expected in future long-read sequencing studies aiming for structural variant identification.

Evolutionary history of the human multigene families reveals widespread gene duplications throughout the history of animals

BACKGROUND

The hypothesis that vertebrates have experienced two ancient, whole genome duplications (WGDs) is of central interest to evolutionary biology and has been implicated in evolution of developmental complexity. Three-way and Four-way paralogy regions in human and other vertebrate genomes are considered as vital evidence to support this hypothesis. Alternatively, it has been proposed that such paralogy regions are created by small-scale duplications that occurred at different intervals over the evolution of life.

RESULTS

To address this debate, the present study investigates the evolutionary history of multigene families with at least three-fold representation on human chromosomes 1, 2, 8 and 20. Phylogenetic analysis and the tree topology comparisons classified the members of 36 multigene families into four distinct co-duplicated groups. Gene families falling within the same co-duplicated group might have duplicated together, whereas genes belong to different co-duplicated groups might have distinct evolutionary origins.

CONCLUSION

Taken together with previous investigations, the current study yielded no proof in favor of WGDs hypothesis. Rather, it appears that the vertebrate genome evolved as a result of small-scale duplication events, that cover the entire span of the animals' history.

Phylogenetic analyses of human 1/2/8/20 paralogons suggest segmental duplications during animal evolution

Susumu Ohno hypothesized that the diversity of vertebrate gene families and genome is due to two rounds of whole genome duplications (also referred as 2R hypothesis). The quadruplicate paralogous blocks present on 1/2/8/20 chromosomes are taken as one of the evidences in favor of the 2R. In this study, we investigated that whether 2R has shaped the vertebrate evolution using gene families residing on chromosomes 1/2/8/20. Evolutionary history of 22 gene families (11 from the current study and 11 from the previous study) was evaluated by the phylogenetic analysis with triplicated or quadruplicated distribution on these human chromosomes 1/2/8/20. The phylogenetic analysis was performed using high-quality whole genomic sequence data of multiple species with neighbor-joining (NJ) and maximum likelihood (ML) methods. The phylogenetic tree topology of these gene families revealed variable duplication time points during invertebrate-vertebrate evolution. Topology comparison approach categorized 22 gene families into three groups. Tree topologies of ten gene families fell into Group 1 (duplications prior to invertebrate-vertebrate split), four in Group 2 (i.e., (AB) (C) (D), topology incongruent with 2R) and eight in Group 3 (((AB) (CD)), 2R congruent topology). Therefore, taken together the current and previous data of 1/2/8/20 paralogons, we propose that, in addition to whole genome duplications events, current developmental, morphological and genomic complexity of the vertebrate genomes may also have originated through segmental duplications occurring at varying time points during the course of animal evolution.

Multi-platform discovery of haplotype-resolved structural variation in human genomes

The incomplete identification of structural variants (SVs) from whole-genome sequencing data limits studies of human genetic diversity and disease association. Here, we apply a suite of long-read, short-read, strand-specific sequencing technologies, optical mapping, and variant discovery algorithms to comprehensively analyze three trios to define the full spectrum of human genetic variation in a haplotype-resolved manner. We identify 818,054 indel variants (<50 bp) and 27,622 SVs (≥50 bp) per genome. We also discover 156 inversions per genome and 58 of the inversions intersect with the critical regions of recurrent microdeletion and microduplication syndromes. Taken together, our SV callsets represent a three to sevenfold increase in SV detection compared to most standard high-throughput sequencing studies, including those from the 1000 Genomes Project. The methods and the dataset presented serve as a gold standard for the scientific community allowing us to make recommendations for maximizing structural variation sensitivity for future genome sequencing studies.

Multi-platform discovery of haplotype-resolved structural variation in human genomes

The incomplete identification of structural variants (SVs) from whole-genome sequencing data limits studies of human genetic diversity and disease association. Here, we apply a suite of long-read, short-read, strand-specific sequencing technologies, optical mapping, and variant discovery algorithms to comprehensively analyze three trios to define the full spectrum of human genetic variation in a haplotype-resolved manner. We identify 818,054 indel variants (<50 bp) and 27,622 SVs (≥50 bp) per genome. We also discover 156 inversions per genome and 58 of the inversions intersect with the critical regions of recurrent microdeletion and microduplication syndromes. Taken together, our SV callsets represent a three to sevenfold increase in SV detection compared to most standard high-throughput sequencing studies, including those from the 1000 Genomes Project. The methods and the dataset presented serve as a gold standard for the scientific community allowing us to make recommendations for maximizing structural variation sensitivity for future genome sequencing studies.

Fine-Scale Characterization of Genomic Structural Variation in the Human Genome Reveals Adaptive and Biomedically Relevant Hotspots

Genomic structural variants (SVs) are distributed nonrandomly across the human genome. The "hotspots" of SVs have been implicated in evolutionary innovations, as well as medical conditions. However, the evolutionary and biomedical features of these hotspots remain incompletely understood. Here, we analyzed data from 2,504 genomes to construct a refined map of 1,148 SV hotspots in human genomes. We confirmed that segmental duplication-related nonallelic homologous recombination is an important mechanistic driver of SV hotspot formation. However, to our surprise, we also found that a majority of SVs in hotspots do not form through such recombination-based mechanisms, suggesting diverse mechanistic and selective forces shaping hotspots. Indeed, our evolutionary analyses showed that the majority of SV hotspots are within gene-poor regions and evolve under relaxed negative selection or neutrality. However, we still found a small subset of SV hotspots harboring genes that are enriched for anthropologically crucial functions and evolve under geography-specific and balancing adaptive forces. These include two independent hotspots on different chromosomes affecting alpha and beta hemoglobin gene clusters. Biomedically, we found that the SV hotspots coincide with breakpoints of clinically relevant, large de novo SVs, significantly more often than genome-wide expectations. For example, we showed that the breakpoints of multiple large SVs, which lead to idiopathic short stature, coincide with SV hotspots. Therefore, the mutational instability in SV hotpots likely enables chromosomal breaks that lead to pathogenic structural variation formations. Overall, our study contributes to a better understanding of the mutational and adaptive landscape of the genome.

GenTree, an integrated resource for analyzing the evolution and function of primate-specific coding genes

The origination of new genes contributes to phenotypic evolution in humans. Two major challenges in the study of new genes are the inference of gene ages and annotation of their protein-coding potential. To tackle these challenges, we created GenTree, an integrated online database that compiles age inferences from three major methods together with functional genomic data for new genes. Genome-wide comparison of the age inference methods revealed that the synteny-based pipeline (SBP) is most suited for recently duplicated genes, whereas the protein-family–based methods are useful for ancient genes. For SBP-dated primate-specific protein-coding genes (PSGs), we performed manual evaluation based on published PSG lists and showed that SBP generated a conservative data set of PSGs by masking less reliable syntenic regions. After assessing the coding potential based on evolutionary constraint and peptide evidence from proteomic data, we curated a list of 254 PSGs with different levels of protein evidence. This list also includes 41 candidate misannotated pseudogenes that encode primate-specific short proteins. Coexpression analysis showed that PSGs are preferentially recruited into organs with rapidly evolving pathways such as spermatogenesis, immune response, mother–fetus interaction, and brain development. For brain development, primate-specific KRAB zinc-finger proteins (KZNFs) are specifically up-regulated in the mid-fetal stage, which may have contributed to the evolution of this critical stage. Altogether, hundreds of PSGs are either recruited to processes under strong selection pressure or to processes supporting an evolving novel organ.

Genomic inversions and GOLGA core duplicons underlie disease instability at the 15q25 locus

Human chromosome 15q25 is involved in several disease-associated structural rearrangements, including microdeletions and chromosomal markers with inverted duplications. Using comparative fluorescence in situ hybridization, strand-sequencing, single-molecule, real-time sequencing and Bionano optical mapping analyses, we investigated the organization of the 15q25 region in human and nonhuman primates. We found that two independent inversions occurred in this region after the fission event that gave rise to phylogenetic chromosomes XIV and XV in humans and great apes. One of these inversions is still polymorphic in the human population today and may confer differential susceptibility to 15q25 microdeletions and inverted duplications. The inversion breakpoints map within segmental duplications containing core duplicons of the GOLGA gene family and correspond to the site of an ancestral centromere, which became inactivated about 25 million years ago. The inactivation of this centromere likely released segmental duplications from recombination repression typical of centromeric regions. We hypothesize that this increased the frequency of ectopic recombination creating a hotspot of hominid inversions where dispersed GOLGA core elements now predispose this region to recurrent genomic rearrangements associated with disease.

Evolution of gene regulation in ruminants differs between evolutionary breakpoint regions and homologous synteny blocks

The role of chromosome rearrangements in driving evolution has been a long-standing question of evolutionary biology. Here we focused on ruminants as a model to assess how rearrangements may have contributed to the evolution of gene regulation. Using reconstructed ancestral karyotypes of Cetartiodactyls, Ruminants, Pecorans, and Bovids, we traced patterns of gross chromosome changes. We found that the lineage leading to the ruminant ancestor after the split from other cetartiodactyls was characterized by mostly intrachromosomal changes, whereas the lineage leading to the pecoran ancestor (including all livestock ruminants) included multiple interchromosomal changes. We observed that the liver cell putative enhancers in the ruminant evolutionary breakpoint regions are highly enriched for DNA sequences under selective constraint acting on lineage-specific transposable elements (TEs) and a set of 25 specific transcription factor (TF) binding motifs associated with recently active TEs. Coupled with gene expression data, we found that genes near ruminant breakpoint regions exhibit more divergent expression profiles among species, particularly in cattle, which is consistent with the phylogenetic origin of these breakpoint regions. This divergence was significantly greater in genes with enhancers that contain at least one of the 25 specific TF binding motifs and located near bovidae-to-cattle lineage breakpoint regions. Taken together, by combining ancestral karyotype reconstructions with analysis of cis regulatory element and gene expression evolution, our work demonstrated that lineage-specific regulatory elements colocalized with gross chromosome rearrangements may have provided valuable functional modifications that helped to shape ruminant evolution.

Mobile elements contribute to the uniqueness of human genome with 15,000 human-specific insertions and 14 Mbp sequence increase

Mobile elements (MEs) collectively contribute to at least 50% of the human genome. Due to their past incremental accumulation and ongoing DNA transposition, MEs serve as a significant source for both inter- and intra-species genetic and phenotypic diversity during primate and human evolution. By making use of the most recent genome sequences for human and many other closely related primates and robust multi-way comparative genomic approach, we identified a total of 14,870 human-specific MEs (HS-MEs) with more than 8,000 being newly identified. Collectively, these HS-MEs contribute to a total of 14.2 Mbp net genome sequence increase. Several new observations were made based on these HS-MEs, including the finding of Y chromosome as a strikingly hot target for HS-MEs and a strong mutual preference for SINE-R/VNTR/Alu (SVAs). Furthermore, ∼8,000 of these HS-MEs were found to locate in the vicinity of ∼4,900 genes, and collectively they contribute to ∼84 kb sequences in the human reference transcriptome in association with over 300 genes, including protein-coding sequences for 40 genes. In conclusion, our results demonstrate that MEs made a significant contribution to the evolution of human genome by participating in gene function in a human-specific fashion.

SV-plaudit: A cloud-based framework for manually curating thousands of structural variants

SV-plaudit is a framework for rapidly curating structural variant (SV) predictions. For each SV, we generate an image that visualizes the coverage and alignment signals from a set of samples. Images are uploaded to our cloud framework where users assess the quality of each image using a client-side web application. Reports can then be generated as a tab-delimited file or annotated Variant Call Format (VCF) file. As a proof of principle, nine researchers collaborated for 1 hour to evaluate 1,350 SVs each. We anticipate that SV-plaudit will become a standard step in variant calling pipelines and the crowd-sourced curation of other biological results.Code available at https://github.com/jbelyeu/SV-plauditDemonstration video available at https://www.youtube.com/watch?v=ono8kHMKxDs.

Ongoing human chromosome end extension revealed by analysis of BioNano and nanopore data

The majority of human chromosome ends remain incompletely assembled due to their highly repetitive structure. In this study, we use BioNano data to anchor and extend chromosome ends from two European trios as well as two unrelated Asian genomes. At least 11 BioNano assembled chromosome ends are structurally divergent from the reference genome, including both missing sequence and extensions. These extensions are heritable and in some cases divergent between Asian and European samples. Six out of nine predicted extension sequences from NA12878 can be confirmed and filled by nanopore data. We identify two multi-kilobase sequence families both enriched more than 100-fold in extension sequence (p-values < 1e-5) whose origins can be traced to interstitial sequence on ancestral primate chromosome 7. Extensive sub-telomeric duplication of these families has occurred in the human lineage subsequent to divergence from chimpanzees.

Deciphering highly similar multigene family transcripts from Iso-Seq data with IsoCon

A significant portion of genes in vertebrate genomes belongs to multigene families, with each family containing several gene copies whose presence/absence, as well as isoform structure, can be highly variable across individuals. Existing de novo techniques for assaying the sequences of such highly-similar gene families fall short of reconstructing end-to-end transcripts with nucleotide-level precision or assigning alternatively spliced transcripts to their respective gene copies. We present IsoCon, a high-precision method using long PacBio Iso-Seq reads to tackle this challenge. We apply IsoCon to nine Y chromosome ampliconic gene families and show that it outperforms existing methods on both experimental and simulated data. IsoCon has allowed us to detect an unprecedented number of novel isoforms and has opened the door for unraveling the structure of many multigene families and gaining a deeper understanding of genome evolution and human diseases.

Effect of Collapsed Duplications on Diversity Estimates: What to Expect

The study of segmental duplications (SDs) and copy-number variants (CNVs) is of great importance in the fields of genomics and evolution. However, SDs and CNVs are usually excluded from genome-wide scans for natural selection. Because of high identity between copies, SDs and CNVs that are not included in reference genomes are prone to be collapsed-that is, mistakenly aligned to the same region-when aligning sequence data from single individuals to the reference. Such collapsed duplications are additionally challenging because concerted evolution between duplications alters their site frequency spectrum and linkage disequilibrium patterns. To investigate the potential effect of collapsed duplications upon natural selection scans we obtained expectations for four summary statistics from simulations of duplications evolving under a range of interlocus gene conversion and crossover rates. We confirm that summary statistics traditionally used to detect the action of natural selection on DNA sequences cannot be applied to SDs and CNVs since in some cases values for known duplications mimic selective signatures. As a proof of concept of the pervasiveness of collapsed duplications, we analyzed data from the 1,000 Genomes Project. We find that, within regions identified as variable in copy number, diversity between individuals with the duplication is consistently higher than between individuals without the duplication. Furthermore, the frequency of single nucleotide variants (SNVs) deviating from Hardy-Weinberg Equilibrium is higher in individuals with the duplication, which strongly suggests that higher diversity is a consequence of collapsed duplications and incorrect evaluation of SNVs within these CNV regions.

A Two-State Model of Tree Evolution and Its Applications to Alu Retrotransposition

Models of tree evolution have mostly focused on capturing the cladogenesis processes behind speciation. Processes that derive the evolution of genomic elements, such as repeats, are not necessarily captured by these existing models. In this article, we design a model of tree evolution that we call the dual-birth model, and we show how it can be useful in studying the evolution of short Alu repeats found in the human genome in abundance. The dual-birth model extends the traditional birth-only model to have two rates of propagation, one for active nodes that propagate often, and another for inactive nodes, that with a lower rate, activate and start propagating. Adjusting the ratio of the rates controls the expected tree balance. We present several theoretical results under the dual-birth model, introduce parameter estimation techniques, and study the properties of the model in simulations. We then use the dual-birth model to estimate the number of active Alu elements and their rates of propagation and activation in the human genome based on a large phylogenetic tree that we build from close to one million Alu sequences.

Unravelling the genetic architecture of autosomal recessive epilepsy in the genomic era

The technological advancement of next-generation sequencing has greatly accelerated the pace of variant discovery in epilepsy. Despite an initial focus on autosomal dominant epilepsy due to the tractable nature of variant discovery with trios under a de novo model, more and more variants are being reported in families with epilepsies consistent with autosomal recessive (AR) inheritance. In this review, we touch on the classical AR epilepsy variants such as the inborn errors of metabolism and malformations of cortical development. However, we also highlight recently reported genes that are being identified by next-generation sequencing approaches and online 'matchmaking' platforms. Syndromes mainly characterized by seizures and complex neurodevelopmental disorders comorbid with epilepsy are discussed as an example of the wide phenotypic spectrum associated with the AR epilepsies. We conclude with a foray into the future, from the application of whole-genome sequencing to identify elusive epilepsy variants, to the promise of precision medicine initiatives to provide novel targeted therapeutics specific to the individual based on their clinical genetic testing.

npInv: accurate detection and genotyping of inversions using long read sub-alignment

Background

Detection of genomic inversions remains challenging. Many existing methods primarily target inzversions with a non repetitive breakpoint, leaving inverted repeat (IR) mediated non-allelic homologous recombination (NAHR) inversions largely unexplored.

Result

We present npInv, a novel tool specifically for detecting and genotyping NAHR inversion using long read sub-alignment of long read sequencing data. We benchmark npInv with other tools in both simulation and real data. We use npInv to generate a whole-genome inversion map for NA12878 consisting of 30 NAHR inversions (of which 15 are novel), including all previously known NAHR mediated inversions in NA12878 with flanking IR less than 7kb. Our genotyping accuracy on this dataset was 94%. We used PCR to confirm the presence of two of these novel inversions. We show that there is a near linear relationship between the length of flanking IR and the minimum inversion size, without inverted repeats.

Conclusion

The application of npInv shows high accuracy in both simulation and real data. The results give deeper insight into understanding inversion.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2252-9) contains supplementary material, which is available to authorized users.

Human-Specific NOTCH2NL Genes Affect Notch Signaling and Cortical Neurogenesis

Genetic changes causing brain size expansion in human evolution have remained elusive. Notch signaling is essential for radial glia stem cell proliferation and is a determinant of neuronal number in the mammalian cortex. We find that three paralogs of human-specific NOTCH2NL are highly expressed in radial glia. Functional analysis reveals that different alleles of NOTCH2NL have varying potencies to enhance Notch signaling by interacting directly with NOTCH receptors. Consistent with a role in Notch signaling, NOTCH2NL ectopic expression delays differentiation of neuronal progenitors, while deletion accelerates differentiation into cortical neurons. Furthermore, NOTCH2NL genes provide the breakpoints in 1q21.1 distal deletion/duplication syndrome, where duplications are associated with macrocephaly and autism and deletions with microcephaly and schizophrenia. Thus, the emergence of human-specific NOTCH2NL genes may have contributed to the rapid evolution of the larger human neocortex, accompanied by loss of genomic stability at the 1q21.1 locus and resulting recurrent neurodevelopmental disorders.

Detection and analysis of ancient segmental duplications in mammalian genomes

Although segmental duplications (SDs) represent hotbeds for genomic rearrangements and emergence of new genes, there are still no easy-to-use tools for identifying SDs. Moreover, while most previous studies focused on recently emerged SDs, detection of ancient SDs remains an open problem. We developed an SDquest algorithm for SD finding and applied it to analyzing SDs in human, gorilla, and mouse genomes. Our results demonstrate that previous studies missed many SDs in these genomes and show that SDs account for at least 6.05% of the human genome (version hg19), a 17% increase as compared to the previous estimate. Moreover, SDquest classified 6.42% of the latest GRCh38 version of the human genome as SDs, a large increase as compared to previous studies. We thus propose to re-evaluate evolution of SDs based on their accurate representation across multiple genomes. Toward this goal, we analyzed the complex mosaic structure of SDs and decomposed mosaic SDs into elementary SDs, a prerequisite for follow-up evolutionary analysis. We also introduced the concept of the breakpoint graph of mosaic SDs that revealed SD hotspots and suggested that some SDs may have originated from circular extrachromosomal DNA (ecDNA), not unlike ecDNA that contributes to accelerated evolution in cancer.

Complex evolution of the GSTM gene family involves sharing of GSTM1 deletion polymorphism in humans and chimpanzees

Background

The common deletion of the glutathione S-transferase Mu 1 (GSTM1) gene in humans has been shown to be involved in xenobiotic metabolism and associated with bladder cancer. However, the evolution of this deletion has not been investigated.

Results

In this study, we conducted comparative analyses of primate genomes. We demonstrated that the GSTM gene family has evolved through multiple structural variations, involving gene duplications, losses, large inversions and gene conversions. We further showed experimentally that the GSTM1 was polymorphically deleted in both humans and also in chimpanzees, through independent deletion events. To generalize our results, we searched for genic deletions that are polymorphic in both humans and chimpanzees. Consequently, we found only two such deletions among the thousands that we have searched, one of them being the GSTM1 deletion and the other surprisingly being another metabolizing gene, the UGT2B17.

Conclusions

Overall, our results support the emerging notion that metabolizing gene families, such as the GSTM, NAT, UGT and CYP, have been evolving rapidly through gene duplication and deletion events in primates, leading to complex structural variation within and among species with unknown evolutionary consequences.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4676-z) contains supplementary material, which is available to authorized users.

Involvement of SPATA31 copy number variable genes in human lifespan

The SPATA31 (alias FAM75A) gene family belongs to the core duplicon families that are thought to have contributed significantly to hominoid evolution. It is also among the gene families with the strongest signal of positive selection in hominoids. It has acquired new protein domains in the primate lineage and a previous study has suggested that the gene family has expanded its function into UV response and DNA repair. Here we show that over-expression of SPATA31A1 in fibroblast cells leads to premature senescence due to interference with aging-related transcription pathways. We show that there are considerable copy number differences for this gene family in human populations and we ask whether this could influence mutation rates and longevity in humans. We find no evidence for an influence on germline mutation rates, but an analysis of long-lived individuals (> 96 years) shows that they carry significantly fewer SPATA31 copies in their genomes than younger individuals in a control group. We propose that the evolution of SPATA31 copy number is an example for antagonistic pleiotropy by providing a fitness benefit during the reproductive phase of life, but negatively influencing the overall life span.

CNVs are associated with genomic architecture in a songbird

Background

Understanding variation in genome structure is essential to understand phenotypic differences within populations and the evolutionary history of species. A promising form of this structural variation is copy number variation (CNV). CNVs can be generated by different recombination mechanisms, such as non-allelic homologous recombination, that rely on specific characteristics of the genome architecture. These structural variants can therefore be more abundant at particular genes ultimately leading to variation in phenotypes under selection. Detailed characterization of CNVs therefore can reveal evolutionary footprints of selection and provide insight in their contribution to phenotypic variation in wild populations.

Results

Here we use genotypic data from a long-term population of great tits (Parus major), a widely studied passerine bird in ecology and evolution, to detect CNVs and identify genomic features prevailing within these regions. We used allele intensities and frequencies from high-density SNP array data from 2,175 birds. We detected 41,029 CNVs concatenated into 8,008 distinct CNV regions (CNVRs). We successfully validated 93.75% of the CNVs tested by qPCR, which were sampled at different frequencies and sizes. A mother-daughter family structure allowed for the evaluation of the inheritance of a number of these CNVs. Thereby, only CNVs with 40 probes or more display segregation in accordance with Mendelian inheritance, suggesting a high rate of false negative calls for smaller CNVs. As CNVRs are a coarse-grained map of CNV loci, we also inferred the frequency of coincident CNV start and end breakpoints. We observed frequency-dependent enrichment of these breakpoints at homologous regions, CpG sites and AT-rich intervals. A gene ontology enrichment analyses showed that CNVs are enriched in genes underpinning neural, cardiac and ion transport pathways.

Conclusion

Great tit CNVs are present in almost half of the genes and prominent at repetitive-homologous and regulatory regions. Although overlapping genes under selection, the high number of false negatives make neutrality or association tests on CNVs detected here difficult. Therefore, CNVs should be further addressed in the light of their false negative rate and architecture to improve the comprehension of their association with phenotypes and evolutionary history.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4577-1) contains supplementary material, which is available to authorized users.

Whole-genome analysis of structural variations between Xiang pigs with larger litter sizes and those with smaller litter sizes

To gain a better knowledge of structural variations (SVs) in Xiang pig, we used next-generation sequencing to analyze the Xiang pigs with larger (XL) or smaller litter sizes (XS). Our analysis yielded 28,040 putative SVs in the Xiang pig. These SVs distributed throughout all of chromosomes. Some functional regions including exons and untranslated regions were less varied than introns and intergenic regions. We detected 4637 and 4119 specific SVs, which contained 1697 and 1582 genes in XL and XS group, respectively. These genes were mainly enriched in the well-known pathways involved in development and reproduction processes. Population validation was carried out on 50 SVs candidates using PCR method in 144 Xiang pig crowds. All of 50 SVs were confirmed by PCR method and 14 SVs were associated with the litter size of Xiang pigs. These results may be helpful for the elucidation of growth and reproduction regulation in Xiang pig.

Copy Number Variation in Fungi and Its Implications for Wine Yeast Genetic Diversity and Adaptation

In recent years, copy number (CN) variation has emerged as a new and significant source of genetic polymorphisms contributing to the phenotypic diversity of populations. CN variants are defined as genetic loci that, due to duplication and deletion, vary in their number of copies across individuals in a population. CN variants range in size from 50 base pairs to whole chromosomes, can influence gene activity, and are associated with a wide range of phenotypes in diverse organisms, including the budding yeast Saccharomyces cerevisiae. In this review, we introduce CN variation, discuss the genetic and molecular mechanisms implicated in its generation, how they can contribute to genetic and phenotypic diversity in fungal populations, and consider how CN variants may influence wine yeast adaptation in fermentation-related processes. In particular, we focus on reviewing recent work investigating the contribution of changes in CN of fermentation-related genes in yeast wine strains and offer notable illustrations of such changes, including the high levels of CN variation among the CUP genes, which confer resistance to copper, a metal with fungicidal properties, and the preferential deletion and duplication of the MAL1 and MAL3 loci, respectively, which are responsible for metabolizing maltose and sucrose. Based on the available data, we propose that CN variation is a substantial dimension of yeast genetic diversity that occurs largely independent of single nucleotide polymorphisms. As such, CN variation harbors considerable potential for understanding and manipulating yeast strains in the wine fermentation environment and beyond.

Mammalian Comparative Genomics Reveals Genetic and Epigenetic Features Associated with Genome Reshuffling in Rodentia

Understanding how mammalian genomes have been reshuffled through structural changes is fundamental to the dynamics of its composition, evolutionary relationships between species and, in the long run, speciation. In this work, we reveal the evolutionary genomic landscape in Rodentia, the most diverse and speciose mammalian order, by whole-genome comparisons of six rodent species and six representative outgroup mammalian species. The reconstruction of the evolutionary breakpoint regions across rodent phylogeny shows an increased rate of genome reshuffling that is approximately two orders of magnitude greater than in other mammalian species here considered. We identified novel lineage and clade-specific breakpoint regions within Rodentia and analyzed their gene content, recombination rates and their relationship with constitutive lamina genomic associated domains, DNase I hypersensitivity sites and chromatin modifications. We detected an accumulation of protein-coding genes in evolutionary breakpoint regions, especially genes implicated in reproduction and pheromone detection and mating. Moreover, we found an association of the evolutionary breakpoint regions with active chromatin state landscapes, most probably related to gene enrichment. Our results have two important implications for understanding the mechanisms that govern and constrain mammalian genome evolution. The first is that the presence of genes related to species-specific phenotypes in evolutionary breakpoint regions reinforces the adaptive value of genome reshuffling. Second, that chromatin conformation, an aspect that has been often overlooked in comparative genomic studies, might play a role in modeling the genomic distribution of evolutionary breakpoints.

Novel gene function and regulation in neocortex expansion

The expansion of the neocortex during human evolution is due to changes in our genome that result in increased and prolonged proliferation of neural stem and progenitor cells during neocortex development. Three principal types of such genomic changes can be distinguished, first, novel gene regulation in human, second, novel function in human of genes existing in both human and non-human species, and third, novel, human-specific genes. The latter comprise both, increases in the copy number of genes existing also in non-human species, and the emergence of genes giving rise to unique, human-specific gene products. Examples of all these types of changes in the human genome have been identified, with ARHGAP11B constituting a paradigmatic example of a unique, human-specific protein.

Living Organisms Author Their Read-Write Genomes in Evolution

Evolutionary variations generating phenotypic adaptations and novel taxa resulted from complex cellular activities altering genome content and expression: (i) Symbiogenetic cell mergers producing the mitochondrion-bearing ancestor of eukaryotes and chloroplast-bearing ancestors of photosynthetic eukaryotes; (ii) interspecific hybridizations and genome doublings generating new species and adaptive radiations of higher plants and animals; and, (iii) interspecific horizontal DNA transfer encoding virtually all of the cellular functions between organisms and their viruses in all domains of life. Consequently, assuming that evolutionary processes occur in isolated genomes of individual species has become an unrealistic abstraction. Adaptive variations also involved natural genetic engineering of mobile DNA elements to rewire regulatory networks. In the most highly evolved organisms, biological complexity scales with "non-coding" DNA content more closely than with protein-coding capacity. Coincidentally, we have learned how so-called "non-coding" RNAs that are rich in repetitive mobile DNA sequences are key regulators of complex phenotypes. Both biotic and abiotic ecological challenges serve as triggers for episodes of elevated genome change. The intersections of cell activities, biosphere interactions, horizontal DNA transfers, and non-random Read-Write genome modifications by natural genetic engineering provide a rich molecular and biological foundation for understanding how ecological disruptions can stimulate productive, often abrupt, evolutionary transformations.

Current Research on Non-Coding Ribonucleic Acid (RNA)

Non-coding ribonucleic acid (RNA) has without a doubt captured the interest of biomedical researchers. The ability to screen the entire human genome with high-throughput sequencing technology has greatly enhanced the identification, annotation and prediction of the functionality of non-coding RNAs. In this review, we discuss the current landscape of non-coding RNA research and quantitative analysis. Non-coding RNA will be categorized into two major groups by size: long non-coding RNAs and small RNAs. In long non-coding RNA, we discuss regular long non-coding RNA, pseudogenes and circular RNA. In small RNA, we discuss miRNA, transfer RNA, piwi-interacting RNA, small nucleolar RNA, small nuclear RNA, Y RNA, single recognition particle RNA, and 7SK RNA. We elaborate on the origin, detection method, and potential association with disease, putative functional mechanisms, and public resources for these non-coding RNAs. We aim to provide readers with a complete overview of non-coding RNAs and incite additional interest in non-coding RNA research.

Copy Number Variation Disorders

PURPOSE OF REVIEW

Copy number variation (CNV) disorders arise from the dosage imbalance of one or more gene(s), resulting from deletions, duplications or other genomic rearrangements that lead to the loss or gain of genetic material. Several disorders, characterized by multiple birth defects and neurodevelopmental abnormalities, have been associated with relatively large (>1 Mb) and often recurrent CNVs. CNVs have also been implicated in the etiology of neuropsychiatric disorders including autism and schizophrenia as well as other common complex diseases. Thus, CNVs have a significant impact on human health and disease.

RECENT FINDINGS

The use of increasingly higher resolution, genomewide analysis has greatly enhanced the detection of genetic variation, including CNVs. Furthermore, the availability of comprehensive genetic variation data from large cohorts of healthy controls has the potential to greatly improve the identification of disease associated genetic variants in patient samples.

SUMMARY

This review discusses the current knowledge about CNV disorders, including the mechanisms underlying their formation and phenotypic outcomes, and the advantages and limitations of current methods of detection and disease association.

Centromere Destiny in Dicentric Chromosomes: New Insights from the Evolution of Human Chromosome 2 Ancestral Centromeric Region

Dicentric chromosomes are products of genomic rearrangements that place two centromeres on the same chromosome. Due to the presence of two primary constrictions, they are inherently unstable and overcome their instability by epigenetically inactivating and/or deleting one of the two centromeres, thus resulting in functionally monocentric chromosomes that segregate normally during cell division. Our understanding to date of dicentric chromosome formation, behavior and fate has been largely inferred from observational studies in plants and humans as well as artificially produced de novo dicentrics in yeast and in human cells. We investigate the most recent product of a chromosome fusion event fixed in the human lineage, human chromosome 2, whose stability was acquired by the suppression of one centromere, resulting in a unique difference in chromosome number between humans (46 chromosomes) and our most closely related ape relatives (48 chromosomes). Using molecular cytogenetics, sequencing, and comparative sequence data, we deeply characterize the relicts of the chromosome 2q ancestral centromere and its flanking regions, gaining insight into the ancestral organization that can be easily broadened to all acrocentric chromosome centromeres. Moreover, our analyses offered the opportunity to trace the evolutionary history of rDNA and satellite III sequences among great apes, thus suggesting a new hypothesis for the preferential inactivation of some human centromeres, including IIq. Our results suggest two possible centromere inactivation models to explain the evolutionarily stabilization of human chromosome 2 over the last 5-6 million years. Our results strongly favor centromere excision through a one-step process.

Structural Variation Shapes the Landscape of Recombination in Mouse

Meiotic recombination is an essential feature of sexual reproduction that ensures faithful segregation of chromosomes and redistributes genetic variants in populations. Multiparent populations such as the Diversity Outbred (DO) mouse stock accumulate large numbers of crossover (CO) events between founder haplotypes, and thus present a unique opportunity to study the role of genetic variation in shaping the recombination landscape. We obtained high-density genotype data from [Formula: see text] DO mice, and localized 2.2 million CO events to intervals with a median size of 28 kb. The resulting sex-averaged genetic map of the DO population is highly concordant with large-scale (order 10 Mb) features of previously reported genetic maps for mouse. To examine fine-scale (order 10 kb) patterns of recombination in the DO, we overlaid putative recombination hotspots onto our CO intervals. We found that CO intervals are enriched in hotspots compared to the genomic background. However, as many as [Formula: see text] of CO intervals do not overlap any putative hotspots, suggesting that our understanding of hotspots is incomplete. We also identified coldspots encompassing 329 Mb, or [Formula: see text] of observable genome, in which there is little or no recombination. In contrast to hotspots, which are a few kilobases in size, and widely scattered throughout the genome, coldspots have a median size of 2.1 Mb and are spatially clustered. Coldspots are strongly associated with copy-number variant (CNV) regions, especially multi-allelic clusters, identified from whole-genome sequencing of 228 DO mice. Genes in these regions have reduced expression, and epigenetic features of closed chromatin in male germ cells, which suggests that CNVs may repress recombination by altering chromatin structure in meiosis. Our findings demonstrate how multiparent populations, by bridging the gap between large-scale and fine-scale genetic mapping, can reveal new features of the recombination landscape.

Segmental duplications: evolution and impact among the current Lepidoptera genomes

Background

Structural variation among genomes is now viewed to be as important as single nucleoid polymorphisms in influencing the phenotype and evolution of a species. Segmental duplication (SD) is defined as segments of DNA with homologous sequence.

Results

Here, we performed a systematic analysis of segmental duplications (SDs) among five lepidopteran reference genomes (Plutella xylostella, Danaus plexippus, Bombyx mori, Manduca sexta and Heliconius melpomene) to understand their potential impact on the evolution of these species. We find that the SDs content differed substantially among species, ranging from 1.2% of the genome in B. mori to 15.2% in H. melpomene. Most SDs formed very high identity (similarity higher than 90%) blocks but had very few large blocks. Comparative analysis showed that most of the SDs arose after the divergence of each linage and we found that P. xylostella and H. melpomene showed more duplications than other species, suggesting they might be able to tolerate extensive levels of variation in their genomes. Conserved ancestral and species specific SD events were assessed, revealing multiple examples of the gain, loss or maintenance of SDs over time. SDs content analysis showed that most of the genes embedded in SDs regions belonged to species-specific SDs (“Unique” SDs). Functional analysis of these genes suggested their potential roles in the lineage-specific evolution. SDs and flanking regions often contained transposable elements (TEs) and this association suggested some involvement in SDs formation. Further studies on comparison of gene expression level between SDs and non-SDs showed that the expression level of genes embedded in SDs was significantly lower, suggesting that structure changes in the genomes are involved in gene expression differences in species.

Conclusions

The results showed that most of the SDs were “unique SDs”, which originated after species formation. Functional analysis suggested that SDs might play different roles in different species. Our results provide a valuable resource beyond the genetic mutation to explore the genome structure for future Lepidoptera research.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-1007-y) contains supplementary material, which is available to authorized users.

Genetic Approaches to Understanding Psychiatric Disease

Human genetic studies have been the driving force in bringing to light the underlying biology of psychiatric conditions. As these studies fill in the gaps in our knowledge of the mechanisms at play, we will be better equipped to design therapies in rational and targeted ways, or repurpose existing therapies in previously unanticipated ways. This review is intended for those unfamiliar with psychiatric genetics as a field and provides a primer on different modes of genetic variation, the technologies currently used to probe them, and concepts that provide context for interpreting the gene-phenotype relationship. Like other subfields in human genetics, psychiatric genetics is moving from microarray technology to sequencing-based approaches as barriers of cost and expertise are removed, and the ramifications of this transition are discussed here. A summary is then given of recent genetic discoveries in a number of neuropsychiatric conditions, with particular emphasis on neurodevelopmental conditions. The general impact of genetics on drug development has been to underscore the extensive etiological heterogeneity in seemingly cohesive diagnostic categories. Consequently, the path forward is not in therapies hoping to reach large swaths of patients sharing a clinically defined diagnosis, but rather in targeting patients belonging to specific "biotypes" defined through a combination of objective, quantifiable data, including genotype.

Does copy number variation of APOL1 gene affect the susceptibility to focal segmental glomerulosclerosis?

Background: APOL1 risk variants (G1 and G2) are associated with increased susceptibility to focal segmental glomerulosclerosis (FSGS) in African population. However, the two risk mutations were not found in Chinese FSGS patients. In this study, we explored the association between the copy number variation (CNV) of APOL1 gene and FSGS.

Methods: APOL1 copy number variations were detected by quantitative real-time PCR with TaqMan probes and compared between 133 FSGS patients and 123 controls. The association between CNV of APOL1 gene and clinical parameters was also investigated.

Results: The distribution of APOL1 CNV did not show significant difference between FSGS patients and controls. The creatinine and proteinuria in the high copy number group (CN ≥ 3) were higher than the other two groups, but the difference was not significant (p > .05). The FSGS pathological types were different among the three groups.

Conclusion: There was no significant difference in the distribution of APOL1 gene copy variants between FSGS patients and normal controls, and there was no significant correlation between the APOL1 gene CNV and the FSGS patients’ clinical manifestations. APOL1 CNVs may be not associated with susceptibility to FSGS.

Ruminant-specific multiple duplication events of PRDM9 before speciation

Background

Understanding the genetic and evolutionary mechanisms of speciation genes in sexually reproducing organisms would provide important insights into mammalian reproduction and fitness. PRDM9, a widely known speciation gene, has recently gained attention for its important role in meiotic recombination and hybrid incompatibility. Despite the fact that PRDM9 is a key regulator of recombination and plays a dominant role in hybrid incompatibility, little is known about the underlying genetic and evolutionary mechanisms that generated multiple copies of PRDM9 in many metazoan lineages.

Results

The present study reports (1) evidence of ruminant-specific multiple gene duplication events, which likely have had occurred after the ancestral ruminant population diverged from its most recent common ancestor and before the ruminant speciation events, (2) presence of three copies of PRDM9, one copy (lineages I) in chromosome 1 (chr1) and two copies (lineages II & III) in chromosome X (chrX), thus indicating the possibility of ancient inter- and intra-chromosomal unequal crossing over and gene conversion events, (3) while lineages I and II are characterized by the presence of variable tandemly repeated C2H2 zinc finger (ZF) arrays, lineage III lost these arrays, and (4) C2H2 ZFs of lineages I and II, particularly the amino acid residues located at positions −1, 3, and 6 have evolved under strong positive selection.

Conclusions

Our results demonstrated two gene duplication events of PRDM9 in ruminants: an inter-chromosomal duplication that occurred between chr1 and chrX, and an intra-chromosomal X-linked duplication, which resulted in two additional copies of PRDM9 in ruminants. The observation of such duplication between chrX and chr1 is rare and may possibly have happened due to unequal crossing-over millions of years ago when sex chromosomes were independently derived from a pair of ancestral autosomes. Two copies (lineages I & II) are characterized by the presence of variable sized tandem-repeated C2H2 ZFs and evolved under strong positive selection and concerted evolution, supporting the notion of well-established Red Queen hypothesis. Collectively, gene duplication, concerted evolution, and positive selection are the likely driving forces for the expansion of ruminant PRDM9 sub-family.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0892-4) contains supplementary material, which is available to authorized users.

The birth of a human-specific neural gene by incomplete duplication and gene fusion

BACKGROUND

Gene innovation by duplication is a fundamental evolutionary process but is difficult to study in humans due to the large size, high sequence identity, and mosaic nature of segmental duplication blocks. The human-specific gene hydrocephalus-inducing 2, HYDIN2, was generated by a 364 kbp duplication of 79 internal exons of the large ciliary gene HYDIN from chromosome 16q22.2 to chromosome 1q21.1. Because the HYDIN2 locus lacks the ancestral promoter and seven terminal exons of the progenitor gene, we sought to characterize transcription at this locus by coupling reverse transcription polymerase chain reaction and long-read sequencing.

RESULTS

5' RACE indicates a transcription start site for HYDIN2 outside of the duplication and we observe fusion transcripts spanning both the 5' and 3' breakpoints. We observe extensive splicing diversity leading to the formation of altered open reading frames (ORFs) that appear to be under relaxed selection. We show that HYDIN2 adopted a new promoter that drives an altered pattern of expression, with highest levels in neural tissues. We estimate that the HYDIN duplication occurred ~3.2 million years ago and find that it is nearly fixed (99.9%) for diploid copy number in contemporary humans. Examination of 73 chromosome 1q21 rearrangement patients reveals that HYDIN2 is deleted or duplicated in most cases.

CONCLUSIONS

Together, these data support a model of rapid gene innovation by fusion of incomplete segmental duplications, altered tissue expression, and potential subfunctionalization or neofunctionalization of HYDIN2 early in the evolution of the Homo lineage.

Segmental duplications and evolutionary acquisition of UV damage response in the SPATA31 gene family of primates and humans

Background

Segmental duplications are an abundant source for novel gene functions and evolutionary adaptations. This mechanism of generating novelty was very active during the evolution of primates particularly in the human lineage. Here, we characterize the evolution and function of the SPATA31 gene family (former designation FAM75A), which was previously shown to be among the gene families with the strongest signal of positive selection in hominoids. The mouse homologue for this gene family is a single copy gene expressed during spermatogenesis.

Results

We show that in primates, the SPATA31 gene duplicated into SPATA31A and SPATA31C types and broadened the expression into many tissues. Each type became further segmentally duplicated in the line towards humans with the largest number of full-length copies found for SPATA31A in humans. Copy number estimates of SPATA31A based on digital PCR show an average of 7.5 with a range of 5–11 copies per diploid genome among human individuals. The primate SPATA31 genes also acquired new protein domains that suggest an involvement in UV response and DNA repair. We generated antibodies and show that the protein is re-localized from the nucleolus to the whole nucleus upon UV-irradiation suggesting a UV damage response. We used CRISPR/Cas mediated mutagenesis to knockout copies of the gene in human primary fibroblast cells. We find that cell lines with reduced functional copies as well as naturally occurring low copy number HFF cells show enhanced sensitivity towards UV-irradiation.

Conclusion

The acquisition of new SPATA31 protein functions and its broadening of expression may be related to the evolution of the diurnal life style in primates that required a higher UV tolerance. The increased segmental duplications in hominoids as well as its fast evolution suggest the acquisition of further specific functions particularly in humans.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3595-8) contains supplementary material, which is available to authorized users.

The evolution and population diversity of human-specific segmental duplications

Segmental duplications contribute to human evolution, adaptation and genomic instability but are often poorly characterized. We investigate the evolution, genetic variation and coding potential of human-specific segmental duplications (HSDs). We identify 218 HSDs based on analysis of 322 deeply sequenced archaic and contemporary hominid genomes. We sequence 550 human and nonhuman primate genomic clones to reconstruct the evolution of the largest, most complex regions with protein-coding potential (n=80 genes/33 gene families). We show that HSDs are non-randomly organized, associate preferentially with ancestral ape duplications termed “core duplicons”, and evolved primarily in an interspersed inverted orientation. In addition to Homo sapiens-specific gene expansions (e.g., TCAF1/2), we highlight ten gene families (e.g., ARHGAP11B and SRGAP2C) where copy number never returns to the ancestral state, there is evidence of mRNA splicing, and no common gene-disruptive mutations are observed in the general population. Such duplicates are candidates for the evolution of human-specific adaptive traits.

Probing the Association between Early Evolutionary Markers and Schizophrenia

Schizophrenia is suggested to be a by-product of the evolution in humans, a compromise for our language, creative thinking and cognitive abilities, and thus, essentially, a human disorder. The time of its origin during the course of human evolution remains unclear. Here we investigate several markers of early human evolution and their relationship to the genetic risk of schizophrenia. We tested the schizophrenia evolutionary hypothesis by analyzing genome-wide association studies of schizophrenia and other human phenotypes in a statistical framework suited for polygenic architectures. We analyzed evolutionary proxy measures: human accelerated regions, segmental duplications, and ohnologs, representing various time periods of human evolution for overlap with the human genomic loci associated with schizophrenia. Polygenic enrichment plots suggest a higher prevalence of schizophrenia associations in human accelerated regions, segmental duplications and ohnologs. However, the enrichment is mostly accounted for by linkage disequilibrium, especially with functional elements like introns and untranslated regions. Our results did not provide clear evidence that markers of early human evolution are more likely associated with schizophrenia. While SNPs associated with schizophrenia are enriched in HAR, Ohno and SD regions, the enrichment seems to be mediated by affiliation to known genomic enrichment categories. Taken together with previous results, these findings suggest that schizophrenia risk may have mainly developed more recently in human evolution.

Doublet-Mediated DNA Rearrangement-A Novel and Potentially Underestimated Mechanism for the Formation of Recurrent Pathogenic Deletions

Deletions and duplications of genomic DNA contribute to evolution, phenotypic diversity, and human disease. The underlying mechanisms are incompletely understood. We identified deletions of exon 10 of the SPAST gene in two unrelated families with hereditary spastic paraplegia. We excluded a founder event, but observed that the breakpoints map to identical repeat regions. These regions likely represent an intragenic "doublet," that is, an enigmatic class of local duplications. The fusion sequences for both deletions are compatible with recombination-based as well as with replication-based mechanisms. Searching the literature, we identified a partial SLC24A4 deletion that involved two copies of another doublet, and was likely formed in an analogous way. Comparing the SPAST and the SLC24A4 doublets with doublets identified previously suggested that many additional doublets have a high potential for triggering rearrangements. Considering that doublets are still being formed in the human genome, and that they likely create high local instability, we suggest that a two-step mechanism consisting of doublet generation and subsequent doublet-mediated deletion/duplication may underlie certain copy-number changes for which other mechanisms are currently assumed. Further studies are necessary to delineate the significance of the thus-far understudied doublets for the formation of copy-number variation.

On Computing Breakpoint Distances for Genomes with Duplicate Genes

A fundamental problem in comparative genomics is to compute the distance between two genomes in terms of its higher level organization (given by genes or syntenic blocks). For two genomes without duplicate genes, we can easily define (and almost always efficiently compute) a variety of distance measures, but the problem is NP-hard under most models when genomes contain duplicate genes. To tackle duplicate genes, three formulations (exemplar, maximum matching, and any matching) have been proposed, all of which aim to build a matching between homologous genes so as to minimize some distance measure. Of the many distance measures, the breakpoint distance (the number of nonconserved adjacencies) was the first one to be studied and remains of significant interest because of its simplicity and model-free property. The three breakpoint distance problems corresponding to the three formulations have been widely studied. Although we provided last year a solution for the exemplar problem that runs very fast on full genomes, computing optimal solutions for the other two problems has remained challenging. In this article, we describe very fast, exact algorithms for these two problems. Our algorithms rely on a compact integer-linear program that we further simplify by developing an algorithm to remove variables, based on new results on the structure of adjacencies and matchings. Through extensive experiments using both simulations and biological data sets, we show that our algorithms run very fast (in seconds) on mammalian genomes and scale well beyond. We also apply these algorithms (as well as the classic orthology tool MSOAR) to create orthology assignment, then compare their quality in terms of both accuracy and coverage. We find that our algorithm for the "any matching" formulation significantly outperforms other methods in terms of accuracy while achieving nearly maximum coverage.

UGT2B28 genomic variation is associated with hepatitis B e-antigen seroconversion in response to antiviral therapy

Seroconversion of hepatitis B virus (HBV) e-antigen (HBeAg) is a critical but often-missed therapeutic goal in standard antiviral treatments. An extreme-phenotype genome-wide association study was performed, comparing untreated spontaneous recoverers (with seroconversion of HBV surface antigen) versus entecavir-treated patients failing to achieve HBeAg seroconversion. A single-nucleotide-polymorphism rs2132039 on the UGT2B28 gene, alongside an adjacent copy number polymorphism (CNP605), manifested the strongest clinical associations (P = 3.4 × 10-8 and 0.001, respectively). Multivariate analysis showed that rs2132039-TT genotypes, but not CNP605 copy numbers, remained associated to spontaneous recoverers (P = 0.009). The clinical association of rs2132039 was validated successfully in an independent cohort (n = 302; P = 0.002). Longitudinal case-only analyses revealed that the rs2132039-TT genotype predicted shorter time-to-HBeAg-seroconversion in all antiviral-treated patients (n = 380, P = 0.012), as well as the peginterferon-treated subgroup (n = 123; P = 0.024, Hazard ratio [HR] = 2.104, Confidence interval [CI] = 1.105-4.007). In the entecavir-treated subgroup, the predictive effect was restricted by pretreatment alanine aminotransferase (ALT) and aspartate aminotransferase (AST) levels, with effective prediction observed in patients with ALT < 200 IU/ml and ALT/AST ratio <2 (n = 132; P = 0.013, HR = 10.538, CI = 1.420-78.196).

The Evolutionary Fates of a Large Segmental Duplication in Mouse

Gene duplication and loss are major sources of genetic polymorphism in populations, and are important forces shaping the evolution of genome content and organization. We have reconstructed the origin and history of a 127-kbp segmental duplication, R2d, in the house mouse (Mus musculus). R2d contains a single protein-coding gene, Cwc22 De novo assembly of both the ancestral (R2d1) and the derived (R2d2) copies reveals that they have been subject to nonallelic gene conversion events spanning tens of kilobases. R2d2 is also a hotspot for structural variation: its diploid copy number ranges from zero in the mouse reference genome to >80 in wild mice sampled from around the globe. Hemizygosity for high copy-number alleles of R2d2 is associated in cis with meiotic drive; suppression of meiotic crossovers; and copy-number instability, with a mutation rate in excess of 1 per 100 transmissions in some laboratory populations. Our results provide a striking example of allelic diversity generated by duplication and demonstrate the value of de novo assembly in a phylogenetic context for understanding the mutational processes affecting duplicate genes.

Achilles' heel of pluripotent stem cells: genetic, genomic and epigenetic variations during prolonged culture

Pluripotent stem cells differentiate into almost any specialized adult cell type of an organism. PSCs can be derived either from the inner cell mass of a blastocyst-giving rise to embryonic stem cells-or after reprogramming of somatic terminally differentiated cells to obtain ES-like cells, named induced pluripotent stem cells. The potential use of these cells in the clinic, for investigating in vitro early embryonic development or for screening the effects of new drugs or xenobiotics, depends on capability to maintain their genome integrity during prolonged culture and differentiation. Both human and mouse PSCs are prone to genomic and (epi)genetic instability during in vitro culture, a feature that seriously limits their real potential use. Culture-induced variations of specific chromosomes or genes, are almost all unpredictable and, as a whole, differ among independent cell lines. They may arise at different culture passages, suggesting the absence of a safe passage number maintaining genome integrity and rendering the control of genomic stability mandatory since the very early culture passages. The present review highlights the urgency for further studies on the mechanisms involved in determining (epi)genetic and chromosome instability, exploiting the knowledge acquired earlier on other cell types.

Monkey-based research on human disease: the implications of genetic differences

Assertions that the use of monkeys to investigate human diseases is valid scientifically are frequently based on a reported 90-93% genetic similarity between the species. Critical analyses of the relevance of monkey studies to human biology, however, indicate that this genetic similarity does not result in sufficient physiological similarity for monkeys to constitute good models for research, and that monkey data do not translate well to progress in clinical practice for humans. Salient examples include the failure of new drugs in clinical trials, the highly different infectivity and pathology of SIV/HIV, and poor extrapolation of research on Alzheimer's disease, Parkinson's disease and stroke. The major molecular differences underlying these inter-species phenotypic disparities have been revealed by comparative genomics and molecular biology - there are key differences in all aspects of gene expression and protein function, from chromosome and chromatin structure to post-translational modification. The collective effects of these differences are striking, extensive and widespread, and they show that the superficial similarity between human and monkey genetic sequences is of little benefit for biomedical research. The extrapolation of biomedical data from monkeys to humans is therefore highly unreliable, and the use of monkeys must be considered of questionable value, particularly given the breadth and potential of alternative methods of enquiry that are currently available to scientists.

Marsupials as models for understanding the role of chromosome rearrangements in evolution and disease

Chromosome rearrangements have been implicated in diseases, such as cancer, and speciation, but it remains unclear whether rearrangements are causal or merely a consequence of these processes. Two marsupial families with very different rates of karyotype evolution provide excellent models in which to study the role of chromosome rearrangements in a disease and evolutionary context. The speciose family Dasyuridae displays remarkable karyotypic conservation, with all species examined to date possessing nearly identical karyotypes. Despite the seemingly high degree of chromosome stability within this family, they appear prone to developing tumours, including transmissible devil facial tumours. In contrast, chromosome rearrangements have been frequent in the evolution of the species-rich family Macropodidae, which displays a high level of karyotypic diversity. In particular, the genus Petrogale (rock-wallabies) displays an extraordinary level of chromosome rearrangement among species. For six parapatric Petrogale species, it appears that speciation has essentially been caught in the act, providing an opportunity to determine whether chromosomal rearrangements are a cause or consequence of speciation in this system. This review highlights the reasons that these two marsupial families are excellent models for testing hypotheses for hotspots of chromosome rearrangement and deciphering the role of chromosome rearrangements in disease and speciation.