Recent Human Genome Research Articles

ELOA3: A primate-specific RNA polymerase II elongation factor encoded by a tandem repeat gene cluster.

The biological role of the repetitive DNA sequences in the human genome remains an outstanding question. Recent long-read human genome assemblies have allowed us to identify a function for one of these repetitive regions. We have uncovered a tandem array of conserved primate-specific retrogenes encoding the protein Elongin A3 (ELOA3), a homolog of the RNA polymerase II (RNAPII) elongation factor Elongin A (ELOA). Our genomic analysis shows that the ELOA3 gene cluster is conserved among primates and the number of ELOA3 gene repeats is variable in the human population and across primate species. Moreover, the gene cluster has undergone concerted evolution and homogenization within primates. Our biochemical studies show that ELOA3 functions as a promoter-associated RNAPII pause-release elongation factor with distinct biochemical and functional features from its ancestral homolog, ELOA. We propose that the ELOA3 gene cluster has evolved to fulfil a transcriptional regulatory function unique to the primate lineage that can be targeted to regulate cellular hyperproliferation.

AlcoR: alignment-free simulation, mapping, and visualization of low-complexity regions in biological data.

Low-complexity data analysis is the area that addresses the search and quantification of regions in sequences of elements that contain low-complexity or repetitive elements. For example, these can be tandem repeats, inverted repeats, homopolymer tails, GC-biased regions, similar genes, and hairpins, among many others. Identifying these regions is crucial because of their association with regulatory and structural characteristics. Moreover, their identification provides positional and quantity information where standard assembly methodologies face significant difficulties because of substantial higher depth coverage (mountains), ambiguous read mapping, or where sequencing or reconstruction defects may occur. However, the capability to distinguish low-complexity regions (LCRs) in genomic and proteomic sequences is a challenge that depends on the model's ability to find them automatically. Low-complexity patterns can be implicit through specific or combined sources, such as algorithmic or probabilistic, and recurring to different spatial distances-namely, local, medium, or distant associations. This article addresses the challenge of automatically modeling and distinguishing LCRs, providing a new method and tool (AlcoR) for efficient and accurate segmentation and visualization of these regions in genomic and proteomic sequences. The method enables the use of models with different memories, providing the ability to distinguish local from distant low-complexity patterns. The method is reference and alignment free, providing additional methodologies for testing, including a highly flexible simulation method for generating biological sequences (DNA or protein) with different complexity levels, sequence masking, and a visualization tool for automatic computation of the LCR maps into an ideogram style. We provide illustrative demonstrations using synthetic, nearly synthetic, and natural sequences showing the high efficiency and accuracy of AlcoR. As large-scale results, we use AlcoR to unprecedentedly provide a whole-chromosome low-complexity map of a recent complete human genome and the haplotype-resolved chromosome pairs of a heterozygous diploid African cassava cultivar. The AlcoR method provides the ability of fast sequence characterization through data complexity analysis, ideally for scenarios entangling the presence of new or unknown sequences. AlcoR is implemented in C language using multithreading to increase the computational speed, is flexible for multiple applications, and does not contain external dependencies. The tool accepts any sequence in FASTA format. The source code is freely provided at https://github.com/cobilab/alcor.

Discovery of a MUC3B gene reconstructs the membrane mucin gene cluster on human chromosome 7.

Human tissue surfaces are coated with mucins, a family of macromolecular sugar-laden proteins serving diverse functions from lubrication to the formation of selective biochemical barriers against harmful microorganisms and molecules. Membrane mucins are a distinct group of mucins that are attached to epithelial cell surfaces where they create a dense glycocalyx facing the extracellular environment. All mucin proteins carry long stretches of tandemly repeated sequences that undergo extensive O-linked glycosylation to form linear mucin domains. However, the repetitive nature of mucin domains makes them prone to recombination and renders their genetic sequences particularly difficult to read with standard sequencing technologies. As a result, human mucin genes suffer from significant sequence gaps that have hampered the investigation of gene function in health and disease. Here we leveraged a recent human genome assembly to characterize a previously unmapped MUC3B gene located at the q22 locus on chromosome 7, within a cluster of four structurally related membrane mucin genes that we name the MUC3 cluster. We found that MUC3B shares high sequence identity with the known MUC3A gene and that the two genes are governed by evolutionarily conserved regulatory elements. Furthermore, we show that MUC3A, MUC3B, MUC12, and MUC17 in the human MUC3 cluster are expressed in intestinal epithelial cells (IECs). Our results complete existing genetic gaps in the MUC3 cluster which is a conserved genetic unit in vertebrates. We anticipate our results to be the starting point for the detection of disease-associated polymorphisms in the human MUC3 cluster. Moreover, our study provides the basis for the exploration of intestinal mucin gene function in widely used experimental models such as human intestinal organoids and genetic mouse models.

Potential Therapeutic Approaches through Modulating the Autophagy Process for Skin Barrier Dysfunction.

Autophagy is an attractive process to researchers who are seeking novel potential treatments for various diseases. Autophagy plays a critical role in degrading damaged cellular organelles, supporting normal cell development, and maintaining cellular homeostasis. Because of the various effects of autophagy, recent human genome research has focused on evaluating the relationship between autophagy and a wide variety of diseases, such as autoimmune diseases, cancers, and inflammatory diseases. The skin is the largest organ in the body and provides the first line of defense against environmental hazards, including UV damage, chemical toxins, injuries, oxidative stress, and microorganisms. Autophagy takes part in endogenous defense mechanisms by controlling skin homeostasis. In this manner, regulating autophagy might contribute to the treatment of skin barrier dysfunctions. Various studies are ongoing to elucidate the association between autophagy and skin-related diseases in order to find potential therapeutic approaches. However, little evidence has been gathered about the relationship between autophagy and the skin. In this review, we highlight the previous findings of autophagy and skin barrier disorders and suggest potential therapeutic strategies. The recent research regarding autophagy in acne and skin aging is also discussed.

Next-Generation Sequencing Enables Spatiotemporal Resolution of Human Centromere Replication Timing.

Centromeres serve a critical function in preserving genome integrity across sequential cell divisions, by mediating symmetric chromosome segregation. The repetitive, heterochromatic nature of centromeres is thought to be inhibitory to DNA replication, but has also led to their underrepresentation in human reference genome assemblies. Consequently, centromeres have been excluded from genomic replication timing analyses, leaving their time of replication unresolved. However, the most recent human reference genome, hg38, included models of centromere sequences. To establish the experimental requirements for achieving replication timing profiles for centromeres, we sequenced G1- and S-phase cells from five human cell lines, and aligned the sequence reads to hg38. We were able to infer DNA replication timing profiles for the centromeres in each of the five cell lines, which showed that centromere replication occurs in mid-to-late S phase. Furthermore, we found that replication timing was more variable between cell lines in the centromere regions than expected, given the distribution of variation in replication timing genome-wide. These results suggest the potential of these, and future, sequence models to enable high-resolution studies of replication in centromeres and other heterochromatic regions.

Similarities and differences between variants called with human reference genome HG19 or HG38

BackgroundReference genome selection is a prerequisite for successful analysis of next generation sequencing (NGS) data. Current practice employs one of the two most recent human reference genome versions: HG19 or HG38. To date, the impact of genome version on SNV identification has not been rigorously assessed.MethodsWe conducted analysis comparing the SNVs identified based on HG19 vs HG38, leveraging whole genome sequencing (WGS) data from the genome-in-a-bottle (GIAB) project. First, SNVs were called using 26 different bioinformatics pipelines with either HG19 or HG38. Next, two tools were used to convert the called SNVs between HG19 and HG38. Lastly we calculated conversion rates, analyzed discordant rates between SNVs called with HG19 or HG38, and characterized the discordant SNVs.ResultsThe conversion rates from HG38 to HG19 (average 95%) were lower than the conversion rates from HG19 to HG38 (average 99%). The conversion rates varied slightly among the various calling pipelines. Around 1.5% SNVs were discordantly converted between HG19 or HG38. The conversions from HG38 to HG19 had more SNVs which failed conversion and more discordant SNVs than the opposite conversion (HG19 to HG38). Most of the discordant SNVs had low read depth, were low confidence SNVs as defined by GIAB, and/or were predominated by G/C alleles (52% observed versus 42% expected).ConclusionA significant number of SNVs could not be converted between HG19 and HG38. Based on careful review of our comparisons, we recommend HG38 (the newer version) for NGS SNV analysis. To summarize, our findings suggest caution when translating identified SNVs between different versions of the human reference genome.

Systematic evaluation of the early access applied biosystems precision ID Globalfiler mixture ID and Globalfiler NGS STR panels for the ion S5 system

Systematic evaluation of the early access applied biosystems precision ID Globalfiler mixture ID and Globalfiler NGS STR panels for the ion S5 system

SNPnexus: assessing the functional relevance of genetic variation to facilitate the promise of precision medicine.

Broader functional annotation of genetic variation is a valuable means for prioritising phenotypically-important variants in further disease studies and large-scale genotyping projects. We developed SNPnexus to meet this need by assessing the potential significance of known and novel SNPs on the major transcriptome, proteome, regulatory and structural variation models. Since its previous release in 2012, we have made significant improvements to the annotation categories and updated the query and data viewing systems. The most notable changes include broader functional annotation of noncoding variants and expanding annotations to the most recent human genome assembly GRCh38/hg38. SNPnexus has now integrated rich resources from ENCODE and Roadmap Epigenomics Consortium to map and annotate the noncoding variants onto different classes of regulatory regions and noncoding RNAs as well as providing their predicted functional impact from eight popular non-coding variant scoring algorithms and computational methods. A novel functionality offered now is the support for neo-epitope predictions from leading tools to facilitate its use in immunotherapeutic applications. These updates to SNPnexus are in preparation for its future expansion towards a fully comprehensive computational workflow for disease-associated variant prioritization from sequencing data, placing its users at the forefront of translational research. SNPnexus is freely available at http://www.snp-nexus.org.

SWI/SNF factors required for cellular resistance to DNA damage include ARID1A and ARID1B and show interdependent protein stability.

The SWI/SNF chromatin-remodeling family contains various protein complexes, which regulate gene expression during cellular development and influence DNA damage response in an ATP- and complex-dependent manner, of which details remain elusive. Recent human genome sequencing of various cancer cells revealed frequent mutations in SWI/SNF factors, especially ARID1A, a variant subunit in the BRG1-associated factor (BAF) complex of the SWI/SNF family. We combined live-cell analysis and gene-suppression experiments to show that suppression of either ARID1A or its paralog ARID1B led to reduced nonhomologous end joining activity of DNA double-strand breaks (DSB), decreased accumulation of KU70/KU80 proteins at DSB, and sensitivity to ionizing radiation, as well as to cisplatin and UV. Thus, in contrast to transcriptional regulation, both ARID1 proteins are required for cellular resistance to various types of DNA damage, including DSB. The suppression of other SWI/SNF factors, namely SNF5, BAF60a, BAF60c, BAF155, or BAF170, exhibits a similar phenotype. Of these factors, ARID1A, ARID1B, SNF5, and BAF60c are necessary for the immediate recruitment of the ATPase subunit of the SWI/SNF complex to DSB, arguing that both ARID1 proteins facilitate the damage response of the complex. Finally, we found interdependent protein stability among the SWI/SNF factors, suggesting their direct interaction within the complex and the reason why multiple factors are frequently lost in parallel in cancer cells. Taken together, we show that cancer cells lacking in the expression of certain SWI/SNF factors, including ARID1A, are deficient in DNA repair and potentially vulnerable to DNA damage.

HIV infection reveals widespread expansion of novel centromeric human endogenous retroviruses.

Human endogenous retroviruses (HERVs) make up 8% of the human genome. The HERV-K (HML-2) family is the most recent group of these viruses to have inserted into the genome, and we have detected the activation of HERV-K (HML-2) proviruses in the blood of patients with HIV-1 infection. We report that HIV-1 infection activates expression of a novel HERV-K (HML-2) provirus, termed K111, present in multiple copies in the centromeres of chromosomes throughout the human genome yet not annotated in the most recent human genome assembly. Infection with HIV-1 or stimulation with the HIV-1 Tat protein leads to the activation of K111 proviruses. K111 is present as a single copy in the genome of the chimpanzee, yet K111 is not found in the genomes of other primates. Remarkably, K111 proviruses appear in the genomes of the extinct Neanderthal and Denisovan, while modern humans have at least 100 K111 proviruses spread across the centromeres of 15 chromosomes. Our studies suggest that the progenitor K111 integrated before the Homo-Pan divergence and expanded in copy number during the evolution of hominins, perhaps by recombination. The expansion of K111 provides sequence evidence suggesting that recombination between the centromeres of various chromosomes took place during the evolution of humans. K111 proviruses show significant sequence variations in each individual centromere, which may serve as markers in future efforts to annotate human centromere sequences. Further, this work is an example of the potential to discover previously unknown genomic sequences through the analysis of nucleic acids found in the blood of patients.

Single Nucleotide Variant Increases the open Probability of Human ENaC

The epithelial Na+ channel (ENaC) mediates Na+ entry into specific epithelial cells and has an essential role in the regulation of blood pressure and airway surface liquid volume. Rare ENaC mutations have been reported in inherited disorders associated with hypertension or hypotension, and ENaC gene variants have been associated with several human disorders. Recent human genome sequencing projects have revealed a large number of ENaC gene variations in both exons and introns. However, the functional consequences of most variants are unknown. We examined the functional properties of the human ENaC variant γL511Q in the Xenopus oocyte expression system. Oocytes expressing αβγL511Q exhibited five-fold greater amiloride-sensitive currents than cells expressing WT channels. Mutant and WT channels had similar levels of surface expression. Single channel recordings with a cell-attached patch revealed that mutant channels had a four-fold higher open probability than WT, but similar unitary currents. The mutant had a significantly reduced Na+ self-inhibition response, reflecting less reduction of ENaC open probability by extracellular Na+. Interestingly, the mutant diminished the activating effect of external Zn2+ and essentially converted Zn2+ from a high-affinity ENaC activator to a low-affinity inhibitor. Furthermore, γL511Q exhibited blunted activation by chymotrypsin. We conclude that γL511Q is a gain-of-function human ENaC variant that is characterized by increased open probability and suppressed responses to extracellular Na+, Zn2+ and chymotrypsin.

Gain-of-function variant of the human epithelial sodium channel

Epithelial Na(+) channel (ENaC) mutations are associated with several human disorders, underscoring the importance of these channels in human health. Recent human genome sequencing projects have revealed a large number of ENaC gene variations, several of which have been found in individuals with salt-sensitive hypertension, cystic fibrosis, and other disorders. However, the functional consequences of most variants are unknown. In this study, we used the Xenopus oocyte expression system to examine the functional properties of a human ENaC variant. Oocytes expressing αβγL511Q human ENaCs showed 4.6-fold greater amiloride-sensitive currents than cells expressing wild-type channels. The γL511Q variant did not significantly alter channel surface expression. Single channel recordings revealed that the variant had fourfold higher open probability than wild type. In addition, γL511Q largely eliminated the Na(+) self-inhibition response, which reflects a downregulation of ENaC open probability by extracellular Na(+). Moreover, γL511Q diminished chymotrypsin-induced activation of the mutant channel. We conclude that γL511Q is a gain-of-function human ENaC variant. Our results suggest that γL511Q enhances ENaC activity by increasing channel open probability and dampens channel regulation by extracellular Na(+) and proteases.

Abstract 9368: Plasma Membrane Calcium ATPase Isoform 1 Protects the Heart from Cardiomyopathy, Arrhythmias and Hypertension

With 1 in 3 people in the USA dying as a result of cardiovascular disease, and the incidence of these diseases rapidly increasing in the developing world, it is essential that we continue to provide mechanistic insights into their development and progression. Several recent human genome wide association studies have identified single nucleotide polymorphisms in the plasma membrane calcium ATPase isoform 1 (PMCA1) gene as the strongest association with hypertension and blood pressure variance; hypertension along with cardiomyopathy and arrhythmias are central features of heart failure. In this study we provide evidence that PMCA1, which is a transmembrane calcium transporter known to eject calcium from the cell, is involved in all of these key features of heart failure. To elucidate the role of PMCA1 in cardiac function we generated mice carrying a cardiomyocyte-specific gene deletion of PMCA1 (PMCA1 cko ). In response to pressure overload hypertrophic stress induced by transverse aortic constriction (TAC) PMCA1 cko mice displayed accelerated cardiac dysfunction shown by significantly reduced fractional shortening compared to controls (FS; PMCA1 cko /TAC: 21±2%, PMCA1 f/f /TAC 32±3%, PMCA1 cko /Sham and [alpha]-MHC-Cre/TAC 35±2%, p <0.05, n=6). This was accompanied by increased lung weight/body weight in PMCA1 cko /TAC (11.5±1 vs 7.8±0.8 mg/g in PMCA1 f/f /TAC, p <0.05, n=6), greater cardiomyocyte cross-sectional area and significantly elevated interstitial fibrosis. Our data shows that in the myocardium PMCA1 has an additional role in stabilising cardiac rhythm. At 3 months of age PMCA1 cko displayed episodes of spontaneous ventricular tachycardia whilst exhibiting an 80% downregulation in the expression of the voltage gated sodium channel Na v 1.5 and a significant reduction in the corresponding current. To investigate the relationship between PMCA1 and blood pressure we generated PMCA1 heterozygous knockout mice (PMCA1 Ht ) in which both systolic and diastolic carotid pressure was significantly increased in PMCA1 Ht in comparison to wild type littermates. In conclusion, this data addresses the mechanistic basis for the association of PMCA1 with hypertension, cardiomyopathy and arrhythmias which are key contributors to heart failure.

Mutation nomenclature in practice: Findings and recommendations from the cystic fibrosis external quality assessment scheme

Currently, two nomenclature systems are in use to describe sequence variants for cystic fibrosis: the established traditional nomenclature system and the more recent Human Genome Variation Society (HGVS) nomenclature system. We have evaluated the use of both systems in the laboratory reports of 217 participants in the cystic fibrosis external quality assessment scheme of 2009. The mutation c.1521_1523delCTT (p.Phe508del, F508del) was described by traditional and HGVS nomenclature by 32 of 216 (15%) laboratories that correctly identified the mutation, whereas 171 (79%) laboratories used traditional nomenclature only and 13 (6%) laboratories used HGVS nomenclature only. Overall, 29 of 631 (5%) reports used nomenclature that was evaluated as being seriously incorrect and/or misleading and 136 (22%) reports contained attempts at HGVS coding, of which 104 (76%) contained no coding errors; just 33 (24%) mentioned the correct cDNA name and cited the nucleotide reference sequence. We recognized an urgent need for more consistent and correct usage of nomenclature. We recommended that cystic fibrosis transmembrane conductance regulator testing reports should include a description of the identified sequence variants in both HGVS and traditional nomenclature and provided basic recommendations and other guidance.

Localization of genes for V+LDL plasma cholesterol levels on two diets in the opossum Monodelphis domestica

Plasma cholesterol levels among individuals vary considerably in response to diet. However, the genes that influence this response are largely unknown. Non-HDL (V+LDL) cholesterol levels vary dramatically among gray, short-tailed opossums fed an atherogenic diet, and we previously reported that two quantitative trait loci (QTLs) influenced V+LDL cholesterol on two diets. We used hypothesis-free, genome-wide linkage analyses on data from 325 pedigreed opossums and located one QTL for V+LDL cholesterol on the basal diet on opossum chromosome 1q [logarithm of the odds (LOD) = 3.11, genomic P = 0.019] and another QTL for V+LDL on the atherogenic diet (i.e., high levels of cholesterol and fat) on chromosome 8 (LOD = 9.88, genomic P = 5 x 10(-9)). We then employed a novel strategy involving combined analyses of genomic resources, expression analysis, sequencing, and genotyping to identify candidate genes for the chromosome 8 QTL. A polymorphism in ABCB4 was strongly associated (P = 9 x 10(-14)) with the plasma V+LDL cholesterol concentrations on the high-cholesterol, high-fat diet. The results of this study indicate that genetic variation in ABCB4, or closely linked genes, is responsible for the dramatic differences among opossums in their V+LDL cholesterol response to an atherogenic diet.

Pathogenesis of Hepatitis E Virus and Hepatitis C Virus in Chimpanzees: Similarities and Differences

The chimpanzee is the only animal model for investigating the pathogenesis of viral hepatitis types A through E in humans. Studies of the host response, including microarray analyses, have relied on the close relationship between these two primate species: chimpanzee samples are commonly tested with human-based reagents. In this study, the host responses to two dissimilar viruses, hepatitis E virus (HEV) and hepatitis C virus (HCV), were compared in multiple experimentally infected chimpanzees. Affymetrix U133+2.0 human microarray chips were used to assess the entire transcriptome in serial liver biopsies obtained over the course of the infections. Respecting the limitations of microarray probes designed for human target transcripts to effectively assay chimpanzee transcripts, we conducted probe-level analysis of the microarray data in conjunction with a custom mapping of the probe sequences to the most recent human and chimpanzee genome sequences. Time points for statistical comparison were chosen based on independently measured viremia levels. Regardless of the viral infection, the alignment of differentially expressed genes to the human genome sequence resulted in a larger number of genes being identified when compared with alignment to the chimpanzee genome sequence. This probably reflects the lesser refinement of gene annotation for chimpanzees. In general, the two viruses demonstrated very distinct temporal changes in host response genes, although both RNA viruses induced genes that were involved in many of the same biological systems, including interferon-induced genes. The host response to HCV infection was more robust in the magnitude and number of differentially expressed genes compared to HEV infection.

Evolution of vertebrate opioid receptors: evidence from cloning and bioinformatics

The proteins that mediate the analgesic and other effects of opioid drugs and endogenous opioid peptides are known as opioid receptors. Opioid receptors consist of a family of four closely‐related proteins belonging to the large superfamily of G‐protein coupled receptors. The three types of opioid receptors shown unequivocally to mediate analgesia in animal models are the mu (MOR), delta (DOR), and kappa (KOR) opioid receptor proteins. The role of the fourth and newest member of the opioid receptor family, the nociceptin or orphanin FQ receptor (ORL), is not as clear as hyperalgesia, analgesia, and no effect was reported after administration of ORL agonists. There are now cDNA sequences for all four types of opioid receptors that are expressed in the brain of six species from three different classes of vertebrates. This paper presents a comparative analysis of vertebrate opioid receptors using cloning studies from our lab, bioinformatic analyses, and data from recent human genome studies. Results indicate that it is likely that opioid receptors arose by two whole genome duplication events early in vertebrate evolution, that there is a vector of opioid receptor divergence such that receptors in mammals and humans exhibit a higher degree of opioid type‐selectivity, and that the MOR protein shows evidence of the most rapid adaptive evolution. Supported by a National Institutes of Health grant, DA0124482.

Development and evaluation of new mask protocols for gene expression profiling in humans and chimpanzees

BackgroundCross-species gene expression analyses using oligonucleotide microarrays designed to evaluate a single species can provide spurious results due to mismatches between the interrogated transcriptome and arrayed probes. Based on the most recent human and chimpanzee genome assemblies, we developed updated and accessible probe masking methods that allow human Affymetrix oligonucleotide microarrays to be used for robust genome-wide expression analyses in both species. In this process, only data from oligonucleotide probes predicted to have robust hybridization sensitivity and specificity for both transcriptomes are retained for analysis.ResultsTo characterize the utility of this resource, we applied our mask protocols to existing expression data from brains, livers, hearts, testes, and kidneys derived from both species and determined the effects probe numbers have on expression scores of specific transcripts. In all five tissues, probe sets with decreasing numbers of probes showed non-linear trends towards increased variation in expression scores. The relationships between expression variation and probe number in brain data closely matched those observed in simulated expression data sets subjected to random probe masking. However, there is evidence that additional factors affect the observed relationships between gene expression scores and probe number in tissues such as liver and kidney. In parallel, we observed that decreasing the number of probes within probe sets lead to linear increases in both gained and lost inferences of differential cross-species expression in all five tissues, which will affect the interpretation of expression data subject to masking.ConclusionWe introduce a readily implemented and updated resource for human and chimpanzee transcriptome analysis through a commonly used microarray platform. Based on empirical observations derived from the analysis of five distinct data sets, we provide novel guidelines for the interpretation of masked data that take the number of probes present in a given probe set into consideration. These guidelines are applicable to other customized applications that involve masking data from specific subsets of probes.

The evolution of vertebrate opioid receptors.

The proteins that mediate the analgesic and other effects of opioid drugs and endogenous opioid peptides are known as opioid receptors. Opioid receptors consist of a family of four closely-related proteins belonging to the large superfamily of G-protein coupled receptors. The three types of opioid receptors shown unequivocally to mediate analgesia in animal models are the mu (MOR), delta (DOR), and kappa (KOR) opioid receptor proteins. The role of the fourth member of the opioid receptor family, the nociceptin or orphanin FQ receptor (ORL), is not as clear as hyperalgesia, analgesia, and no effect was reported after administration of ORL agonists. There are now cDNA sequences for all four types of opioid receptors that are expressed in the brain of six species from three different classes of vertebrates. This review presents a comparative analysis of vertebrate opioid receptors using bioinformatics and data from recent human genome studies. Results indicate that opioid receptors arose by gene duplication, that there is a vector of opioid receptor divergence, and that MOR shows evidence of rapid evolution.

The population genetics of structural variation

Population genetics is central to our understanding of human variation, and by linking medical and evolutionary themes, it enables us to understand the origins and impacts of our genomic differences. Despite current limitations in our knowledge of the locations, sizes and mutational origins of structural variants, our characterization of their population genetics is developing apace, bringing new insights into recent human adaptation, genome biology and disease. We summarize recent dramatic advances, describe the diverse mutational origins of chromosomal rearrangements and argue that their complexity necessitates a re-evaluation of existing population genetic methods.