Chromosome Inversions,Genomic Differentiation and Speciation in the African Malaria Mosquito Anopheles gambiae

The African malaria vector, Anopheles gambiae, is characterized by multiple polymorphic chromosomal inversions and has become widely studied as a system for exploring models of speciation. Near complete reproductive isolation between different inversion types, known as chromosomal forms, has led to the suggestion that A. gambiae is in early stages of speciation, with divergence evolving in the face of considerable gene flow. We compared the standard chromosomal arrangement (Savanna form) with genomes homozygous for j, b, c, and u inversions (Bamako form) in order to identify regions of genomic divergence with respect to inversion polymorphism. We found levels of divergence between the two sub-taxa within some of these inversions (2Rj and 2Rb), but at a level lower than expected and confined near the inversion breakpoints, consistent with a gene flux model. Unexpectedly, we found that the majority of diverged regions were located on the X chromosome, which contained half of all significantly diverged regions, with much of this divergence located within exons. This is surprising given that the Bamako and Savanna chromosomal forms are both within the S molecular form that is defined by a locus near centromere of X chromosome. Two X-linked genes (a heat shock protein and P450 encoding genes) involved in reproductive isolation between the M and S molecular forms of A. gambiae were also significantly diverged between the two chromosomal forms. These results suggest that genes mediating reproductive isolation are likely located on the X chromosome, as is thought to be the case for the M and S molecular forms. We conclude that genes located on the sex chromosome may be the major force driving speciation between these chromosomal forms of A. gambiae.     

The evolutionary history of Drosophila buzzatii. XXXV. Inversion polymorphism and nucleotide variability in different regions of the second chromosome

Inversions are portions of a chromosome where the gene order is reversed relative to a standard reference orientation. Because of reduced levels of recombination in heterokaryotypes, inversions have a potentially important effect on patterns of nucleotide variability in those genomic regions close to, or included in, the inverted fragments. Here we report sequence variation at three anonymous regions (STSs) located at different positions in relation to second-chromosome inversion breakpoints in 29 isochromosomal lines derived from an Argentinean population of Drosophila buzzatii. In agreement with previous findings in Drosophila, gene flux (crossing over and/or gene conversion) between arrangements seems to appreciably increase as we approach the middle sections of inversion 2j, and patterns of nucleotide variability within, as well as genetic differentiation between chromosome arrangements, are comparable to those observed at the molecular marker outside the inverted fragments. On the other hand, nucleotide diversity near the proximal breakpoint of inversion 2j is reduced when contrasted with that found at the other regions, particularly in the case of derived inverted chromosomes. Using the data from the breakpoint, we estimate that the inversion polymorphism is approximately 1.63 N generations old, where N is the effective population size. An excess of low-frequency segregating polymorphisms is detected; mostly in the ancestral 2st arrangement and probably indicating a population expansion that predates the coalescent time of inversion 2j. Heterogeneity in mutation rates between the markers linked to the inversions may be sufficient to explain the different levels of nucleotide diversity observed. When considered in the context of other studies on patterns of variation relative to physical distance to inversion breakpoints, our data appear to be consistent with the conclusion that inversions are unlikely to be "long-lived" balanced polymorphisms.     

Sex-specific Trans-regulatory Variation on the Drosophila melanogaster X Chromosome

The X chromosome constitutes a unique genomic environment because it is present in one copy in males, but two copies in females. This simple fact has motivated several theoretical predictions with respect to how standing genetic variation on the X chromosome should differ from the autosomes. Unmasked expression of deleterious mutations in males and a lower census size are expected to reduce variation, while allelic variants with sexually antagonistic effects, and potentially those with a sex-specific effect, could accumulate on the X chromosome and contribute to increased genetic variation. In addition, incomplete dosage compensation of the X chromosome could potentially dampen the male-specific effects of random mutations, and promote the accumulation of X-linked alleles with sexually dimorphic phenotypic effects. Here we test both the amount and the type of genetic variation on the X chromosome within a population of Drosophila melanogaster, by comparing the proportion of X linked and autosomal trans-regulatory SNPs with a sexually concordant and discordant effect on gene expression. We find that the X chromosome is depleted for SNPs with a sexually concordant effect, but hosts comparatively more SNPs with a sexually discordant effect. Interestingly, the contrasting results for SNPs with sexually concordant and discordant effects are driven by SNPs with a larger influence on expression in females than expression in males. Furthermore, the distribution of these SNPs is shifted towards regions where dosage compensation is predicted to be less complete. These results suggest that intrinsic properties of dosage compensation influence either the accumulation of different types of trans-factors and/or their propensity to accumulate mutations. Our findings document a potential mechanistic basis for sex-specific genetic variation, and identify the X as a reservoir for sexually dimorphic phenotypic variation. These results have general implications for X chromosome evolution, as well as the genetic basis of sex-specific evolutionary change.     

Structure and decay of a proto-Y region in Tilapia,Oreochromis niloticus

Background

Sex-determination genes drive the evolution of adjacent chromosomal regions. Sexually antagonistic selection favors the accumulation of inversions that reduce recombination in regions adjacent to the sex-determination gene. Once established, the clonal inheritance of sex-linked inversions leads to the accumulation of deleterious alleles, repetitive elements and a gradual decay of sex-linked genes. This in turn creates selective pressures for the evolution of mechanisms that compensate for the unequal dosage of gene expression. Here we use whole genome sequencing to characterize the structure of a young sex chromosome and quantify sex-specific gene expression in the developing gonad.

Results

We found an 8.8 Mb block of strong differentiation between males and females that corresponds to the location of a previously mapped sex-determiner on linkage group 1 of Oreochromis niloticus. Putatively disruptive mutations are found in many of the genes within this region. We also found a significant female-bias in the expression of genes within the block of differentiation compared to those outside the block of differentiation. Eight candidate sex-determination genes were identified within this region.

Conclusions

This study demonstrates a block of differentiation on linkage group 1, suggestive of an 8.8 Mb inversion encompassing the sex-determining locus. The enrichment of female-biased gene expression inside the proposed inversion suggests incomplete dosage compensation. This study helps establish a model for studying the early-to-intermediate stages of sex chromosome evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-975) contains supplementary material, which is available to authorized users.     

Ability of a Bacterial Chromosome Segment to Invert Is Dictated by Included Material Rather than Flanking Sequence

Homologous recombination between sequences present in inverse order within the same chromosome can result in inversion formation. We have previously shown that inverse order sequences at some sites (permissive) recombine to generate the expected inversion; no inversions are found when the same inverse order sequences flank other (nonpermissive) regions of the chromosome. In hopes of defining how permissive and nonpermissive intervals are determined, we have constructed a strain that carries a large chromosomal inversion. Using this inversion mutant as the parent strain, we have determined the "permissivity" of a series of chromosomal sites for secondary inversions. For the set of intervals tested, permissivity seems to be dictated by the nature of the genetic material present within the chromosomal interval being tested rather than the flanking sequences or orientation of this material in the chromosome. Almost all permissive intervals include the origin or terminus of replication. We suggest that the rules for recovery of inversions reflect mechanistic restrictions on the occurrence of inversions rather than lethal consequences of the completed rearrangement.     

Balanced translocation t(3;18)(p13;q22.3) and points mutation in the ZNF407 gene detected in patients with both moderate non-syndromic intellectual disability and autism

Intellectual disability (ID) is a common disease. While the etiology remains incompletely understood, genetic defects are a major contributor, which include mutations in genes encoding zinc finger proteins. These proteins modulate gene expression via binding to DNA. Consistent with this knowledge, we report here the identification of mutations in the ZNF407 gene in ID/autistic patients. In our study of an ID patient with autism, a reciprocal translocation 46,XY,t(3;18)(p13;q22.3) was detected. By using FISH and long-range PCR approaches, we have precisely mapped the breakpoints associated with this translocation in a gene-free region in chromosome 3 and in the third intron of the ZNF407 gene in chromosome18. The latter reduces ZNF407 expression. Consistent with this observation, in our subsequent investigation of 105 ID/autism patients with similar clinical presentations, two missense mutations Y460C and P1195A were identified. These mutations cause non-conservative amino acid substitutions in the linker regions between individual finger structures. In line with the linker regions being critical for the integrity of zinc finger motifs, both mutations may result in loss of ZNF407 function. Taken together, we demonstrate that mutations in the ZNF407 gene contribute to the pathogenesis of a group of ID patients with autism.     

The Question of the Total Gene Number in DROSOPHILA MELANOGASTER

A statistical analysis has been carried out on the distribution and allelism of nearly 500 sex-linked, X-ray-induced, cytologically normal and rearranged lethal mutations in Drosophila melanogaster that were obtained by G. Lefevre. The mutations were induced in four different regions of the X chromosome: (1) 1A1-3E8, (2) 6D1-8A5, (3) 9E1-11A7 and (4) 19A1-20F4, which together comprise more than one-third of the entire chromosome.--The analysis shows that the number of alleles found at different loci does not fit a Poisson distribution, even when the proper procedures are taken to accommodate the truncated nature of the data. However, the allele distribution fits a truncated negative binomial distribution quite well, with cytologically normal mutations fitting better than rearrangement mutations. This indicates that genes are not equimutable, as required for the data to fit a Poisson distribution.--Using the negative binomial parameters to estimate the number of genes that did not produce a detectable lethal mutation in our experiment (n0) gave a larger number than that derived from the use of the Poisson parameter. Unfortunately, we cannot estimate the total numbers of nonvital loci, loci with undetectable phenotypes and loci having extremely low mutabilities. In any event, our estimate of the total vital gene number was far short of the total number of bands in the analyzed regions; yet, in several short intervals, we have found more vital genes than bands; in other intervals, fewer. We conclude that the one-band, one-gene hypothesis, in its literal sense, is not true; furthermore, it is difficult to support, even approximately.--The question of the total gene number in Drosophila will, not doubt, eventually be solved by molecular analyses, not by statistical analysis of mutation data or saturation studies.     

Complex Patterns of Chromosome 11 Aberrations in Myeloid Malignancies Target CBL,MLL, DDB1 and LMO2

Exome sequencing of primary tumors identifies complex somatic mutation patterns. Assignment of relevance of individual somatic mutations is difficult and poses the next challenge for interpretation of next generation sequencing data. Here we present an approach how exome sequencing in combination with SNP microarray data may identify targets of chromosomal aberrations in myeloid malignancies. The rationale of this approach is that hotspots of chromosomal aberrations might also harbor point mutations in the target genes of deletions, gains or uniparental disomies (UPDs). Chromosome 11 is a frequent target of lesions in myeloid malignancies. Therefore, we studied chromosome 11 in a total of 813 samples from 773 individual patients with different myeloid malignancies by SNP microarrays and complemented the data with exome sequencing in selected cases exhibiting chromosome 11 defects. We found gains, losses and UPDs of chromosome 11 in 52 of the 813 samples (6.4%). Chromosome 11q UPDs frequently associated with mutations of CBL. In one patient the 11qUPD amplified somatic mutations in both CBL and the DNA repair gene DDB1. A duplication within MLL exon 3 was detected in another patient with 11qUPD. We identified several common deleted regions (CDR) on chromosome 11. One of the CDRs associated with de novo acute myeloid leukemia (P=0.013). One patient with a deletion at the LMO2 locus harbored an additional point mutation on the other allele indicating that LMO2 might be a tumor suppressor frequently targeted by 11p deletions. Our chromosome-centered analysis indicates that chromosome 11 contains a number of tumor suppressor genes and that the role of this chromosome in myeloid malignancies is more complex than previously recognized.     

Two new balancer chromosomes on mouse chromosome 4 to facilitate functional annotation of human chromosome 1p

To facilitate genetic screens to identify and maintain recessive mutations that map to the short arm of human chromosome 1, we have utilized chromosome engineering to generate two mouse strains that carry large inversions on the distal region of mouse chromosome 4. The inversion intervals are 16 and 22 cM in size together they cover approximately half of chromosome 4. Since recombination between the wild-type and inversion chromosomes does not occur within these inversion intervals, mutant alleles of genes mapping to this region can be identified and maintained. Therefore, these inversion chromosomes work as balancer chromosomes. These inversions have the additional advantage that they are tagged with genes encoding the visible coat color markers tyrosinase and agouti, and therefore the dosage of the inversion chromosome (+/+, Inv/+, Inv/Inv) can be visually recognized. These inversion strains will be extremely useful for mutagenesis screens that focus on functional annotation of human chromosome 1p.     

Chromosome-wide linkage disequilibrium caused by an inversion polymorphism in the white-throated sparrow (Zonotrichia albicollis)

Chromosomal inversions have been of long-standing interest to geneticists because they are capable of suppressing recombination and facilitating the formation of adaptive gene complexes. An exceptional inversion polymorphism (ZAL2^m) in the white-throated sparrow (Zonotrichia albicollis) is linked to variation in plumage, social behavior and mate choice, and is maintained in the population by negative assortative mating. The ZAL2^m polymorphism is a complex inversion spanning >100 Mb and has been proposed to be a strong suppressor of recombination, as well as a potential model for studying neo-sex chromosome evolution. To quantify and evaluate these features of the ZAL2^m polymorphism, we generated sequence from 8 ZAL2^m and 16 ZAL2 chromosomes at 58 loci inside and 4 loci outside the inversion. Inside the inversion we found that recombination was completely suppressed between ZAL2 and ZAL2^m, resulting in uniformly high levels of genetic differentiation (F_ST=0.94), the formation of two distinct haplotype groups representing the alternate chromosome arrangements and extensive linkage disequilibrium spanning ∼104 Mb within the inversion, whereas gene flow was not suppressed outside the inversion. Finally, although ZAL2^m homozygotes are exceedingly rare in the population, occurring at a frequency of <1%, we detected evidence of historical recombination between ZAL2^m chromosomes inside the inversion, refuting its potential status as a non-recombining autosome.     

STRUCTURE AND POPULATION GENETICS OF THE BREAKPOINTS OF A POLYMORPHIC INVERSION IN DROSOPHILA SUBOBSCURA

Drosophila subobscura is a paleartic species of the obscura group with a rich chromosomal polymorphism. To further our understanding on the origin of inversions and on how they regain variation, we have identified and sequenced the two breakpoints of a polymorphic inversion of D. subobscura—inversion 3 of the O chromosome—in a population sample. The breakpoints could be identified as two rather short fragments (～300 bp and 60 bp long) with no similarity to any known transposable element family or repetitive sequence. The presence of the ～300‐bp fragment at the two breakpoints of inverted chromosomes implies its duplication, an indication of the inversion origin via staggered double‐strand breaks. Present results and previous findings support that the mode of origin of inversions is neither related to the inversion age nor species‐group specific. The breakpoint regions do not consistently exhibit the lower level of variation within and stronger genetic differentiation between arrangements than more internal regions that would be expected, even in moderately small inversions, if gene conversion were greatly restricted at inversion breakpoints. Comparison of the proximal breakpoint region in species of the obscura group shows that this breakpoint lies in a small high‐turnover fragment within a long collinear region (～300 kb).     

Lethal mutations defining 112 complementation groups in a 4.5 Mb sequenced region of Caenorhabditis elegans chromosome III

The central gene cluster of chromosome III was one of the first regions to be sequenced by the Caenorhabditis elegans genome project. We have performed an essential gene analysis on the left part of this cluster, in the region around dpy-17III balanced by the duplication sDp3. We isolated 151 essential gene mutations and characterized them with regard to their arrest stages. To facilitate positioning of these mutations, we generated six new deficiencies that, together with preexisting chromosomal rearrangements, subdivide the region into 14 zones. The 151 mutations were mapped into these zones. They define 112 genes, of which 110 were previously unidentified. Thirteen of the zones have been anchored to the physical sequence by polymerase chain reaction deficiency mapping. Of the 112 essential genes mapped, 105 are within these 13 zones. They span 4.2?Mb of nucleotide sequence. From the nucleotide sequence data, 920 genes are predicted. From a Poisson distribution of our mutations, we predict that 234 of the genes will be essential genes. Thus, the 105 genes constitute 45% of the estimated number of essential genes in the physically defined zones and between 2 and 5% of all essential genes in C. elegans.     

An HMM-Based Comparative Genomic Framework for Detecting Introgression in Eukaryotes

One outcome of interspecific hybridization and subsequent effects of evolutionary forces is introgression, which is the integration of genetic material from one species into the genome of an individual in another species. The evolution of several groups of eukaryotic species has involved hybridization, and cases of adaptation through introgression have been already established. In this work, we report on PhyloNet-HMM—a new comparative genomic framework for detecting introgression in genomes. PhyloNet-HMM combines phylogenetic networks with hidden Markov models (HMMs) to simultaneously capture the (potentially reticulate) evolutionary history of the genomes and dependencies within genomes. A novel aspect of our work is that it also accounts for incomplete lineage sorting and dependence across loci. Application of our model to variation data from chromosome 7 in the mouse (Mus musculus domesticus) genome detected a recently reported adaptive introgression event involving the rodent poison resistance gene Vkorc1, in addition to other newly detected introgressed genomic regions. Based on our analysis, it is estimated that about 9% of all sites within chromosome 7 are of introgressive origin (these cover about 13 Mbp of chromosome 7, and over 300 genes). Further, our model detected no introgression in a negative control data set. We also found that our model accurately detected introgression and other evolutionary processes from synthetic data sets simulated under the coalescent model with recombination, isolation, and migration. Our work provides a powerful framework for systematic analysis of introgression while simultaneously accounting for dependence across sites, point mutations, recombination, and ancestral polymorphism.     

Sequence signatures of a recent chromosomal rearrangement in <Emphasis Type="Italic">Drosophila mojavensis</Emphasis>

The X-chromosome inversion, Xe, distinguishes Drosophila mojavensis and D. arizonae. Earlier work mapped the breakpoints of this inversion to large intervals and provided hypotheses for the locations of the breakpoints within 3000-bp intergenic regions on the D. mojavensis genome sequence assembly. Here, we sequenced these regions directly in the putatively ancestral D. arizonae X-chromosome. We find that the two inversion breakpoints are near an inverted gene duplication and a common repetitive element, respectively, and these features were likely present in the non-inverted ancestral chromosome on the D. mojavensis lineage. Contrary to an earlier hypothesis, the inverted gene duplication appears to predate the inversion. We find no sequence similarity between the breakpoint regions in the D. mojavensis ancestor, excluding an ectopic-exchange model of chromosome rearrangements. We also found no evidence that staggered single-strand breaks caused the inversion. We suggest these features may have contributed to the chromosomal breakages resulting in this inversion. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Screening <Emphasis Type="Italic">in silico</Emphasis> predicted remotely acting <Emphasis Type="Italic">NF1</Emphasis>gene regulatory elements for mutations in patients with neurofibromatosis type 1

Neurofibromatosis type 1 (NF1), a neuroectodermal disorder, is caused by germline mutations in the NF1 gene. NF1 affects approximately 1/3,000 individuals worldwide, with about 50% of cases representing de novo mutations. Although the NF1 gene was identified in 1990, the underlying gene mutations still remain undetected in a small but obdurate minority of NF1 patients. We postulated that in these patients, hitherto undetected pathogenic mutations might occur in regulatory elements far upstream of the NF1 gene. In an attempt to identify such remotely acting regulatory elements, we reasoned that some of them might reside within DNA sequences that (1) have the potential to interact at distance with the NF1 gene and (2) lie within a histone H3K27ac-enriched region, a characteristic of active enhancers. Combining Hi-C data, obtained by means of the chromosome conformation capture technique, with data on the location and level of histone H3K27ac enrichment upstream of the NF1 gene, we predicted in silico the presence of two remotely acting regulatory regions, located, respectively, approximately 600 kb and approximately 42 kb upstream of the NF1 gene. These regions were then sequenced in 47 NF1 patients in whom no mutations had been found in either the NF1 or SPRED1 gene regions. Five patients were found to harbour DNA sequence variants in the distal H3K27ac-enriched region. Although these variants are of uncertain pathological significance and still remain to be functionally characterized, this approach promises to be of general utility for the detection of mutations underlying other inherited disorders that may be caused by mutations in remotely acting regulatory elements.     

Mapping regions within cosmopolitan inversion In(3R)Payne associated with natural variation in body size in Drosophila melanogaster

Associations between genotypes for inversions and quantitative traits have been reported in several organisms, but little has been done to localize regions within inversions controlling variation in these traits. Here, we use an association mapping technique to identify genomic regions controlling variation in wing size within the cosmopolitan inversion In(3R)Payne in Drosophila melanogaster. Previous studies have shown that this inversion strongly influences variation in wing size, a trait highly correlated with body size. We found three alleles from two separate regions within In(3R)Payne with significant additive effects on wing size after the additional effect of the inversion itself had been taken into account. There were also several alleles with significant genotype-by-inversion interaction effects on wing size. None of the alleles tested had a significant additive effect on development time, suggesting different genes control these traits and that clinal patterns in them have therefore arisen independently. The presence of multiple regions within In(3R)Payne controlling size is consistent with the idea that inversions persist in populations because they contain multiple sets of locally adapted alleles, but more work needs to be done to test if they are indeed coadapted.     

Approaches to Half-Tetrad Analysis in Bacteria: Recombination between Repeated, Inverse-Order Chromosomal Sequences

In standard bacterial recombination assays, a linear fragment of DNA is transferred to a recipient cell and, at most, a single selected recombinant type is recovered from each merozygote. This contrasts with fungal systems, for which tetrads allow recovery of all meiotic products, including both ultimate recombinant products of an apparent single act of recombination. We have developed a bacterial recombination system in which two recombining sequences are placed in inverse order at widely separated sites in the circular chromosome of Salmonella typhimurium. Recombination can reassort markers between these repeated sequences (double recombination and apparent gene conversion), or can exchange flanking sequences, leading to inversion of the chromosome segment between the recombining sequences. Since two recombinant products remain in the chromosome of a recombinant with an inversion, one can, in principle, approach the capability of tetrad analysis. Using this system, the following observations have been made. (a) When long sequences (40 kb) recombine, conversion frequently accompanies exchange of flanking sequences. (b) When short sequences (5 kb) recombine, conversion rarely accompanies exchange of flanks. (c) Both recA and recB mutations eliminate inversion formation. (d) The frequency of exchanges between short repeats is more sensitive to the distance separating the recombining sequences in the chromosome. The results are presented with the assumption that inversions occur by simple interaction of two sequences in the same circular chromosome. In an appendix we discuss mechanistically more complex possibilities, some of which could also apply to standard fungal systems.     

An Excess of Gene Expression Divergence on the X Chromosome in Drosophila Embryos: Implications for the Faster-X Hypothesis

The X chromosome is present as a single copy in the heterogametic sex, and this hemizygosity is expected to drive unusual patterns of evolution on the X relative to the autosomes. For example, the hemizgosity of the X may lead to a lower chromosomal effective population size compared to the autosomes, suggesting that the X might be more strongly affected by genetic drift. However, the X may also experience stronger positive selection than the autosomes, because recessive beneficial mutations will be more visible to selection on the X where they will spend less time being masked by the dominant, less beneficial allele—a proposal known as the faster-X hypothesis. Thus, empirical studies demonstrating increased genetic divergence on the X chromosome could be indicative of either adaptive or non-adaptive evolution. We measured gene expression in Drosophila species and in D. melanogaster inbred strains for both embryos and adults. In the embryos we found that expression divergence is on average more than 20% higher for genes on the X chromosome relative to the autosomes; but in contrast, in the inbred strains, gene expression variation is significantly lower on the X chromosome. Furthermore, expression divergence of genes on Muller''s D element is significantly greater along the branch leading to the obscura sub-group, in which this element segregates as a neo-X chromosome. In the adults, divergence is greatest on the X chromosome for males, but not for females, yet in both sexes inbred strains harbour the lowest level of gene expression variation on the X chromosome. We consider different explanations for our results and conclude that they are most consistent within the framework of the faster-X hypothesis.