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.

A wealth of data and comprehensive reviews exist on pancreas development in mammals, primarily mice, and other vertebrates. By contrast, human pancreatic development has been less comprehensively reviewed. Here, we draw together those studies conducted directly in human embryonic and fetal tissue to provide an overview of what is known about human pancreatic development. We discuss the relevance of this work to manufacturing insulin-secreting ß-cells from pluripotent stem cells and to different aspects of diabetes, especially permanent neonatal diabetes, and its underlying causes

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.

Failure to predict hepatotoxic drugs in pre-clinical 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 judgement 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. Typical biomarkers showed that adult hepatocytes functioned equally well in sandwich or air-liquid interface culture. Fetal cells, however, were viable over longer culture periods and their function was enhanced over time 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. The clustering also shows that primary cells 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 hepatocytes in culture and provide a new protein ‘toolkit’ for determining human hepatocyte maturity in cultured cells.

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.

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.

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.

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.

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.

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

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.

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.

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.

Sequence analysis is widely used to infer function from one protein sequence to another. One of the remaining hurdles in pushing further the limits of efficient usage is the modular and mobile architecture of proteins. Although many resources provide precompiled signatures, they are not yet used in conjunction with pairwise comparisons to investigate the modular architecture and phylogeny. We present a program, doMosaic, which combines existing domain definitions with an adapted sequence alignment algorithm to analyse the possible phylogeny of domain architecture in proteins. It is based on the Smith-Waterman scheme, can be combined with trees derived from the domains and provides a user-friendly graphical interface. The method enables fast and efficient analysis of domain duplication events within one or in-between two sequences. Therefore, it enables refined functional annotation. (pp. 15-20)