Associations Among Parental Caregiving Quality, Cannabinoid Receptor 1 Expression-Based Polygenic Scores, and Infant-Parent Attachment: Evidence for Differential Genetic Susceptibility?

Attachment is a biological evolutionary system contributing to infant survival. When primary caregivers/parents are sensitive and responsive to their infants' needs, infants develop a sense of security. Secure infant attachment has been linked to healthy brain and organ-system development. Belsky and colleagues proposed the term differential susceptibility to describe context-dependent associations between genetic variations and behavioral outcomes as a function of parenting environments. Variations in the Cannabinoid Receptor Gene 1 (CNR1) are associated with memory, mood, and reward and connote differential susceptibility to more and less optimal parental caregiving quality in predicting children's behavioral problems. AIM: To determine if parental caregiving quality interacts with children's expression-based polygenic risk score (ePRS) for the CNR1 gene networks in the prefrontal cortex, striatum, and hippocampus in predicting the probability of attachment security and disorganized attachment. DESIGN: Prospective correlational methods examined maternal-infant pairs (n = 142) from which infants provided DNA samples at 3 months. Parental caregiving quality was assessed via the Child Adult Relationship Experiment (CARE)-index at 6 months, and attachment security via the Strange Situation Procedure at a mean age of 22 months. The CNR1 ePRSs include genes co-expressed with the CNR1 genes in the prefrontal cortex, striatum, or hippocampus, and were calculated using the effect size of the association between the individual single nucleotide polymorphisms from those genes and region-specific gene expression (GTEx). Logistic regression was employed (alpha < 0.05, two-tailed) to examine the main and interaction effects between parental caregiving quality and ePRSs in predicting attachment patterns. Interpretation of results was aided by analyses that distinguished between differential susceptibility and diathesis-stress. RESULTS: Significant interactions were observed between (1) maternal sensitivity and ePRS in the striatum in predicting attachment security, (2) maternal unresponsiveness with the ePRS in the hippocampus in predicting disorganization, and (3) maternal controlling with the ePRS in the hippocampus in predicting disorganization. CONCLUSION: These findings offer support for genetic differential susceptibility to the quality of maternal sensitivity in the context of the ePRS in the striatum. However, the significant interactions between hippocampal ePRS and maternal unresponsiveness and controlling in predicting the probability of disorganization were more suggestive of the diathesis-stress model.

Highly Elevated Prevalence of Spinobulbar Muscular Atrophy in Indigenous Communities in Canada Due to a Founder Effect.

OBJECTIVE: Spinobulbar muscular atrophy (SBMA) is an X-linked adult-onset neuromuscular disorder that causes progressive weakness and androgen insensitivity in hemizygous males. This condition is reported to be extremely rare, but has higher prevalence in certain populations due to multiple founder effects. Anecdotal observations of a higher prevalence of SBMA in patients of Indigenous descent in Saskatchewan led us to perform this study, to estimate the disease prevalence, and to attempt to identify a founder effect. METHODS: For our prevalence estimation, we identified patients with confirmed SBMA diagnosis from the Saskatoon neuromuscular clinic database for comparison with population data available from Statistics Canada. For our haplotype analysis, participants with SBMA were recruited from 2 neuromuscular clinics, as well as 5 control participants. Clinical data were collected, as well as a DNA sample using saliva kits. We performed targeted quantification of DXS1194, DXS1111, DXS135, and DXS1125 microsatellite repeats and the AR GGC repeat to attempt to identify a disease haplotype and compare it with prior studies. RESULTS: We estimate the prevalence of SBMA among persons of Indigenous descent in Saskatchewan as 14.7 per 100,000 population. Although we believe that this is an underestimate, this still appears to be the highest population prevalence for SBMA in the world. A total of 21 participants were recruited for the haplotype study, and we identified a unique haplotype that was shared among 13 participants with Indigenous ancestry. A second shared haplotype was identified in 2 participants, which may represent a second founder haplotype, but this would need to be confirmed with future studies. CONCLUSIONS: We describe a very high prevalence of SBMA in western Canadians of Indigenous descent, which appears to predominantly be due to a founder effect. This necessitates further studies of SBMA in these populations to comprehensively ascertain the disease prevalence and allow appropriate allocation of resources to support individuals living with this chronic disease.

Characterization of a novel variant in the HR1 domain of MFN2 in a patient with ataxia, optic atrophy and sensorineural hearing loss.

Background: Pathogenic variants in MFN2 cause Charcot-Marie-Tooth disease (CMT) type 2A (CMT2A) and are the leading cause of the axonal subtypes of CMT. CMT2A is characterized by predominantly distal motor weakness and muscle atrophy, with highly variable severity and onset age. Notably, some MFN2 variants can also lead to other phenotypes such as optic atrophy, hearing loss and lipodystrophy. Despite the clear link between MFN2 and CMT2A, our mechanistic understanding of how dysfunction of the MFN2 protein causes human disease pathologies remains incomplete. This lack of understanding is due in part to the multiple cellular roles of MFN2. Though initially characterized for its role in mediating mitochondrial fusion, MFN2 also plays important roles in mediating interactions between mitochondria and other organelles, such as the endoplasmic reticulum and lipid droplets. Additionally, MFN2 is also important for mitochondrial transport, mitochondrial autophagy, and has even been implicated in lipid transfer. Though over 100 pathogenic MFN2 variants have been described to date, only a few have been characterized functionally, and even then, often only for one or two functions. Method: Several MFN2-mediated functions were characterized in fibroblast cells from a patient presenting with cerebellar ataxia, deafness, blindness, and diffuse cerebral and cerebellar atrophy, who harbours a novel homozygous MFN2 variant, D414V, which is found in a region of the HR1 domain of MFN2 where few pathogenic variants occur. Results: We found evidence for impairment of several MFN2-mediated functions. Consistent with reduced mitochondrial fusion, patient fibroblasts exhibited more fragmented mitochondrial networks and had reduced mtDNA copy number. Additionally, patient fibroblasts had reduced oxygen consumption, fewer mitochondrial-ER contacts, and altered lipid droplets that displayed an unusual perinuclear distribution. Conclusion: Overall, this work characterizes D414V as a novel variant in MFN2 and expands the phenotypic presentation of MFN2 variants to include cerebellar ataxia.

Parenting Interacts With Plasticity Genes in Predicting Behavioral Outcomes in Preschoolers.

BACKGROUND: Public health and pediatric nurses typically focus on supporting parenting to reduce the likelihood of children's behavioral problems. Studies have identified interactions between early exposures to stress in caregiving and child genotype in predicting children's behavioral problems, such that certain genotypes connote greater differential susceptibility or plasticity to environmental stressors. We sought to uncover the interaction between observational measures of parent-child relationship quality and genotype in predicting early-onset behavioral problems in 24-month-olds, using prospective methods. METHODS: We conducted a secondary analysis of data collected on a subsample of 176 women and their infants enrolled during pregnancy in the ongoing Alberta Pregnancy Outcomes and Nutrition (APrON) cohort study. Inclusion criteria required mothers to be ≥18 years of age, English speaking and ≤22 weeks gestational age at enrollment. Genetic data were obtained from blood leukocytes and buccal epithelial cell samples, collected from infants at three months of age. For each child, the presence of plasticity alleles was determined for BDNF, CNR1, DRD2/ANKK1, DRD4, DAT1, 5-HTTLPR, and MAOA and an overall index was calculated to summarize the number of plasticity alleles present. Observational assessments of parent-child relationship quality (sensitivity, controlling, and unresponsiveness) were conducted at six months of age. Children's internalizing (e.g., emotionally reactive, anxious/depressed, somatic complaint, withdrawn) and externalizing (e.g., aggression, inattention) behaviors were assessed at 24 months of age. After extracting genetic data, a maximum likelihood method for regressions was employed with Akaike Information Criterion (AIC) for model selection. RESULTS: When parents were less responsive and children possessed more plasticity alleles, children were more likely to be emotionally reactive, anxious/depressed, report somatic complaints, and withdrawn, while when parents were less responsive and children possessed fewer plasticity alleles, children were less likely to display these internalizing behaviors, in a differentially susceptible manner. Furthermore, when parents were more responsive, and children possessed more plasticity alleles, children were less likely to display internalizing behaviors (P = 0.034). Similarly, children who possessed either the CNR1-A plasticity allele (P = 0.010) or DAT1 9-repeat plasticity allele (P = 0.036) and experienced more/less parental control displayed more/fewer externalizing problems, respectively, in a differentially susceptible manner. CONCLUSIONS: The plasticity index score interacted with parental unresponsiveness in predicting anxiety and depressive behavioral problems in children, while individual genetic variants interacted with parental controlling behavior in predicting aggression and inattention in children, suggestive of differential susceptibility to caregiving. Especially in the context of nursing interventions designed to support childrearing and children's development, nurses need to be aware of the interactions between child genotype and parenting in understanding how well interventions will work in promoting optimal child behavior.

Molecular Adaptations for Sensing and Securing Prey and Insight into Amniote Genome Diversity from the Garter Snake Genome.

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.

Allele Age Under Non-Classical Assumptions is Clarified by an Exact Computational Markov Chain Approach.

Determination of the age of an allele based on its population frequency is a well-studied problem in population genetics, for which a variety of approximations have been proposed. We present a new result that, surprisingly, allows the expectation and variance of allele age to be computed exactly (within machine precision) for any finite absorbing Markov chain model in a matter of seconds. This approach makes none of the classical assumptions (e.g., weak selection, reversibility, infinite sites), exploits modern sparse linear algebra techniques, integrates over all sample paths, and is rapidly computable for Wright-Fisher populations up to N e = 100,000. With this approach, we study the joint effect of recurrent mutation, dominance, and selection, and demonstrate new examples of "selective strolls" where the classical symmetry of allele age with respect to selection is violated by weakly selected alleles that are older than neutral alleles at the same frequency. We also show evidence for a strong age imbalance, where rare deleterious alleles are expected to be substantially older than advantageous alleles observed at the same frequency when population-scaled mutation rates are large. These results highlight the under-appreciated utility of computational methods for the direct analysis of Markov chain models in population genetics.

EzMap: a simple pipeline for reproducible analysis of the human virome.

SUMMARY: In solid-organ transplant recipients, a delicate balance between immunosuppression and immunocompetence must be achieved, which can be difficult to monitor in real-time. Shotgun sequencing of cell-free DNA (cfDNA) has been recently proposed as a new way to indirectly assess immune function in transplant recipients through analysis of the status of the human virome. To facilitate exploration of the utility of the human virome as an indicator of immune status, and to enable rapid, straightforward analyses by clinicians, we developed a fully automated computational pipeline, EzMap, for performing metagenomic analysis of the human virome. EzMap combines a number of tools to clean, filter, and subtract WGS reads by mapping to a reference human assembly. The relative abundance of each virus present is estimated using a maximum likelihood approach that accounts for genome size, and results are presented with interactive visualizations and taxonomy-based summaries that enable rapid insights. The pipeline is automated to run on both workstations and computing clusters for all steps. EzMap automates an otherwise tedious and time-consuming protocol and aims to facilitate rapid and reproducible insights from cfDNA. AVAILABILITY AND IMPLEMENTATION: EzMap is freely available at https://github.com/dekoning-lab/ezmap. CONTACT: jason.dekoning@ucalgary.ca. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Wright-Fisher exact solver (WFES): scalable analysis of population genetic models without simulation or diffusion theory.

Motivation: The simplifying assumptions that are used widely in theoretical population genetics may not always be appropriate for empirical population genetics. General computational approaches that do not require the assumptions of classical theory are therefore quite desirable. One such general approach is provided by the theory of absorbing Markov chains, which can be used to obtain exact results by directly analyzing population genetic Markov models, such as the classic bi-allelic Wright-Fisher model. Although these approaches are sometimes used, they are usually forgone in favor of simulation methods, due to the perception that they are too computationally burdensome. Here we show that, surprisingly, direct analysis of virtually any Markov chain model in population genetics can be made quite efficient by exploiting transition matrix sparsity and by solving restricted systems of linear equations, allowing a wide variety of exact calculations (within machine precision) to be easily and rapidly made on modern workstation computers. Results: We introduce Wright-Fisher Exact Solver (WFES), a fast and scalable method for direct analysis of Markov chain models in population genetics. WFES can rapidly solve for both long-term and transient behaviours including fixation and extinction probabilities, expected times to fixation or extinction, sojourn times, expected allele age and variance, and others. Our implementation requires only seconds to minutes of runtime on modern workstations and scales to biological population sizes ranging from humans to model organisms. Availability and Implementation: The code is available at https://github.com/dekoning-lab/wfes. Contact: jason.dekoning@ucalgary.ca. Supplementary information: Supplementary data are available at Bioinformatics online.

The king cobra genome reveals dynamic gene evolution and adaptation in the snake venom system.

Snakes are limbless predators, and many species use venom to help overpower relatively large, agile prey. Snake venoms are complex protein mixtures encoded by several multilocus gene families that function synergistically to cause incapacitation. To examine venom evolution, we sequenced and interrogated the genome of a venomous snake, the king cobra (Ophiophagus hannah), and compared it, together with our unique transcriptome, microRNA, and proteome datasets from this species, with data from other vertebrates. In contrast to the platypus, the only other venomous vertebrate with a sequenced genome, we find that snake toxin genes evolve through several distinct co-option mechanisms and exhibit surprisingly variable levels of gene duplication and directional selection that correlate with their functional importance in prey capture. The enigmatic accessory venom gland shows a very different pattern of toxin gene expression from the main venom gland and seems to have recruited toxin-like lectin genes repeatedly for new nontoxic functions. In addition, tissue-specific microRNA analyses suggested the co-option of core genetic regulatory components of the venom secretory system from a pancreatic origin. Although the king cobra is limbless, we recovered coding sequences for all Hox genes involved in amniote limb development, with the exception of Hoxd12. Our results provide a unique view of the origin and evolution of snake venom and reveal multiple genome-level adaptive responses to natural selection in this complex biological weapon system. More generally, they provide insight into mechanisms of protein evolution under strong selection.

The Burmese python genome reveals the molecular basis for extreme adaptation in snakes.

Snakes possess many extreme morphological and physiological adaptations. Identification of the molecular basis of these traits can provide novel understanding for vertebrate biology and medicine. Here, we study snake biology using the genome sequence of the Burmese python (Python molurus bivittatus), a model of extreme physiological and metabolic adaptation. We compare the python and king cobra genomes along with genomic samples from other snakes and perform transcriptome analysis to gain insights into the extreme phenotypes of the python. We discovered rapid and massive transcriptional responses in multiple organ systems that occur on feeding and coordinate major changes in organ size and function. Intriguingly, the homologs of these genes in humans are associated with metabolism, development, and pathology. We also found that many snake metabolic genes have undergone positive selection, which together with the rapid evolution of mitochondrial proteins, provides evidence for extensive adaptive redesign of snake metabolic pathways. Additional evidence for molecular adaptation and gene family expansions and contractions is associated with major physiological and phenotypic adaptations in snakes; genes involved are related to cell cycle, development, lungs, eyes, heart, intestine, and skeletal structure, including GRB2-associated binding protein 1, SSH, WNT16, and bone morphogenetic protein 7. Finally, changes in repetitive DNA content, guanine-cytosine isochore structure, and nucleotide substitution rates indicate major shifts in the structure and evolution of snake genomes compared with other amniotes. Phenotypic and physiological novelty in snakes seems to be driven by system-wide coordination of protein adaptation, gene expression, and changes in the structure of the genome.

The western painted turtle genome, a model for the evolution of extreme physiological adaptations in a slowly evolving lineage.

BACKGROUND: We describe the genome of the western painted turtle, Chrysemys picta bellii, one of the most widespread, abundant, and well-studied turtles. We place the genome into a comparative evolutionary context, and focus on genomic features associated with tooth loss, immune function, longevity, sex differentiation and determination, and the species' physiological capacities to withstand extreme anoxia and tissue freezing. RESULTS: Our phylogenetic analyses confirm that turtles are the sister group to living archosaurs, and demonstrate an extraordinarily slow rate of sequence evolution in the painted turtle. The ability of the painted turtle to withstand complete anoxia and partial freezing appears to be associated with common vertebrate gene networks, and we identify candidate genes for future functional analyses. Tooth loss shares a common pattern of pseudogenization and degradation of tooth-specific genes with birds, although the rate of accumulation of mutations is much slower in the painted turtle. Genes associated with sex differentiation generally reflect phylogeny rather than convergence in sex determination functionality. Among gene families that demonstrate exceptional expansions or show signatures of strong natural selection, immune function and musculoskeletal patterning genes are consistently over-represented. CONCLUSIONS: Our comparative genomic analyses indicate that common vertebrate regulatory networks, some of which have analogs in human diseases, are often involved in the western painted turtle's extraordinary physiological capacities. As these regulatory pathways are analyzed at the functional level, the painted turtle may offer important insights into the management of a number of human health disorders.

Thousands of microsatellite loci from the venomous coralsnake Micrurus fulvius and variability of select loci across populations and related species.

Studies of population genetics increasingly use next-generation DNA sequencing to identify microsatellite loci in nonmodel organisms. There are, however, relatively few studies that validate the feasibility of transitioning from marker development to experimental application across populations and species. North American coralsnakes of the Micrurus fulvius species complex occur in the United States and Mexico, and little is known about their population structure and phylogenetic relationships. This absence of information and population genetics markers is particularly concerning because they are highly venomous and have important implications on human health. To alleviate this problem in coralsnakes, we investigated the feasibility of using 454 shotgun sequences for microsatellite marker development. First, a genomic shotgun library from a single individual was sequenced (approximately 7.74 megabases; 26,831 reads) to identify potentially amplifiable microsatellite loci (PALs). We then hierarchically sampled 76 individuals from throughout the geographic distribution of the species complex and examined whether PALs were amplifiable and polymorphic. Approximately half of the loci tested were readily amplifiable from all individuals, and 80% of the loci tested for variation were variable and thus informative as population genetic markers. To evaluate the repetitive landscape characteristics across multiple snakes, we also compared microsatellite content between the coralsnake and two other previously sampled snakes, the venomous copperhead (Agkistrodon contortrix) and Burmese python (Python molurus).

Phylogenetics, likelihood, evolution and complexity.

SUMMARY: Phylogenetics, likelihood, evolution and complexity (PLEX) is a flexible and fast Bayesian Markov chain Monte Carlo software program for large-scale analysis of nucleotide and amino acid data using complex evolutionary models in a phylogenetic framework. The program gains large speed improvements over standard approaches by implementing 'partial sampling of substitution histories', a data augmentation approach that can reduce data analysis times from months to minutes on large comparative datasets. A variety of nucleotide and amino acid substitution models are currently implemented, including non-reversible and site-heterogeneous mixture models. Due to efficient algorithms that scale well with data size and model complexity, PLEX can be used to make inferences from hundreds to thousands of taxa in only minutes on a desktop computer. It also performs probabilistic ancestral sequence reconstruction. Future versions will support detection of co-evolutionary interactions between sites, probabilistic tests of convergent evolution and rigorous testing of evolutionary hypotheses in a Bayesian framework. AVAILABILITY AND IMPLEMENTATION: PLEX v1.0 is licensed under GPL. Source code and documentation will be available for download at www.evolutionarygenomics.com/ProgramsData/PLEX. PLEX is implemented in C++ and supported on Linux, Mac OS X and other platforms supporting standard C++ compilers. Example data, control files, documentation and accessory Perl scripts are available from the website. CONTACT: David.Pollock@UCDenver.edu. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The interface of protein structure, protein biophysics, and molecular evolution.

Abstract The interface of protein structural biology, protein biophysics, molecular evolution, and molecular population genetics forms the foundations for a mechanistic understanding of many aspects of protein biochemistry. Current efforts in interdisciplinary protein modeling are in their infancy and the state-of-the art of such models is described. Beyond the relationship between amino acid substitution and static protein structure, protein function, and corresponding organismal fitness, other considerations are also discussed. More complex mutational processes such as insertion and deletion and domain rearrangements and even circular permutations should be evaluated. The role of intrinsically disordered proteins is still controversial, but may be increasingly important to consider. Protein geometry and protein dynamics as a deviation from static considerations of protein structure are also important. Protein expression level is known to be a major determinant of evolutionary rate and several considerations including selection at the mRNA level and the role of interaction specificity are discussed. Lastly, the relationship between modeling and needed high-throughput experimental data as well as experimental examination of protein evolution using ancestral sequence resurrection and in vitro biochemistry are presented, towards an aim of ultimately generating better models for biological inference and prediction.

Rapid microsatellite identification from Illumina paired-end genomic sequencing in two birds and a snake.

Identification of microsatellites, or simple sequence repeats (SSRs), can be a time-consuming and costly investment requiring enrichment, cloning, and sequencing of candidate loci. Recently, however, high throughput sequencing (with or without prior enrichment for specific SSR loci) has been utilized to identify SSR loci. The direct "Seq-to-SSR" approach has an advantage over enrichment-based strategies in that it does not require a priori selection of particular motifs, or prior knowledge of genomic SSR content. It has been more expensive per SSR locus recovered, however, particularly for genomes with few SSR loci, such as bird genomes. The longer but relatively more expensive 454 reads have been preferred over less expensive Illumina reads. Here, we use Illumina paired-end sequence data to identify potentially amplifiable SSR loci (PALs) from a snake (the Burmese python, Python molurus bivittatus), and directly compare these results to those from 454 data. We also compare the python results to results from Illumina sequencing of two bird genomes (Gunnison Sage-grouse, Centrocercus minimus, and Clark's Nutcracker, Nucifraga columbiana), which have considerably fewer SSRs than the python. We show that direct Illumina Seq-to-SSR can identify and characterize thousands of potentially amplifiable SSR loci for as little as $10 per sample--a fraction of the cost of 454 sequencing. Given that Illumina Seq-to-SSR is effective, inexpensive, and reliable even for species such as birds that have few SSR loci, it seems that there are now few situations for which prior hybridization is justifiable.

Repetitive elements may comprise over two-thirds of the human genome.

Transposable elements (TEs) are conventionally identified in eukaryotic genomes by alignment to consensus element sequences. Using this approach, about half of the human genome has been previously identified as TEs and low-complexity repeats. We recently developed a highly sensitive alternative de novo strategy, P-clouds, that instead searches for clusters of high-abundance oligonucleotides that are related in sequence space (oligo "clouds"). We show here that P-clouds predicts >840 Mbp of additional repetitive sequences in the human genome, thus suggesting that 66%-69% of the human genome is repetitive or repeat-derived. To investigate this remarkable difference, we conducted detailed analyses of the ability of both P-clouds and a commonly used conventional approach, RepeatMasker (RM), to detect different sized fragments of the highly abundant human Alu and MIR SINEs. RM can have surprisingly low sensitivity for even moderately long fragments, in contrast to P-clouds, which has good sensitivity down to small fragment sizes (∼25 bp). Although short fragments have a high intrinsic probability of being false positives, we performed a probabilistic annotation that reflects this fact. We further developed "element-specific" P-clouds (ESPs) to identify novel Alu and MIR SINE elements, and using it we identified ∼100 Mb of previously unannotated human elements. ESP estimates of new MIR sequences are in good agreement with RM-based predictions of the amount that RM missed. These results highlight the need for combined, probabilistic genome annotation approaches and suggest that the human genome consists of substantially more repetitive sequence than previously believed.

Bayesian analysis of high-throughput quantitative measurement of protein-DNA interactions.

Transcriptional regulation depends upon the binding of transcription factor (TF) proteins to DNA in a sequence-dependent manner. Although many experimental methods address the interaction between DNA and proteins, they generally do not comprehensively and accurately assess the full binding repertoire (the complete set of sequences that might be bound with at least moderate strength). Here, we develop and evaluate through simulation an experimental approach that allows simultaneous high-throughput quantitative analysis of TF binding affinity to thousands of potential DNA ligands. Tens of thousands of putative binding targets can be mixed with a TF, and both the pre-bound and bound target pools sequenced. A hierarchical Bayesian Markov chain Monte Carlo approach determines posterior estimates for the dissociation constants, sequence-specific binding energies, and free TF concentrations. A unique feature of our approach is that dissociation constants are jointly estimated from their inferred degree of binding and from a model of binding energetics, depending on how many sequence reads are available and the explanatory power of the energy model. Careful experimental design is necessary to obtain accurate results over a wide range of dissociation constants. This approach, which we call Simultaneous Ultra high-throughput Ligand Dissociation EXperiment (SULDEX), is theoretically capable of rapid and accurate elucidation of an entire TF-binding repertoire.

Discovery of highly divergent repeat landscapes in snake genomes using high-throughput sequencing.

We conducted a comprehensive assessment of genomic repeat content in two snake genomes, the venomous copperhead (Agkistrodon contortrix) and the Burmese python (Python molurus bivittatus). These two genomes are both relatively small (∼1.4 Gb) but have surprisingly extensive differences in the abundance and expansion histories of their repeat elements. In the python, the readily identifiable repeat element content is low (21%), similar to bird genomes, whereas that of the copperhead is higher (45%), similar to mammalian genomes. The copperhead's greater repeat content arises from the recent expansion of many different microsatellites and transposable element (TE) families, and the copperhead had 23-fold greater levels of TE-related transcripts than the python. This suggests the possibility that greater TE activity in the copperhead is ongoing. Expansion of CR1 LINEs in the copperhead genome has resulted in TE-mediated microsatellite expansion ("microsatellite seeding") at a scale several orders of magnitude greater than previously observed in vertebrates. Snakes also appear to be prone to horizontal transfer of TEs, particularly in the copperhead lineage. The reason that the copperhead has such a small genome in the face of so much recent expansion of repeat elements remains an open question, although selective pressure related to extreme metabolic performance is an obvious candidate. TE activity can affect gene regulation as well as rates of recombination and gene duplication, and it is therefore possible that TE activity played a role in the evolution of major adaptations in snakes; some evidence suggests this may include the evolution of venom repertoires.

Adaptive molecular convergence: Molecular evolution versus molecular phylogenetics.

Definitive identification of convergent evolution, the acquisition of the same biological trait in unrelated lineages, provides one of the most compelling sources of evidence for natural selection. Although numerous examples of convergent morphological evolution are well known (such as the independent development of wings in birds and mammals), cases of convergent evolution at the molecular-genetic level appear to be quite rare. We recently discovered a remarkable case of convergent molecular evolution involving more than 100 parallel amino-acid changes across all 13 mitochondrially-encoded proteins of snakes and agamid lizards. Just a few of these convergent substitutions were sufficient to positively mislead the inference of phylogeny, even with thousands of sites providing latent support for the correct underlying relationships. Since this example demonstrates that molecular convergence can happen en masse in nature, affecting multiple genes, it is important to consider the threat this poses to molecular systematics, and careful genome-wide assays for convergent molecular evolution are warranted. This result implies that the protein adaptive landscape is sometimes highly constrained.