Molecular phylogenetics of western North American frogs of the Rana boylii species group

Phylogenetic relationships among frogs of the genus Rana from western North America are investigated using 2013 aligned bases of mitochondrial DNA sequence from the genes encoding ND1 (subunit one of NADH dehydrogenase), tRNA(Ile), tRNA(Gln), tRNA(Met), ND2, tRNA(Trp), tRNA(Ala), tRNA(Asn), tRNA(Cys), tRNA(Tyr), and COI (subunit I of cytochrome c oxidase), plus the origin for light-strand replication (O(L)) between the tRNA(Asn) and tRNA(Cys) genes. The aligned sequences contain 401 phylogenetically informative characters. A well-resolved phylogenetic hypothesis in which the Rana boylii species group (R. aurora, R. boylii, R. cascadae, R. muscosa, and R. pretiosa) is monophyletic is obtained. Molecular sequence divergence suggests that the R. boylii species group is approximately 8 million years old. The traditional hypothesis showing monophyly of the yellow-legged frogs (R. boylii and R. muscosa) is statistically rejected in favor of a hypothesis in which R. aurora, R. cascadae, and R. muscosa form a clade. Reanalyses of published nuclear ribosomal DNA restriction-site data and allozymic data support a monophyletic R. boylii group, but do not effectively resolve relationships among species within this group. Eight populations of R. muscosa form two major clades separated by a biogeographic break in the Sierra Nevada of California. This biogeographic break is broadly concordant with breaks found in four other amphibian and reptilian taxa. The two major clades within R. muscosa are estimated to have diverged approximately 2.2 million years before present. Each of these major clades contains two subgroups showing approximately 1.5 million years divergence, implicating climatic effects of Pleistocene glaciation in vicariance. The four distinct subgroups of R. muscosa separated by at least 1.4 million years of evolutionary divergence are suggested as potential units for conservation.     

Technical Advance: a novel technique for the sensitive quantification of acyl CoA esters from plant tissues

We report a novel, highly sensitive and selective method for the extraction and quantification of acyl CoA esters from plant tissues. The method detects acyl CoA esters with acyl chain lengths from C4 to C20 down to concentrations as low as 6 fmol in extracts. Acyl CoA esters from standard solutions or plant extracts were derived to their fluorescent acyl etheno CoA esters in the presence of chloroacetaldehyde, separated by ion-paired reversed-phase high-performance liquid chromatography, and detected fluorometrically. This derivitization procedure circumvents the selectivity problems associated with previously published enzymatic methods, and methods that rely on acyl chain or thiol group modification for acyl CoA ester detection. The formation of acyl etheno CoA esters was verified by mass spectrometry, which was also used to identify unknown peaks from chromatograms of plant extracts. Using this method, we report the composition and concentration of the acyl CoA pool during lipid synthesis in maturing Brassica napus seeds and during storage lipid breakdown in 2-day-old Arabidopsis thaliana seedlings. The concentrations measured were in the 3--6 microM range for both tissue types. We also demonstrate the utility of acyl CoA profiling in a transgenic B. napus line that has high levels of lauric acid. To our knowledge, this is the first time that reliable estimates of acyl CoA ester concentrations have been made for higher plants, and the ability to profile these metabolites provides a valuable new tool for the investigation of gene function.     

What can phylogeographic breaks in the livebearer Poeciliopsis turrubarensis tell about fish community boundaries in western Costa Rica?

This study asks if the geographic boundary delineating two fish communities in western Costa Rica is congruent with a phylogeographic break in a single widespread fish species Poeciliopsis turrubarensis (Poeciliidae) that spans this area. Such congruence would suggest that a common historical event (e.g. geological or climatic) could be responsible for both patterns. It was found that there was a shared break across a region in central Costa Rica suggesting a common cause may be responsible for both the abrupt shift in fish community composition and the genetic break in P. turrubarensis.     

Molecular Evolution across Mouse Spermatogenesis

Genes involved in spermatogenesis tend to evolve rapidly, but we lack a clear understanding of how protein sequences and patterns of gene expression evolve across this complex developmental process. We used fluorescence-activated cell sorting (FACS) to generate expression data for early (meiotic) and late (postmeiotic) cell types across 13 inbred strains of mice (Mus) spanning ∼7 My of evolution. We used these comparative developmental data to investigate the evolution of lineage-specific expression, protein-coding sequences, and expression levels. We found increased lineage specificity and more rapid protein-coding and expression divergence during late spermatogenesis, suggesting that signatures of rapid testis molecular evolution are punctuated across sperm development. Despite strong overall developmental parallels in these components of molecular evolution, protein and expression divergences were only weakly correlated across genes. We detected more rapid protein evolution on the X chromosome relative to the autosomes, whereas X-linked gene expression tended to be relatively more conserved likely reflecting chromosome-specific regulatory constraints. Using allele-specific FACS expression data from crosses between four strains, we found that the relative contributions of different regulatory mechanisms also differed between cell types. Genes showing cis-regulatory changes were more common late in spermatogenesis, and tended to be associated with larger differences in expression levels and greater expression divergence between species. In contrast, genes with trans-acting changes were more common early and tended to be more conserved across species. Our findings advance understanding of gene evolution across spermatogenesis and underscore the fundamental importance of developmental context in molecular evolutionary studies.     

Molecular pathogenesis of a new glycogenosis caused by a glycogenin-1 mutation

Glycogenin-1 initiates the glycogen synthesis in skeletal muscle by the autocatalytic formation of a short oligosaccharide at tyrosine 195. Glycogenin-1 catalyzes both the glucose-O-tyrosine linkage and the α1,4 glucosidic bonds linking the glucose molecules in the oligosaccharide. We recently described a patient with glycogen depletion in skeletal muscle as a result of a non-functional glycogenin-1. The patient carried a Thr83Met substitution in glycogenin-1. In this study we have investigated the importance of threonine 83 for the catalytic activity of glycogenin-1. Non-glucosylated glycogenin-1 constructs, with various amino acid substitutions in position 83 and 195, were expressed in a cell-free expression system and autoglucosylated in vitro. The autoglucosylation was analyzed by gel-shift on western blot, incorporation of radiolabeled UDP-(14)C-glucose and nano-liquid chromatography with tandem mass spectrometry (LC/MS/MS). We demonstrate that glycogenin-1 with the Thr83Met substitution is unable to form the glucose-O-tyrosine linkage at tyrosine 195 unless co-expressed with the catalytically active Tyr195Phe glycogenin-1. Our results explain the glycogen depletion in the patient expressing only Thr83Met glycogenin-1 and why heterozygous carriers without clinical symptoms show a small proportion of unglucosylated glycogenin-1.     

Visualization of proteins in acrylamide gels using ultraviolet illumination.

Proteins in polyacrylamide gels can be rapidly visualized by soaking in trichloroacetic acid or chloroform followed by illumination with UV light. The UV-light-driven reaction of tryptophan in the presence of trichlorocompounds yields products that emit sufficiently in the visible region to identify the location of the protein bands on the gel. This method can be used to rapidly identify protein bands on a gel in less than 20 min. On thin polyacrylamide gels, 1.0 microg of protein can easily be detected for proteins with typical tryptophan percentages.     

Mechanistic Phenotypes: An Aggregative Phenotyping Strategy to Identify Disease Mechanisms Using GWAS Data

A single mutation can alter cellular and global homeostatic mechanisms and give rise to multiple clinical diseases. We hypothesized that these disease mechanisms could be identified using low minor allele frequency (MAF<0.1) non-synonymous SNPs (nsSNPs) associated with “mechanistic phenotypes”, comprised of collections of related diagnoses. We studied two mechanistic phenotypes: (1) thrombosis, evaluated in a population of 1,655 African Americans; and (2) four groupings of cancer diagnoses, evaluated in 3,009 white European Americans. We tested associations between nsSNPs represented on GWAS platforms and mechanistic phenotypes ascertained from electronic medical records (EMRs), and sought enrichment in functional ontologies across the top-ranked associations. We used a two-step analytic approach whereby nsSNPs were first sorted by the strength of their association with a phenotype. We tested associations using two reverse genetic models and standard additive and recessive models. In the second step, we employed a hypothesis-free ontological enrichment analysis using the sorted nsSNPs to identify functional mechanisms underlying the diagnoses comprising the mechanistic phenotypes. The thrombosis phenotype was solely associated with ontologies related to blood coagulation (Fisher''s p = 0.0001, FDR p = 0.03), driven by the F5, P2RY12 and F2RL2 genes. For the cancer phenotypes, the reverse genetics models were enriched in DNA repair functions (p = 2×10−5, FDR p = 0.03) (POLG/FANCI, SLX4/FANCP, XRCC1, BRCA1, FANCA, CHD1L) while the additive model showed enrichment related to chromatid segregation (p = 4×10−6, FDR p = 0.005) (KIF25, PINX1). We were able to replicate nsSNP associations for POLG/FANCI, BRCA1, FANCA and CHD1L in independent data sets. Mechanism-oriented phenotyping using collections of EMR-derived diagnoses can elucidate fundamental disease mechanisms.     

ANATOMICAL CHANGES DURING LEAF ONTOGENY IN POPULUS DELTOIDES

The leaf plastochron index (LPI) was used to interpret the anatomical changes during leaf ontogeny in the developing leaf zone of young cottonwood trees and to relate leaf anatomical structure to physiological function. The lamina tip matured precociously with respect to both structure and function. Below the lamina tip the intercellular spaces, stomates, and minor veins matured basipetally, while the major veins developed acropetally. Ontogenetically, maturation progressed from LPI –1.0, which was anatomically immature except for its lamina tip, to the first fully expanded leaf at LPI 6.0, which was anatomically mature. Physiological maturity also occurred at LPI 6.0, thus signifying a transition with respect to both structure and function. By evaluating the anatomical observations in conjunction with physiological data collected at comparable LPI's in other experiments, it could be demonstrated that anatomical development was a limiting factor in photosynthesis and translocation of assimilates.     

The future is now: Amplicon sequencing and sequence capture usher in the conservation genomics era

The genomics revolution has initiated a new era of population genetics where genome‐wide data are frequently used to understand complex patterns of population structure and selection. However, the application of genomic tools to inform management and conservation has been somewhat rare outside a few well studied species. Fortunately, two recently developed approaches, amplicon sequencing and sequence capture, have the potential to significantly advance the field of conservation genomics. Here, amplicon sequencing refers to highly multiplexed PCR followed by high‐throughput sequencing (e.g., GTseq), and sequence capture refers to using capture probes to isolate loci from reduced‐representation libraries (e.g., Rapture). Both approaches allow sequencing of thousands of individuals at relatively low costs, do not require any specialized equipment for library preparation, and generate data that can be analyzed without sophisticated computational infrastructure. Here, we discuss the advantages and disadvantages of each method and provide a decision framework for geneticists who are looking to integrate these methods into their research programme. While it will always be important to consider the specifics of the biological question and system, we believe that amplicon sequencing is best suited for projects aiming to genotype <500 loci on many individuals (>1,500) or for species where continued monitoring is anticipated (e.g., long‐term pedigrees). Sequence capture, on the other hand, is best applied to projects including fewer individuals or where >500 loci are required. Both of these techniques should smooth the transition from traditional genetic techniques to genomics, helping to usher in the conservation genomics era.