Evolution of zinc finger-Y and zinc finger-X genes in oryzomyne-akodontine rodents (Cricetidae)

Summary Zinc finger-Y (Zfy) and zinc finger-X (Zfx) genes were analyzed by Southern blotting in male and female specimens of 10 species belonging to the oryzomyne-akodontine stock of Cricetidae rodents. DNA fragments were used as characters to construct a parsimony tree of the genes. Zfx and Zfy trees in general coincide with the evolutionary history of the taxa. Both trees show Oryzomys longicaudatus genes as the outgroup whereas Akodon xanthorrhinus genes are also distant from those of the other species. Oxymycterus rufus and Bolomys obscurus share related sequences, while genes from the other six Akodon species form a group of their own. It was found that 9 out of the 10 species analyzed show Zfy amplification in a range varying from 2 to 24 copies and with a pattern that is clade specific. The estimation of the average changes per character strongly suggests that Zfy has evolved more rapidly than Zfx; our estimates of the rate of nucleotide substitution are 4.6 times higher for Zfy than for Zfx. Offprint requests to: N.O. Bianchi     

Retroposon compensatory mechanism hypothesis not supported: Zfa knockout mice are fertile

It is hypothesized that autosomal retroposons compensate for the loss of their inactivated essential X-chromosome progenitors during spermatogenesis. Here we test this Retroposon Compensatory Mechanism (RCM) hypothesis using the Zfy gene family. The mouse autosomal retroposon Zfa is expressed in testes at the same developmental time points at which Zfx levels decline, which correspond to the time of male sex chromosome inactivation, suggesting that Zfa may compensate for the loss of Zfx during spermatogenesis. We examined the effect of Zfa-targeted mutagenesis on spermatogenesis in three genetically distinct mouse strains. Surprisingly, Zfa knockout mice showed no detectable fertility, sperm count, or testes morphology defects. We therefore conclude that Zfa is not an essential gene for spermatogenesis and fertility. This surprising finding now challenges the RCM hypothesis at least for the Zfy gene family. It also forces us to reevaluate the original data underpinning the RCM hypothesis for this family and to propose alternative hypotheses.     

Likelihood ratio tests for detecting positive selection and application to primate lysozyme evolution

An excess of nonsynonymous substitutions over synonymous ones is an important indicator of positive selection at the molecular level. A lineage that underwent Darwinian selection may have a nonsynonymous/synonymous rate ratio (dN/dS) that is different from those of other lineages or greater than one. In this paper, several codon-based likelihood models that allow for variable dN/dS ratios among lineages were developed. They were then used to construct likelihood ratio tests to examine whether the dN/dS ratio is variable among evolutionary lineages, whether the ratio for a few lineages of interest is different from the background ratio for other lineages in the phylogeny, and whether the dN/dS ratio for the lineages of interest is greater than one. The tests were applied to the lysozyme genes of 24 primate species. The dN/dS ratios were found to differ significantly among lineages, indicating that the evolution of primate lysozymes is episodic, which is incompatible with the neutral theory. Maximum- likelihood estimates of parameters suggested that about nine nonsynonymous and zero synonymous nucleotide substitutions occurred in the lineage leading to hominoids, and the dN/dS ratio for that lineage is significantly greater than one. The corresponding estimates for the lineage ancestral to colobine monkeys were nine and one, and the dN/dS ratio for the lineage is not significantly greater than one, although it is significantly higher than the background ratio. The likelihood analysis thus confirmed most, but not all, conclusions Messier and Stewart reached using reconstructed ancestral sequences to estimate synonymous and nonsynonymous rates for different lineages.      

Zfa is an expressed retroposon derived from an alternative transcript of the Zfx gene.

ZFY, a gene on the Y chromosome encoding a zinc finger protein, has been proposed as a candidate for the human testis determining gene. Sequences related to ZFY, called ZFX, are present on the X chromosome of a wide range of placental mammals. Unlike most mammals the mouse has four genes homologous to ZFY; two on the Y chromosome, Zfy-1 and Zfy-2, an X-linked gene, Zfx, and an autosomal gene, Zfa. We show here that Zfa has arisen recently by retroposition of one of at least three alternatively spliced mRNAs transcribed from the Zfx gene. Zfa is an unusual retroposon in that it has retained an open reading frame and is expressed, although its function may be limited or altered by the presence of a potentially inactivating mutation in the third of its zinc fingers. This mutation must have occurred at the same time or soon after the retroposition event as it is also present in the Zfa gene of Mus spretus. Interestingly the third finger of the M. musculus musculus Zfy-2 gene has also sustained a mutation suggesting that this gene family may be rapidly evolving in mice.     

Localization of Murine X and Autosomal Sequences Homologous to the Human Y Located Testis-Determining Region

Recently a candidate gene for the primary testis-determining factor (TDF) encoding a zinc finger protein (ZFY) has been cloned from the human Y chromosome. A highly homologous X-linked copy has also been identified. Using this human sequence it is possible to identify two Y loci, an X and an autosomal locus in the mouse (Zfy-1, Zfy-2, Zfx and Zfa, respectively). Suprisingly ZFY is more homologous to the mouse X and autosomal sequences than it is to either of the Y-linked loci. Both Zfy-1 and Zfy-2 are present in the Sxr region of the Y but Zfy-2 is absent in the Sxr deletion variant Sxrb (or Sxr") suggesting it is not necessary for male determination. Extensive backcross analyses map Zfa to mouse chromosome 10 and Zfx to a 5-cM interval between anonymous X probe MDXS120 and the tabby locus (Ta). We also show that the mouse androgen receptor locus (m-AR) believed to underlie the testicular feminization mutation (Tfm) shows complete linkage to Zfx. Comparative mapping indicates that in man these genes lie in separate conserved DNA segments.     

Preponderance of slightly deleterious polymorphism in mitochondrial DNA: nonsynonymous/synonymous rate ratio is much higher within species than between species

We estimated synonymous (dN) and nonsynonymous (dS) substitution rates for protein-coding genes of the mitochondrial genome from two individuals each of the species human, chimpanzee, and gorilla. The genes were analyzed both separately and in a combined data set. Pairwise sequence comparisons suggest that the dN/dS rate ratios are about 5-10 times higher in within-species comparisons than in between-species comparisons. This result is confirmed by a more rigorous likelihood ratio test, which rejected the null hypothesis that the dN/dS rate ratios are identical within and between species. The likelihood models account for the genetic code structure, transition/transversion rate ratio, and codon usage bias and are expected to produce more reliable results than the commonly used contingency test. Separate analyses of different genes show that the dN/dS rate ratios are higher within species than between species for all 13 mitochondrial genes, with the difference being statistically significant for all except three small or slowly evolving genes. Furthermore, in conserved genes, nonsynonymous rates within species tend to be higher than the between-species rates by a greater proportion than in fast-changing genes. Our findings confirm and extend earlier results obtained from smaller data sets and suggest the operation of slightly deleterious mutations throughout the mitochondrial genome in the hominoids. Implications of the results for evolutionary studies and, in particular, for studies of the origin of modern humans, are discussed.     

Lineage‐specific sequence evolution and exon edge conservation partially explain the relationship between evolutionary rate and expression level in A. thaliana

Rapidly evolving proteins can aid the identification of genes underlying phenotypic adaptation across taxa, but functional and structural elements of genes can also affect evolutionary rates. In plants, the ‘edges’ of exons, flanking intron junctions, are known to contain splice enhancers and to have a higher degree of conservation compared to the remainder of the coding region. However, the extent to which these regions may be masking indicators of positive selection or account for the relationship between dN/dS and other genomic parameters is unclear. We investigate the effects of exon edge conservation on the relationship of dN/dS to various sequence characteristics and gene expression parameters in the model plant Arabidopsis thaliana. We also obtain lineage‐specific dN/dS estimates, making use of the recently sequenced genome of Thellungiella parvula, the second closest sequenced relative after the sister species Arabidopsis lyrata. Overall, we find that the effect of exon edge conservation, as well as the use of lineage‐specific substitution estimates, upon dN/dS ratios partly explains the relationship between the rates of protein evolution and expression level. Furthermore, the removal of exon edges shifts dN/dS estimates upwards, increasing the proportion of genes potentially under adaptive selection. We conclude that lineage‐specific substitutions and exon edge conservation have an important effect on dN/dS ratios and should be considered when assessing their relationship with other genomic parameters.     

Regional assignments of the zinc finger Y-linked gene (ZFY) and related sequences on human and mouse chromosomes

Recent chromosome walking experiments have identified a candidate gene (ZFY) for the testis-determining factor on the human Y chromosome (Page et al., 1987). We report here the regional assignments of the ZFY gene and related sequences in the human and the mouse. By in situ hybridization, we assigned ZFX and ZFY to human chromosome bands Xp21 and Yp11.3, respectively. Although the mouse harbors two Zfy genes, only one site at band A1 of its Y chromosome was significantly labeled. The mouse Zfx gene and the Zfa gene on chromosome 10 were assigned to bands XD and 10B5, respectively. These assignments of the ZFX gene in human and mouse add another marker to the conserved syntenic group for evaluating the evolutionary relationship of the human and mouse X chromosomes.     

Phylogenetic relationships of the Cobitoidea (Teleostei: Cypriniformes) inferred from mitochondrial and nuclear genes with analyses of gene evolution

The superfamily Cobitoidea of the order Cypriniformes is a diverse group of fishes, inhabiting freshwater ecosystems across Eurasia and North Africa. The phylogenetic relationships of this well-corroborated natural group and diverse clade are critical to not only informing scientific communities of the phylogeny of the order Cypriniformes, the world's largest freshwater fish order, but are key to every area of comparative biology examining the evolution of traits, functional structures, and breeding behaviors to their biogeographic histories, speciation, anagenetic divergence, and divergence time estimates. In the present study, two mitochondrial gene sequences (COI, ND4+5) and four single-copy nuclear gene segments (RH1, RAG1, EGR2B, IRBP) were used to infer the phylogenetic relationships of the Cobitoidea as reconstructed from maximum likelihood (ML) and partitioned Bayesian Analysis (BA). Analyses of the combined mitochondrial/nuclear gene datasets revealed five strongly supported monophyletic Cobitoidea families and their sister-group relationships: Botiidae+(Vaillantellidae+(Cobitidae+(Nemacheilidae+Balitoridae))). These recovered relationships are in agreement with previous systematic studies on the order Cypriniformes and/or those focusing on the superfamily Cobitoidea. Using these relationships, our analyses revealed pattern lineage- or ecological-group-specific evolution of these genes for the Cobitoidea. These observations and results corroborate the hypothesis that these group-specific-ancestral ecological characters have contributed in the diversification and/or adaptations within these groups. Positive selections were detected in RH1 of nemacheilids and in RAG1 of nemacheilids and genus Vaillantella, which indicated that evolution of RH1 (related to eye's optic sense) and RAG1 (related to immunity) genes appeared to be important for the diversification of these groups. The balitorid lineage (those species inhabiting fast-flowing riverine habitats) had, as compared with other cobitoid lineages, significantly different dN/dS, dN and dS values for ND4 and IRBP genes. These significant differences are usually indicative of weaker selection pressure, and lineage-specific evolution on genes along the balitorid lineage. Furthermore, within Cobitoidea, excluding balitorids, species living in subtropics had significantly higher dN/dS values in RAG1 and IRBP genes than those living in temperate and tropical zones. Among tropical cobitoids, genes COI, ND5, EGR2B, IRBP and RH1, had a significantly higher mean dS value than those species in subtropical and temperate groups. These findings suggest that the evolution of these genes could also be ecological-group-specific and may have played an important role in the adaptive evolution and diversification of these groups. Thus, we hypothesize that the genes included in the present study were actively involved in lineage- and/or ecological-group-specific evolutionary processes of the highly diverse Cobitoidea. These two evolutionary patterns, both subject to further testing, are hypothesized as integral in the diversification with this major clade of the world's most diverse group of freshwater fishes.     

Patterns of Y and X chromosome DNA sequence divergence during the Felidae radiation.

The 37 species of modern cats have evolved from approximately eight phylogenetic lineages within the past 10 to 15 million years. The Felidae family has been described with multiple measures of morphologic and molecular evolutionary methods that serve as a framework for tracking gene divergence during brief evolutionary periods. In this report, we compare the mode and tempo of evolution of noncoding sequences of a large intron within Zfy (783 bp) and Zfx (854 bp), homologous genes located on the felid Y and X chromosomes, respectively. Zfy sequence variation evolves at about twice the rate of Zfx, and both gene intron sequences track feline hierarchical topologies accurately. As homoplasies are infrequent in patterns of nucleotide substitution, the Y chromosome sequence displays a remarkable degree of phylogenetic consistency among cat species and provides a highly informative glimpse of divergence of sex chromosome sequences in Felidae.     

ZFX基因同源序列在黄鳝基因组中的检出及其染色体定位

以大熊猫锌指蛋白基因Zfx为探针 ,在黄鳝基因组DNA中检测到一条长约 9 5kb的杂交带。依据哺乳类和爬行类动物锌指蛋白基因 (ZFX/Zfc)编码第 7～ 13个锌指结构的DNA序列保守性设计引物 ,在黄鳝基因组DNA中仅扩增到一条 5 12bp的DNA片段。将此片段克隆至载体 pBS中 ,从雌性、雄性个体中分别挑选 4个含有插入片段的白色克隆进行测序。测序结果表明 ,这些克隆中插入片段的核苷酸序列一致。该DNA片段在核苷酸水平上与人类ZFX和ZFY分别具有 88%和 87%同源性 ,但其与美洲鳄鱼Zfc的同源性可达 90 % ,而在氨基酸水平上则分别存在 95 9%、95 9%和 93 5 %的同源性 (170个氨基酸 )。该基因命名为黄鳝锌指蛋白基因Zfa ,并运用FISH将其定位于黄鳝 1号染色体 ,距离着丝粒的相对位置为 6 0 1± 0 38。通过进一步研究证明 ,黄鳝 1号染色体上存在有真兽类哺乳动物X染色质同源的保守片段 ,该保守片段有可能就是哺乳动物X染色体起源和进化的原始物质基础之一。应用哺乳动物X染色体连锁的其他基因在鱼类开展染色体比较定位研究 ,将有望促进脊椎动物性染色体进化的深入研究     

Independent effects of alternative splicing and structural constraint on the evolution of mammalian coding exons

Alternative splicing (AS) is known to significantly affect exon-level protein evolutionary rates in mammals. Particularly, alternatively spliced exons (ASEs) have a higher nonsynonymous-to-synonymous substitution rate (dN/dS) ratio than constitutively spliced exons (CSEs), possibly because the former are required only occasionally for normal biological functions. Meanwhile, intrinsically disordered regions (IDRs), the protein regions lacking fixed 3D structures, are also reported to have an increased evolutionary rate due to lack of structural constraint. Interestingly, IDRs tend to be located in alternative protein regions. Yet which of these two factors is the major determinant of the increased dN/dS in mammalian ASEs remains unclear. By comparing human-macaque and human-mouse one-to-one orthologous genes, we demonstrate that AS and protein structural disorder have independent effects on mammalian exon evolution. We performed analyses of covariance to demonstrate that the slopes of the (dN/dS-percentage of IDR) regression lines differ significantly between CSEs and ASEs. In other words, the dN/dS ratios of both ASEs and CSEs increase with the proportion of IDR (PIDR), whereas ASEs have higher dN/dS ratios than CSEs when they have similar PIDRs. Since ASEs and IDRs may less frequently overlap with protein domains (which also affect dN/dS), we also examined the correlations between dN/dS ratio and exon type/PIDR by controlling for the density of protein domain. We found that the effects of exon type and PIDR on dN/dS are both independent of domain density. Our results imply that nature can select for different biological features with regard to ASEs and IDRs, even though the two biological features tend to be localized in the same protein regions.     

The molecular evolution of ZFY-related genes in birds and mammals

Summary We report the isolation and nucleotide sequence determination of clones derived from five ZFY-related zinc-finger genes from birds and mammals. These sequences are analyzed with reference to the previously published human genes, ZFX and ZFY, and mouse genes, Zfx, Zfa, Zfy-1, and Zfy-2. The analysis indicates that ZFY-related genes are highly conserved in birds and mammals, and that the rate of nucleotide substitution in the Y-linked genes is not as high as predicted. However, the mouse Zfy-1 and Zfy-2 genes are markedly divergent members of the ZFY gene family; we suggest this relates to X-inactivation of the mouse gene Zfx.     

Maximum likelihood estimation of ancestral codon usage bias parameters in Drosophila

We present a likelihood method for estimating codon usage bias parameters along the lineages of a phylogeny. The method is an extension of the classical codon-based models used for estimating dN/dS ratios along the lineages of a phylogeny. However, we add one extra parameter for each lineage: the selection coefficient for optimal codon usage (S), allowing joint maximum likelihood estimation of S and the dN/dS ratio. We apply the method to previously published data from Drosophila melanogaster, Drosophila simulans, and Drosophila yakuba and show, in accordance with previous results, that the D. melanogaster lineage has experienced a reduction in the selection for optimal codon usage. However, the D. melanogaster lineage has also experienced a change in the biological mutation rates relative to D. simulans, in particular, a relative reduction in the mutation rate from A to G and an increase in the mutation rate from C to T. However, neither a reduction in the strength of selection nor a change in the mutational pattern can alone explain all of the data observed in the D. melanogaster lineage. For example, we also confirm previous results showing that the Notch locus has experienced positive selection for previously classified unpreferred mutations.     

Selection on protein-coding genes of natural cyanobacterial populations

We examined the distribution of synonymous and non-synonymous changes in 12 protein-coding genes of natural populations of cyanobacteria to infer changes in gene functionality. By comparing mutation distributions within and across species using the McDonald-Kreitman test, we found data sets to contain evidence for purifying selection (hetR of Trichodesmium, nifH of Cylindrospermopsis raceborskii and rpoC1 of Anabaena lemmermannii) and positive selection (kaiC of Microcoleus chthonoplastes and rbcX of Anabaena and Aphanizomenon sp.). Other genes from the same set of clonal isolates (petB and rbcL in M. chthonoplastes and Anabaena/Aphanizomenon, respectively) did not harbour evidence for either form of selection. The results of branch models of codon evolution agreed fully with the results of the McDonald-Kreitman test in terms of significance and absolute value of the dN/dS estimates. The high frequency of gene-specific mutation patterns and their association with branches that separate closely related cyanobacterial genera suggest that evolutionary tests are suited to uncover gene-specific selective differentiation in cyanobacterial genomes. At the same time, given the lack of information about the history of cyanobacteria, analysis of larger numbers of protein-coding genes of clonal cyanobacterial isolates will produce more detailed pictures of the effects of natural selection.     

Dolphin genome provides evidence for adaptive evolution of nervous system genes and a molecular rate slowdown

Cetaceans (dolphins and whales) have undergone a radical transformation from the original mammalian bodyplan. In addition, some cetaceans have evolved large brains and complex cognitive capacities. We compared approximately 10 000 protein-coding genes culled from the bottlenose dolphin genome with nine other genomes to reveal molecular correlates of the remarkable phenotypic features of these aquatic mammals. Evolutionary analyses demonstrated that the overall synonymous substitution rate in dolphins has slowed compared with other studied mammals, and is within the range of primates and elephants. We also discovered 228 genes potentially under positive selection (dN/dS > 1) in the dolphin lineage. Twenty-seven of these genes are associated with the nervous system, including those related to human intellectual disabilities, synaptic plasticity and sleep. In addition, genes expressed in the mitochondrion have a significantly higher mean dN/dS ratio in the dolphin lineage than others examined, indicating evolution in energy metabolism. We encountered selection in other genes potentially related to cetacean adaptations such as glucose and lipid metabolism, dermal and lung development, and the cardiovascular system. This study underlines the parallel molecular trajectory of cetaceans with other mammalian groups possessing large brains.     

Mitochondrial DNA as a tool for reconstructing past life‐history traits in mammals

Reconstructing the ancestral characteristics of species is a major goal in evolutionary and comparative biology. Unfortunately, fossils are not always available and sufficiently informative, and phylogenetic methods based on models of character evolution can be unsatisfactory. Genomic data offer a new opportunity to estimate ancestral character states, through (i) the correlation between DNA evolutionary processes and species life‐history traits and (ii) available reliable methods for ancestral sequence inference. Here, we assess the relevance of mitochondrial DNA – the most popular molecular marker in animals – as a predictor of ancestral life‐history traits in mammals, using the order of Cetartiodactyla as a benchmark. Using the complete set of 13 mitochondrial protein‐coding genes, we show that the lineage‐specific nonsynonymous over synonymous substitution rate ratio (dN/dS) is closely correlated with the species body mass, longevity and age of sexual maturity in Cetartiodactyla and can be used as a marker of ancestral traits provided that the noise introduced by short branches is appropriately dealt with. Based on ancestral dN/dS estimates, we predict that the first cetartiodactyls were relatively small animals (around 20 kg). This finding is in accordance with Cope's rule and the fossil record but could not be recovered via continuous character evolution methods.     

Differentiation of Xylella fastidiosa strains via multilocus sequence analysis of environmentally mediated genes (MLSA-E)

Isolates of the plant pathogen Xylella fastidiosa are genetically very similar, but studies on their biological traits have indicated differences in virulence and infection symptomatology. Taxonomic analyses have identified several subspecies, and phylogenetic analyses of housekeeping genes have shown broad host-based genetic differences; however, results are still inconclusive for genetic differentiation of isolates within subspecies. This study employs multilocus sequence analysis of environmentally mediated genes (MLSA-E; genes influenced by environmental factors) to investigate X. fastidiosa relationships and differentiate isolates with low genetic variability. Potential environmentally mediated genes, including host colonization and survival genes related to infection establishment, were identified a priori. The ratio of the rate of nonsynonymous substitutions to the rate of synonymous substitutions (dN/dS) was calculated to select genes that may be under increased positive selection compared to previously studied housekeeping genes. Nine genes were sequenced from 54 X. fastidiosa isolates infecting different host plants across the United States. Results of maximum likelihood (ML) and Bayesian phylogenetic (BP) analyses are in agreement with known X. fastidiosa subspecies clades but show novel within-subspecies differentiation, including geographic differentiation, and provide additional information regarding host-based isolate variation and specificity. dN/dS ratios of environmentally mediated genes, though <1 due to high sequence similarity, are significantly greater than housekeeping gene dN/dS ratios and correlate with increased sequence variability. MLSA-E can more precisely resolve relationships between closely related bacterial strains with low genetic variability, such as X. fastidiosa isolates. Discovering the genetic relationships between X. fastidiosa isolates will provide new insights into the epidemiology of populations of X. fastidiosa, allowing improved disease management in economically important crops.