Genetics and Lineage-Specific Evolution of a Lethal Hybrid Incompatibility Between Drosophila mauritiana and Its Sibling Species

The Dobzhansky–Muller model posits that intrinsic postzygotic reproductive isolation—the sterility or lethality of species hybrids—results from the evolution of incompatible epistatic interactions between species: favorable or neutral alleles that become fixed in the genetic background of one species can cause sterility or lethality in the genetic background of another species. The kind of hybrid incompatibility that evolves between two species, however, depends on the particular evolutionary history of the causative substitutions. An allele that is functionally derived in one species can be incompatible with an allele that is functionally derived in the other species (a derived-derived hybrid incompatibility). But an allele that is functionally derived in one species can also be incompatible with an allele that has retained the ancestral state in the other species (a derived-ancestral hybrid incompatibility). The relative abundance of such derived-derived vs. derived-ancestral hybrid incompatibilities is unknown. Here, we characterize the genetics and evolutionary history of a lethal hybrid incompatibility between Drosophila mauritiana and its two sibling species, D. sechellia and D. simulans. We show that a hybrid lethality factor(s) in the pericentric heterochromatin of the D. mauritiana X chromosome, hybrid lethal on the X (hlx), is incompatible with a factor(s) in the same small autosomal region from both D. sechellia and D. simulans, Suppressor of hlx [Su(hlx)]. By combining genetic and phylogenetic information, we infer that hlx-Su(hlx) hybrid lethality is likely caused by a derived-ancestral incompatibility, a hypothesis that can be tested directly when the genes are identified.     

Estimates of natural selection due to protein tertiary structure inform the ancestry of biallelic loci

We consider the inference of which of two alleles is ancestral when the alleles have a single nonsynonymous difference and when natural selection acts via protein tertiary structure. Whereas the probability that an allele is ancestral under neutrality is equal to its frequency, under selection this probability depends on allele frequency and on the magnitude and direction of selection pressure. Although allele frequencies can be well estimated from intraspecific data, small fitness differences have a large evolutionary impact but can be difficult to estimate with only intraspecific data. Methods for predicting aspects of phenotype from genotype can supplement intraspecific sequence data. Recently developed statistical techniques can assess effects of phenotypes, such as protein tertiary structure on molecular evolution. While these techniques were initially designed for comparing protein-coding genes from different species, the resulting interspecific inferences can be assigned population genetic interpretations to assess the effect of selection pressure, and we use them here along with intraspecific allele frequency data to estimate the probability that an allele is ancestral. We focus on 140 nonsynonymous single nucleotide polymorphisms of humans that are in proteins with known tertiary structures. We find that our technique for employing protein tertiary structure information yields some biologically plausible results but that it does not substantially improve the inference of ancestral human allele types.     

Ancestral processes for non-neutral models of complex diseases

We consider non-neutral models for unlinked loci, where the fitness of a chromosome or individual is not multiplicative across loci. Such models are suitable for many complex diseases, where there are gene-interactions. We derive a genealogical process for such models, called the complex selection graph (CSG). This coalescent-type process is related to the ancestral selection graph, and is derived from the ancestral influence graph by considering the limit as the recombination rate between loci gets large. We analyse the CSG both theoretically and via simulation. The main results are that the gene-interactions do not produce linkage disequilibrium, but do produce dependencies in allele frequencies between loci. For small selection rates, the distributions of the genealogy and the allele frequencies at a single locus are well-approximated by their distributions under a single locus model, where the fitness of each allele is the average of the true fitnesses of that allele with respect to the distribution of alleles at other loci.     

基于染色体计数和ITS序列初步探讨横断山区柴胡属植物(伞形科)的系统发育

横断山区是中国柴胡属Bupleurum植物的分布中心。本文对横断山区6个种2变种进行了染色体记数报道,其中4个种2变种是首次报道。对横断山区的10个种4个变种、中国北方（河北和黑龙江）的3个种的nrDNA ITS进行测序,同时从GenBank里面下载同属的来自非洲和地中海西部的16个nrDNA ITS序列数据,结合染色体数目变化结果,初步探讨了横断山区柴胡属植物的系统发育。结果表明横断山区可能是现代柴胡属植物的频度中心和多样分布中心之一。它们的祖先种可能是非洲北部的木本柴胡属植物B.fruticosum,或者是地中海西部的柴胡属植物,推测是通过中东和高加索扩散而形成的,其中与非洲南部特有种B.mundtii的亲缘关系也较近;染色体基数演化趋势是：8是较原始基数,6和7是次生基数,其染色体异基数变异和多倍化可能是物种形成、进化以及向外扩散的主要方式;在ITS系统发育树中,中国柴胡属植物染色体基数为8的种类聚为一支,染色体基数为6和7的种类聚为了一支,不支持舒璞等（1998）关于中国柴胡属的属下分类系统。结合已有的形态学、细胞学、孢粉学证据和ITS系统发育树,建议窄竹叶柴胡B.marginatum var.stenophyllum独立成种。     

Phylogenetic inference of the genus Bupleurum (Apiaceae) in Hengduan Mountains based on chromosome counts and nuclear ribosomal DNA ITS sequences

pleurum of Apiaceae in China. This paper reports chromosome numbers of six species and two varieties of Bupleurum, and for four species and two varieties their chromosome numbers are reported for the first time. The phylogeny of Bupleurum was investigated based on the ITS region of the nuclear ribosomal DNA (nrDNA) of 14 taxa from the Hengduan Mountains, 3 taxa from the North China (Hebei and Heilongjiang), and 16 taxa from Africa and the Mediterranean region. Varia-tions in chromosome numbers and the ITS sequences were used to infer phylogenetic relationships between Bupleurum species in Hengduan Mountains. The results showed that the Hengduan Mountains might represent one of the frequency and diversity centers for Bupleurum. The ancestors of Bupleurum species in the Hengduan Mountains may be related to the woody B. fruticosum in North Africa, or the species in the western Mediterranean region. It is postulated that the ancestral population migrated into Hengduan Mountains through the Middle East and the Caucasus. Furthermore, the neo-endemic B. mundtii in South Africa appeared to be a close relative of the species in the Hengduan Mountains. In the trend of basic chromosome number evolution, x = 8 should be regarded as the ancestral basic number, while x = 6, 7 as the derived ones. The Bupleurum species in the Hengduan Moun-tains have been undergoing changes in the basic chromosome numbers or the ploidy level. The ITS phylogenetic tree showed that the Chinese species were divided into two clades: one with the basic chromosome number x = 8, and the other with x = 6, 7. The results rejected the previous infrageneric classification of Bupleurum in China. We further suggested to raise B. marginatum var. stenophyllum to species rank based on the combined evidence from morphology, karyology, pollen morphology, and the ITS phylogenetic tree.     

A worldwide polymorphism in aldehyde dehydrogenase in Drosophila melanogaster: evidence for selection mediated by dietary ethanol

Clinally varying traits in Drosophila melanogaster provide good opportunities for elucidating the genetic basis of adaptation. Resistance to ethanol, a natural component of D. melanogaster's breeding sites, increases with latitude on multiple continents, indicating that the trait is under selection. Although the well-studied Alcohol dehydrogenase (Adh) polymorphism makes a contribution to the clines, it accounts for only a small proportion of the phenotypic variation. We describe an amino acid replacement polymorphism in Aldehyde dehydrogenase (Aldh), the gene encoding the second enzyme in the ethanol degradation pathway, that shows hallmarks of also contributing to the clines. The derived Aldh allele, like the Adh-Fast allele, increases in frequency in laboratory populations selected for ethanol resistance, and increases in frequency with latitude in wild populations. Moreover, strains with the derived allele have significantly higher ALDH enzyme activity with acetaldehyde (the breakdown product of ethanol) as a substrate than strains with the ancestral allele. As is the case with the Adh-Fast allele, chromosomes with the derived Aldh allele show markedly reduced molecular variation in the vicinity of the replacement polymorphism compared to those with the ancestral allele, suggesting a single, relatively recent origin. Nonetheless, the Aldh polymorphism differs from the Adh polymorphism in that the ethanol-associated allele remains in relatively low frequency in most populations. We present evidence that this is likely to be the result of a trade-off in catalytic activity, with the advantage of the derived allele in acetaldehyde detoxification being offset by a disadvantage in detoxification of other aldehydes.     

Evidence consistent with human L1 retrotransposition in maternal meiosis I

We have used a unique polymorphic 3' transduction to show that a human L1, or LINE-1 (long interspersed nucleotide element-1), retrotransposition event most likely occurred in the maternal primary oocyte during meiosis I. We characterized a truncated L1 retrotransposon with a 3' transduction that was inserted, in a Dutch male patient, into the X-linked gene CYBB, thereby causing chronic granulomatous disease. We used the unique flanking sequence to localize the precursor L1 locus, LRE3, to chromosome 2q24.1. In a cell culture assay, the retrotransposition frequency of LRE3 is greater than that for any other element that has been tested to date. The patient's mother had two LRE3 alleles that differed slightly in the 3'-flanking genomic DNA. The patient had a single LRE3 allele that was identical to one of the maternal alleles; however, the patient's insertion matched the maternal LRE3 allele that he did not inherit. Other data indicate that there is only a small chance that the father (unavailable for analysis) carries the precursor LRE3 allele. In addition, paternal origin of the insertion would have required that an LRE3 mRNA transcribed before meiosis II be carried separately from its precursor LRE3 allele in the fertilizing sperm. Since the mother carries a potential precursor allele and the insertion was on the patient's maternal X chromosome, it is highly likely that the insertion originated during maternal meiosis I.     

Comparative genomics: tracking chromosome evolution in the family ursidae using reciprocal chromosome painting

The Ursidae family includes eight species, the karyotype of which diverges somewhat, in both chromosome number and morphology, from that of other families in the order Carnivora. The combination of consensus molecular phylogeny and high-resolution trypsin G-banded karyotype analysis has suggested that ancestral chromosomal fissions and at least two fusion events are associated with the development of the different ursid species. Here, we revisit this hypothesis by hybridizing reciprocal chromosome painting probes derived from the giant panda (Ailuropoda melanoleuca), domestic cat (Felis catus), and man (Homo sapiens) to representative bear species karyotypes. Comparative analysis of the different chromosome segment homologies allowed reconstruction of the genomic composition of a putative ancestral bear karyotype based upon the recognition of 39 chromosome segments defined by painting as the smallest conserved evolutionary unit segments (pSCEUS) among these species. The different pSCEUS combinations occurring among modern bear species support and extend the postulated sequence of chromosomal rearrangements and provide a framework to propose patterns of genome reorganization among carnivores and other mammal radiations.     

The pattern of phylogenomic evolution of the Canidae

Canidae species fall into two categories with respect to their chromosome composition: those with high numbered largely acrocentric karyotypes and others with a low numbered principally metacentric karyotype. Those species with low numbered metacentric karyotypes are derived from multiple independent fusions of chromosome segments found as acrocentric chromosomes in the high numbered species. Extensive chromosome homology is apparent among acrocentric chromosome arms within Canidae species; however, little chromosome arm homology exists between Canidae species and those from other Carnivore families. Here we use Zoo-FISH (fluorescent in situ hybridization, also called chromosomal painting) probes from flow-sorted chromosomes of the Japanese raccoon dog (Nyctereutes procyonoides) to examine two phylogenetically divergent canids, the arctic fox (Alopex lagopus) and the crab-eating fox (Cerdocyon thous). The results affirm intra-canid chromosome homologies, also implicated by G-banding. In addition, painting probes from domestic cat (Felis catus), representative of the ancestral carnivore karyotype (ACK), and giant panda (Ailuropoda melanoleuca) were used to define primitive homologous segments apparent between canids and other carnivore families. Canid chromosomes seem unique among carnivores in that many canid chromosome arms are mosaics of two to four homology segments of the ACK chromosome arms. The mosaic pattern apparently preceded the divergence of modern canid species since conserved homology segments among different canid species are common, even though those segments are rearranged relative to the ancestral carnivore genome arrangement. The results indicate an ancestral episode of extensive centric fission leading to an ancestral canid genome organization that was subsequently reorganized by multiple chromosome fusion events in some but not all Canidae lineages.     

"Bar-coding" primate chromosomes: molecular cytogenetic screening for the ancestral hominoid karyotype

Two recently introduced multicolor FISH approaches, cross-species color banding (also termed Rx-FISH) and multiplex FISH using painting probes derived from somatic cell hybrids retaining fragments of human chromosomes, were applied in a comparative molecular cytogenetic study of higher primates. We analyzed these "chromosome bar code" patterns to obtain an overview of chromosomal rearrangements that occurred during higher primate evolution. The objective was to reconstruct the ancestral genome organization of hominoids using the macaque as outgroup species. Approximately 160 individual and discernible molecular cytogenetic markers were assigned in these species. Resulting comparative maps allowed us to identify numerous intra-chromosomal rearrangements, to discriminate them from previous contradicting chromosome banding interpretations and to propose an ancestral karyotype for hominoids. From 25 different chromosome forms in an ancestral karyotype for all hominoids of 2N=48 we propose 21. Probes for chromosomes 2p, 4, 9 and Y were not informative in the present experiments. The orangutan karyotype was very similar to the proposed ancestral organization and conserved 19 of the 21 ancestral forms; thus most chromosomes were already present in early hominoid evolution, while African apes and human show various derived changes.     

A New Test for Selection Applied to Codon Usage in Drosophila simulans and D. mauritiana

In many organisms, synonymous codon usage is biased by a history of natural selection. However, codon bias, itself, does not indicate that selection is ongoing; it may be a vestige of past selection. Simple statistical tests have been devised to infer ongoing selection on codon usage by comparing the derived state frequency spectra at polymorphic sites segregating either derived preferred codons or derived unpreferred codons; if selection is effective, the frequency of derived states should be higher in the former. We propose a new test that uses the inferred degree of preference, essentially calculating the correlation of derived state frequency and the difference in preference between the derived and the ancestral states; the correlation should be positive if selection is effective. When implementing the test, derived and ancestral states can be assigned by parsimony or on the basis of relative probability. In either case, statistical significance is estimated by a simple permutation test. We explored the statistical power of the test by sampling polymorphism data from 14 loci in 16 strains of D. simulans, finding that the test retains 80% power even when quite a few of the data are discarded. The power of the test likely reflects better use of multiple features of the data, combining population frequencies of polymorphic variants and quantitative estimates of codon preferences. We also applied this novel test to 14 newly sequenced loci in five strains of D. mauritiana, showing for the first time ongoing selection on codon usage in this species.     

Microsatellite variation and evolutionary history of PCDHX/Y gene pair within the Xq21.3/Yp11.2 hominid-specific homology block

To better understand the evolutionary dynamics of repetitive sequences in human sex chromosomes, we have analyzed seven new X/Y homologous microsatellites located within PCDHX/Y, one of the two recently described gene pairs in the Xq21.3/Yp11.2 hominid-specific homology block, in samples from Portugal and Mozambique. Sharp differences were observed on X/Y allele distributions, concerning both the presence of private alleles and a different modal repeat length for X-linked and Y-linked markers, and this difference was statistically significant. Higher diversity was found in X-linked microsatellites than in their Y chromosome counterparts; when comparing populations, Mozambicans showed more allele diversity for the X chromosome, but the contrary was true for the Y chromosome microsatellites. Evolutionary patterns, relying on intragenic PCDHX/Y SNPs, also revealed distinct scenarios for X and Y chromosomes. Greater microsatellite diversity was displayed by African X chromosomes within the most common haplotypes shared by both populations, whereas higher microsatellite diversity was found in Portugal for the ancestral Y chromosome haplotype. The most frequent PCDHY haplotype in Portuguese was the derived one, and it was not found in Mozambicans. TMRCA estimated by the rho parameter resulted in 13,700 years (7,500-20,000 years), which is consistent with a recent, post-Out-of-Africa origin for this haplotype. In conclusion, the newly described microsatellite loci generally displayed greater X-linked to Y-linked diversity and this pattern was also detected with slower evolving markers, with a remarkable differentiation between populations observed for Y chromosome haplotypes and, thus, greater divergence among Y chromosomes in human populations.     

Sexual Antagonism,Temporally Fluctuating Selection,and Variable Dominance Affect a Regulatory Polymorphism in Drosophila melanogaster

Understanding how genetic variation is maintained within species is a major goal of evolutionary genetics that can shed light on the preservation of biodiversity. Here, we examined the maintenance of a regulatory single-nucleotide polymorphism (SNP) of the X-linked Drosophila melanogaster gene fezzik. The derived variant at this site is at intermediate frequency in many worldwide populations but absent in populations from the ancestral species range in sub-Saharan Africa. We collected and genotyped wild-caught individuals from a single European population biannually over a period of 5 years, which revealed an overall difference in allele frequency between the sexes and a consistent change in allele frequency across seasons in females but not in males. Modeling based on the observed allele and genotype frequencies suggested that both sexually antagonistic and temporally fluctuating selection may help maintain variation at this site. The derived variant is predicted to be female-beneficial and mostly recessive; however, there was uncertainty surrounding our dominance estimates and long-term modeling projections suggest that it is more likely to be dominant. By examining gene expression phenotypes, we found that phenotypic dominance was variable and dependent upon developmental stage and genetic background, suggesting that dominance may be variable at this locus. We further determined that fezzik expression and genotype are associated with starvation resistance in a sex-dependent manner, suggesting a potential phenotypic target of selection. By characterizing the mechanisms of selection acting on this SNP, our results improve our understanding of how selection maintains genetic and phenotypic variation in natural populations.     

Speciation and evolutionary dynamics of asymmetric mating preference

Asymmetric mating preferences occur in two closely related species, if females of one species are highly selective against males of the second, while females of the second show less selection against males of the first species. It has been suggested that such asymmetry is an indicator of common ancestry between the two species, but actual observations are contradictory and inconclusive. We developed a scenario of speciation history and asymmetric mating preference, incorporating invasion dynamicsvia frequency-dependent interspecific sexual competition. A newly isolated (derived) species may form at the periphery of the ancestral species’ distribution by invading a new range. Only a few closely related species would be expected in the new area, while many related species are expected to coexist with the ancestral species. In a peripherally derived species, female mating preferences should be relaxed through sexual character release, owing to a lack of sympatric species and a scarcity of intraspecific mating opportunities. Secondary contacts may then happen as: 1. repeated invasions, i.e. subsequent invasion by the ancestral species into the new range or, 2. backward invasions, i.e. derived species incursions into the ancestral range. Repeated invasions could lead to the coexistence of both the derived species and the newly invading ancestor. Backward invasions by the derived species can succeed only when the derived females develop a strict mating discrimination against the ancestral males. We then expect strong character displacement in the derived species. Thus, peripheral isolation and repeated invasions lead to the relaxed female mating preferences in the derived species and backward invasions lead to stronger female mating preferences in the derived species. This agrees withDrosophila data from Hawaii and the continents. Experimental data of theDrosophila arizonaemojavensis species cluster also support the hypothesis.     

Chromosomal homologies between Cebus and Ateles (primates) based on ZOO-FISH and G-banding comparisons

ZOO-FISH (Fluorescent "in vitro" hybridization) was used to establish the chromosomal homology between humans (HSA) and Cebus nigrivitatus (CNI) and Ateles belzebuth hybridus (ABH). These two species belong to different New World monkey families (Cebidae and Atelidae, respectively) which differ greatly in chromosome number and in chromosome morphology. The molecular results were followed by a detailed banding analysis. The ancestral karyotype of Cebus was then determined by a comparison of in situ hybridization results, as well as chromosomal morphology and banding in other Platyrrhini species. The karyotypes of the four species belonging to the genus Cebus differ from each other by three inversions and one fusion as well as in the location and amounts of heterochromatin. Results obtained by ZOO-FISH in ABH are in general agreement with previous gene-mapping and in situ hybridization data in Ateles, which show that spider monkeys have highly derived genomes. The chromosomal rearrangements detected between HSA and ABH on a band-to-band basis were 27 fusions/fissions, 12 centromeric shifts, and six pericentric inversions. The ancestral karyotype of Cebus was then compared with that of Ateles. The rearrangements detected were 20 fusions/fissions, nine centromeric shifts, and five inversions. Atelidae species are linked by a fragmentation of chromosome 4 into three segments forming an association of 4/15, while Ateles species are linked by 13 derived associations. The results also helped clarify the content of the ancestral platyrrhine karyotype and the mode of chromosomal evolution in these primates. In particular, associations 2/16 and 5/7 should be included in the ancestral karyotype of New World monkeys.     

Origin of agouti-melanistic polymorphism in Wild Black Rats (Rattus rattus) inferred from Mc1r gene sequences

We examined nucleotide changes that underlie coat color variation in Black Rats (the Rattus rattus species complex), which show polymorphism in dorsal fur color, including either grayish brown (agouti) or black (melanistic) forms. We examined the full coding sequence of a gene known to produce melanism in other vertebrates-melanocortin-1-receptor gene Mc1r (954 bp) -using samples of both R. rattus (with 2n = 38) and its close relative Asian Black Rat (R. tanezumi; 2n = 42). We used 61 specimens from Japan with karyotype-known individuals and four samples from Pakistan. We found 11 allele sequences and constructed a network tree that shows two distinct clusters, with allelic segregation according to karyotype and by inference, representing the two species. We found that a nucleotide substitution from G to A at site 280, producing an amino acid change from glutamic acid to lysine, was associated with the dominant trait of the melanistic form of the coat color in R. rattus. Notably, the derived SNP 280A was found in a single allele, with the ancestral SNP 280G present in seven alleles. By contrast, all three alleles for R. tanezumi retain the ancestral SNP 280G. These results suggest a possible recent origin of melanism in R. rattus.     

Patterns of replication in the neo-sex chromosomes ofDrosophila nasuta albomicans

Drosophila nasuta albomicans (with 2n = 6), contains a pair of metacentric neo-sex chromosomes. Phylogenetically these are products of centric fusion between ancestral sex (X, Y) chromosomes and an autosome (chromosome 3). The polytene chromosome complement of males with a neo-X- and neo-Y-chromosomes has revealed asynchrony in replication between the two arms of the neo-sex chromosomes. The arm which represents the ancestral X-chromosome is faster replicating than the arm which represents ancestral autosome. The latter arm of the neo-sex chromosome is synchronous with other autosomes of the complement. We conclude that one arm of the neo-X/Y is still mimicking the features of an autosome while the other arm has the features of a classical X/Y-chromosome. This X-autosome translocation differs from the other evolutionary X-autosome translocations known in certain species ofDrosophila.     

A Bayesian method for jointly estimating allele age and selection intensity.

The problem of jointly estimating the intensity of past selection affecting an allele and the allele's age is formulated in a Bayesian framework. The prior distribution of allele age given its frequency is obtained from existing population genetics theory. The prior distribution of selection intensity is assumed to reflect the fact that positive selection on a new mutant is more likely to be weak than strong. The general approach is illustrated by the development of an importance sampling method applicable to low-frequency alleles. This method can be used either when the haplotypes of closely linked marker loci are known or when the lengths of linked ancestral chromosomal segments can be inferred. The method is illustrated with an application to the A-allele of G6PD in Africa. Because changes in allele frequency and recombination are both intrinsically stochastic, there are limits to the accuracy achievable with any method.