The contribution of common regulatory and protein-coding TYR variants to the genetic architecture of albinism.

Genetic diseases have been historically segregated into rare Mendelian disorders and common complex conditions. Large-scale studies using genome sequencing are eroding this distinction and are gradually unmasking the underlying complexity of human traits. Here, we analysed data from the Genomics England 100,000 Genomes Project and from a cohort of 1313 individuals with albinism aiming to gain insights into the genetic architecture of this archetypal rare disorder. We investigated the contribution of protein-coding and regulatory variants both rare and common. We focused on TYR, the gene encoding tyrosinase, and found that a high-frequency promoter variant, TYR c.-301C>T [rs4547091], modulates the penetrance of a prevalent, albinism-associated missense change, TYR c.1205G>A (p.Arg402Gln) [rs1126809]. We also found that homozygosity for a haplotype formed by three common, functionally-relevant variants, TYR c.[-301C;575C>A;1205G>A], is associated with a high probability of receiving an albinism diagnosis (OR>82). This genotype is also associated with reduced visual acuity and with increased central retinal thickness in UK Biobank participants. Finally, we report how the combined analysis of rare and common variants can increase diagnostic yield and can help inform genetic counselling in families with albinism.

Dynamic changes in the epigenomic landscape regulate human organogenesis and link to developmental disorders.

How the genome activates or silences transcriptional programmes governs organ formation. Little is known in human embryos undermining our ability to benchmark the fidelity of stem cell differentiation or cell programming, or interpret the pathogenicity of noncoding variation. Here, we study histone modifications across thirteen tissues during human organogenesis. We integrate the data with transcription to build an overview of how the human genome differentially regulates alternative organ fates including by repression. Promoters from nearly 20,000 genes partition into discrete states. Key developmental gene sets are actively repressed outside of the appropriate organ without obvious bivalency. Candidate enhancers, functional in zebrafish, allow imputation of tissue-specific and shared patterns of transcription factor binding. Overlaying more than 700 noncoding mutations from patients with developmental disorders allows correlation to unanticipated target genes. Taken together, the data provide a comprehensive genomic framework for investigating normal and abnormal human development.

An integrative transcriptomic atlas of organogenesis in human embryos.

Human organogenesis is when severe developmental abnormalities commonly originate. However, understanding this critical embryonic phase has relied upon inference from patient phenotypes and assumptions from in vitro stem cell models and non-human vertebrates. We report an integrated transcriptomic atlas of human organogenesis. By lineage-guided principal components analysis, we uncover novel relatedness of particular developmental genes across different organs and tissues and identified unique transcriptional codes which correctly predicted the cause of many congenital disorders. By inference, our model pinpoints co-enriched genes as new causes of developmental disorders such as cleft palate and congenital heart disease. The data revealed more than 6000 novel transcripts, over 90% of which fulfil criteria as long non-coding RNAs correlated with the protein-coding genome over megabase distances. Taken together, we have uncovered cryptic transcriptional programs used by the human embryo and established a new resource for the molecular understanding of human organogenesis and its associated disorders.

Phenotypic and functional analyses show stem cell-derived hepatocyte-like cells better mimic fetal rather than adult hepatocytes.

BACKGROUND & AIMS: Hepatocyte-like cells (HLCs), differentiated from pluripotent stem cells by the use of soluble factors, can model human liver function and toxicity. However, at present HLC maturity and whether any deficit represents a true fetal state or aberrant differentiation is unclear and compounded by comparison to potentially deteriorated adult hepatocytes. Therefore, we generated HLCs from multiple lineages, using two different protocols, for direct comparison with fresh fetal and adult hepatocytes. METHODS: Protocols were developed for robust differentiation. Multiple transcript, protein and functional analyses compared HLCs to fresh human fetal and adult hepatocytes. RESULTS: HLCs were comparable to those of other laboratories by multiple parameters. Transcriptional changes during differentiation mimicked human embryogenesis and showed more similarity to pericentral than periportal hepatocytes. Unbiased proteomics demonstrated greater proximity to liver than 30 other human organs or tissues. However, by comparison to fresh material, HLC maturity was proven by transcript, protein and function to be fetal-like and short of the adult phenotype. The expression of 81% phase 1 enzymes in HLCs was significantly upregulated and half were statistically not different from fetal hepatocytes. HLCs secreted albumin and metabolized testosterone (CYP3A) and dextrorphan (CYP2D6) like fetal hepatocytes. In seven bespoke tests, devised by principal components analysis to distinguish fetal from adult hepatocytes, HLCs from two different source laboratories consistently demonstrated fetal characteristics. CONCLUSIONS: HLCs from different sources are broadly comparable with unbiased proteomic evidence for faithful differentiation down the liver lineage. This current phenotype mimics human fetal rather than adult hepatocytes.

Genome-Wide Responses of Female Fruit Flies Subjected to Divergent Mating Regimes.

Elevated rates of mating and reproduction cause decreased female survival and lifetime reproductive success across a wide range of taxa from flies to humans. These costs are fundamentally important to the evolution of life histories. Here we investigate the potential mechanistic basis of this classic life history component. We conducted 4 independent replicated experiments in which female Drosophila melanogaster were subjected to 'high' and 'low' mating regimes, resulting in highly significant differences in lifespan. We sampled females for transcriptomic analysis at day 10 of life, before the visible onset of ageing, and used Tiling expression arrays to detect differential gene expression in two body parts (abdomen versus head+thorax). The divergent mating regimes were associated with significant differential expression in a network of genes showing evidence for interactions with ecdysone receptor. Preliminary experimental manipulation of two genes in that network with roles in post-transcriptional modification (CG11486, eyegone) tended to enhance sensitivity to mating costs. However, the subtle nature of those effects suggests substantial functional redundancy or parallelism in this gene network, which could buffer females against excessive responses. There was also evidence for differential expression in genes involved in germline maintenance, cell proliferation and in gustation / odorant reception. Interestingly, we detected differential expression in three specific genes (EcR, keap1, lbk1) and one class of genes (gustation / odorant receptors) with previously reported roles in determining lifespan. Our results suggest that high and low mating regimes that lead to divergence in lifespan are associated with changes in the expression of genes such as reproductive hormones, that influence resource allocation to the germ line, and that may modify post-translational gene expression. This predicts that the correct signalling of nutrient levels to the reproductive system is important for maintaining organismal integrity.

Proteome-wide analyses of human hepatocytes during differentiation and dedifferentiation.

UNLABELLED: Failure to predict hepatotoxic drugs in preclinical testing makes it imperative to develop better liver models with a stable phenotype in culture. Stem cell-derived models offer promise, with differentiated hepatocyte-like cells currently considered to be "fetal-like" in their maturity. However, this judgment is based on limited biomarkers or transcripts and lacks the required proteomic datasets that directly compare fetal and adult hepatocytes. Here, we quantitatively compare the proteomes of human fetal liver, adult hepatocytes, and the HepG2 cell line. In addition, we investigate the proteome changes in human fetal and adult hepatocytes when cultured in a new air-liquid interface format compared to conventional submerged extracellular matrix sandwich culture. From albumin and urea secretion, and luciferase-based cytochrome P450 activity, adult hepatocytes were viable in either culture model over 2 weeks. The function of fetal cells was better maintained in the air-liquid interface system. Strikingly, the proteome was qualitatively similar across all samples but hierarchical clustering showed that each sample type had a distinct quantitative profile. HepG2 cells more closely resembled fetal than adult hepatocytes. Furthermore, clustering showed that primary adult hepatocytes cultured at the air-liquid interface retained a proteome that more closely mimicked their fresh counterparts than conventional culture, which acquired myofibroblast features. Principal component analysis extended these findings and identified a simple set of proteins, including cytochrome P450 2A6, glutathione S transferase P, and alcohol dehydrogenases as specialized indicators of hepatocyte differentiation. CONCLUSION: Our quantitative datasets are the first that directly compare multiple human liver cells, define a model for enhanced maintenance of the hepatocyte proteome in culture, and provide a new protein "toolkit" for determining human hepatocyte maturity in cultured cells.

An excess of gene expression divergence on the X chromosome in Drosophila embryos: implications for the faster-X hypothesis.

The X chromosome is present as a single copy in the heterogametic sex, and this hemizygosity is expected to drive unusual patterns of evolution on the X relative to the autosomes. For example, the hemizgosity of the X may lead to a lower chromosomal effective population size compared to the autosomes, suggesting that the X might be more strongly affected by genetic drift. However, the X may also experience stronger positive selection than the autosomes, because recessive beneficial mutations will be more visible to selection on the X where they will spend less time being masked by the dominant, less beneficial allele--a proposal known as the faster-X hypothesis. Thus, empirical studies demonstrating increased genetic divergence on the X chromosome could be indicative of either adaptive or non-adaptive evolution. We measured gene expression in Drosophila species and in D. melanogaster inbred strains for both embryos and adults. In the embryos we found that expression divergence is on average more than 20% higher for genes on the X chromosome relative to the autosomes; but in contrast, in the inbred strains, gene expression variation is significantly lower on the X chromosome. Furthermore, expression divergence of genes on Muller's D element is significantly greater along the branch leading to the obscura sub-group, in which this element segregates as a neo-X chromosome. In the adults, divergence is greatest on the X chromosome for males, but not for females, yet in both sexes inbred strains harbour the lowest level of gene expression variation on the X chromosome. We consider different explanations for our results and conclude that they are most consistent within the framework of the faster-X hypothesis.

REDfly v3.0: toward a comprehensive database of transcriptional regulatory elements in Drosophila.

The REDfly database of Drosophila transcriptional cis-regulatory elements provides the broadest and most comprehensive available resource for experimentally validated cis-regulatory modules and transcription factor binding sites among the metazoa. The third major release of the database extends the utility of REDfly as a powerful tool for both computational and experimental studies of transcription regulation. REDfly v3.0 includes the introduction of new data classes to expand the types of regulatory elements annotated in the database along with a roughly 40% increase in the number of records. A completely redesigned interface improves access for casual and power users alike; among other features it now automatically provides graphical views of the genome, displays images of reporter gene expression and implements improved capabilities for database searching and results filtering. REDfly is freely accessible at http://redfly.ccr.buffalo.edu.

Quantifying the life-history response to increased male exposure in female Drosophila melanogaster.

Precise estimates of costs and benefits, the fitness economics, of mating are of key importance in understanding how selection shapes the coevolution of male and female mating traits. However, fitness is difficult to define and quantify. Here, we used a novel application of an established analytical technique to calculate individual- and population-based estimates of fitness-including those sensitive to the timing of reproduction-to measure the effects on females of increased exposure to males. Drosophila melanogaster females were exposed to high and low frequencies of contact with males, and life-history traits for each individual female were recorded. We then compared different fitness estimates to determine which of them best described the changes in life histories. We predicted that rate-sensitive estimates would be more accurate, as mating influences the rate of offspring production in this species. The results supported this prediction. Increased exposure to males led to significantly decreased fitness within declining but not stable or increasing populations. There was a net benefit of increased male exposure in expanding populations, despite a significant decrease in lifespan. The study shows how a more accurate description of fitness, and new insights can be achieved by considering individual life-history strategies within the context of population growth.

Gene expression divergence recapitulates the developmental hourglass model.

The observation that animal morphology tends to be conserved during the embryonic phylotypic period (a period of maximal similarity between the species within each animal phylum) led to the proposition that embryogenesis diverges more extensively early and late than in the middle, known as the hourglass model. This pattern of conservation is thought to reflect a major constraint on the evolution of animal body plans. Despite a wealth of morphological data confirming that there is often remarkable divergence in the early and late embryos of species from the same phylum, it is not yet known to what extent gene expression evolution, which has a central role in the elaboration of different animal forms, underpins the morphological hourglass pattern. Here we address this question using species-specific microarrays designed from six sequenced Drosophila species separated by up to 40 million years. We quantify divergence at different times during embryogenesis, and show that expression is maximally conserved during the arthropod phylotypic period. By fitting different evolutionary models to each gene, we show that at each time point more than 80% of genes fit best to models incorporating stabilizing selection, and that for genes whose evolutionarily optimal expression level is the same across all species, selective constraint is maximized during the phylotypic period. The genes that conform most to the hourglass pattern are involved in key developmental processes. These results indicate that natural selection acts to conserve patterns of gene expression during mid-embryogenesis, and provide a genome-wide insight into the molecular basis of the hourglass pattern of developmental evolution.

CEACAM3: an innate immune receptor directed against human-restricted bacterial pathogens.

Carcinoembryonic antigen-related cell adhesion molecule 3 (CEACAM3) is an immunoglobulin-related glycoprotein exclusively expressed on granulocytes. In contrast to other members of the CEACAM family, CEACAM3 does not support cell-cell adhesion, but rather mediates the opsonin-independent recognition and elimination of a restricted set of human-specific Gram-negative bacterial pathogens including Neisseria gonorrhoeae, Haemophilus influenzae, and Moraxella catarrhalis. Within the last 4 years, molecular determinants of CEACAM3 function and CEACAM3-initiated signaling pathways have been elucidated. Sequence comparison between CEACAM3 and other CEACAM family members points to a chimeric origin of this receptor with the bacteria-binding extracellular domain and the function-promoting intracellular domain derived from different genes. This review summarizes the current knowledge about the structure-function relationship of CEACAM3 and tries to combine these molecular aspects with a plausible scenario concerning the evolutionary origin of this phagocyte receptor in the light of host-pathogen adaptation.

Positive selection and gene conversion in SPP120, a fertilization-related gene, during the East African cichlid fish radiation.

The ability to infer historical natural selection from sequence data aides in finding genes that might be important in adaptation and the formation of new species. As the fastest evolving and largest known vertebrate radiation, the cichlid fish of the African Great Lakes exhibit a wide range of recent morphological diversification. We used DNA databases, mostly of expressed sequence tags, to find candidate orthologous coding sequences from 2 tribes of cichlids and, using an automated procedure, scanned these sequence pairs for high dN/dS, the signal of positive selection and protein adaptation. The results included vertebrate genes commonly found to be under selection (e.g., major histocompatibility complex [MHC] loci) as well as genes known to be important specifically in the cichlid radiation (e.g., long-wave-sensitive opsins). Further investigation focused on a gene encoding a fertilization-related protein, SPP120, which was previously known only from cichlids. Using maximum likelihood analysis on novel SPP120 cDNA sequences from a range of African cichlids, we demonstrate the influence of positive selection in a specific subregion of the protein. We also show that SPP120 is a tandemly arranged, multicopy gene evolving with occasional interlocus gene conversion. A search of the Medaka genome database also revealed a tandem arrangement of multiple SPP120 copies and evolutionary rate differences between Medaka gene subregions mirroring those found for cichlids. Combined, these results suggest that SPP120 has been under repeated diversifying selection for over 100 Myr.

Positive and negative selection on mammalian Y chromosomes.

Y chromosomes are genetically degenerate in most organisms studied. The loss of genes from Y chromosomes is thought to be due to the inefficiency of purifying selection in nonrecombining regions, which leads to the accumulation of deleterious mutations via the processes of hitchhiking, background selection, and Muller's ratchet. As the severity of these processes depends on the number of functional genes linked together on the nonrecombining Y, it is not clear whether these processes are still at work on the old, gene-poor mammalian Y chromosomes. If purifying selection is indeed less efficient in the Y-linked, compared to the X-linked genes, deleterious nonsynonymous substitutions are expected to accumulate faster on the Y chromosome. However, positive selection on Y-linked genes could also increase the rate of amino acid-changing substitutions. Thus, the previous reports of an elevated nonsynonymous substitution rate in Y-linked genes are still open to interpretation. Here, we report evidence for positive selection in two out of three studied mammalian Y-linked genes, suggesting that adaptive Darwinian evolution may be common on mammalian Y chromosomes. Taking positive selection into account, we demonstrate that purifying selection is less efficient in mammalian Y-linked genes compared to their X-linked homologues, suggesting that these genes continue to degenerate.

High mutation rates in human and ape pseudoautosomal genes.

It has been suggested that recombination may be mutagenic, which, if true, would inflate intraspecies diversity and interspecies silent divergence in regions of high recombination. Here, we test this hypothesis comparing human/orangutan genome-wide non-coding divergence (K) to that in the pseudoautosomal genes which were reported to recombine much more frequently than the rest of the genome. We demonstrate that, compared to the average human/orangutan non-coding divergence (K=3%), the substitution rate is significantly elevated in the introns of SHOX (K=5.7%), PPP2R3L (K=8.7%) and ASMT (K=6.5%) genes located in the human and orangutan Xp/Yp pseudoautosomal region (p-PAR), where recombination is over 20-fold higher than the genomic average. On the other hand, human/orangutan non-coding divergence at the Xp/Yp pseudoautosomal boundary (K=3.5%) and in the SYBL1 gene (K=2.7%), located in the human Xq/Yq pseudoautosomal region (q-PAR), where recombination is known to be less frequent than in p-PAR, was not significantly higher than the genome average. The data are consistent with the hypothesis that recombination may be mutagenic.