Evolutionary Genetic Studies Research Articles

Mutation patterns of mtDNA: Empirical inferences for the coding region

BackgroundHuman mitochondrial DNA (mtDNA) has been extensively used in population and evolutionary genetics studies. Thus, a valid estimate of human mtDNA evolutionary rate is important in many research fields. The small number of estimations performed for the coding region of the molecule, showed important differences between phylogenetic and empirical approaches. We analyzed a portion of the coding region of mtDNA (tRNALeu, ND1 and tRNAIle genes), using individuals belonging to extended families from the Azores Islands (Portugal) with the main aim of providing empirical estimations of the mutation rate of the coding region of mtDNA under different assumptions, and hence to better understand the mtDNA evolutionary process.ResultsHeteroplasmy was detected in 6.5% (3/46) of the families analyzed. In all of the families the presence of mtDNA heteroplasmy resulted from three new point mutations, and no cases of insertions or deletions were identified. Major differences were found in the proportion and type of heteroplasmy found in the genes studied when compared to those obtained in a previous report for the D-loop. Our empirical estimation of mtDNA coding region mutation rate, calculated taking into account the sex of individuals carrying new mutations, the probability of intra-individual fixation of mutations present in heteroplasmy and, to the possible extent, the effect of selection, is similar to that obtained using phylogenetic approaches.ConclusionBased on our results, the discrepancy previously reported between the human mtDNA coding region mutation rates observed along evolutionary timescales and estimations obtained using family pedigrees can be resolved when correcting for the previously cited factors.

Observation of non-nucleated erythrocytes in the peripheral blood of medaka, Oryzias latipes

The medaka, Oryzias latipes is a useful animal model for the study of vertebrate developmental genetics. Using May-Grünwald-Giemsa stain, we found non-nucleated erythrocytes in the peripheral blood of medaka. Eleven of 50 fish occasionally showed non-nucleated erythrocytes in their peripheral blood. We expect that this observation will be useful in future studies involving screening for hematologic mutants.

Molecular basis for skeletal variation: insights from developmental genetic studies in mice.

Skeletal variations are common in humans, and potentially are caused by genetic as well as environmental factors. We here review molecular principles in skeletal development to develop a knowledge base of possible alterations that could explain variations in skeletal element number, shape or size. Environmental agents that induce variations, such as teratogens, likely interact with the molecular pathways that regulate skeletal development.

Development of Y Chromosome Intraspecific Polymorphic Markers in the Felidae

Y chromosome haplotyping based on microsatellites and single nucleotide polymorphisms (SNPs) has proved to be a powerful tool for population genetic studies of humans. However, the promise of the approach is hampered in the majority of nonhuman mammals by the lack of Y-specific polymorphic markers. We were able to identify new male-specific polymorphisms in the domestic cat Felis catus and 6 additional Felidae species with a combination of molecular genetic and cytogenetic approaches including 1) identifying domestic cat male-specific microsatellites from markers generated from a male cat microsatellite-enriched genomic library, a flow-sorted Y cosmid library, or a Y-specific cat bacteria artificial chromosome (BAC) clone, (2) constructing microsatellite-enriched libraries from flow-sorted Y chromosomes isolated directly from focal wildcat species, and (3) screening Y chromosome conserved anchored tagged sequences primers in Felidae species. Forty-one male-specific microsatellites were identified, but only 6 were single-copy loci, consistent with the repetitive nature of the Y chromosome. Nucleotide diversity (pi) of Y-linked intron sequences (2.1 kbp) was in the range of 0 (tiger) to 9.95 x 10(-4) (marbled cat), and the number of SNPs ranged from none in the tiger to 7 in the Asian leopard cat. The Y haplotyping system described here, consisting of 4 introns (SMCY3, SMCY7, UTY11, and DBY7) and 1 polymorphic microsatellite (SMCY-STR), represents the first available markers for tracking intraspecific male lineage polymorphisms in Felidae species and promises to provide significant insights to evolutionary and population genetic studies of the species.

Isolation and characterization of polymorphic microsatellite markers in the endangered Nothotsuga longibracteata (Pinaceae)

AbstractTwelve polymorphic microsatellite loci were isolated and characterized from an AC‐enriched genomic library of Nothotsuga longibracteata. The average allele number of these microsatellites was 8.3 per locus, ranging from two to 13. The ranges of observed and expected heterozygosities were 0.03–0.97 and 0.09–0.88, respectively. These polymorphic markers provide useful tools for the study of evolutionary history and conservation genetics of N. longibracteata.

Abundant sequence variation around the mitochondrial origin of replication in the human opportunistic yeast pathogen Candida albicans from a tropical island in China

Due to its haploid nature and its predominantly uniparental mode of inheritance, the mitochondrial genome has been analyzed extensively in population and evolutionary genetic studies of eukaryotes. Among the mitochondrial DNA (mtDNA) fragments, the region surrounding the origin of replication is the most commonly studied. However, most of the studies have focused on animals and little is known about the extent and patterns of sequence variation around the mtDNA origin of replication (mtOri) in fungi. In this study, we found abundant variation in a 597 bp fragment surrounding the mtOri for 53 isolates of the pathogenic yeast Candida albicans obtained from a diverse group of hosts in Hainan, a tropical island of China. Within this DNA fragment, a total of 17 haplotypes were found for the 53 isolates. The extent of sequence variation for this group of strains was similar to that for 24 strains that represented the global nuclear genotype diversity. In contrast to those in animals where there were significant biases in favor of transitional mutations (e.g. the transition to transversion ratio is about 20 in human mtDNA), our data showed a transition to transversion ratio of ∼0.5 around mtOri in C. albicans. Our analysis revealed no apparent geographic pattern of sequence variation based on the birthplaces of the analyzed hosts. However, the sample from patients had a lower genotypic diversity than those from healthy hosts borne either in Hainan or elsewhere in China. Our results suggest that C. albicans mtDNA has a base substitution pattern different from its nuclear genome and that sequences from the mtOri region could enhance our understanding of C. albicans genome evolution and population structuring.

Morphing the hyomandibular skeleton in development and evolution

How might changes in developmental regulatory pathways underlie evolutionary changes in morphology? Here we focus on a particular pathway regulated by a secreted, signaling peptide, Endothelin1 (Edn1). Developmental genetic analyses show the Edn1-pathway to be crucial for hyomandibular patterning, and we discuss our work with zebrafish suggesting how the signal may function in regulating numbers of skeletal elements, their sizes and their shapes. We then review a broader collection of comparative studies that examine morphological evolution of a subset of the same skeletal elements-the opercular-branchiostegal series of bones of the hyoid arch. We find that phenotypic changes in zebrafish mutants copy evolutionary changes that recur along many actinopterygian lineages. Hence the developmental genetic studies are informative for providing candidate pathways for macroevolution of facial morphology, as well as for our understanding of how these pathways work.

Simple sequence repeats in zebra finch (Taeniopygia guttata) expressed sequence tags: a new resource for evolutionary genetic studies of passerines

BackgroundPasserines (perching birds) are widely studied across many biological disciplines including ecology, population biology, neurobiology, behavioural ecology and evolutionary biology. However, understanding the molecular basis of relevant traits is hampered by the paucity of passerine genomics tools. Efforts to address this problem are underway, and the zebra finch (Taeniopygia guttata) will be the first passerine to have its genome sequenced. Here we describe a bioinformatic analysis of zebra finch expressed sequence tag (EST) Genbank entries.ResultsA total of 48,862 ESTs were downloaded from GenBank and assembled into contigs, representing an estimated 17,404 unique sequences. The unique sequence set contained 638 simple sequence repeats (SSRs) or microsatellites of length ≥20 bp and purity ≥90% and 144 simple sequence repeats of length ≥30 bp. A chromosomal location for the majority of SSRs was predicted by BLASTing against assembly 2.1 of the chicken genome sequence. The relative exonic location (5' untranslated region, coding region or 3' untranslated region) was predicted for 218 of the SSRs, by BLAST search against the ENSEMBL chicken peptide database. Ten loci were examined for polymorphism in two zebra finch populations and two populations of a distantly related passerine, the house sparrow Passer domesticus. Linkage was confirmed for four loci that were predicted to reside on the passerine homologue of chicken chromosome 7.ConclusionWe show that SSRs are abundant within zebra finch ESTs, and that their genomic location can be predicted from sequence similarity with the assembled chicken genome sequence. We demonstrate that a useful proportion of zebra finch EST-SSRs are likely to be polymorphic, and that they can be used to build a linkage map. Finally, we show that many zebra finch EST-SSRs are likely to be useful in evolutionary genetic studies of other passerines.

Challenges and prospects of population genetic studies in terns (Charadriiformes, Aves)

Little information is available about the population structure of communally nesting terns (Sternidae) and skimmers (Rynchopidae) throughout the world. In order to fill this gap, a survey of molecular markers was carried out for six species of terns (Anous stolidus, Sterna hirundinacea, S. fuscata, S. superciliaris, Thalasseus maximus and Phaetusa simplex) and one species of skimmer (Rynchops niger). First, we describe the results of the construction of genomic DNA libraries and document problems encountered during this procedure. Secondly, we tested the cross-amplification of 18 microsatellite loci previously described for related species (the number of polymorphic loci ranged from three to seven). Thirdly, we tested the usefulness of mtDNA (control region, ND2, Cytochrome b and ATPase 6/8) for phylogeographic studies in this group of birds. The occurrence of nuclear copies of the mitochondrial control region is reported. Nucleotide divergence in the mtDNA genes analyzed ranged from 0.0 to 0.006. Despite the difficulties associated with the selection of variable markers in this group of seabirds, we were able to select polymorphic markers for each species tested and we anticipate these results will help the development of genetic studies concerning important biological questions in terns.

Isolation and characterization of microsatellite markers in pangolins (Mammalia, Pholidota, Manis spp.)

AbstractThirty‐four polymorphic dinucleotide microsatellite loci were developed in the Malayan pangolin Manis javanica. Of the 34 markers, 32 and 18 were also amplified, respectively, in the Chinese pangolin (Manis pentadactyla) and the African tree pangolin (Manis tricuspis). Analysis of 24 Malayan, 12 Chinese and 2 African tree pangolins showed high levels of variability (heterozygosity ranging from 0.321 to 0.708). These are the first available microsatellite markers in Pholidota and will be an invaluable tool for evolutionary and conservation genetic studies in pangolins.

Inflammatory caspases and inflammasomes: master switches of inflammation

Fifteen years have passed since the cloning and characterization of the interleukin-1beta-converting enzyme (ICE/caspase-1), the first identified member of a family of proteases currently known as caspases. Caspase-1 is the prototypical member of a subclass of caspases involved in cytokine maturation termed inflammatory caspases that also include caspase-4 caspase -5, caspase -11 and caspase -12. Efforts to elucidate the molecular mechanisms involved in the activation of these proteases have uncovered an important role for the NLR family members, NALPs, NAIP and IPAF. These proteins promote the assembly of multiprotein complexes termed inflammasomes, which are required for activation of inflammatory caspases. This article will review some evolutionary aspects, biochemical evidences and genetic studies, underlining the role of inflammasomes and inflammatory caspases in innate immunity against pathogens, autoinflammatory syndromes and in the biology of reproduction.

POWSIM: a computer program for assessing statistical power when testing for genetic differentiation

AbstractKnowledge of statistical power is essential for sampling design and data evaluation when testing for genetic differentiation. Yet, such information is typically missing in studies of conservation and evolutionary genetics, most likely because of complex interactions between the many factors that affect power. powsim is a 32‐bit Windows/DOS simulation‐based computer program that estimates power (and α error) for chi‐square and Fisher's exact tests when evaluating the hypothesis of genetic homogeneity. Optional combinations include the number of samples, sample sizes, number of loci and alleles, allele frequencies, and degree of differentiation (quantified as FST). powsim is available at http://www.zoologi.su.se/~ryman.

JADE: a distributed Java application for deleterious genomic mutation (DGM) estimation

The characterization of deleterious genomic mutation (DGM) is of central significance for evolutionary biology and genetic studies. Fitness moment method has been developed to efficiently characterize DGM from natural population directly. In order to enable researchers to employ this method for theoretical and empirical research on characterizing DGM, we here present a distributed Java Application for DGM Estimation (JADE). http://orclinux.creighton.edu/DGM/index.htm.

MECHANISM OF LEAF-SHAPE DETERMINATION

Biodiversity of plant shape is mainly attributable to biodiversity of leaf shape and the shape of floral organs, the modified leaves. However, the exact mechanisms of leaf-shape determination remain unclear due to the complexity of flat-structure organogenesis that includes the simultaneous cell cycling and cell enlargement in primordia. Recent studies in developmental and molecular genetics have revealed several important aspects of leaf-shape control mechanisms. For example, understanding of polar control in leaf-blade expansion has advanced greatly. A curious phenomenon called "compensated cell enlargement" found in leaf organogenesis studies should also provide interesting clues regarding the mechanisms of multicellular organ development. This paper reviews recent research findings with a focus on leaf development in Arabidopsis thaliana.

Novel Relationships Among Ten Fish Model Species Revealed Based on a Phylogenomic Analysis Using ESTs

The power of comparative phylogenomic analyses also depends on the amount of data that are included in such studies. We used expressed sequence tags (ESTs) from fish model species as a proof of principle approach in order to test the reliability of using ESTs for phylogenetic inference. As expected, the robustness increases with the amount of sequences. Although some progress has been made in the elucidation of the phylogeny of teleosts, relationships among the main lineages of the derived fish (Euteleostei) remain poorly defined and are still debated. We performed a phylogenomic analysis of a set of 42 of orthologous genes from 10 available fish model systems from seven different orders (Salmoniformes, Siluriformes, Cypriniformes, Tetraodontiformes, Cyprinodontiformes, Beloniformes, and Perciformes) of euteleostean fish to estimate divergence times and evolutionary relationships among those lineages. All 10 fish species serve as models for developmental, aquaculture, genomic, and comparative genetic studies. The phylogenetic signal and the strength of the contribution of each of the 42 orthologous genes were estimated with randomly chosen data subsets. Our study revealed a molecular phylogeny of higher-level relationships of derived teleosts, which indicates that the use of multiple genes produces robust phylogenies, a finding that is expected to apply to other phylogenetic issues among distantly related taxa. Our phylogenomic analyses confirm that the euteleostean superorders Ostariophysi and Acanthopterygii are monophyletic and the Protacanthopterygii and Ostariophysi are sister clades. In addition, and contrary to the traditional phylogenetic hypothesis, our analyses determine that killifish (Cyprinodontiformes), medaka (Beloniformes), and cichlids (Perciformes) appear to be more closely related to each other than either of them is to pufferfish (Tetraodontiformes). All 10 lineages split before or during the fragmentation of the supercontinent Pangea in the Jurassic.

Expression profiling the developing mammalian enteric nervous system identifies marker and candidate Hirschsprung disease genes

The enteric nervous system (ENS) is composed of neurons and glial cells, organized as interconnected ganglia within the gut wall, which controls peristalsis of the gut wall and secretions from its glands. The Ret receptor tyrosine kinase is expressed throughout enteric neurogenesis and is required for normal ENS development; humans with mutations in the RET locus have Hirschsprung disease (HSCR, an absence of ganglia in the colon), and mice lacking Ret have total intestinal aganglionosis. The Ret mutant mouse provides a tool for identifying genes implicated in development of the ENS. By using RNA from WT and Ret mutant (aganglionic) gut tissue and DNA microarrays, we have conducted a differential screen for ENS-expressed genes and have identified hundreds of candidate ENS-expressed genes. Forty-seven genes were selected for further analysis, representing diverse functional classes. We show that all of the analyzed genes are expressed in the ENS and that the screen was sensitive enough to identify genes marking only subpopulations of ENS cells. Our screen, therefore, was reliable and sensitive and has identified many previously undescribed genes for studying ENS development. Moreover, two of the genes identified in our screen Arhgef3 and Ctnnal1, have human homologues that map to previously identified HSCR susceptibility loci, thus representing excellent candidates for HSCR genes. This comprehensive profile of ENS gene expression refines our understanding of ENS development and serves as a resource for future developmental, biochemical, and human genetic studies.

A Comparison of Phasing Algorithms for Trios and Unrelated Individuals

Knowledge of haplotype phase is valuable for many analysis methods in the study of disease, population, and evolutionary genetics. Considerable research effort has been devoted to the development of statistical and computational methods that infer haplotype phase from genotype data. Although a substantial number of such methods have been developed, they have focused principally on inference from unrelated individuals, and comparisons between methods have been rather limited. Here, we describe the extension of five leading algorithms for phase inference for handling father-mother-child trios. We performed a comprehensive assessment of the methods applied to both trios and to unrelated individuals, with a focus on genomic-scale problems, using both simulated data and data from the HapMap project. The most accurate algorithm was PHASE (v2.1). For this method, the percentages of genotypes whose phase was incorrectly inferred were 0.12%, 0.05%, and 0.16% for trios from simulated data, HapMap Centre d'Etude du Polymorphisme Humain (CEPH) trios, and HapMap Yoruban trios, respectively, and 5.2% and 5.9% for unrelated individuals in simulated data and the HapMap CEPH data, respectively. The other methods considered in this work had comparable but slightly worse error rates. The error rates for trios are similar to the levels of genotyping error and missing data expected. We thus conclude that all the methods considered will provide highly accurate estimates of haplotypes when applied to trio data sets. Running times differ substantially between methods. Although it is one of the slowest methods, PHASE (v2.1) was used to infer haplotypes for the 1 million-SNP HapMap data set. Finally, we evaluated methods of estimating the value of r(2) between a pair of SNPs and concluded that all methods estimated r(2) well when the estimated value was >or=0.8.

A Hyperspace Model to Decipher the Genetic Architecture of Developmental Processes: Allometry Meets Ontogeny

To better utilize limited resources for their survival and reproduction, all organisms undergo developmental changes in both body size and shape during ontogeny. The genetic analysis of size change with increasing age, i.e., growth, has received considerable attention in quantitative developmental genetic studies, but the genetic architecture of ontogenetic changes in body shape and its associated allometry have been poorly understood partly due to the lack of analytical tools. In this article, we attempt to construct a multivariate statistical framework for studying the genetic regulation of ontogenetic growth and shape. We have integrated biologically meaningful mathematical functions of growth curves and developmental allometry into the estimation process of genetic mapping aimed at identifying individual quantitative trait loci (QTL) for phenotypic variation. This model defined with high dimensions can characterize the ontogenetic patterns of genetic effects of QTL over the lifetime of an organism and assess the interplay between genetic actions/interactions and phenotypic integration. The closed forms for the residual covariance matrix and its determinant and inverse were derived to overcome the computational complexity typical of our high-dimensional model. We used a worked example to validate the utility of this model. The implications of this model for genetic research of evo-devo are discussed.

Fundamentals of Fungal Molecular Population Genetic Analyses

The last two decades have seen tremendous growth in the development and application of molecular methods in the analyses of fungal species and populations. In this paper, I provide an overview of the molecular techniques and the basic analytical tools used to address various fundamental population and evolutionary genetic questions in fungi. With increasing availability and decreasing cost, DNA sequencing is becoming a mainstream data acquisition method in fungal evolutionary genetic studies. However, other methods, especially those based on the polymerase chain reaction, remain powerful in addressing specific questions for certain groups of taxa. These developments are bringing fungal population and evolutionary genetics into mainstream ecology and evolutionary biology.

S100 Expression in the Lateral Line System of Zebrafish from Juvenile to Adult Stages

The lateral line system (LLS) is a sensory organ present in fishes and in the larval stage of aquatic amphibians, which provides the animal with information about its surrounding environment. Mechanoreceptive sensory cells in the LLS of fishes are grouped into neuromasts that are considered the morphological unit of this system. Structurally the neuromast consists of a central cluster of hair cells with a ring of support cells with some lateral mantle cells. The LLS is localized in both the body surface (superficial neuromast or pit organs) and in a system of canals distributed along the head, trunk and tail (deep neuromast) (Cernuda‐Cernuda et al., 1996). The study of hair cells in the neuromast of fishes is becoming of a great interest, since they are similar in structure and function to the sensory cells of the inner ear and show some identical immunohistochemical properties (Foster et al., 1993). Recently, we have demonstrated the occurrence and the presence of S100 protein in the neuromast hair cells of some teleostean species (Abbate et al., 2002; Germanà et al., 2004). The zebrafish has become an ideal model for the study of developmental genetics, and the neuromast has been used as a model for the study of vertebrate auditory system (Withfield, 2002), so the need for cell selective markers has increased dramatically. Therefore, in this study we have investigated by tissue and whole‐mount immunohistochemistry the presence of S100 protein in both deep and superficial neuromasts in the zebrafish larval, juvenile and adult stages. The results of this work demonstrated that the localization of S100 protein was restricted only to the hair cells of the neuromast at all the examined ages. The mantle cells and supporting cells within the neuromast were unreactive. Nevertheless, heterogeneity has been detected in the S100 protein expression by hair cells, independent of their localization within the neuromast. Our data suggest that the S100 protein is a specific marker for sensory cells in the neuromast of zebrafish at post‐embryonic different stages of development, and it could be important for the analysis of development, regeneration and turnover of hair cells.