A genome-wide survey of structural variation between human and chimpanzee

Chromosomal Instability in Genome Evolution: From Cancer to Macroevolution

The integrity of the genome is crucial for the survival of all living organisms. However, genomes need to adapt to survive certain pressures, and for this purpose use several mechanisms to diversify. Chromosomal instability (CIN) is one of the main mechanisms leading to the creation of genomic heterogeneity by altering the number of chromosomes and changing their structures. In this review, we will discuss the different chromosomal patterns and changes observed in speciation, in evolutional biology as well as during tumor progression. By nature, the human genome shows an induction of diversity during gametogenesis but as well during tumorigenesis that can conclude in drastic changes such as the whole genome doubling to more discrete changes as the complex chromosomal rearrangement chromothripsis. More importantly, changes observed during speciation are strikingly similar to the genomic evolution observed during tumor progression and resistance to therapy. The different origins of CIN will be treated as the importance of double-strand breaks (DSBs) or the consequences of micronuclei. We will also explain the mechanisms behind the controlled DSBs, and recombination of homologous chromosomes observed during meiosis, to explain how errors lead to similar patterns observed during tumorigenesis. Then, we will also list several diseases associated with CIN, resulting in fertility issues, miscarriage, rare genetic diseases, and cancer. Understanding better chromosomal instability as a whole is primordial for the understanding of mechanisms leading to tumor progression.

Galactomyces Ferment Filtrate Potentiates an Anti-Inflammaging System in Keratinocytes

Skincare products play a crucial role in preventing the dry skin induced by various causes. Certain ingredients can help to improve the efficacy of skincare products. Galactomyces ferment filtrate (GFF) is such a functional ingredient. Its use originated from the empirical observation that the hands of sake brewers who deal with yeast fermentation retain a beautiful and youthful appearance. Consequently, skincare products based on GFF are widely used throughout the world. Recent studies have demonstrated that GFF activates an aryl hydrocarbon receptor (AHR) and upregulates the expression of filaggrin, a pivotal endogenous source of natural moisturizing factors, in epidermal keratinocytes. It also activates nuclear factor erythroid-2-related factor 2 (NRF2), the antioxidative master transcription factor, and exhibits potent antioxidative activity against oxidative stress induced by ultraviolet irradiation and proinflammatory cytokines, which also accelerate inflammaging. GFF-mediated NRF2 activation downregulates the expression of CDKN2A, which is known to be overexpressed in senescent keratinocytes. Moreover, GFF enhances epidermal terminal differentiation by upregulating the expression of caspase-14, claudin-1, and claudin-4. It also promotes the synthesis of the antiinflammatory cytokine IL-37 and downregulates the expression of proallergic cytokine IL-33 in keratinocytes. In addition, GFF downregulates the expression of the CXCL14 and IL6R genes, which are involved in inflammaging. These beneficial properties might underpin the potent barrier-protecting and anti-inflammaging effects of GFF-containing skin formulae.

A survey of current methods to detect and genotype inversions

Polymorphic inversions are ubiquitous in humans, and they have been linked to both adaptation and disease. Following their discovery in Drosophila more than a century ago, inversions have proved to be more elusive than other structural variants. A wide variety of methods for the detection and genotyping of inversions have recently been developed: multiple techniques based on selective amplification by PCR, short- and long-read sequencing approaches, principal component analysis of small variant haplotypes, template strand sequencing, optical mapping, and various genome assembly methods. Many methods apply complex wet lab protocols or increasingly refined bioinformatic analyses. This review is an attempt to provide a practical summary and comparison of the methods that are in current use, with a focus on metrics such as the maximum size of segmental duplications at inversion breakpoints that each method can tolerate, the size range of inversions that they recover, their throughput, and whether the locations of putative inversions must be known beforehand. This article is protected by copyright. All rights reserved.

The emerging role of Interleukin 37 in bone homeostasis and inflammatory bone diseases

Interleukin 37 (IL-37) is a newly identified cytokine that belongs to the IL-1 family. Unlike other members of the IL-1 family, it has been demonstrated that IL-37 possesses anti-inflammatory characteristics in both innate and acquired immune responses. Recently, significant progress has been made in understanding the role of IL-37 in inflammatory signaling pathways. Meanwhile, IL-37 has also attracted more and more attention in bone homeostasis and inflammatory bone diseases. The latest studies have revealed that IL-37 palys an essential role in the regulation of osteoclastogenesis and osteoblastogenesis. The levels of IL-37 are abnormal in patients with inflammatory bone diseases such as rheumatoid arthritis (RA), osteoarthritis (OA), ankylosing spondylitis (AS), and periodontitis. In addition, in vivo studies have further confirmed that recombinant IL-37 treatment displayed therapeutic potential in these diseases. The present review article aims to provide an overview describing the biological functions of IL-37 in bone homeostasis and inflammatory bone diseases, thus shedding new light on a novel therapeutic strategy in the future.

The Cumulative Indel Model: Fast and Accurate Statistical Evolutionary Alignment

Sequence alignment is essential for phylogenetic and molecular evolution inference, as well as in many other areas of bioinformatics and evolutionary biology. Inaccurate alignments can lead to severe biases in most downstream statistical analyses. Statistical alignment based on probabilistic models of sequence evolution addresses these issues by replacing heuristic score functions with evolutionary model-based probabilities. However, score-based aligners and fixed-alignment phylogenetic approaches are still more prevalent than methods based on evolutionary indel models, mostly due to computational convenience. Here, I present new techniques for improving the accuracy and speed of statistical evolutionary alignment. The "cumulative indel model" approximates realistic evolutionary indel dynamics using differential equations. "Adaptive banding" reduces the computational demand of most alignment algorithms without requiring prior knowledge of divergence levels or pseudo-optimal alignments. Using simulations, I show that these methods lead to fast and accurate pairwise alignment inference. Also, I show that it is possible, with these methods, to align and infer evolutionary parameters from a single long synteny block ($\approx$530 kbp) between the human and chimp genomes. The cumulative indel model and adaptive banding can therefore improve the performance of alignment and phylogenetic methods. [Evolutionary alignment; pairHMM; sequence evolution; statistical alignment; statistical genetics.].

Single-cell strand sequencing of a macaque genome reveals multiple nested inversions and breakpoint reuse during primate evolution

Rhesus macaque is an Old World monkey that shared a common ancestor with human ∼25 Myr ago and is an important animal model for human disease studies. A deep understanding of its genetics is therefore required for both biomedical and evolutionary studies. Among structural variants, inversions represent a driving force in speciation and play an important role in disease predisposition. Here we generated a genome-wide map of inversions between human and macaque, combining single-cell strand sequencing with cytogenetics. We identified 375 total inversions between 859 bp and 92 Mbp, increasing by eightfold the number of previously reported inversions. Among these, 19 inversions flanked by segmental duplications overlap with recurrent copy number variants associated with neurocognitive disorders. Evolutionary analyses show that in 17 out of 19 cases, the Hominidae orientation of these disease-associated regions is always derived. This suggests that duplicated sequences likely played a fundamental role in generating inversions in humans and great apes, creating architectures that nowadays predispose these regions to disease-associated genetic instability. Finally, we identified 861 genes mapping at 156 inversions breakpoints, with some showing evidence of differential expression in human and macaque cell lines, thus highlighting candidates that might have contributed to the evolution of species-specific features. This study depicts the most accurate fine-scale map of inversions between human and macaque using a two-pronged integrative approach, such as single-cell strand sequencing and cytogenetics, and represents a valuable resource toward understanding of the biology and evolution of primate species.

Recurrent inversion toggling and great ape genome evolution

Inversions play an important role in disease and evolution but are difficult to characterize because their breakpoints map to large repeats. We increased by sixfold the number (n = 1,069) of previously reported great ape inversions by using single-cell DNA template strand and long-read sequencing. We find that the X chromosome is most enriched (2.5-fold) for inversions, on the basis of its size and duplication content. There is an excess of differentially expressed primate genes near the breakpoints of large (>100 kilobases (kb)) inversions but not smaller events. We show that when great ape lineage-specific duplications emerge, they preferentially (approximately 75%) occur in an inverted orientation compared to that at their ancestral locus. We construct megabase-pair scale haplotypes for individual chromosomes and identify 23 genomic regions that have recurrently toggled between a direct and an inverted state over 15 million years. The direct orientation is most frequently the derived state for human polymorphisms that predispose to recurrent copy number variants associated with neurodevelopmental disease.

A new target for the treatment of inflammatory bowel disease: Interleukin-37

Interleukin (IL)-37 belongs to the IL-1 cytokine family. It has anti-inflammatory effects on numerous autoimmune diseases such as asthma, psoriasis, inflammatory bowel disease (IBD), systemic lupus erythematosus (SLE), multiple sclerosis (MS) and rheumatoid arthritis (RA). Mechanistically, IL-37 plays an anti-inflammatory role by regulating the expression of inflammatory factors in two ways: binding extracellular receptors IL-18R or transferring into the nucleus with Smad3. IBD is a kind of idiopathic intestinal inflammatory disease with unknown etiology and pathogenesis. Recent researches had proved that IL-37 is negatively involved in the pathogenesis and development of IBD. Among various inflammatory diseases, IL-37 has been shown to regulate inflammatory development by acting on various immune cells such as neutrophils, macrophages (Mϕ), dendritic cells (DCs), T cells and intestinal epithelial cells. This review summarizes the biological role of IL-37, and its immunoregulatory effects on the immune cells, especially anti-inflammatory function in both human and experimental models of IBD.

Identification of Structural Variation in Chimpanzees Using Optical Mapping and Nanopore Sequencing

Recent efforts to comprehensively characterize great ape genetic diversity using short-read sequencing and single-nucleotide variants have led to important discoveries related to selection within species, demographic history, and lineage-specific traits. Structural variants (SVs), including deletions and inversions, comprise a larger proportion of genetic differences between and within species, making them an important yet understudied source of trait divergence. Here, we used a combination of long-read and -range sequencing approaches to characterize the structural variant landscape of two additional Pan troglodytes verus individuals, one of whom carries 13% admixture from Pan troglodytes troglodytes. We performed optical mapping of both individuals followed by nanopore sequencing of one individual. Filtering for larger variants (>10 kbp) and combined with genotyping of SVs using short-read data from the Great Ape Genome Project, we identified 425 deletions and 59 inversions, of which 88 and 36, respectively, were novel. Compared with gene expression in humans, we found a significant enrichment of chimpanzee genes with differential expression in lymphoblastoid cell lines and induced pluripotent stem cells, both within deletions and near inversion breakpoints. We examined chromatin-conformation maps from human and chimpanzee using these same cell types and observed alterations in genomic interactions at SV breakpoints. Finally, we focused on 56 genes impacted by SVs in >90% of chimpanzees and absent in humans and gorillas, which may contribute to chimpanzee-specific features. Sequencing a greater set of individuals from diverse subspecies will be critical to establish the complete landscape of genetic variation in chimpanzees.

New genome assembly of the barn owl (Tyto alba alba)

New genomic tools open doors to study ecology, evolution, and population genomics of wild animals. For the Barn owl species complex, a cosmopolitan nocturnal raptor, a very fragmented draft genome was assembled for the American species (Tyto furcata pratincola) (Jarvis et al. 2014). To improve the genome, we assembled de novo Illumina and Pacific Biosciences (PacBio) long reads sequences of its European counterpart (Tyto alba alba). This genome assembly of 1.219 Gbp comprises 21,509 scaffolds and results in a N50 of 4,615,526 bp. BUSCO (Universal Single-Copy Orthologs) analysis revealed an assembly completeness of 94.8% with only 1.8% of the genes missing out of 4,915 avian orthologs searched, a proportion similar to that found in the genomes of the zebra finch (Taeniopygia guttata) or the collared flycatcher (Ficedula albicollis). By mapping the reads of the female American barn owl to the male European barn owl reads, we detected several structural variants and identified 70 Mbp of the Z chromosome. The barn owl scaffolds were further mapped to the chromosomes of the zebra finch. In addition, the completeness of the European barn owl genome is demonstrated with 94 of 128 proteins missing in the chicken genome retrieved in the European barn owl transcripts. This improved genome will help future barn owl population genomic investigations.

Chromoanagenesis: a piece of the macroevolution scenario

Over the last decade, new types of massive and complex chromosomal rearrangements based on the chaotic shattering and restructuring of chromosomes have been identified in cancer cells as well as in patients with congenital diseases and healthy individuals. These unanticipated phenomena are named chromothripsis, chromoanasynthesis and chromoplexy, and are grouped under the term of chromoanagenesis. As mechanisms for rapid and profound genome modifications in germlines and early development, these processes can be regarded as credible pathways for genomic evolution and speciation process. Their discovery confirms the importance of genome-centric investigations to fully understand organismal evolution.

Because they oppose the model of progressive acquisition of driver mutations or rearrangements, these phenomena conceptually give support to the concept of macroevolution, known through the models of “Hopeful Monsters” and the “Punctuated Equilibrium”. In this review, we summarize mechanisms underlying chromoanagenesis processes and we show that numerous cases of chromosomal speciation and short-term adaptation could be correlated to chromoanagenesis-related mechanisms.

In the frame of a modern and integrative analysis of eukaryote evolutionary processes, it seems important to consider the unexpected chromoanagenesis phenomena.

Insight into interleukin-37: The potential therapeutic target in allergic diseases

Allergic diseases are ubiquitous diseases with detrimental effects on the quality of life of people worldwide. Common allergic diseases include asthma, allergic rhinitis (AR) and allergic dermatitis (AD). Recently, studies have shown that interleukin (IL)-37, a novel cytokine in the IL-1 family, exhibits broad protective properties in various diseases, such as autoimmune diseases and cancer. IL-37 displays its anti-inflammatory effect on diseases by curbing innate and acquired immunity as well as inflammatory reactions. IL-37 functions by forming a complex with IL-18Rα and IL-1R8 extracellularly and can be translocated to the nucleus upon forming a complex with mothers against decapentaplegic homolog 3 (Smad3) intracellularly, thereby affecting gene transcription and signaling pathway activation. In addition, increasing evidence confirms that IL-37 expression is aberrant in asthma, AR and AD, which indicates that IL-37 may also play essential roles in allergic diseases. Furthermore, accumulating data obtained from recombinant IL-37 (rIL-37)-treated mice and from IL-37 transgenic (IL-37tg) mice suggest a protective role for IL-37. This review will detail the role of IL-37 in the occurrence and development of allergic diseases and discuss the potential of IL-37 as a therapeutic target in allergic diseases.

Evolutionary Changes in Transcriptional Regulation: Insights into Human Behavior and Neurological Conditions

Understanding the biological basis for human-specific cognitive traits presents both immense challenges and unique opportunities. Although the question of what makes us human has been investigated with several different methods, the rise of comparative genomics, epigenomics, and medical genetics has provided tools to help narrow down and functionally assess the regions of the genome that seem evolutionarily relevant along the human lineage. In this review, we focus on how medical genetic cases have provided compelling functional evidence for genes and loci that appear to have interesting evolutionary signatures in humans. Furthermore, we examine a special class of noncoding regions, human accelerated regions (HARs), that have been suggested to show human-lineage-specific divergence, and how the use of clinical and population data has started to provide functional information to examine these regions. Finally, we outline methods that provide new insights into functional noncoding sequences in evolution.

Cynomolgus macaque IL37 polymorphism and control of SIV infection

The association between gene polymorphisms and plasma virus load at the set point (SP-PVL) was investigated in Mauritian macaques inoculated with SIV. Among 44 macaques inoculated with 50 AID50, six individuals were selected: three with SP-PVL among the highest and three with SP-PVL among the lowest. The exons of 390 candidate genes of these six animals were sequenced. Twelve non-synonymous single nucleotide polymorphisms (NS-SNPs) lying in nine genes potentially associated with PVL were genotyped in 23 animals. Three NS-SNPs with probabilities of association with PVL less than 0.05 were genotyped in a total of 44 animals. One NS-SNP lying in exon 1 of the IL37 gene displayed a significant association (p = 3.33 × 10⁻⁴) and a strong odds ratio (19.52). Multiple linear regression modeling revealed three significant predictors of SP-PVL, including the IL37 exon 1 NS-SNP (p = 0.0004) and the MHC Class IB haplotypes M2 (p = 0.0007) and M6 (p = 0.0013). These three factors in conjunction explained 48% of the PVL variance (p = 4.8 × 10⁻⁶). The potential role of IL37 in the control of SIV infection is discussed.

InDel markers: An extended marker resource for molecular breeding in chickpea

Chickpea is one of the most important food legumes that holds the key to meet rising global food and nutritional demand. In order to deploy molecular breeding approaches in crop improvement programs, user friendly and cost effective marker resources remain prerequisite. The advent of next generation sequencing (NGS) technology has resulted in the generation of several thousands of markers as part of several large scale genome sequencing and re-sequencing initiatives. Very recently, PCR based Insertion-deletions (InDels) are becoming a popular gel based genotyping solution because of their co-dominant, inexpensive, and highly polymorphic nature. With an objective to expand marker resources for genomics assisted breeding (GAB) in chickpea, whole genome re-sequencing data generated on five parental lines of one interspecific (ICC 4958 × PI 489777) and two intra-specific (ICC 283 × ICC 8261 and ICC 4958 × ICC 1882) mapping populations, were used for identification of InDels. A total of 231,658 InDels were identified using Dindel software with default parameters. Further, a total of 8,307 InDels with ≥20 bp size were selected for development of gel based markers, of which primers could be designed for 7,523 (90.56%) markers. On average, markers appeared at a frequency of 1,038 InDels/LG with a maximum number of markers on CaLG04 (1,952 InDels) and minimum on CaLG08 (360 InDels). In order to validate these InDels, a total of 423 primer pairs were randomly selected and tested on the selected parental lines. A high amplification rate of 80% was observed ranging from 46.06 to 58.01% polymorphism rate across parents on 3% agarose gel. This study clearly reflects the usefulness of available sequence data for the development of genome-wide InDels in chickpea that can further contribute and accelerate a wide range of genetic and molecular breeding activities in chickpea.

The molecular characterization of fixed inversions breakpoints unveils the ancestral character of the Drosophila guanche chromosomal arrangements

Cytological studies revealed that the number of chromosomes and their organization varies across species. The increasing availability of whole genome sequences of multiple species across specific phylogenies has confirmed and greatly extended these cytological observations. In the Drosophila genus, the ancestral karyotype consists of five rod-like acrocentric chromosomes (Muller elements A to E) and one dot-like chromosome (element F), each exhibiting a generally conserved gene content. Chromosomal fusions and paracentric inversions are thus the major contributors, respectively, to chromosome number variation among species and to gene order variation within chromosomal element. The subobscura cluster of Drosophila consists in three species that retain the genus ancestral karyotype and differ by a reduced number of fixed inversions. Here, we have used cytological information and the D. guanche genome sequence to identify and molecularly characterize the breakpoints of inversions that became fixed since the D. guanche-D. subobscura split. Our results have led us to propose a modified version of the D. guanche cytological map of its X chromosome, and to establish that (i) most inversions became fixed in the D. subobscura lineage and (ii) the order in which the four X chromosome overlapping inversions occurred and became fixed.

Structural variation during dog domestication: insights from gray wolf and dhole genomes

Several processes like phenotypic evolution, disease susceptibility and environmental adaptations, which fashion the domestication of animals, are largely attributable to structural variations (SVs) in the genome. Here, we present high-quality draft genomes of the gray wolf (Canis lupus) and dhole (Cuon alpinus) with scaffold N50 of 6.04 Mb and 3.96 Mb, respectively. Sequence alignment comprising genomes of three canid species reveals SVs specific to the dog, particularly 16 315 insertions, 2565 deletions, 443 repeats, 16 inversions and 15 translocations. Functional annotation of the dog SVs associated with genes indicates their enrichments in energy metabolisms, neurological processes and immune systems. Interestingly, we identify and verify at population level an insertion fully covering a copy of the AKR1B1 (Aldo-Keto Reductase Family 1 Member B) transcript. Transcriptome analysis reveals a high level of expression of the new AKR1B1 copy in the small intestine and liver, implying an increase in de novo fatty acid synthesis and antioxidant ability in dog compared to gray wolf, likely in response to dietary shifts during the agricultural revolution. For the first time, we report a comprehensive analysis of the evolutionary dynamics of SVs during the domestication step of dogs. Our findings demonstrate that retroposition can birth new genes to facilitate domestication, and affirm the importance of large-scale genomic variants in domestication studies.

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.

Suppression of inflammation and acquired immunity by IL-37

IL-37 is a unique member of the IL-1 family of cytokines, which functions as a natural suppressor of inflammatory and immune responses. Immune and non-immune cells produce IL-37 precursor following pro-inflammatory stimuli. Following activating cleavage by caspase-1, mature IL-37 translocates to the nucleus, where it suppresses transcription of pro-inflammatory genes. Both precursor and mature IL-37 are also secreted in the extracellular space, where they bind IL-18Rα and recruit the IL-1R8 (formerly TIR8 or SIGIRR), which transduces anti-inflammatory signals by suppressing NF-kB and MAPK and by activating Mer-PTEN-DOK pathways. During inflammation, IL-37 restores the metabolism of the cell by reducing succinate, inhibiting mTOR, and activating AMPK. Transgenic mice expressing human IL-37 and wild type mice treated with recombinant human IL-37 are protected from several experimental models of inflammation, including endotoxin shock, colitis, lung and spinal cord injury, coronary artery disease, arthritis and inflammation-induced fatigue, while also exhibiting reduced adaptive immune responses. In humans, IL-37 likely functions to limit excessive inflammation: accordingly, IL-37 levels are abnormal in patients with inflammatory and autoimmune diseases. In this review, we provide an overview of the discovery and biology of IL-37, and discuss the potential for development of this cytokine as a therapeutic agent.

High-resolution comparative analysis of great ape genomes

Genetic studies of human evolution require high-quality contiguous ape genome assemblies that are not guided by the human reference. We coupled long-read sequence assembly and full-length complementary DNA sequencing with a multiplatform scaffolding approach to produce ab initio chimpanzee and orangutan genome assemblies. By comparing these with two long-read de novo human genome assemblies and a gorilla genome assembly, we characterized lineage-specific and shared great ape genetic variation ranging from single- to mega-base pair-sized variants. We identified ~17,000 fixed human-specific structural variants identifying genic and putative regulatory changes that have emerged in humans since divergence from nonhuman apes. Interestingly, these variants are enriched near genes that are down-regulated in human compared to chimpanzee cerebral organoids, particularly in cells analogous to radial glial neural progenitors.

Inversion variants in human and primate genomes

For many years, inversions have been proposed to be a direct driving force in speciation since they suppress recombination when heterozygous. Inversions are the most common large-scale differences among humans and great apes. Nevertheless, they represent large events easily distinguishable by classical cytogenetics, whose resolution, however, is limited. Here, we performed a genome-wide comparison between human, great ape, and macaque genomes using the net alignments for the most recent releases of genome assemblies. We identified a total of 156 putative inversions, between 103 kb and 91 Mb, corresponding to 136 human loci. Combining literature, sequence, and experimental analyses, we analyzed 109 of these loci and found 67 regions inverted in one or multiple primates, including 28 newly identified inversions. These events overlap with 81 human genes at their breakpoints, and seven correspond to sites of recurrent rearrangements associated with human disease. This work doubles the number of validated primate inversions larger than 100 kb, beyond what was previously documented. We identified 74 sites of errors, where the sequence has been assembled in the wrong orientation, in the reference genomes analyzed. Our data serve two purposes: First, we generated a map of evolutionary inversions in these genomes representing a resource for interrogating differences among these species at a functional level; second, we provide a list of misassembled regions in these primate genomes, involving over 300 Mb of DNA and 1978 human genes. Accurately annotating these regions in the genome references has immediate applications for evolutionary and biomedical studies on primates.

A comparative transcriptional landscape of maize and sorghum obtained by single-molecule sequencing

Maize and sorghum are both important crops with similar overall plant architectures, but they have key differences, especially in regard to their inflorescences. To better understand these two organisms at the molecular level, we compared expression profiles of both protein-coding and noncoding transcripts in 11 matched tissues using single-molecule, long-read, deep RNA sequencing. This comparative analysis revealed large numbers of novel isoforms in both species. Evolutionarily young genes were likely to be generated in reproductive tissues and usually had fewer isoforms than old genes. We also observed similarities and differences in alternative splicing patterns and activities, both among tissues and between species. The maize subgenomes exhibited no bias in isoform generation; however, genes in the B genome were more highly expressed in pollen tissue, whereas genes in the A genome were more highly expressed in endosperm. We also identified a number of splicing events conserved between maize and sorghum. In addition, we generated comprehensive and high-resolution maps of poly(A) sites, revealing similarities and differences in mRNA cleavage between the two species. Overall, our results reveal considerable splicing and expression diversity between sorghum and maize, well beyond what was reported in previous studies, likely reflecting the differences in architecture between these two species.

Chromothripsis and the Macroevolution Theory

The recent discovery of a new class of massive chromosomal rearrangements, occurring during one unique cellular event and baptized "chromothripsis," deeply modifies our perception on the genesis of complex genomic rearrangements, but also it raises the question of the potential driving role of chromothripsis in species evolution. The occurrence of chromothripsis appears to be in good agreement with macroevolution models proposed as a complement to phyletic gradualism. The emergence of this unexpected phenomenon may help to demonstrate the contribution of chromosome rearrangements to speciation process.

The potential of interleukin-37 as an effective therapeutic agent in asthma

Interleukin (IL)-37 belongs to the IL-1 cytokine family. It binds to IL-18Rα and recruits the orphan decoy IL-1R8. Emerging evidence shows that IL-37 is a key player in the regulation of inflammation, cellular differentiation, and proliferation. Altered IL-37 expression has been demonstrated in many inflammatory disease conditions, including asthma. In rheumatoid arthritis, IL-37 is involved in the regulation of proliferation, production of inflammatory mediators, and activation of inflammatory cells. Furthermore, this cytokine acts as a negative regulator of inflammation in inflammatory bowel disease. Similarly, IL-37 also appears to suppress allergic inflammation in asthma. In a murine model of asthma, local administration of IL-37 markedly reduced the degree of inflammatory cell infiltration and airway hyper-responsiveness. IL-37 has also been shown to be involved in a number of aspects of allergic inflammation, such as eosinophil and neutrophil recruitment, as well as inhibition of Th1/Th2/Th17 inflammatory mediators. However, the exact molecular mechanisms underlying the function of IL-37 in human asthma have yet to be fully elucidated. This review describes the current evidence regarding the role of IL-37 in the pathophysiology of asthma and evaluates both the potential of IL-37 as a biomarker for airway inflammation and a therapeutic target for the treatment of asthma.

Suppression of innate inflammation and immunity by interleukin-37

IL-37 is unique in the IL-1 family in that unlike other members of the family, IL-37 broadly suppresses innate immunity. IL-37 can be elevated in humans with inflammatory and autoimmune diseases where it likely functions to limit inflammation. Transgenic mice expressing human IL-37 (IL37-tg) exhibit less severe inflammation in models of endotoxin shock, colitis, myocardial infarction, lung, and spinal cord injury. IL37-tg mice have reduced antigen-specific responses and dendritic cells (DCs) from these mice exhibit characteristics of tolerogenic DCs. Compared to aging wild-type (WT) mice, aging IL37-tg mice are protected against B-cell leukemogenesis and heart failure. Treatment of WT mice with recombinant human IL-37 has been shown to be protective in several models of inflammation and injury. IL-37 binds to the IL-18 receptor but then recruits the orphan IL-1R8 (formerly TIR8 or SIGIRR) in order to function as an inhibitor. Here, we review the discovery of IL-37, its production, release, and mechanisms by which IL-37 reduces inflammation and suppresses immune responses. The data reviewed here suggest a therapeutic potential for IL-37.

Single-Cell-Based Platform for Copy Number Variation Profiling through Digital Counting of Amplified Genomic DNA Fragments

We develop a novel single-cell-based platform through digital counting of amplified genomic DNA fragments, named multifraction amplification (mfA), to detect the copy number variations (CNVs) in a single cell. Amplification is required to acquire genomic information from a single cell, while introducing unavoidable bias. Unlike prevalent methods that directly infer CNV profiles from the pattern of sequencing depth, our mfA platform denatures and separates the DNA molecules from a single cell into multiple fractions of a reaction mix before amplification. By examining the sequencing result of each fraction for a specific fragment and applying a segment-merge maximum likelihood algorithm to the calculation of copy number, we digitize the sequencing-depth-based CNV identification and thus provide a method that is less sensitive to the amplification bias. In this paper, we demonstrate a mfA platform through multiple displacement amplification (MDA) chemistry. When performing the mfA platform, the noise of MDA is reduced; therefore, the resolution of single-cell CNV identification can be improved to 100 kb. We can also determine the genomic region free of allelic drop-out with mfA platform, which is impossible for conventional single-cell amplification methods.

Interleukin 37 expression in mice alters sleep responses to inflammatory agents and influenza virus infection

Multiple interactions between the immune system and sleep are known, including the effects of microbial challenge on sleep or the effects of sleep loss on facets of the immune response. Cytokines regulate, in part, sleep and immune responses. Here we examine the role of an anti-inflammatory cytokine, interleukin-37 (IL-37) on sleep in a mouse strain that expresses human IL-37b (IL37tg mice). Constitutive expression of the IL-37 gene in the brains of these mice under resting conditions is low; however, upon an inflammatory stimulus, expression increases dramatically. We measured sleep in three conditions; (a) under baseline conditions and after 6 h of sleep loss, (b) after bolus intraperitoneal administration of lipopolysaccharide (LPS) or IL-1β and (c) after intranasal influenza virus challenge. Under baseline conditions, the IL37tg mice had 7% more spontaneous non-rapid eye movement sleep (NREMS) during the light period than wild-type (WT) mice. After sleep deprivation both WT mice and IL37tg mice slept an extra 21% and 12%, respectively, during the first 6 h of recovery. NREMS responses after sleep deprivation did not significantly differ between WT mice and IL37tg mice. However, in response to either IL-1β or LPS, the increases in time spent in NREMS were about four-fold greater in the WT mice than in the IL37tg mice. In contrast, in response to a low dose of mouse-adapted H1N1 influenza virus, sleep responses developed slowly over the 6 day recording period. By day 6, NREMS increased by 10% and REMS increased by 18% in the IL37tg mice compared to the WT mice. Further, by day 4 IL37tg mice lost less weight, remained more active, and retained their body temperatures closer to baseline values than WT mice. We conclude that conditions that promote IL-37 expression attenuate morbidity to severe inflammatory challenge.

•

Sleep responses to mild acute sleep deprivation are similar in mice transgenic for interleukin-37 (IL37tg) IL37tg and wild type (WT) mice.

•

Sleep responses induced by either IL-β or LPS are greatly attenuated in IL37tg mice compared to WT mice.

•

After influenza virus challenge, IL37tg mice have reduced morbidities and enhanced sleep responses.

Common genetic heterogeneity of human interleukin-37 leads to functional variance

Interleukin-37 (IL-37) is an inhibitory member of the IL-1 family of cytokines. We previously found that balanced selection maintains common variations of the human IL37 gene. However, the functional consequences of this selection have yet to be validated. Here, using cells expressing exogenous IL-37 variants, including IL-37 Ref and IL-37 Var1 and Var2, we found that the three variants of IL-37 exhibited different immunoregulatory potencies in response to immune stimulation. The protein level of IL-37 Var2 was found to be significantly less than that of IL-37 Ref or Var1, despite the comparable mRNA levels of all three variants. Further study showed that IL-37 Var2 was rapidly degraded by a proteasome-dependent mechanism mediated by enhanced polyubiquitination, leading to a transient upregulation of IL-37 Var2 after immune stimulation. Finally, when ectopically expressed in cells, human IL-37 Var2 exerted less inhibition on proinflammatory cytokine production than did other IL-37 variants. Conversely, purified extracellular IL-37 variant proteins demonstrated comparable inhibitory abilities in vitro. In conclusion, our study reveals that common genetic variants of IL37 lead to different immune-inhibitory potencies, primarily as a result of differences in IL-37 protein stability, suggesting the possible involvement of these variants in various human diseases.

African Indigenous Cattle: Unique Genetic Resources in a Rapidly Changing World

At least 150 indigenous African cattle breeds have been named, but the majority of African cattle populations remain largely uncharacterized. As cattle breeds and populations in Africa adapted to various local environmental conditions, they acquired unique features. We know now that the history of African cattle was particularly complex and while several of its episodes remain debated, there is no doubt that African cattle population evolved dramatically over time. Today, we find a mosaic of genetically diverse population from the purest Bos taurus to the nearly pure Bos indicus. African cattle are now found all across the continent, with the exception of the Sahara and the river Congo basin. They are found on the rift valley highlands as well as below sea level in the Afar depression. These unique livestock genetic resources are in danger to disappear rapidly following uncontrolled crossbreeding and breed replacements with exotic breeds. Breeding improvement programs of African indigenous livestock remain too few while paradoxically the demand of livestock products is continually increasing. Many African indigenous breeds are endangered now, and their unique adaptive traits may be lost forever. This paper reviews the unique known characteristics of indigenous African cattle populations while describing the opportunities, the necessity and urgency to understand and utilize these resources to respond to the needs of the people of the continent and to the benefit of African farmers.

Copy number variation in the speciation of pigs: a possible prominent role for olfactory receptors

Background

Unraveling the genetic mechanisms associated with reduced gene flow between genetically differentiated populations is key to understand speciation. Different types of structural variations (SVs) have been found as a source of genetic diversity in a wide range of species. Previous studies provided detailed knowledge on the potential evolutionary role of SVs, especially copy number variations (CNVs), between well diverged species of e.g. primates. However, our understanding of their significance during ongoing speciation processes is limited due to the lack of CNV data from closely related species. The genus Sus (pig and its close relatives) which started to diverge ~4 Mya presents an excellent model for studying the role of CNVs during ongoing speciation.

Results

In this study, we identified 1408 CNV regions (CNVRs) across the genus Sus. These CNVRs encompass 624 genes and were found to evolve ~2.5 times faster than single nucleotide polymorphisms (SNPs). The majority of these copy number variable genes are olfactory receptors (ORs) known to play a prominent role in food foraging and mate recognition in Sus. Phylogenetic analyses, including novel Bayesian analysis, based on CNVRs that overlap ORs retain the well-accepted topology of the genus Sus whereas CNVRs overlapping genes other than ORs show evidence for random drift and/or admixture.

Conclusion

We hypothesize that inter-specific variation in copy number of ORs provided the means for rapid adaptation to different environments during the diversification of the genus Sus in the Pliocene. Furthermore, these regions might have acted as barriers preventing massive gene flow between these species during the multiple hybridization events that took place later in the Pleistocene suggesting a possible prominent role of ORs in the ongoing Sus speciation.

Electronic supplementary material

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

IL-37 requires IL-18Rα and SIGIRR/IL-1R8 to diminish allergic airway inflammation in mice

BACKGROUND

Interleukin (IL) 37 has been described as a negative regulator of innate immunity, as it reduces the activation and cytokine production of different innate immune cells. Recently, results from the CLARA childhood asthma cohort suggested an implication of IL-37 for human asthma pathogenesis. This study aimed to investigate the effects of IL-37 on allergic airway inflammation in a mouse model of experimental asthma.

METHODS

Peripheral blood mononuclear cells (PBMCs) of children were cultured for 48 h (anti-CD3/anti-CD28 stimulation or unstimulated), and IL-37 concentrations in supernatants were determined. Wild-type, IL-18Rα-deficient ((-/-) ), and SIGIRR(-/-) C57BL/6 mice were sensitized to ovalbumin (OVA) and challenged with OVA aerosol to induce acute experimental asthma, and IL-37 was applied intranasally prior to each OVA challenge. Airway hyper-responsiveness (AHR), airway inflammation, cytokine levels in broncho-alveolar lavage fluid, and mucus production were determined.

RESULTS

IL-37 production of human PBMCs was significantly lower in allergic asthmatics vs healthy children. In wild-type mice, intranasal administration of IL-37 ablated allergic airway inflammation as well as cytokine production and subsequently diminished the hallmarks of experimental asthma including mucus hyperproduction and AHR. In contrast, local application of IL-37 produced none of these effects in mice lacking either IL18Rα or SIGIRR/IL-1R8.

CONCLUSIONS

This study demonstrates that IL-37 is able to ablate a TH2 cell-directed allergic inflammatory response and the hallmarks of experimental asthma in mice, suggesting that IL-37 may be critical for asthma pathogenesis. Furthermore, these data suggest a mode of action of IL-37 that involves IL18Rα as well as the orphan receptor SIGIRR/IL-1R8.

A molecular perspective on a complex polymorphic inversion system with cytological evidence of multiply reused breakpoints

Genome sequence comparison across the Drosophila genus revealed that some fixed inversion breakpoints had been multiply reused at this long timescale. Cytological studies of Drosophila inversion polymorphism had previously shown that, also at this shorter timescale, some breakpoints had been multiply reused. The paucity of molecularly characterized polymorphic inversion breakpoints has so far precluded contrasting whether cytologically shared breakpoints of these relatively young inversions are actually reused at the molecular level. The E chromosome of Drosophila subobscura stands out because it presents several inversion complexes. This is the case of the E1+2+9+3 arrangement that originated from the ancestral Est arrangement through the sequential accumulation of four inversions (E1, E2, E9 and E3) sharing some breakpoints. We recently identified the breakpoints of inversions E1 and E2, which allowed establishing reuse at the molecular level of the cytologically shared breakpoint of these inversions. Here, we identified and sequenced the breakpoints of inversions E9 and E3, because they share breakpoints at sections 58D and 64C with those of inversions E1 and E2. This has allowed establishing that E9 and E3 originated through the staggered-break mechanism. Most importantly, sequence comparison has revealed the multiple reuse at the molecular level of the proximal breakpoint (section 58D), which would have been used at least by inversions E2, E9 and E3. In contrast, the distal breakpoint (section 64C) might have been only reused once by inversions E1 and E2, because the distal E3 breakpoint is displaced >70 kb from the other breakpoint limits.

Interleukin-37 gene variants segregated anciently coexist during hominid evolution

IL37 is a member of IL-1 cytokine family but conveys anti-inflammatory functions. The biological characteristic and genetic heterogeneity of IL37 are not fully understood yet. Here using the whole-genome sequencing data from 1000 Genomes Project, we performed population and evolutionary genetic analysis of human IL37 gene. First, 2184 IL37 gene sequences from different human populations were retrieved. The IL37 protein sequences were inferred from the coding DNA sequences and multiple species alignment was made. Then, the phylogenetic tree of IL37 was built and dN/dS ratios were calculated for each evolutionary branch, the classic McDonald and Kreitman test was also performed. Next, we conducted intraspecific evolutionary genetic analysis and built the genealogy network of 116 unique IL37 haplotypes through median-joining network analysis. Finally, we compared IL37 sequences between the modern and archaic humans. Our results for the first time provide solid evidence that common IL37 variants other than NCBI reference sequence are present worldwide. Our data also supports that IL37 variants are shaped and maintained by selection instead of neutral evolution. We further identified that human IL37 variants consist of two major haplogroups and their presence in archaic humans corroborates its ancient origin in hominid evolution. In conclusion, these data indicate that common IL37 variants are maintained among human populations by selective force, suggesting their potential involvements in immune regulation and human diseases. In addition, the ancient history of IL37 variants reveals interesting insight into the complicated human evolutionary history.

Global patterns of apparent copy number variation in birds revealed by cross-species comparative genomic hybridization

There is a growing interest in copy number variation (CNV) and the recognition of its importance in phenotype, disease, adaptation and speciation. CNV data is usually ascertained by array-CGH within-species, but similar inter-species comparisons have also been made in primates, mice and domestic mammals. Here, we conducted a broad appraisal of putative cross-species CNVs in birds, 16 species in all, using the standard array-CGH approach. Using a chicken oligonucleotide microarray, we detected 790 apparent CNVs within 135 unique regions and developed a bioinformatic tool 'CNV Analyser' for analysing and visualising cross-species data sets. We successfully addressed four hypotheses as follows: (a) Cross-species CNVs (compared to chicken) are, as suggested from preliminary evidence, smaller and fewer in number than in mammals; this 'dogma' was rejected in the light of the new evidence. (b) CNVs in birds are likely to have a functional effect through an association with genes; a large proportion of detected regions (70 %) were indeed associated with genes (suggesting functional significance), however, not necessarily more so than in mammals. (c) There are more CNVs in birds with more rearranged karyotypes; this hypothesis was rejected. Indeed, Falco species contained fewer than most with relatively standard (chicken-like) karyotypes. (d) There are more CNVs per megabase on micro-chromosomes than macrochromosomes; this hypothesis was accepted. Indeed, in species with rearranged karyotypes characterised by chromosomal fusions, the fused former microchromosomes still 'behaved' as though they were their microchromosomal ancestors. Gene ontology analysis of CNVRs revealed enrichment in immune response and antigen presentation genes and five CNVRs were perfectly correlated with the unique loss of sexual dichromatism in one Galliformes species.

Novel H3K4me3 marks are enriched at human- and chimpanzee-specific cytogenetic structures

Human and chimpanzee genomes are 98.8% identical within comparable sequences. However, they differ structurally in nine pericentric inversions, one fusion that originated human chromosome 2, and content and localization of heterochromatin and lineage-specific segmental duplications. The possible functional consequences of these cytogenetic and structural differences are not fully understood and their possible involvement in speciation remains unclear. We show that subtelomeric regions—regions that have a species-specific organization, are more divergent in sequence, and are enriched in genes and recombination hotspots—are significantly enriched for species-specific histone modifications that decorate transcription start sites in different tissues in both human and chimpanzee. The human lineage-specific chromosome 2 fusion point and ancestral centromere locus as well as chromosome 1 and 18 pericentric inversion breakpoints showed enrichment of human-specific H3K4me3 peaks in the prefrontal cortex. Our results reveal an association between plastic regions and potential novel regulatory elements.

A fine-scale recombination map of the human-chimpanzee ancestor reveals faster change in humans than in chimpanzees and a strong impact of GC-biased gene conversion

Recombination is a major determinant of adaptive and nonadaptive evolution. Understanding how the recombination landscape has evolved in humans is thus key to the interpretation of human genomic evolution. Comparison of fine-scale recombination maps of human and chimpanzee has revealed large changes at fine genomic scales and conservation over large scales. Here we demonstrate how a fine-scale recombination map can be derived for the ancestor of human and chimpanzee, allowing us to study the changes that have occurred in human and chimpanzee since these species diverged. The map is produced from more than one million accurately determined recombination events. We find that this new recombination map is intermediate to the maps of human and chimpanzee but that the recombination landscape has evolved more rapidly in the human lineage than in the chimpanzee lineage. We use the map to show that recombination rate, through the effect of GC-biased gene conversion, is an even stronger determinant of base composition evolution than previously reported.

Sequence assembly demystified

Advances in sequencing technologies and increased access to sequencing services have led to renewed interest in sequence and genome assembly. Concurrently, new applications for sequencing have emerged, including gene expression analysis, discovery of genomic variants and metagenomics, and each of these has different needs and challenges in terms of assembly. We survey the theoretical foundations that underlie modern assembly and highlight the options and practical trade-offs that need to be considered, focusing on how individual features address the needs of specific applications. We also review key software and the interplay between experimental design and efficacy of assembly.

Molecular trajectories leading to the alternative fates of duplicate genes

Gene duplication generates extra gene copies in which mutations can accumulate without risking the function of pre-existing genes. Such mutations modify duplicates and contribute to evolutionary novelties. However, the vast majority of duplicates appear to be short-lived and experience duplicate silencing within a few million years. Little is known about the molecular mechanisms leading to these alternative fates. Here we delineate differing molecular trajectories of a relatively recent duplication event between humans and chimpanzees by investigating molecular properties of a single duplicate: DNA sequences, gene expression and promoter activities. The inverted duplication of the Glutathione S-transferase Theta 2 (GSTT2) gene had occurred at least 7 million years ago in the common ancestor of African great apes and is preserved in chimpanzees (Pan troglodytes), whereas a deletion polymorphism is prevalent in humans. The alternative fates are associated with expression divergence between these species, and reduced expression in humans is regulated by silencing mutations that have been propagated between duplicates by gene conversion. In contrast, selective constraint preserved duplicate divergence in chimpanzees. The difference in evolutionary processes left a unique DNA footprint in which dying duplicates are significantly more similar to each other (99.4%) than preserved ones. Such molecular trajectories could provide insights for the mechanisms underlying duplicate life and death in extant genomes.

What have studies of genomic disorders taught us about our genome?

The elucidation of genomic disorders began with molecular technologies that enabled detection of genomic changes which were (a) smaller than those resolved by traditional cytogenetics (less than 5 Mb) and (b) larger than what could be determined by conventional gel electrophoresis. Methods such as pulsed field gel electrophoresis (PFGE) and fluorescent in situ hybridization (FISH) could resolve such changes but were limited to locus-specific studies. The study of genomic disorders has rapidly advanced with the development of array-based techniques. These enabled examination of the entire human genome at a higher level of resolution, thus allowing elucidation of the basis of many new disorders, mechanisms that result in genomic changes that can result in copy number variation (CNV), and most importantly, a deeper understanding of the characteristics, features, and plasticity of our genome. In this chapter, we focus on the structural and architectural features of the genome, which can potentially result in genomic instability, delineate how mechanisms, such as NAHR, NHEJ, and FoSTeS/MMBIR lead to disease-causing rearrangements, and briefly describe the relationship between the leading methods presently used in studying genomic disorders. We end with a discussion on our new understanding about our genome including: the contribution of new mutation CNV to disease, the abundance of mosaicism, the extent of subtelomeric rearrangements, the frequency of de novo rearrangements associated with sporadic birth defects, the occurrence of balanced and unbalanced translocations, the increasing discovery of insertional translocations, the exploration of complex rearrangements and exonic CNVs. In the postgenomic era, our understanding of the genome has advanced very rapidly as the level of technical resolution has become higher. This leads to a greater understanding of the effects of rearrangements present both in healthy subjects and individuals with clinically relevant phenotypes.

A fine-scale chimpanzee genetic map from population sequencing

To study the evolution of recombination rates in apes, we developed methodology to construct a fine-scale genetic map from high-throughput sequence data from 10 Western chimpanzees, Pan troglodytes verus. Compared to the human genetic map, broad-scale recombination rates tend to be conserved, but with exceptions, particularly in regions of chromosomal rearrangements and around the site of ancestral fusion in human chromosome 2. At fine scales, chimpanzee recombination is dominated by hotspots, which show no overlap with those of humans even though rates are similarly elevated around CpG islands and decreased within genes. The hotspot-specifying protein PRDM9 shows extensive variation among Western chimpanzees, and there is little evidence that any sequence motifs are enriched in hotspots. The contrasting locations of hotspots provide a natural experiment, which demonstrates the impact of recombination on base composition.

Advances in BAC-based physical mapping and map integration strategies in plants

In the advent of next-generation sequencing (NGS) platforms, map-based sequencing strategy has been recently suppressed being too expensive and laborious. The detailed studies on NGS drafts alone indicated these assemblies remain far from gold standard reference quality, especially when applied on complex genomes. In this context the conventional BAC-based physical mapping has been identified as an important intermediate layer in current hybrid sequencing strategy. BAC-based physical map construction and its integration with high-density genetic maps have benefited from NGS and high-throughput array platforms. This paper addresses the current advancements of BAC-based physical mapping and high-throughput map integration strategies to obtain densely anchored well-ordered physical maps. The resulted maps are of immediate utility while providing a template to harness the maximum benefits of the current NGS platforms.

Generation of long insert pairs using a Cre-LoxP Inverse PCR approach

Large insert mate pair reads have a major impact on the overall success of de novo assembly and the discovery of inherited and acquired structural variants. The positional information of mate pair reads generally improves genome assembly by resolving repeat elements and/or ordering contigs. Currently available methods for building such libraries have one or more of limitations, such as relatively small insert size; unable to distinguish the junction of two ends; and/or low throughput. We developed a new approach, Cre-LoxP Inverse PCR Paired-End (CLIP-PE), which exploits the advantages of (1) Cre-LoxP recombination system to efficiently circularize large DNA fragments, (2) inverse PCR to enrich for the desired products that contain both ends of the large DNA fragments, and (3) the use of restriction enzymes to introduce a recognizable junction site between ligated fragment ends and to improve the self-ligation efficiency. We have successfully created CLIP-PE libraries up to 22 kb that are rich in informative read pairs and low in small fragment background. These libraries have demonstrated the ability to improve genome assemblies. The CLIP-PE methodology can be implemented with existing and future next-generation sequencing platforms.

Fidgetin-like1 is a strong candidate for a dynamic impairment of male meiosis leading to reduced testis weight in mice

Background

In a previous work, using an interspecific recombinant congenic mouse model, we reported a genomic region of 23 Mb on mouse chromosome 11 implicated in testis weight decrease and moderate teratozoospermia (∼20–30%), a Quantitative Trait Locus (QTL) called Ltw1. The objective of the present study is to identify the gene underlying this phenotype.

Results

In the present study, we refined the QTL position to a 5 Mb fragment encompassing only 11 genes. We showed that the low testis weight phenotype was due to kinetic alterations occurring during the first wave of the spermatogenesis where we could point out to an abnormal lengthening of spermatocyte prophase. We identify Fidgetin-like 1 (Fignl1) as the gene underlying the phenotype, since if fulfilled both the physiological and molecular characteristics required. Indeed, amongst the 11 positional candidates it is the only gene that is expressed during meiosis at the spermatocyte stage, and that presents with non-synonymous coding variations differentiating the two mouse strains at the origin of the cross.

Conclusions

This work prompted us to propose Fignl1 as a novel actor in mammal's male meiosis dynamics which has fundamental interest. Besides, this gene is a new potential candidate for human infertilities caused by teratozoospermia and blockades of spermatogenesis. In addition this study demonstrates that interspecific models may be useful for understanding complex quantitative traits.

Evolutionary patterns of recently emerged animal duplogs

Duplogs, or intraspecies paralogs, constitute the important portion of eukaryote genomes and serve as a major source of functional innovation. We conducted detailed analyses of recently emerged animal duplogs. Genome data of three vertebrate species (Homo sapiens, Mus musculus, and Danio rerio), Caenorhabditis elegans, and two Drosophila species (Drosophila melanogaster and D. pseudoobscura) were used. Duplication events were divided into six age-groups according to the synonymous distance (dS) up to 0.6. Duplogs were classified into four equal-sized classes on physical distances and into three classes on relative orientations. We observed the following shared characteristics among intrachromosomal multiexon duplogs: 1) inverted duplogs account for 20-50%, and about a half of the physically most distant 25%; 2) except for C. elegans, the composition of physical distances, that of relative orientations, and the proportion of inverted duplogs in each physical distance category are more or less uniform; 3) except for C. elegans, the characteristics of the youngest (dS < 0.01) duplogs are similar to the overall characteristics of the entire set. These results suggest that intrachromosomal duplogs with fairly long physical distances were generated at once, rather than resulting from tandem duplications and subsequent genomic rearrangements. This is different from the three well-known modes of gene duplication: tandem duplication, retrotransposition, and genome duplication. We termed this new mode as "drift" duplication. The drift duplication has been producing duplicate copies at paces comparable with tandem duplications since the common ancestor of vertebrates, and it may have already operated in the common ancestor of bilateral animals.

Copy number variation analysis in the great apes reveals species-specific patterns of structural variation

Copy number variants (CNVs) are increasingly acknowledged as an important source of evolutionary novelties in the human lineage. However, our understanding of their significance is still hindered by the lack of primate CNV data. We performed intraspecific comparative genomic hybridizations to identify loci harboring copy number variants in each of the four great apes: bonobos, chimpanzees, gorillas, and orangutans. For the first time, we could analyze differences in CNV location and frequency in these four species, and compare them with human CNVs and primate segmental duplication (SD) maps. In addition, for bonobo and gorilla, patterns of CNV and nucleotide diversity were studied in the same individuals. We show that CNVs have been subject to different selective pressures in different lineages. Evidence for purifying selection is stronger in gorilla CNVs overlapping genes, while positive selection appears to have driven the fixation of structural variants in the orangutan lineage. In contrast, chimpanzees and bonobos present high levels of common structural polymorphism, which is indicative of relaxed purifying selection together with the higher mutation rates induced by the known burst of segmental duplication in the ancestor of the African apes. Indeed, the impact of the duplication burst is noticeable by the fact that bonobo and chimpanzee share more CNVs with gorilla than expected. Finally, we identified a number of interesting genomic regions that present high-frequency CNVs in all great apes, while containing only very rare or even pathogenic structural variants in humans.

Transcriptional variations mediated by an alternative promoter of the FPR3 gene

Formyl peptide receptor 3 (FPR3) is a potential player in innate immunity and appears with FPR2 as a FPR cluster during primate evolution. Comparative genome analyses indicate that a segmental duplication (SD) event upstream of the FPR3 gene after the divergence of New and Old World monkeys led to the emergence of an alternative promoter. In this study we combined computational and experimental approaches to identify a FPR3 gene that is controlled by an alternative promoter derived during a SD event. Its transcriptional activity was detected by quantitative reverse transcription polymerase chain reaction. Human alternative transcripts (FPR3-1 and FPR3-2) showed tissue-specific patterns with strong expressions in lung or uterus, while the FPR3-1 transcript of rhesus macaque is broadly expressed in various tissues. Overall, transcriptional variations of FPR3 occur by an alternative promoter during primate evolution.

Compressing genomic sequence fragments using SlimGene

With the advent of next generation sequencing technologies, the cost of sequencing whole genomes is poised to go below $1000 per human individual in a few years. As more and more genomes are sequenced, analysis methods are undergoing rapid development, making it tempting to store sequencing data for long periods of time so that the data can be re-analyzed with the latest techniques. The challenging open research problems, huge influx of data, and rapidly improving analysis techniques have created the need to store and transfer very large volumes of data. Compression can be achieved at many levels, including trace level (compressing image data), sequence level (compressing a genomic sequence), and fragment-level (compressing a set of short, redundant fragment reads, along with quality-values on the base-calls). We focus on fragment-level compression, which is the pressing need today. Our article makes two contributions, implemented in a tool, SlimGene. First, we introduce a set of domain specific loss-less compression schemes that achieve over 40× compression of fragments, outperforming bzip2 by over 6×. Including quality values, we show a 5× compression using less running time than bzip2. Second, given the discrepancy between the compression factor obtained with and without quality values, we initiate the study of using "lossy" quality values. Specifically, we show that a lossy quality value quantization results in 14× compression but has minimal impact on downstream applications like SNP calling that use the quality values. Discrepancies between SNP calls made between the lossy and loss-less versions of the data are limited to low coverage areas where even the SNP calls made by the loss-less version are marginal.

Gorilla genome structural variation reveals evolutionary parallelisms with chimpanzee

Structural variation has played an important role in the evolutionary restructuring of human and great ape genomes. Recent analyses have suggested that the genomes of chimpanzee and human have been particularly enriched for this form of genetic variation. Here, we set out to assess the extent of structural variation in the gorilla lineage by generating 10-fold genomic sequence coverage from a western lowland gorilla and integrating these data into a physical and cytogenetic framework of structural variation. We discovered and validated over 7665 structural changes within the gorilla lineage, including sequence resolution of inversions, deletions, duplications, and mobile element insertions. A comparison with human and other ape genomes shows that the gorilla genome has been subjected to the highest rate of segmental duplication. We show that both the gorilla and chimpanzee genomes have experienced independent yet convergent patterns of structural mutation that have not occurred in humans, including the formation of subtelomeric heterochromatic caps, the hyperexpansion of segmental duplications, and bursts of retroviral integrations. Our analysis suggests that the chimpanzee and gorilla genomes are structurally more derived than either orangutan or human genomes.

Close 3D proximity of evolutionary breakpoints argues for the notion of spatial synteny

Background

Folding and intermingling of chromosomes has the potential of bringing close to each other loci that are very distant genomically or even on different chromosomes. On the other hand, genomic rearrangements also play a major role in the reorganisation of loci proximities. Whether the same loci are involved in both mechanisms has been studied in the case of somatic rearrangements, but never from an evolutionary standpoint.

Results

In this paper, we analysed the correlation between two datasets: (i) whole-genome chromatin contact data obtained in human cells using the Hi-C protocol; and (ii) a set of breakpoint regions resulting from evolutionary rearrangements which occurred since the split of the human and mouse lineages. Surprisingly, we found that two loci distant in the human genome but adjacent in the mouse genome are significantly more often observed in close proximity in the human nucleus than expected. Importantly, we show that this result holds for loci located on the same chromosome regardless of the genomic distance separating them, and the signal is stronger in gene-rich and open-chromatin regions.

Conclusions

These findings strongly suggest that part of the 3D organisation of chromosomes may be conserved across very large evolutionary distances. To characterise this phenomenon, we propose to use the notion of spatial synteny which generalises the notion of genomic synteny to the 3D case.