Statistical properties of the variation at linked microsatellite loci: implications for the history of human Y chromosomes [published erratum appears in Mol Biol Evol 1997 Mar;14(3):354]

It has recently been suggested that observed levels of variation at microsatellite loci can be used to infer patterns of selection in genomes and to assess demographic history. In order to evaluate the feasibility of these suggestions it is necessary to know something about how levels of variation at microsatellite loci are expected to fluctuate due simply to stochasticity in the processes of mutation and inheritance (genetic sampling). Here we use recently derived properties of the stepwise mutation model to place confidence intervals around the variance in repeat score that is expected at mutation-drift equilibrium and outline a statistical test for whether an observed value differs significantly from expectation. We also develop confidence intervals for the time course of the buildup of variation following a complete elimination of variation, such as might be caused by a selective sweep or an extreme population bottleneck. We apply these methods to the variation observed at human Y-specific microsatellites. Although a number of authors have suggested the possibility of a very recent sweep, our analyses suggest that a sweep or extreme bottleneck is unlikely to have occurred anytime during the last approximately 74,000 years. To generate this result we use a recently estimated mutation rate for microsatellite loci of 5.6 x 10(-4) along with the variation observed at autosomal microsatellite loci to estimate the human effective population size. This estimate is 18,000, implying an effective number of 4,500 Y chromosomes. One important general conclusion to emerge from this study is that in order to reject mutation-drift equilibrium at a set of linked microsatellite loci it is necessary to have an unreasonably large number of loci unless the observed variance is far below that expected at mutation-drift equilibrium.      

The Mitochondrial Genomes of Nuttalliella namaqua (Ixodoidea: Nuttalliellidae) and Argas africolumbae (Ixodoidae: Argasidae): Estimation of Divergence Dates for the Major Tick Lineages and Reconstruction of Ancestral Blood-Feeding Characters

Ixodida are composed of hard (Ixodidae), soft (Argasidae) and the monotypic Nuttalliellidae (Nuttalliella namaqua) tick families. Nuclear 18S rRNA analysis suggested that N. namaqua was the closest extant relative to the last common ancestral tick lineage. The mitochondrial genomes of N. namaqua and Argas africolumbae were determined using next generation sequencing and de novo assembly to investigate this further. The latter was included since previous estimates on the divergence times of argasids lacked data for this major genus. Mitochondrial gene order for both was identical to that of the Argasidae and Prostriata. Bayesian analysis of the COI, Cytb, ND1, ND2 and ND4 genes confirmed the monophyly of ticks, the basal position of N. namaqua to the other tick families and the accepted systematic relationships of the other tick genera. Molecular clock estimates were derived for the divergence of the major tick lineages and supported previous estimates on the origins of ticks in the Carboniferous. N. namaqua larvae fed successfully on lizards and mice in a prolonged manner similar to many argasids and all ixodids. Excess blood meal-derived water was secreted via the salivary glands, similar to ixodids. We propose that this prolonged larval feeding style eventually gave rise to the long feeding periods that typify the single larval, nymphal and adult stages of ixodid ticks and the associated secretion of water via the salivary glands. Ancestral reconstruction of characters involved in blood-feeding indicates that most of the characteristics unique to either hard or soft tick families were present in the ancestral tick lineage.     

Mutational analysis of the pseudoknot in the tRNA-like structure of turnip yellow mosaic virus RNA. Aminoacylation efficiency and RNA pseudoknot stability.

Site-directed mutations were introduced in the connecting loops and one of the two stem regions of the RNA pseudoknot in the tRNA-like structure of turnip yellow mosaic virus RNA. The kinetic parameters of valylation for each mutated RNA were determined in a cell-free extract from wheat germ. Structure mapping was performed on most mutants with enzymic probes, like RNase T1, nuclease S1 and cobra venom ribonuclease. An insertion of four A residues in the four-membered connecting loop L1 that crosses the deep groove of the pseudoknot reduces aminoacylation efficiency. Deletions up to three nucleotides do not affect aminoacylation or RNA pseudoknot formation. Deletion of the entire loop abolishes aminoacylation. Although elimination of the pseudoknot is presumed, this could not be demonstrated. Unlike the mutations in loop L1, all mutations in the three-membered connecting loop L2 that crosses the shallow groove of the RNA pseudoknot decrease the aminoacylation efficiency considerably. Nonetheless, the RNA pseudoknot is still present in most mutated RNAs. These results indicate that a number of mutations can be introduced in both loops without abolishing aminoacylation. Results obtained with the introduction of mismatches and A.U base-pairs in stem S1 of the pseudoknot, containing three G.C base-pairs in wild-type RNA, indicate that the pseudoknot is only marginally stable. Our estimation of the gain of free energy due to the pseudoknot formation is at most 2.0 kcal/mol. The pseudoknot structure can, however, be stabilized upon binding the valyl-tRNA synthetase.     

MetaTM - a consensus method for transmembrane protein topology prediction

Background  

Transmembrane (TM) proteins are proteins that span a biological membrane one or more times. As their 3-D structures are hard to determine, experiments focus on identifying their topology (i. e. which parts of the amino acid sequence are buried in the membrane and which are located on either side of the membrane), but only a few topologies are known. Consequently, various computational TM topology predictors have been developed, but their accuracies are far from perfect. The prediction quality can be improved by applying a consensus approach, which combines results of several predictors to yield a more reliable result.     

Cloning and characterization of 29 tetranucleotide and two dinucleotide polymorphic microsatellite loci from the endangered marbled murrelet (Brachyramphus marmoratus)

We developed 31 novel, polymorphic microsatellite loci in the marbled murrelet (Brachyramphus marmoratus), a critically endangered seabird. Variability was tested on 15 individuals from the Santa Cruz, California population, with each locus characterized by two to 12 alleles. Observed levels of heterozygosity ranged from 0.13 to 0.93. These loci provide a valuable means of assessing the population structure and demographic parameters of this species, which may be critical to its conservation across a fragmented habitat.     

Crystal structure of the murine cytomegalovirus MHC-I homolog m144

Large DNA viruses of the herpesvirus family produce proteins that mimic host MHC-I molecules as part of their immunoevasive strategy. The m144 glycoprotein, expressed by murine cytomegalovirus, is thought to be an MHC-I homolog whose expression prolongs viral survival in vivo by preventing natural killer cell activation. To explore the structural basis of this m144 function, we have determined the three-dimensional structure of an m144/beta2-microglobulin (beta2m) complex at 1.9A resolution. This structure reveals the canonical features of MHC-I molecules including readily identifiable alpha1, alpha2, and alpha3 domains. A unique disulfide bond links the alpha1 helix to the beta-sheet floor, explaining the known thermal stability of m144. Close juxtaposition of the alpha1 and alpha2 helices and the lack of critical residues that normally contribute to anchoring the peptide N and C termini eliminates peptide binding. A region of 13 amino acid residues, corresponding to the amino-terminal portion of the alpha2 helix, is missing in the electron density map, suggesting an area of structural flexibility that may be involved in ligand binding.