Systematics within Gyps vultures: a clade at risk

Background  

Populations of the Oriental White-backed Vulture (Gyps bengalensis) have declined by over 95% within the past decade. This decline is largely due to incidental consumption of the non-steroidal anti-inflammatory veterinary pharmaceutical diclofenac, commonly used to treat domestic livestock. The conservation status of other Gyps vultures in southern Asia is also of immediate concern, given the lack of knowledge regarding status of their populations and the continuing existence of taxonomic uncertainties. In this study, we assess phylogenetic relationships for all recognized species and the majority of subspecies within the genus Gyps. The continuing veterinary use of diclofenac is an unknown but potential risk to related species with similar feeding habits to Gyps bengalensis. Therefore, an accurate assessment of the phylogenetic relationships among Gyps vultures should aid in their conservation by clarifying taxonomic uncertainties, and enabling inference of their respective relatedness to susceptible G. bengalensis.     

Patterns of divergence during evolution of alpha 1-proteinase inhibitors in mammals

alpha 1-Proteinase inhibitor (alpha 1-PI), a member of the serine proteinase inhibitor superfamily, has a primary role in controlling neutrophil elastase activity within the mammalian circulation. Several studies have indicated that the reactive center region of alpha 1-PI, the amino acid sequence of which is critical to recognition of and binding to target proteinases, is highly divergent within and among species. This appears to be a consequence of accelerated rates of evolution that may have been driven by positive Darwinian selection. In order to examine this and other features of alpha 1-PI evolution in more detail, we have isolated and sequenced cDNAs representing alpha 1- PI mRNAs of the mouse species Mus saxicola and Mus minutoides and have compared these with a number of other mammalian alpha 1-PI mRNAs. Relative to other mammalian mRNAs, the extent of nonsynonymous substitution is generally high throughout the alpha 1-PI mRNA molecule, indicating greater overall rates of amino acid substitution. Within and among mouse species, the 5'-half of the mRNA, but not the 3'-half, has been homogenized by concerted evolution. Finally, the reactive center is under diversifying or positive Darwinian selection in murid rodents (rats, mice) and guinea pigs yet is under purifying selection in primates and artiodactyls. The significance of these findings to alpha 1-PI function and the possible selective forces driving evolution of serpins in general are discussed.      

Mitochondrial DNA polymorphisms in subterranean mole-rats of the Spalax ehrenbergi superspecies in Israel, and its peripheral isolates

Patterns of mitochondrial DNA (mtDNA) variation were examined in 133 mole-rats constituting all four chromosomal species (2n = 52, 2n = 54, 2n = 58, and 2n = 60) of the Spalax ehrenbergi superspecies in Israel, as well as the peripheral isolates of 2n = 60. In the main range of the complex, a total of 28 mtDNA haplotypes were found in 64 mole-rats, with most haplotypes being unique to either a single chromosomal species or population. mtDNA divergence increased from low to high diploid number in a north-to-south direction in Israel. Overall levels of mtDNA diversity were unexpectedly the highest in the 2n = 60, the youngest species of the complex. The mtDNA haplotypes can be separated into two major groups, 2n = 52-54 and 2n = 58-60, and a phylogenetic analysis for each group revealed evidence of a few haplotypes not sorted by diploid number. The overall patterns of mtDNA divergence seen within and among the four chromosomal species are consistent with the parapatric mode of speciation as suggested from previous studies of allozyme and DNA hybridization. In a separate data set the patterns of mtDNA variation were examined across the main geographic range and across peripheral semi-isolates and isolates of the 2n = 60 chromosomal species. Fifteen haplotypes were found in 69 mole-rats. High levels of mtDNA diversity characterized the main range, semi-isolated, and even some desert isolated populations. The peripheral isolates contain much mtDNA diversity, including novel haplotypes.      

Charge-charge interactions are key determinants of the pK values of ionizable groups in ribonuclease Sa (pI=3.5) and a basic variant (pI=10.2)

The pK values of the titratable groups in ribonuclease Sa (RNase Sa) (pI=3.5), and a charge-reversed variant with five carboxyl to lysine substitutions, 5K RNase Sa (pI=10.2), have been determined by NMR at 20 degrees C in 0.1M NaCl. In RNase Sa, 18 pK values and in 5K, 11 pK values were measured. The carboxyl group of Asp33, which is buried and forms three intramolecular hydrogen bonds in RNase Sa, has the lowest pK (2.4), whereas Asp79, which is also buried but does not form hydrogen bonds, has the most elevated pK (7.4). These results highlight the importance of desolvation and charge-dipole interactions in perturbing pK values of buried groups. Alkaline titration revealed that the terminal amine of RNase Sa and all eight tyrosine residues have significantly increased pK values relative to model compounds.A primary objective in this study was to investigate the influence of charge-charge interactions on the pK values by comparing results from RNase Sa with those from the 5K variant. The solution structures of the two proteins are very similar as revealed by NMR and other spectroscopic data, with only small changes at the N terminus and in the alpha-helix. Consequently, the ionizable groups will have similar environments in the two variants and desolvation and charge-dipole interactions will have comparable effects on the pK values of both. Their pK differences, therefore, are expected to be chiefly due to the different charge-charge interactions. As anticipated from its higher net charge, all measured pK values in 5K RNase are lowered relative to wild-type RNase Sa, with the largest decrease being 2.2 pH units for Glu14. The pK differences (pK(Sa)-pK(5K)) calculated using a simple model based on Coulomb's Law and a dielectric constant of 45 agree well with the experimental values. This demonstrates that the pK differences between wild-type and 5K RNase Sa are mainly due to changes in the electrostatic interactions between the ionizable groups. pK values calculated using Coulomb's Law also showed a good correlation (R=0.83) with experimental values. The more complex model based on a finite-difference solution to the Poisson-Boltzmann equation, which considers desolvation and charge-dipole interactions in addition to charge-charge interactions, was also used to calculate pK values. Surprisingly, these values are more poorly correlated (R=0.65) with the values from experiment. Taken together, the results are evidence that charge-charge interactions are the chief perturbant of the pK values of ionizable groups on the protein surface, which is where the majority of the ionizable groups are positioned in proteins.     

Phylogenetic utility of the nuclear gene arginine decarboxylase: an example from Brassicaceae

Arginine decarboxylase (ADC) is an important enzyme in the production of putrescine and polyamines in plants. It is encoded by a single or low-copy nuclear gene that lacks introns in sequences studied to date. The rate of Adc amino acid sequence evolution is similar to that of ndhF for the angiosperm family studied. Highly conserved regions provide several target sites for PCR priming and sequencing and aid in nucleotide and amino acid sequence alignment across a range of taxonomic levels, while a variable region provides an increased number of potentially informative characters relative to ndhF for the taxa surveyed. The utility of the Adc gene in plant molecular systematic studies is demonstrated by analysis of its partial nucleotide sequences obtained from 13 representatives of Brassicaceae and 3 outgroup taxa, 2 from the mustard oil clade (order Capparales) and 1 from the related order Malvales. Two copies of the Adc gene, Adc1 and Adc2, are found in all members of the Brassicaceae studied to data except the basal genus Aethionema. The resulting Adc gene tree provides robust phylogenetic data regarding relationships within the complex mustard family, as well as independent support for proposed tribal realignments based on other molecular data sets such as those from chloroplast DNA.      

Isocitrate dehydrogenase 1 R132C mutation occurs exclusively in microsatellite stable colorectal cancers with the CpG island methylator phenotype

The CpG Island Methylator Phenotype (CIMP) is fundamental to an important subset of colorectal cancer; however, its cause is unknown. CIMP is associated with microsatellite instability but is also found in BRAF mutant microsatellite stable cancers that are associated with poor prognosis. The isocitrate dehydrogenase 1 (IDH1) gene causes CIMP in glioma due to an activating mutation that produces the 2-hydroxyglutarate oncometabolite. We therefore examined IDH1 alteration as a potential cause of CIMP in colorectal cancer. The IDH1 mutational hotspot was screened in 86 CIMP-positive and 80 CIMP-negative cancers. The entire coding sequence was examined in 81 CIMP-positive colorectal cancers. Forty-seven cancers varying by CIMP-status and IDH1 mutation status were examined using Illumina 450K DNA methylation microarrays. The R132C IDH1 mutation was detected in 4/166 cancers. All IDH1 mutations were in CIMP cancers that were BRAF mutant and microsatellite stable (4/45, 8.9%). Unsupervised hierarchical cluster analysis identified an IDH1 mutation-like methylation signature in approximately half of the CIMP-positive cancers. IDH1 mutation appears to cause CIMP in a small proportion of BRAF mutant, microsatellite stable colorectal cancers. This study provides a precedent that a single gene mutation may cause CIMP in colorectal cancer, and that this will be associated with a specific epigenetic signature and clinicopathological features.