Search

Bronchogenic carcinomas arise from bronchial epithelial cells (BECs). Inhalation exposure of BECs to nitrosamines in cigarette smoke is an important exogenous risk factor for malignant transformation of BECs. Thus, an important endogenous risk factor is likely to be the capacity of BECs to metabolize nitrosamines. Among the cytochrome P450 enzymes capable of metabolizing nitrosamines, CYP2A6, CYP2E1 and CYP2B6 are expressed in BECs. In this study, we used quantitative RT-PCR to evaluate expression of CYP2A6 and CYP2E1 in primary human BECs from 12 non-smokers and eight smokers. CYP2A6 was expressed in 20/20 cases and quantifiable in 18/20 cases, with a mean level of 580 mRNA/106 β-actin mRNA. CYP2E1 expression was observed in 9/20 cases, but in all cases it was expressed at levels below our limit of quantification (10 mRNA/106 β-actin mRNA). There was significant (P < 0.05) 20 fold inter-individual variation in expression of CYP2A6. Further, the mean level of CYP2A6 among smokers (260 mRNA/106 β-actin mRNA) was significantly lower than among non-smokers (740 mRNA/106 β-actin mRNA). It is hypothesized that: (i) inter-individual variation in CYP2A6 gene expression may contribute to inter-individual variation in risk for bronchogenic carcinoma; (ii) smoking may reduce the level of expression of CYP2A6 in the BECs of some individuals; and (iii) CYP2A6 is more important than CYP2E1 for metabolic activation of nitrosamines in bronchial epithelial cells. [ABSTRACT FROM PUBLISHER]

Much of evolutionary theory is predicated on assumptions about the relative importance of simple additive versus complex epistatic genetic architectures. Previous work suggests traits strongly associated with fitness will lack additive genetic variation, whereas traits less strongly associated with fitness are expected to exhibit more additive genetic variation. We use a quantitative genetics method, line cross analysis, to infer genetic architectures that contribute to trait divergence. By parsing over 1,600 datasets by trait type, clade, and cross divergence, we estimated the relative importance of epistasis across the tree of life. In our comparison between life-history traits and morphological traits, we found greater epistatic contributions to life-history traits. Our comparison between plants and animals showed that animals have more epistatic contribution to trait divergence than plants. In our comparison of within-species versus between-species crosses, we found that only animals exhibit a greater epistatic contribution to trait divergence as divergence increases. While many scientists have argued that epistasis is ultimately of little importance, our results show that epistasis underlies much of trait divergence and must be accounted for in theory and practical applications like domestication, conservation breeding design, and understanding complex diseases., (© The Author(s) 2024. Published by Oxford University Press on behalf of The Society for the Study of Evolution (SSE).)

Chromosomal mutations such as fusions and fissions are often thought to be deleterious, especially in heterozygotes (underdominant), and consequently are unlikely to become fixed. Yet, many models of chromosomal speciation ascribe an important role to chromosomal mutations. When the effective population size (Ne) is small, the efficacy of selection is weakened, and the likelihood of fixing underdominant mutations by genetic drift is greater. Thus, it is possible that ecological and phenotypic transitions that modulate Ne facilitate the fixation of chromosome changes, increasing the rate of karyotype evolution. We synthesize all available chromosome number data in Coleoptera and estimate the impact of traits expected to change Ne on the rate of karyotype evolution in the family Carabidae and 12 disparate clades from across Coleoptera. Our analysis indicates that in Carabidae, wingless clades have faster rates of chromosome number increase. Additionally, our analysis indicates clades exhibiting multiple traits expected to reduce Ne, including strict inbreeding, oligophagy, winglessness, and island endemism, have high rates of karyotype evolution. Our results suggest that chromosome number changes are likely fixed by genetic drift despite an initial fitness cost and that chromosomal speciation models may be important to consider in clades with very small Ne., (© The American Genetic Association. 2024.)

Tribolium castaneum, the red flour beetle, is among the most well-studied eukaryotic genetic model organisms. Tribolium often serves as a comparative bridge from highly derived Drosophila traits to other organisms. Simultaneously, as a member of the most diverse order of metazoans, Coleoptera, Tribolium informs us about innovations that accompany hyper diversity. However, understanding the tempo and mode of evolutionary innovation requires well-resolved, time-calibrated phylogenies, which are not available for Tribolium. The most recent effort to understand Tribolium phylogenetics used two mitochondrial and three nuclear markers. The study concluded that the genus may be paraphyletic and reported a broad range for divergence time estimates. Here we employ recent advances in Bayesian methods to estimate the relationships and divergence times among Tribolium castaneum, T. brevicornis, T. confusum, T. freemani, and Gnatocerus cornutus using 1368 orthologs conserved across all five species and an independent substitution rate estimate. We find that the most basal split within Tribolium occurred ~86 Mya [95% HPD 85.90-87.04 Mya] and that the most recent split was between T. freemani and T. castaneum at ~14 Mya [95% HPD 13.55-14.00]. Our results are consistent with broader phylogenetic analyses of insects and suggest that Cenozoic climate changes played a role in the Tribolium diversification., (Copyright © 2021 Elsevier Inc. All rights reserved.)

Here we use a chromosome-level genome assembly of a prairie rattlesnake ( Crotalus viridis ), together with Hi-C, RNA-seq, and whole-genome resequencing data, to study key features of genome biology and evolution in reptiles. We identify the rattlesnake Z Chromosome, including the recombining pseudoautosomal region, and find evidence for partial dosage compensation driven by an evolutionary accumulation of a female-biased up-regulation mechanism. Comparative analyses with other amniotes provide new insight into the origins, structure, and function of reptile microchromosomes, which we demonstrate have markedly different structure and function compared to macrochromosomes. Snake microchromosomes are also enriched for venom genes, which we show have evolved through multiple tandem duplication events in multiple gene families. By overlaying chromatin structure information and gene expression data, we find evidence for venom gene-specific chromatin contact domains and identify how chromatin structure guides precise expression of multiple venom gene families. Further, we find evidence for venom gland-specific transcription factor activity and characterize a complement of mechanisms underlying venom production and regulation. Our findings reveal novel and fundamental features of reptile genome biology, provide insight into the regulation of snake venom, and broadly highlight the biological insight enabled by chromosome-level genome assemblies., (© 2019 Schield et al.; Published by Cold Spring Harbor Laboratory Press.)

Colubridae represents the most phenotypically diverse and speciose family of snakes, yet no well-assembled and annotated genome exists for this lineage. Here, we report and analyze the genome of the garter snake, Thamnophis sirtalis, a colubrid snake that is an important model species for research in evolutionary biology, physiology, genomics, behavior, and the evolution of toxin resistance. Using the garter snake genome, we show how snakes have evolved numerous adaptations for sensing and securing prey, and identify features of snake genome structure that provide insight into the evolution of amniote genomes. Analyses of the garter snake and other squamate reptile genomes highlight shifts in repeat element abundance and expansion within snakes, uncover evidence of genes under positive selection, and provide revised neutral substitution rate estimates for squamates. Our identification of Z and W sex chromosome-specific scaffolds provides evidence for multiple origins of sex chromosome systems in snakes and demonstrates the value of this genome for studying sex chromosome evolution. Analysis of gene duplication and loss in visual and olfactory gene families supports a dim-light ancestral condition in snakes and indicates that olfactory receptor repertoires underwent an expansion early in snake evolution. Additionally, we provide some of the first links between secreted venom proteins, the genes that encode them, and their evolutionary origins in a rear-fanged colubrid snake, together with new genomic insight into the coevolutionary arms race between garter snakes and highly toxic newt prey that led to toxin resistance in garter snakes.

Many taxa exhibit plastic immune responses initiated after primary microbial exposure that provide increased protection against disease-induced mortality and the fitness costs of infection. In several arthropod species, this protection can even be passed from parents to offspring through a phenomenon called trans-generational immune priming. Here, we first demonstrate that trans-generational priming is a repeatable phenomenon in flour beetles (Tribolium castaneum) primed and infected with Bacillus thuringiensis (Bt). We then quantify the within-host dynamics of microbes and host physiological responses in infected offspring from primed and unprimed mothers by monitoring bacterial density and using mRNA-seq to profile host gene expression, respectively, over the acute infection period. We find that priming increases inducible resistance against Bt around a critical temporal juncture where host septicaemic trajectories, and consequently survival, may be determined in unprimed individuals. Our results identify a highly differentially expressed biomarker of priming, containing an EIF4-e domain, in uninfected individuals, as well as several other candidate genes. Moreover, the induction and decay dynamics of gene expression over time suggest a metabolic shift in primed individuals. The identified bacterial and gene expression dynamics are likely to influence patterns of bacterial fitness and disease transmission in natural populations., (© 2017 John Wiley & Sons Ltd.)

Background: Relatively little is known about the genomic basis and evolution of wood-feeding in beetles. We undertook genome sequencing and annotation, gene expression assays, studies of plant cell wall degrading enzymes, and other functional and comparative studies of the Asian longhorned beetle, Anoplophora glabripennis, a globally significant invasive species capable of inflicting severe feeding damage on many important tree species. Complementary studies of genes encoding enzymes involved in digestion of woody plant tissues or detoxification of plant allelochemicals were undertaken with the genomes of 14 additional insects, including the newly sequenced emerald ash borer and bull-headed dung beetle., Results: The Asian longhorned beetle genome encodes a uniquely diverse arsenal of enzymes that can degrade the main polysaccharide networks in plant cell walls, detoxify plant allelochemicals, and otherwise facilitate feeding on woody plants. It has the metabolic plasticity needed to feed on diverse plant species, contributing to its highly invasive nature. Large expansions of chemosensory genes involved in the reception of pheromones and plant kairomones are consistent with the complexity of chemical cues it uses to find host plants and mates., Conclusions: Amplification and functional divergence of genes associated with specialized feeding on plants, including genes originally obtained via horizontal gene transfer from fungi and bacteria, contributed to the addition, expansion, and enhancement of the metabolic repertoire of the Asian longhorned beetle, certain other phytophagous beetles, and to a lesser degree, other phytophagous insects. Our results thus begin to establish a genomic basis for the evolutionary success of beetles on plants.

There are few patterns in evolution that are as rigidly held as Haldane's rule (HR), which states, "When in the first generation between hybrids between 2 species, 1 sex is absent, rare, or sterile, that sex is always the heterogametic sex." Yet despite considerable attention for almost a century, questions persist as to how many independent examples exist and what is (are) the underlying genetic cause(s). Here, we review recent evidence extending HR to plants, where previously it has only been documented in animals. We also discuss recent comparative analyses that show much more variation in sex-chromosome composition than previously recognized, thus increasing the number of potential independent origins of HR dramatically. Finally, we review the standing of genetic theories proposed to explain HR in light of the new examples and new molecular understanding., (© The American Genetic Association. 2016. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

The pace and direction of evolution in response to selection, drift, and mutation are governed by the genetic architecture that underlies trait variation. Consequently, much of evolutionary theory is predicated on assumptions about whether genes can be considered to act in isolation, or in the context of their genetic background. Evolutionary biologists have disagreed, sometimes heatedly, over which assumptions best describe evolution in nature. Methods for estimating genetic architectures that favor simpler (i.e., additive) models contribute to this debate. Here we address one important source of bias, model selection in line cross analysis (LCA). LCA estimates genetic parameters conditional on the best model chosen from a vast model space using relatively few line means. Current LCA approaches often favor simple models and ignore uncertainty in model choice. To address these issues we introduce Software for Analysis of Genetic Architecture (SAGA), which comprehensively assesses the potential model space, quantifies model selection uncertainty, and uses model weighted averaging to accurately estimate composite genetic effects. Using simulated data and previously published LCA studies, we demonstrate the utility of SAGA to more accurately define the components of complex genetic architectures, and show that traditional approaches have underestimated the importance of epistasis., (© 2015 The Author(s). Evolution © 2015 The Society for the Study of Evolution.)

Loss of the Y-chromosome is a common feature of species with chromosomal sex determination. However, our understanding of why some lineages frequently lose Y-chromosomes while others do not is limited. The fragile Y hypothesis proposes that in species with chiasmatic meiosis the rate of Y-chromosome aneuploidy and the size of the recombining region have a negative correlation. The fragile Y hypothesis provides a number of novel insights not possible under traditional models. Specifically, increased rates of Y aneuploidy may impose positive selection for (i) gene movement off the Y; (ii) translocations and fusions which expand the recombining region; and (iii) alternative meiotic segregation mechanisms (achiasmatic or asynaptic). These insights as well as existing evidence for the frequency of Y-chromosome aneuploidy raise doubt about the prospects for long-term retention of the human Y-chromosome despite recent evidence for stable gene content in older non-recombining regions., (© 2015 WILEY Periodicals, Inc.)

Recent efforts to catalog the diversity of sex chromosome systems coupled with genome sequencing projects are adding a new level of resolution to our understanding of insect sex chromosome origins. Y-chromosome degeneration makes sequencing difficult and may erase homology so rapidly that their origins will often remain enigmatic. X-chromosome origins are better understood, but thus far prove to be remarkably labile, often lacking homology even among close relatives. Furthermore, evidence now suggests that differentiated X or Y-chromosomes may both revert to autosomal inheritance. Data for ZW systems is scarcer, but W and Y-chromosomes seem to share many characteristics. Limited evidence suggests that Z-chromosome homology is more conserved than X counterparts, but broader sampling of both sex chromosome systems is needed., (Copyright © 2015 Elsevier Inc. All rights reserved.)

Chromosomal sex determination is phylogenetically widespread, having arisen independently in many lineages. Decades of theoretical work provide predictions about sex chromosome differentiation that are well supported by observations in both XY and ZW systems. However, the phylogenetic scope of previous work gives us a limited understanding of the pace of sex chromosome gain and loss and why Y or W chromosomes are more often lost in some lineages than others, creating XO or ZO systems. To gain phylogenetic breadth we therefore assembled a database of 4724 beetle species' karyotypes and found substantial variation in sex chromosome systems. We used the data to estimate rates of Y chromosome gain and loss across a phylogeny of 1126 taxa estimated from seven genes. Contrary to our initial expectations, we find that highly degenerated Y chromosomes of many members of the suborder Polyphaga are rarely lost, and that cases of Y chromosome loss are strongly associated with chiasmatic segregation during male meiosis. We propose the "fragile Y" hypothesis, that recurrent selection to reduce recombination between the X and Y chromosome leads to the evolution of a small pseudoautosomal region (PAR), which, in taxa that require XY chiasmata for proper segregation during meiosis, increases the probability of aneuploid gamete production, with Y chromosome loss. This hypothesis predicts that taxa that evolve achiasmatic segregation during male meiosis will rarely lose the Y chromosome. We discuss data from mammals, which are consistent with our prediction., (Copyright © 2014 by the Genetics Society of America.)

Examination of the genetic architecture of hybrid breakdown can provide insight into the genetic mechanisms of commonly observed isolating phenomena such as Haldane's rule. We used line-cross analysis to dissect the genetic architecture of divergence between two plant species that exhibit Haldane's rule for male sterility and rarity, Silene latifolia and Silene diclinis. We made 15 types of crosses, including reciprocal F1, F2, backcrosses, and later-generation crosses, grew the seeds to flowering, and measured the number of viable ovules, proportion of viable pollen, and sex ratio. Typically, Haldane's rule for male rarity in XY animal hybrids is explained by interactions involving recessive X-linked alleles that are deleterious when hemizygous (dominance theory), whereas sterility is explained by rapid evolution of spermatogenesis genes (faster-male evolution). In contrast, we found that the genetic mechanisms underlying Haldane's rule between the two Silene species did not follow these conventions. Dominance theory was sufficient to explain male sterility, but male rarity likely involved faster-male evolution. We also found an effect of the neo-sex chromosomes of S. diclinis on the extreme rarity of some hybrid males. Our findings suggest that the genetic architecture of Haldane's rule in dioecious plants may differ from those commonly found in animals., (© 2013 The Author(s). Evolution © 2013 The Society for the Study of Evolution.)

Gene expression levels correlate with multiple aspects of gene sequence and gene structure in phylogenetically diverse taxa, suggesting an important role of gene expression levels in the evolution of protein-coding genes. Here we present results of a genome-wide study of the influence of gene expression on synonymous codon usage, amino acid composition, and gene structure in the red flour beetle, Tribolium castaneum. Consistent with the action of translational selection, we find that synonymous codon usage bias increases with gene expression. However, the correspondence between tRNA gene copy number and optimal codons is weak. At the amino acid level, translational selection is suggested by the positive correlation between tRNA gene numbers and amino acid usage, which is stronger for highly expressed genes. In addition, there is a clear trend for increased use of metabolically cheaper, less complex amino acids as gene expression increases. tRNA gene numbers also correlate negatively with amino acid size/complexity (S/C) score indicating the coupling between translational selection and selection to minimize the use of large/complex amino acids. Interestingly, the analysis of 10 additional genomes suggests that the correlation between tRNA gene numbers and amino acid S/C score is widespread and might be explained by selection against negative consequences of protein misfolding. At the level of gene structure, three major trends are detected: 1) complete coding region length increases across low and intermediate expression levels but decreases in highly expressed genes; 2) the average intron size shows the opposite trend, first decreasing with expression, followed by a slight increase in highly expressed genes; and 3) intron density remains nearly constant across all expression levels. These changes in gene architecture are only in partial agreement with selection favoring reduced cost of biosynthesis.

During the process of speciation, diverging taxa often hybridize and produce offspring wherein the heterogametic sex (i.e., XY or ZW) is unfit (Haldane's rule). Dominance theory seeks to explain Haldane's rule in terms of the difference in X-linked dominance regimes experienced by the sexes. However, X inactivation in female mammals extends the effects of hemizygosity to both sexes. Here, we highlight where the assumptions of dominance theory are particularly problematic in marsupials, where X inactivation uniformly results in silencing the paternal X. We then present evidence of Haldane's rule for sterility but not for viability in marsupials, as well as the first violations of Haldane's rule for these traits among all mammals. Marsupials represent a large taxonomic group possessing heteromorphic sex chromosomes, where the dominance theory cannot explain Haldane's rule. In this light, we evaluate alternative explanations for the preponderance of male sterility in interspecific hybrids, including faster male evolution, X-Y interactions, and genomic conflict hypotheses.

While parasites and immunity are widely believed to play important roles in the evolution of male ornaments, their potential influence on systems where male weaponry is the object of sexual selection is poorly understood. We experimentally infect larval broad-horned flour beetles with a tapeworm and study the consequent effects on: 1) adult male morphology 2) male-male contests for mating opportunities, and 3) induction of the innate immune system. We find that infection significantly reduces adult male size in ways that are expected to reduce mating opportunities in nature. The sum of our morphological, competition, and immunological data indicate that during a life history stage where no new resources are acquired, males allocate their finite resources in a way that increases future mating potential.

We look at sex-limited chromosome (Y or W) evolution with particular emphasis on the importance of palindromes. Y chromosome palindromes consist of inverted duplicates that allow for local recombination in an otherwise nonrecombining chromosome. Since palindromes enable intrachromosomal gene conversion that can help eliminate deleterious mutations, they are often highlighted as mechanisms to protect against Y degeneration. However, the adaptive significance of recombination resides in its ability to decouple the evolutionary fates of linked mutations, leading to both a decrease in degeneration rate and an increase in adaptation rate. Our paper emphasizes the latter, that palindromes may exist to accelerate adaptation by increasing the potential targets and fixation rates of incoming beneficial mutations. This hypothesis helps reconcile two enigmatic features of the "palindromes as protectors" view: (1) genes that are not located in palindromes have been retained under purifying selection for tens of millions of years, and (2) under models that only consider deleterious mutations, gene conversion benefits duplicate gene maintenance but not initial fixation. We conclude by looking at ways to test the hypothesis that palindromes enhance the rate of adaptive evolution of Y-linked genes and whether this effect can be extended to palindromes on other chromosomes.

Bacterial two-hybrid analysis identified the Staphylococcus aureus RNA degradosome-like complex to include RNase J1, RNase J2, RNase Y, polynucleotide phosphorylase (PNPase), enolase, phosphofructokinase, and a DEAD box RNA helicase. Results also revealed that the recently recognized RNase RnpA interacts with the S. aureus degradosome and that this interaction is conserved in other Gram-positive organisms.

A genome's ability to produce two separate sexually dimorphic phenotypes is an intriguing biological mystery. Microarray-based studies of a handful of model systems suggest that much of the mystery can be explained by sex-biased gene expression evolved in response to sexually antagonistic selection. We present the first whole-genome study of sex-biased expression in the red flour beetle, Tribolium castaneum. Tribolium is a model for the largest eukaryotic order, Coleoptera, and we show that in whole-body adults, approximately 20% of the transcriptome is differentially regulated between the sexes. Among T. castaneum, Drosophila melanogaster, and Anopheles gambiae, we identify 416 1:1:1 orthologs with conserved sex-biased expression. Overrepresented functional categories among sex-biased genes are primarily those involved in gamete production and development. The genomic distribution of sex-biased genes in T. castaneum is distinctly nonrandom, with the strongest deficit of male-biased genes on the X chromosome (9 of 793) of any species studied to date. Tribolium also shows a significant enrichment of X-linked female-biased genes (408 of 793). Our analyses suggest that the extensive female bias of Tribolium X chromosome gene expression is due to hyperexpression of X-linked genes in both males and females. We propose that the overexpression of X chromosomes in females is an evolutionary side effect of the need to dosage compensate in males and that mechanisms to reduce female X chromosome gene expression to autosomal levels are sufficient but imperfect.

Duplicate genes act as a source of genetic material from which new functions arise. They exist in large numbers in every sequenced eukaryotic genome and may be responsible for many differences in phenotypes between species. However, recent work searching for the targets of positive selection in humans has largely ignored duplicated genes due to complications in orthology assignment. Here we find that a high proportion of young gene duplicates in the human, macaque, mouse, and rat genomes have experienced adaptive natural selection. Approximately 10% of all lineage-specific duplicates show evidence for positive selection on their protein sequences, larger than any reported amount of selection among single-copy genes in these lineages using similar methods. We also find that newly duplicated genes that have been transposed to new chromosomal locations are significantly more likely to have undergone positive selection than the ancestral copy. Human-specific duplicates evolving under adaptive natural selection include a surprising number of genes involved in neuronal and cognitive functions. Our results imply that genome scans for selection that ignore duplicated loci are missing a large fraction of all adaptive substitutions. The results are also in agreement with the classical model of evolution by gene duplication, supporting a common role for neofunctionalization in the long-term maintenance of gene duplicates.

One of the unique insights provided by the growing number of fully sequenced genomes is the pervasiveness of gene duplication and gene loss. Indeed, several metrics now suggest that rates of gene birth and death per gene are only 10-40% lower than nucleotide substitutions per site, and that per nucleotide, the consequent lineage-specific expansion and contraction of gene families may play at least as large a role in adaptation as changes in orthologous sequences. While gene family evolution is pervasive, it may be especially important in our own evolution since it appears that the "revolving door" of gene duplication and loss has undergone multiple accelerations in the lineage leading to humans. In this paper, we review current understanding of gene family evolution including: methods for inferring copy number change, evidence for adaptive expansion and adaptive contraction of gene families, the origins of new families and deaths of previously established ones, and finally we conclude with a perspective on challenges and promising directions for future research.

Tribolium castaneum is a member of the most species-rich eukaryotic order, a powerful model organism for the study of generalized insect development, and an important pest of stored agricultural products. We describe its genome sequence here. This omnivorous beetle has evolved the ability to interact with a diverse chemical environment, as shown by large expansions in odorant and gustatory receptors, as well as P450 and other detoxification enzymes. Development in Tribolium is more representative of other insects than is Drosophila, a fact reflected in gene content and function. For example, Tribolium has retained more ancestral genes involved in cell-cell communication than Drosophila, some being expressed in the growth zone crucial for axial elongation in short-germ development. Systemic RNA interference in T. castaneum functions differently from that in Caenorhabditis elegans, but nevertheless offers similar power for the elucidation of gene function and identification of targets for selective insect control.

The molecular changes responsible for the evolution of modern humans have primarily been discussed in terms of individual nucleotide substitutions in regulatory or protein coding sequences. However, rates of nucleotide substitution are slowed in primates, and thus humans and chimpanzees are highly similar at the nucleotide level. We find that a third source of molecular evolution, gene gain and loss, is accelerated in primates relative to other mammals. Using a novel method that allows estimation of rate heterogeneity among lineages, we find that the rate of gene turnover in humans is more than 2.5 times faster than in other mammals and may be due to both mutational and selective forces. By reconciling the gene trees for all of the gene families included in the analysis, we are able to independently verify the numbers of inferred duplications. We also use two methods based on the genome assembly of rhesus macaque to further verify our results. Our analyses identify several gene families that have expanded or contracted more rapidly than is expected even after accounting for an overall rate acceleration in primates, including brain-related families that have more than doubled in size in humans. Many of the families showing large expansions also show evidence for positive selection on their nucleotide sequences, suggesting that selection has been important in shaping copy-number differences among mammals. These findings may help explain why humans and chimpanzees show high similarity between orthologous nucleotides yet great morphological and behavioral differences.

Background: Duchenne muscular dystrophy results from mutation of the dystrophin gene, causing skeletal and cardiac muscle loss of function. The mdx mouse model of Duchenne muscular dystrophy is widely utilized to evaluate the potential of therapeutic regimens to modulate the loss of skeletal muscle function associated with dystrophin mutation. Importantly, progressive loss of diaphragm function is the most consistent striated muscle effect observed in the mdx mouse model, which is the same as in patients suffering from Duchenne muscular dystrophy., Methods: Using the mdx mouse model, we have evaluated the effect that corticotrophin releasing factor 2 receptor (CRF2R) agonist treatment has on diaphragm function, morphology and gene expression., Results: We have observed that treatment with the potent CRF2R-selective agonist PG-873637 prevents the progressive loss of diaphragm specific force observed during aging of mdx mice. In addition, the combination of PG-873637 with glucocorticoids not only prevents the loss of diaphragm specific force over time, but also results in recovery of specific force. Pathological analysis of CRF2R agonist-treated diaphragm muscle demonstrates that treatment reduces fibrosis, immune cell infiltration, and muscle architectural disruption. Gene expression analysis of CRF2R-treated diaphragm muscle showed multiple gene expression changes including globally decreased immune cell-related gene expression, decreased extracellular matrix gene expression, increased metabolism-related gene expression, and, surprisingly, modulation of circadian rhythm gene expression., Conclusion: Together, these data demonstrate that CRF2R activation can prevent the progressive degeneration of diaphragm muscle associated with dystrophin gene mutation.

The rhesus macaque (Macaca mulatta) is an abundant primate species that diverged from the ancestors of Homo sapiens about 25 million years ago. Because they are genetically and physiologically similar to humans, rhesus monkeys are the most widely used nonhuman primate in basic and applied biomedical research. We determined the genome sequence of an Indian-origin Macaca mulatta female and compared the data with chimpanzees and humans to reveal the structure of ancestral primate genomes and to identify evidence for positive selection and lineage-specific expansions and contractions of gene families. A comparison of sequences from individual animals was used to investigate their underlying genetic diversity. The complete description of the macaque genome blueprint enhances the utility of this animal model for biomedical research and improves our understanding of the basic biology of the species.

The heterogametic sex tends to be rare, absent, sterile, or deformed in F1 hybrid crosses between species, a pattern called Haldane's rule (HR). The introgression of single genes or chromosomal regions from one drosophilid species into the genetic background of another have shown that HR is most often associated with fixed genetic differences in inter-specific crosses. However, because such introgression studies have involved species diverged several hundred thousand generations from a common ancestor, it is not clear whether HR attends the speciation process or results from the accumulation of epistatically acting genes postspeciation. We report the first evidence for HR prior to speciation in crosses between two populations of the red flour beetle, Tribolium castaneum, collected 931 km apart in Colombia and Ecuador. In this cross, HR is manifested as an increase in the proportion of deformed males compared to females and the expression of HR is temperature dependent. Neither population, when crossed to a geographically distant population from Japan, exhibits HR at any rearing temperature. Using joint-scaling analysis and additional data from backcrosses and F2's, we find that the hybrid incompatibilities and the emergence of HR are concurrent processes involving interactions between X-linked and autosomal genes. However, we also find many examples of incompatibilities manifest by F2 and backcross hybrids but not by F1 hybrids and most incompatibilities are not sex different in their effects, even when they involve both X-autosomal interactions and genotype-by-environment interactions. We infer that incipient speciation in flour beetles can occur with or without HR and that significant hybrid incompatibilities result from the accumulation of epistatically acting gene differences between populations without differentially affecting the heterogametic sex in F1 hybrids. The temperature dependence of the incompatibilities supports the inference that genotype-by-environment interactions and adaptation to different environments contribute to the genetic divergence important to postzygotic reproductive isolation.

We used joint-scaling analyses in conjunction with rearing temperature variation to investigate the contributions of additive, non-additive, and environmental effects to genetic divergence and incipient speciation among 12 populations of the red flour beetle, Tribolium castaneum, with small levels of pairwise nuclear genetic divergence (0.033 < Nei's D < 0.125). For 15 population pairs we created a full spectrum of line crosses (two parental, two reciprocal F1's, four F2's, and eight backcrosses), reared them at multiple temperatures, and analyzed the numbers and developmental defects of offspring. We assayed a total of 219,388 offspring from 5147 families. Failed crosses occurred predominately in F2's, giving evidence of F2 breakdown within this species. In all cases where a significant model could be fit to the data on offspring number, we observed at least one type of digenic epistasis. We also found maternal and cytoplasmic effects to be common components of divergence among T. castaneum populations. In some cases, the most complex model tested (additive, dominance, epistatic, maternal, and cytoplasmic effects) did not provide a significant fit to the data, suggesting that linkage or higher order epistasis is involved in differentiation between some populations. For the limb deformity data, we observed significant genotype-by-environment interaction in most crosses and pure parent crosses tended to have fewer deformities than hybrid crosses. Complexity of genetic architecture was not correlated with either geographic distance or genetic distance. Our results support the view that genetic incompatibilities responsible for postzygotic isolation, an important component of speciation, may be a natural but serendipitous consequence of nonadditive genetic effects and structured populations.

Population genetic theory predicts that maternal effect genes will evolve differently than genes expressed in both sexes because selection is only half as effective on autosomal genes expressed in one sex but not the other. Here, we use sequences of the tandem gene duplicates, bicoid (bcd) and zerknüllt (zen), to test the prediction that, with similar coefficients of purifying selection, a maternal effect gene evolves more rapidly than a zygotic gene because of this reduction in selective constraint. We find that the maternal effect gene, bcd, is evolving more rapidly than zygotically expressed, zen, providing the first direct confirmation of this prediction of maternal effect theory from molecular evidence. Our results extend current explanations for the accelerated rate of bcd evolution by providing an evolutionary mechanism, relaxed selective constraint, that allows bcd the evolutionary flexibility to escape the typical functional constraints of early developmental genes. We discuss general implications of our findings for the role of maternal effect genes in early developmental patterning.

Gene families are groups of homologous genes that are likely to have highly similar functions. Differences in family size due to lineage-specific gene duplication and gene loss may provide clues to the evolutionary forces that have shaped mammalian genomes. Here we analyze the gene families contained within the whole genomes of human, chimpanzee, mouse, rat, and dog. In total we find that more than half of the 9,990 families present in the mammalian common ancestor have either expanded or contracted along at least one lineage. Additionally, we find that a large number of families are completely lost from one or more mammalian genomes, and a similar number of gene families have arisen subsequent to the mammalian common ancestor. Along the lineage leading to modern humans we infer the gain of 689 genes and the loss of 86 genes since the split from chimpanzees, including changes likely driven by adaptive natural selection. Our results imply that humans and chimpanzees differ by at least 6% (1,418 of 22,000 genes) in their complement of genes, which stands in stark contrast to the oft-cited 1.5% difference between orthologous nucleotide sequences. This genomic "revolving door" of gene gain and loss represents a large number of genetic differences separating humans from our closest relatives.

Summary: We present CAFE (Computational Analysis of gene Family Evolution), a tool for the statistical analysis of the evolution of the size of gene families. It uses a stochastic birth and death process to model the evolution of gene family sizes over a phylogeny. For a specified phylogenetic tree, and given the gene family sizes in the extant species, CAFE can estimate the global birth and death rate of gene families, infer the most likely gene family size at all internal nodes, identify gene families that have accelerated rates of gain and loss (quantified by a p-value) and identify which branches cause the p-value to be small for significant families., Availability: Software is available from http://www.bio.indiana.edu/~hahnlab/Software.html

Objective: Studies have reported that administration of stromal cell-derived factor-1 (SDF-1), the ligand for the G-protein coupled receptor CXCR4, increased collateral blood flow in a mouse model of vascular insufficiency via recruitment of endothelial precursor cells (EPC). The present study investigated the contribution of mature endothelial cells in the actions of SDF-1., Methods: The regulation of SDF-1 and CXCR4 was examined in the rat cornea cauterization (CC) and aortic ring (AR) model. The functional significance of the SDF-1/CXCR4 pathway was explored in cultured endothelial cells, the AR model, and on collateral blood flow in a rat model of vascular insufficiency., Results: In the present study, the CXCR4 transcript was dramatically upregulated in the rat CC and AR explants, systems containing and lacking bone marrow-derived EPCs, respectively. Addition of AMD3100, a selective CXCR4 antagonist, had no effect on vessel growth in the AR alone, but completely inhibited SDF-1 mediated increases in vascular sprouting. In cultured endothelial cells, SDF-1 alone or in combination with vascular endothelial growth factor (VEGF) significantly enhanced cell survival and migration. Finally, systemic administration of SDF-1 in a rat model of arterial insufficiency enhanced collateral blood flow above vehicle control and equal to that of VEGF after 2 weeks of treatment., Conclusion: These studies support activation of the SDF-1/CXCR4 axis as a means to promote blood vessel growth and enhance collateral blood flow, at least in part, via direct effects on vascular endothelial cells.

The origin of evolutionary novelty is believed to involve both positive selection and relaxed developmental constraint. In flies, the redesign of anterior patterning during embryogenesis is a major developmental innovation and the rapidly evolving Hox gene, bicoid (bcd), plays a critical role. We report evidence for relaxation of selective constraint acting on bicoid as a result of its maternal pattern of gene expression. Evolutionary theory predicts 2-fold greater sequence diversity for maternal effect genes than for zygotically expressed genes, because natural selection is only half as effective acting on autosomal genes expressed in one sex as it is on genes expressed in both sexes. We sample an individual from ten populations of Drosophila melanogaster and nine populations of D. simulans for polymorphism in the tandem gene duplicates bcd, which is maternally expressed, and zerknüllt (zen), which is zygotically expressed. In both species, we find the ratio of bcd to zen nucleotide diversity to be two or more in the coding regions but one in the noncoding regions, providing the first quantitative support for the theoretical prediction of relaxed selective constraint on maternal-effect genes resulting from sex-limited expression. Our results suggest that the accelerated rate of evolution observed for bcd is owing, at least partly, to variation generated by relaxed selective constraint., Competing Interests: Competing interests. The authors have declared that no competing interests exist.

We present a quantitative genetic (QG) interpretation of the Bateson-Dobzhansky-Muller (BDM) genetic model of speciation in order to unify the theoretical framework for understanding how the genetic differentiation of populations is associated with the process of speciation. Specifically, we compare the QG theory of joint scaling with the Turelli-Orr mathematical formulation of the BDM model. By formally linking the two models, we show that a wealth of empirical methods from QG can be brought to bear on the study of the genetic architecture of hybrid phenotypes to better understand the connections, if any, between microevolution within populations and macroevolution in the origin of species. By integrating the two theories, we make additional novel predictions that enrich the opportunities for empirically testing speciation genetic theory or facets of it, such as Haldane's rule. We show that the connection between the two theories is simple and straightforward for autosomal genes but not for sex-linked genes. Differences between the two approaches highlight key conceptual issues concerning the relevance of epistasis to evolution within and among lineages and to differences in the process of speciation in hermaphrodites and in organisms with separate sexes.

In a verbal model, Trivers and Willard proposed that, whenever there is sexual selection among males, natural selection should favor mothers that produce sons when in good condition but daughters when in poor condition. The predictions of this model have been the subject of recent debate. We present an explicit population genetic model for the evolution of a maternal-effect gene that biases offspring sex ratio. We show that, like local mate competition, sexual selection favors female-biased sex ratios whenever maternal condition affects the reproductive competitive ability of sons. However, Fisherian sex-ratio selection, which favors a balanced sex ratio, is an opposing force. We show that the evolution of maternal sex-ratio biasing by these opposing selection forces requires a positive covariance across environments between the sex-ratio bias toward sons (b) and the mating success of sons (r). This covariance alone is not a sufficient condition for the evolution of maternal sex-ratio biasing; it must be sufficiently positive to outweigh the opposing sex-ratio selection. To identify the necessary and sufficient conditions, we partition total evolutionary change into three components: (1) maternal sex-ratio bias, (2) sexual selection on sons, and (3) sex-ratio selection. Because the magnitude of the first component asymmetrically affects the strength of the second, biasing broods toward females in a poor environment evolves faster than the same degree of bias toward males in a good environment. Consequently, female-biased sex ratios, rather than male-biased sex ratios, are more likely to evolve. We discuss our findings in the context of the primary sex-ratio biases observed in strongly sexually selected species and indicate how this perspective can assist the experimental study of sex ratio evolution.

Morphological analysis of cytologic samples obtained by fine-needle aspirate (FNA) or bronchoscopy is an important method for diagnosing bronchogenic carcinoma. However, this approach has only about 65 to 80% diagnostic sensitivity. Based on previous studies, the c-myc x E2F-1/p21WAF1/CIP1 (p21 hereafter) gene expression index is highly sensitive and specific for distinguishing normal from malignant bronchial epithelial tissues. In an effort to improve sensitivity of diagnosing lung cancer in cytologic specimens, we used Standardized Reverse Transcriptase Polymerase Chain Reaction (StaRT-PCR) to measure the c-myc x E2F-1/p21 index in cDNA samples from 14 normal lung samples (6 normal lung parenchyma and 8 normal bronchial epithelial cell [NBEC] biopsies), and 16 FNA biopsies from 14 suspected tumors. Based on cytomorphologic criteria, 11 of the 14 suspected tumors were diagnosed as bronchogenic carcinoma and three specimens were non-diagnostic. Subsequent biopsy samples confirmed that the three non-diagnostic samples were derived from lung carcinomas. The index value for each bronchogenic carcinoma was above a cut-off value of 7000 and the index value of all but one normal sample was below 7000. Thus the c-myc x E2F-1/p21 index may augment cytomorphologic diagnosis of bronchogenic carcinoma biopsy samples, particularly those considered non-diagnostic by cytomorphologic criteria.

Radon exposure is associated with increased risk for bronchogenic carcinoma. Mutagenesis analyses have revealed that radon induces mostly multi-locus chromosome deletions. Based on these findings, it was hypothesized that deletion analysis of multiple radon-induced malignant transformants would reveal common mutations in chromosomal regions containing tumor suppressor genes responsible for malignant transformation. This hypothesis was supported by a previous study in which tumorigenic derivatives of the human papillomavirus 18-immortalized human bronchial epithelial cell line BEP2D were established following irradiation with 30 cGy of high linear energy transfer radon-simulated alpha-particles. Herein, we describe the analyses of 10 additional tumorigenic derivative cell lines resulting from the irradiation of five additional independent BEP2D populations. The new transformants have common cytogenetic changes, including the loss of chromosome (ch)Y, one of three copies of ch8, one of two copies of ch11p15-pter and one of three copies of ch14. These changes are the same as those reported previously. Analysis of PCR-amplified short tandem repeats of informative loci confirmed the loss of heterozygosity (LOH) at 12 loci spanning the length of ch8 in cell lines from four of the total of eight irradiation treatments to date and the loss of chY in all cell lines (8 of 8). LOH analysis with a total of 17 informative loci confirmed loss on ch14 in transformants from seven of eight irradiation treatments and indicated a 0.5-1.7 cM region of common involvement centered around locus D14S306. No LOH was detected at any of the informative loci on ch11. The overall results support our stated hypothesis. Further studies are currently in progress to determine whether the ch8 and ch14 regions contain genes with tumor suppressor function in bronchial epithelial cells.

Bronchogenic carcinomas arise from bronchial epithelial cells (BECs). Inhalation exposure of BECs to nitrosamines in cigarette smoke is an important exogenous risk factor for malignant transformation of BECs. Thus, an important endogenous risk factor is likely to be the capacity of BECs to metabolize nitrosamines. Among the cytochrome P450 enzymes capable of metabolizing nitrosamines, CYP2A6, CYP2E1 and CYP2B6 are expressed in BECs. In this study, we used quantitative RT-PCR to evaluate expression of CYP2A6 and CYP2E1 in primary human BECs from 12 non-smokers and eight smokers. CYP2A6 was expressed in 20/20 cases and quantifiable in 18/20 cases, with a mean level of 580 mRNA/10(6) beta-actin mRNA. CYP2E1 expression was observed in 9/20 cases, but in all cases it was expressed at levels below our limit of quantification (10 mRNA/10(6) beta-actin mRNA). There was significant (P < 0.05) 20-fold inter-individual variation in expression of CYP2A6. Further, the mean level of CYP2A6 among smokers (260 mRNA/10(6) beta-actin mRNA) was significantly lower than among non-smokers (740 mRNA/10(6) beta-actin mRNA). It is hypothesized that: (i) inter-individual variation in CYP2A6 gene expression may contribute to inter-individual variation in risk for bronchogenic carcinoma; (ii) smoking may reduce the level of expression of CYP2A6 in the BECs of some individuals; and (iii) CYP2A6 is more important than CYP2E1 for metabolic activation of nitrosamines in bronchial epithelial cells.

Progress toward complete sequencing of all human genes through the Human Genome Project has already resulted in a need for methods that allow quantitative expression measurement of multiple genes simultaneously. It is increasingly recognized that relative measurement of multiple genes will provide more mechanistic information regarding cell pathophysiology than measurement of individual genes one by one or by methods that do not allow direct intergene comparison. In this study, previously described quantitative reverse transcription-polymerase chain reaction methods were modified in an effort to provide a rapid, simple method for this purpose. Internal standard competitive templates (CTs) were prepared for each gene and were combined in a single solution containing CTs for more than 40 genes at defined concentrations relative to one another. Any subsequent dilution of the CT mixture did not alter the relationship of one CT to another. Because the same CT standard solution or a dilution of it was used in all experiments, data obtained from different experiments were easily compared. The use of multiple CT mixtures with different housekeeping gene to target gene ratios provided a linear dynamic range spanning the range of expression of all genes thus far evaluated. CT stock solutions were used to simultaneously quantify the expression of 25 genes relative to beta-actin and glyceraldehyde-3-phosphate dehydrogenase in normal and malignant bronchial epithelial cells. Because the CT concentrations were known, data in the form of both absolute messenger RNA (mRNA) copy number and mRNA relative to housekeeping gene mRNA were obtained. The methods and reagents described will allow rapid, quantitative measurement of multiple genes simultaneously, using inexpensive and widely available equipment. Furthermore, the CT standard solution may be distributed to other investigators for interlaboratory standardization of experimental conditions.

Expression of the small, proline-rich protein (spr1) squamous differentiation marker was measured in five cultured normal and 12 malignant human bronchial epithelial cell (BEC) populations by quantitative reverse transcriptase polymerase chain reaction (RT-PCR). Whereas spr1 expression was quantifiable and inducible in all five cultured normal cell populations, in all 12 carcinoma cell lines evaluated it was neither quantifiable nor inducible. Primers spanning the entire spr1 coding sequence amplified full-length PCR product from genomic DNA; therefore, large deletions in the coding region were not responsible for the loss of expression measurable by RT-PCR. This is the first molecular genetic marker reported that distinguishes all normal from all carcinoma cell populations evaluated. Because the spr1 protein is a component of the crosslinked envelope that forms during the squamous differentiation process, we hypothesize that the apparent loss of spr1 gene expression disrupts mechanisms for terminal squamous differentiation in the bronchial epithelium, thereby contributing to malignant transformation.

Recent methodological developments allow expression measurement of many genes simultaneously, thereby revealing patterns of gene expression that can be related to phenotype. We hypothesized that through the use of such methods we could identify patterns of gene expression associated with the malignant phenotype in human bronchial epithelial cells (BEC). To test this hypothesis, a recently developed quantitative reverse transcriptase polymerase chain reaction method was used to assess simultaneously expression of 15 genes mechanistically associated with cell-cycle control (c-myc, E2F-1, p21, rb, PCNA, cyclin D2, cyclin D3, cyclin E, cdc2, CDK2, CDK4, mad, max p21, max p22, and p53) in normal cell cultures from five individuals and in nine different malignant BEC lines. Relative to the mean expression levels in cultured normal cell populations, expression of c-myc, E2F-1, PCNA, cyclin E, and CDK4 messenger RNA (mRNA) were significantly increased and expression of p21 and p53 mRNA were significantly decreased in one or two, but not all three subtypes (squamous, adenocarcinoma and small cell) of carcinoma cell lines evaluated. No single cell-cycle control gene discriminated all three subtypes from normal cell populations. In contrast, the gene expression index c-myc x E2F-1/p21 separated all carcinoma cell lines from all normal cell populations initially evaluated. This malignancy index was validated in an additional three cultured normal BEC and three carcinoma cell lines, as well as three pairs of matched primary normal bronchial epithelial and primary bronchogenic carcinoma samples, and three pairs of matched primary normal lung parenchyma and primary bronchogenic carcinoma tissue. Again, the c-myc x E2F-1/ p21 index successfully discriminated all cultured and primary normal from malignant samples and thereby had a predictive value of 1 (no false positives and no false negatives). We hypothesize that because of functional mutations in cell-cycle regulatory genes (e.g., p53 and/or rb), cells lose the ability to maintain a pattern of gene expression mechanistically associated with normal, division-limited homeostatic equilibrium. Because the c-myc x E2F-1/p21 gene expression index has high specificity for malignant tissue, it will allow confirmation that there is a significant amount of tumor tissue present in small (e.g., fine-needle) biopsy specimens prior to evaluating them for expression of other genes, such as those involved in chemoresistance or radioresistance. In addition, the goal of most gene therapy efforts is to alter levels of gene expression quantitatively. This index and others derived in a similar manner may better define potential gene therapy targets as well as response of targeted genes to therapy.