Identifying the genetic variants that increase the risk of type 2 diabetes (T2D) in humans has been a formidable challenge. Adopting a genome-wide association strategy, we genotyped 1161 Finnish T2D cases and 1174 Finnish normal glucose-tolerant (NGT) controls with >315,000 single-nucleotide polymorphisms (SNPs) and imputed genotypes for an additional >2 million autosomal SNPs. We carried out association analysis with these SNPs to identify genetic variants that predispose to T2D, compared our T2D association results with the results of two similar studies, and genotyped 80 SNPs in an additional 1215 Finnish T2D cases and 1258 Finnish NGT controls. We identify T2D-associated variants in an intergenic region of chromosome 11p12, contribute to the identification of T2D-associated variants near the genes IGF2BP2 and CDKAL1 and the region of CDKN2A and CDKN2B, and confirm that variants near TCF7L2, SLC30A8, HHEX, FTO, PPARG, and KCNJ11 are associated with T2D risk. This brings the number of T2D loci now confidently identified to at least 10.

BACKGROUND

We tested the hypothesis that the level of circulating tumor cells can predict survival in metastatic breast cancer.

METHODS

In a prospective, multicenter study, we tested 177 patients with measurable metastatic breast cancer for levels of circulating tumor cells both before the patients were to start a new line of treatment and at the first follow-up visit. The progression of the disease or the response to treatment was determined with the use of standard imaging studies at the participating centers.

RESULTS

Outcomes were assessed according to levels of circulating tumor cells at baseline, before the patients started a new treatment for metastatic disease. Patients in a training set with levels of circulating tumor cells equal to or higher than 5 per 7.5 ml of whole blood, as compared with the group with fewer than 5 circulating tumor cells per 7.5 ml, had a shorter median progression-free survival (2.7 months vs. 7.0 months, P<0.001) and shorter overall survival (10.1 months vs. >18 months, P<0.001). At the first follow-up visit after the initiation of therapy, this difference between the groups persisted (progression-free survival, 2.1 months vs. 7.0 months; P<0.001; overall survival, 8.2 months vs. >18 months; P<0.001), and the reduced proportion of patients (from 49 percent to 30 percent) in the group with an unfavorable prognosis suggested that there was a benefit from therapy. The multivariate Cox proportional-hazards regression showed that, of all the variables in the statistical model, the levels of circulating tumor cells at baseline and at the first follow-up visit were the most significant predictors of progression-free and overall survival.

CONCLUSIONS

The number of circulating tumor cells before treatment is an independent predictor of progression-free survival and overall survival in patients with metastatic breast cancer.

Background: Coronavirus disease 2019 (Covid-19) is associated with diffuse lung damage. Glucocorticoids may modulate inflammation-mediated lung injury and thereby reduce progression to respiratory failure and death.

Methods: In this controlled, open-label trial comparing a range of possible treatments in patients who were hospitalized with Covid-19, we randomly assigned patients to receive oral or intravenous dexamethasone (at a dose of 6 mg once daily) for up to 10 days or to receive usual care alone. The primary outcome was 28-day mortality. Here, we report the preliminary results of this comparison.

Results: A total of 2104 patients were assigned to receive dexamethasone and 4321 to receive usual care. Overall, 482 patients (22.9%) in the dexamethasone group and 1110 patients (25.7%) in the usual care group died within 28 days after randomization (age-adjusted rate ratio, 0.83; 95% confidence interval [CI], 0.75 to 0.93; P<0.001). The proportional and absolute between-group differences in mortality varied considerably according to the level of respiratory support that the patients were receiving at the time of randomization. In the dexamethasone group, the incidence of death was lower than that in the usual care group among patients receiving invasive mechanical ventilation (29.3% vs. 41.4%; rate ratio, 0.64; 95% CI, 0.51 to 0.81) and among those receiving oxygen without invasive mechanical ventilation (23.3% vs. 26.2%; rate ratio, 0.82; 95% CI, 0.72 to 0.94) but not among those who were receiving no respiratory support at randomization (17.8% vs. 14.0%; rate ratio, 1.19; 95% CI, 0.91 to 1.55).

Conclusions: In patients hospitalized with Covid-19, the use of dexamethasone resulted in lower 28-day mortality among those who were receiving either invasive mechanical ventilation or oxygen alone at randomization but not among those receiving no respiratory support. (Funded by the Medical Research Council and National Institute for Health Research and others; RECOVERY ClinicalTrials.gov number, NCT04381936; ISRCTN number, 50189673.).

There are many sources of systematic variation in cDNA microarray experiments which affect the measured gene expression levels (e.g. differences in labeling efficiency between the two fluorescent dyes). The term normalization refers to the process of removing such variation. A constant adjustment is often used to force the distribution of the intensity log ratios to have a median of zero for each slide. However, such global normalization approaches are not adequate in situations where dye biases can depend on spot overall intensity and/or spatial location within the array. This article proposes normalization methods that are based on robust local regression and account for intensity and spatial dependence in dye biases for different types of cDNA microarray experiments. The selection of appropriate controls for normalization is discussed and a novel set of controls (microarray sample pool, MSP) is introduced to aid in intensity-dependent normalization. Lastly, to allow for comparisons of expression levels across slides, a robust method based on maximum likelihood estimation is proposed to adjust for scale differences among slides.

We present GOseq, an application for performing Gene Ontology (GO) analysis on RNA-seq data. GO analysis is widely used to reduce complexity and highlight biological processes in genome-wide expression studies, but standard methods give biased results on RNA-seq data due to over-detection of differential expression for long and highly expressed transcripts. Application of GOseq to a prostate cancer data set shows that GOseq dramatically changes the results, highlighting categories more consistent with the known biology.

Cells continually experience stress and damage from exogenous and endogenous sources, and their responses range from complete recovery to cell death. Proliferating cells can initiate an additional response by adopting a state of permanent cell-cycle arrest that is termed cellular senescence. Understanding the causes and consequences of cellular senescence has provided novel insights into how cells react to stress, especially genotoxic stress, and how this cellular response can affect complex organismal processes such as the development of cancer and ageing.

It is hypothesized that collective efficacy, defined as social cohesion among neighbors combined with their willingness to intervene on behalf of the common good, is linked to reduced violence. This hypothesis was tested on a 1995 survey of 8782 residents of 343 neighborhoods in Chicago, Illinois. Multilevel analyses showed that a measure of collective efficacy yields a high between-neighborhood reliability and is negatively associated with variations in violence, when individual-level characteristics, measurement error, and prior violence are controlled. Associations of concentrated disadvantage and residential instability with violence are largely mediated by collective efficacy.

The terminus of a DNA helix has been called its Achilles' heel. Thus to prevent possible incomplete replication and instability of the termini of linear DNA, eukaryotic chromosomes end in characteristic repetitive DNA sequences within specialized structures called telomeres. In immortal cells, loss of telomeric DNA due to degradation or incomplete replication is apparently balanced by telomere elongation, which may involve de novo synthesis of additional repeats by novel DNA polymerase called telomerase. Such a polymerase has been recently detected in HeLa cells. It has been proposed that the finite doubling capacity of normal mammalian cells is due to a loss of telomeric DNA and eventual deletion of essential sequences. In yeast, the est1 mutation causes gradual loss of telomeric DNA and eventual cell death mimicking senescence in higher eukaryotic cells. Here, we show that the amount and length of telomeric DNA in human fibroblasts does in fact decrease as a function of serial passage during ageing in vitro and possibly in vivo. It is not known whether this loss of DNA has a causal role in senescence.

BACKGROUND

Usually the researchers performing meta-analysis of continuous outcomes from clinical trials need their mean value and the variance (or standard deviation) in order to pool data. However, sometimes the published reports of clinical trials only report the median, range and the size of the trial.

METHODS

In this article we use simple and elementary inequalities and approximations in order to estimate the mean and the variance for such trials. Our estimation is distribution-free, i.e., it makes no assumption on the distribution of the underlying data.

RESULTS

We found two simple formulas that estimate the mean using the values of the median (m), low and high end of the range (a and b, respectively), and n (the sample size). Using simulations, we show that median can be used to estimate mean when the sample size is larger than 25. For smaller samples our new formula, devised in this paper, should be used. We also estimated the variance of an unknown sample using the median, low and high end of the range, and the sample size. Our estimate is performing as the best estimate in our simulations for very small samples (n < or = 15). For moderately sized samples (15 < n < or = 70), our simulations show that the formula range/4 is the best estimator for the standard deviation (variance). For large samples (n>> 70), the formula range/6 gives the best estimator for the standard deviation (variance). We also include an illustrative example of the potential value of our method using reports from the Cochrane review on the role of erythropoietin in anemia due to malignancy.

CONCLUSIONS

Using these formulas, we hope to help meta-analysts use clinical trials in their analysis even when not all of the information is available and/or reported.

There is persuasive clinical and experimental evidence that macrophages promote cancer initiation and malignant progression. During tumor initiation, they create an inflammatory environment that is mutagenic and promotes growth. As tumors progress to malignancy, macrophages stimulate angiogenesis, enhance tumor cell migration and invasion, and suppress antitumor immunity. At metastatic sites, macrophages prepare the target tissue for arrival of tumor cells, and then a different subpopulation of macrophages promotes tumor cell extravasation, survival, and subsequent growth. Specialized subpopulations of macrophages may represent important new therapeutic targets.

OBJECTIVE

To establish a unified working diagnostic tool for the metabolic syndrome (MetS) that is convenient to use in clinical practice and that can be used world-wide so that data from different countries can be compared. An additional aim was to highlight areas where more research into the MetS is needed.

METHODS

The International Diabetes Federation (IDF) convened a workshop held 12-14 May 2004 in London, UK. The 21 participants included experts in the fields of diabetes, public health, epidemiology, lipidology, genetics, metabolism, nutrition and cardiology. There were participants from each of the five continents as well as from the World Health Organization (WHO) and the National Cholesterol Education Program-Third Adult Treatment Panel (ATP III). The workshop was sponsored by an educational grant from AstraZeneca Pharmaceuticals.

METHODS

The consensus statement emerged following detailed discussions at the IDF workshop. After the workshop, a writing group produced a consensus statement which was reviewed and approved by all participants.

CONCLUSIONS

The IDF has produced a new set of criteria for use both epidemiologically and in clinical practice world-wide with the aim of identifying people with the MetS to clarify the nature of the syndrome and to focus therapeutic strategies to reduce the long-term risk of cardiovascular disease. Guidance is included on how to compensate for differences in waist circumference and in regional adipose tissue distribution between different populations. The IDF has also produced recommendations for additional criteria that should be included when studying the MetS for research purposes. Finally, the IDF has identified areas where more studies are currently needed; these include research into the aetiology of the syndrome.

New theoretical and practical concepts are presented for considerably enhancing the performance of magnetic resonance imaging (MRI) by means of arrays of multiple receiver coils. Sensitivity encoding (SENSE) is based on the fact that receiver sensitivity generally has an encoding effect complementary to Fourier preparation by linear field gradients. Thus, by using multiple receiver coils in parallel scan time in Fourier imaging can be considerably reduced. The problem of image reconstruction from sensitivity encoded data is formulated in a general fashion and solved for arbitrary coil configurations and k-space sampling patterns. Special attention is given to the currently most practical case, namely, sampling a common Cartesian grid with reduced density. For this case the feasibility of the proposed methods was verified both in vitro and in vivo. Scan time was reduced to one-half using a two-coil array in brain imaging. With an array of five coils double-oblique heart images were obtained in one-third of conventional scan time. Magn Reson Med 42:952-962, 1999.

New information regarding neuronal circuits that control food intake and their hormonal regulation has extended our understanding of energy homeostasis, the process whereby energy intake is matched to energy expenditure over time. The profound obesity that results in rodents (and in the rare human case as well) from mutation of key signalling molecules involved in this regulatory system highlights its importance to human health. Although each new signalling pathway discovered in the hypothalamus is a potential target for drug development in the treatment of obesity, the growing number of such signalling molecules indicates that food intake is controlled by a highly complex process. To better understand how energy homeostasis can be achieved, we describe a model that delineates the roles of individual hormonal and neuropeptide signalling pathways in the control of food intake and the means by which obesity can arise from inherited or acquired defects in their function.

BACKGROUND

Trastuzumab, a recombinant monoclonal antibody against HER2, has clinical activity in advanced breast cancer that overexpresses HER2. We investigated its efficacy and safety after excision of early-stage breast cancer and completion of chemotherapy.

METHODS

This international, multicenter, randomized trial compared one or two years of trastuzumab given every three weeks with observation in patients with HER2-positive and either node-negative or node-positive breast cancer who had completed locoregional therapy and at least four cycles of neoadjuvant or adjuvant chemotherapy.

RESULTS

Data were available for 1694 women randomly assigned to two years of treatment with trastuzumab, 1694 women assigned to one year of trastuzumab, and 1693 women assigned to observation. We report here the results only of treatment with trastuzumab for one year or observation. At the first planned interim analysis (median follow-up of one year), 347 events (recurrence of breast cancer, contralateral breast cancer, second nonbreast malignant disease, or death) were observed: 127 events in the trastuzumab group and 220 in the observation group. The unadjusted hazard ratio for an event in the trastuzumab group, as compared with the observation group, was 0.54 (95 percent confidence interval, 0.43 to 0.67; P<0.0001 by the log-rank test, crossing the interim analysis boundary), representing an absolute benefit in terms of disease-free survival at two years of 8.4 percentage points. Overall survival in the two groups was not significantly different (29 deaths with trastuzumab vs. 37 with observation). Severe cardiotoxicity developed in 0.5 percent of the women who were treated with trastuzumab.

CONCLUSIONS

One year of treatment with trastuzumab after adjuvant chemotherapy significantly improves disease-free survival among women with HER2-positive breast cancer. (ClinicalTrials.gov number, NCT00045032.)

NCBI's Conserved Domain Database (CDD) is a resource for the annotation of protein sequences with the location of conserved domain footprints, and functional sites inferred from these footprints. CDD includes manually curated domain models that make use of protein 3D structure to refine domain models and provide insights into sequence/structure/function relationships. Manually curated models are organized hierarchically if they describe domain families that are clearly related by common descent. As CDD also imports domain family models from a variety of external sources, it is a partially redundant collection. To simplify protein annotation, redundant models and models describing homologous families are clustered into superfamilies. By default, domain footprints are annotated with the corresponding superfamily designation, on top of which specific annotation may indicate high-confidence assignment of family membership. Pre-computed domain annotation is available for proteins in the Entrez/Protein dataset, and a novel interface, Batch CD-Search, allows the computation and download of annotation for large sets of protein queries. CDD can be accessed via http://www.ncbi.nlm.nih.gov/Structure/cdd/cdd.shtml.

Proteins often differ in amino-acid sequence across species. This difference has evolved by the accumulation of neutral mutations by random drift, the fixation of adaptive mutations by selection, or a mixture of the two. Here we propose a simple statistical test of the neutral protein evolution hypothesis based on a comparison of the number of amino-acid replacement substitutions to synonymous substitutions in the coding region of a locus. If the observed substitutions are neutral, the ratio of replacement to synonymous fixed differences between species should be the same as the ratio of replacement to synonymous polymorphisms within species. DNA sequence data on the Adh locus (encoding alcohol dehydrogenase, EC 1.1.1.1) in three species in the Drosophila melanogaster species subgroup do not fit this expectation; instead, there are more fixed replacement differences between species than expected. We suggest that these excess replacement substitutions result from adaptive fixation of selectively advantageous mutations.

The innate immune system senses viral infection by recognizing a variety of viral components (including double-stranded (ds)RNA) and triggers antiviral responses. The cytoplasmic helicase proteins RIG-I (retinoic-acid-inducible protein I, also known as Ddx58) and MDA5 (melanoma-differentiation-associated gene 5, also known as Ifih1 or Helicard) have been implicated in viral dsRNA recognition. In vitro studies suggest that both RIG-I and MDA5 detect RNA viruses and polyinosine-polycytidylic acid (poly(I:C)), a synthetic dsRNA analogue. Although a critical role for RIG-I in the recognition of several RNA viruses has been clarified, the functional role of MDA5 and the relationship between these dsRNA detectors in vivo are yet to be determined. Here we use mice deficient in MDA5 (MDA5-/-) to show that MDA5 and RIG-I recognize different types of dsRNAs: MDA5 recognizes poly(I:C), and RIG-I detects in vitro transcribed dsRNAs. RNA viruses are also differentially recognized by RIG-I and MDA5. We find that RIG-I is essential for the production of interferons in response to RNA viruses including paramyxoviruses, influenza virus and Japanese encephalitis virus, whereas MDA5 is critical for picornavirus detection. Furthermore, RIG-I-/- and MDA5-/- mice are highly susceptible to infection with these respective RNA viruses compared to control mice. Together, our data show that RIG-I and MDA5 distinguish different RNA viruses and are critical for host antiviral responses.

OBJECTIVE

To improve clinical recognition and provide research diagnostic criteria for three clinical syndromes associated with frontotemporal lobar degeneration.

METHODS

Consensus criteria for the three prototypic syndromes-frontotemporal dementia, progressive nonfluent aphasia, and semantic dementia-were developed by members of an international workshop on frontotemporal lobar degeneration. These criteria build on earlier published clinical diagnostic guidelines for frontotemporal dementia produced by some of the workshop members.

RESULTS

The consensus criteria specify core and supportive features for each of the three prototypic clinical syndromes and provide broad inclusion and exclusion criteria for the generic entity of frontotemporal lobar degeneration. The criteria are presented in lists, and operational definitions for features are provided in the text.

CONCLUSIONS

The criteria ought to provide the foundation for research work into the neuropsychology, neuropathology, genetics, molecular biology, and epidemiology of these important clinical disorders that account for a substantial proportion of cases of primary degenerative dementia occurring before the age of 65 years.

In observational studies, investigators have no control over the treatment assignment. The treated and non-treated (that is, control) groups may have large differences on their observed covariates, and these differences can lead to biased estimates of treatment effects. Even traditional covariance analysis adjustments may be inadequate to eliminate this bias. The propensity score, defined as the conditional probability of being treated given the covariates, can be used to balance the covariates in the two groups, and therefore reduce this bias. In order to estimate the propensity score, one must model the distribution of the treatment indicator variable given the observed covariates. Once estimated the propensity score can be used to reduce bias through matching, stratification (subclassification), regression adjustment, or some combination of all three. In this tutorial we discuss the uses of propensity score methods for bias reduction, give references to the literature and illustrate the uses through applied examples.

We report the draft genome of the black cottonwood tree, Populus trichocarpa. Integration of shotgun sequence assembly with genetic mapping enabled chromosome-scale reconstruction of the genome. More than 45,000 putative protein-coding genes were identified. Analysis of the assembled genome revealed a whole-genome duplication event; about 8000 pairs of duplicated genes from that event survived in the Populus genome. A second, older duplication event is indistinguishably coincident with the divergence of the Populus and Arabidopsis lineages. Nucleotide substitution, tandem gene duplication, and gross chromosomal rearrangement appear to proceed substantially more slowly in Populus than in Arabidopsis. Populus has more protein-coding genes than Arabidopsis, ranging on average from 1.4 to 1.6 putative Populus homologs for each Arabidopsis gene. However, the relative frequency of protein domains in the two genomes is similar. Overrepresented exceptions in Populus include genes associated with lignocellulosic wall biosynthesis, meristem development, disease resistance, and metabolite transport.

The central dogma of gene expression is that DNA is transcribed into messenger RNAs, which in turn serve as the template for protein synthesis. The discovery of extensive transcription of large RNA transcripts that do not code for proteins, termed long noncoding RNAs (lncRNAs), provides an important new perspective on the centrality of RNA in gene regulation. Here, we discuss genome-scale strategies to discover and characterize lncRNAs. An emerging theme from multiple model systems is that lncRNAs form extensive networks of ribonucleoprotein (RNP) complexes with numerous chromatin regulators and then target these enzymatic activities to appropriate locations in the genome. Consistent with this notion, lncRNAs can function as modular scaffolds to specify higher-order organization in RNP complexes and in chromatin states. The importance of these modes of regulation is underscored by the newly recognized roles of long RNAs for proper gene control across all kingdoms of life.

Rice, one of the world's most important food plants, has important syntenic relationships with the other cereal species and is a model plant for the grasses. Here we present a map-based, finished quality sequence that covers 95% of the 389 Mb genome, including virtually all of the euchromatin and two complete centromeres. A total of 37,544 non-transposable-element-related protein-coding genes were identified, of which 71% had a putative homologue in Arabidopsis. In a reciprocal analysis, 90% of the Arabidopsis proteins had a putative homologue in the predicted rice proteome. Twenty-nine per cent of the 37,544 predicted genes appear in clustered gene families. The number and classes of transposable elements found in the rice genome are consistent with the expansion of syntenic regions in the maize and sorghum genomes. We find evidence for widespread and recurrent gene transfer from the organelles to the nuclear chromosomes. The map-based sequence has proven useful for the identification of genes underlying agronomic traits. The additional single-nucleotide polymorphisms and simple sequence repeats identified in our study should accelerate improvements in rice production.

BACKGROUND

Artemisinin-based combination therapies are the recommended first-line treatments of falciparum malaria in all countries with endemic disease. There are recent concerns that the efficacy of such therapies has declined on the Thai-Cambodian border, historically a site of emerging antimalarial-drug resistance.

METHODS

In two open-label, randomized trials, we compared the efficacies of two treatments for uncomplicated falciparum malaria in Pailin, western Cambodia, and Wang Pha, northwestern Thailand: oral artesunate given at a dose of 2 mg per kilogram of body weight per day, for 7 days, and artesunate given at a dose of 4 mg per kilogram per day, for 3 days, followed by mefloquine at two doses totaling 25 mg per kilogram. We assessed in vitro and in vivo Plasmodium falciparum susceptibility, artesunate pharmacokinetics, and molecular markers of resistance.

RESULTS

We studied 40 patients in each of the two locations. The overall median parasite clearance times were 84 hours (interquartile range, 60 to 96) in Pailin and 48 hours (interquartile range, 36 to 66) in Wang Pha (P<0.001). Recrudescence confirmed by means of polymerase-chain-reaction assay occurred in 6 of 20 patients (30%) receiving artesunate monotherapy and 1 of 20 (5%) receiving artesunate-mefloquine therapy in Pailin, as compared with 2 of 20 (10%) and 1 of 20 (5%), respectively, in Wang Pha (P=0.31). These markedly different parasitologic responses were not explained by differences in age, artesunate or dihydroartemisinin pharmacokinetics, results of isotopic in vitro sensitivity tests, or putative molecular correlates of P. falciparum drug resistance (mutations or amplifications of the gene encoding a multidrug resistance protein [PfMDR1] or mutations in the gene encoding sarco-endoplasmic reticulum calcium ATPase6 [PfSERCA]). Adverse events were mild and did not differ significantly between the two treatment groups.

CONCLUSIONS

P. falciparum has reduced in vivo susceptibility to artesunate in western Cambodia as compared with northwestern Thailand. Resistance is characterized by slow parasite clearance in vivo without corresponding reductions on conventional in vitro susceptibility testing. Containment measures are urgently needed. (ClinicalTrials.gov number, NCT00493363, and Current Controlled Trials number, ISRCTN64835265.)

The origin of organismal complexity is generally thought to be tightly coupled to the evolution of new gene functions arising subsequent to gene duplication. Under the classical model for the evolution of duplicate genes, one member of the duplicated pair usually degenerates within a few million years by accumulating deleterious mutations, while the other duplicate retains the original function. This model further predicts that on rare occasions, one duplicate may acquire a new adaptive function, resulting in the preservation of both members of the pair, one with the new function and the other retaining the old. However, empirical data suggest that a much greater proportion of gene duplicates is preserved than predicted by the classical model. Here we present a new conceptual framework for understanding the evolution of duplicate genes that may help explain this conundrum. Focusing on the regulatory complexity of eukaryotic genes, we show how complementary degenerative mutations in different regulatory elements of duplicated genes can facilitate the preservation of both duplicates, thereby increasing long-term opportunities for the evolution of new gene functions. The duplication-degeneration-complementation (DDC) model predicts that (1) degenerative mutations in regulatory elements can increase rather than reduce the probability of duplicate gene preservation and (2) the usual mechanism of duplicate gene preservation is the partitioning of ancestral functions rather than the evolution of new functions. We present several examples (including analysis of a new engrailed gene in zebrafish) that appear to be consistent with the DDC model, and we suggest several analytical and experimental approaches for determining whether the complementary loss of gene subfunctions or the acquisition of novel functions are likely to be the primary mechanisms for the preservation of gene duplicates. For a newly duplicated paralog, survival depends on the outcome of the race between entropic decay and chance acquisition of an advantageous regulatory mutation. Sidow 1996(p. 717) On one hand, it may fix an advantageous allele giving it a slightly different, and selectable, function from its original copy. This initial fixation provides substantial protection against future fixation of null mutations, allowing additional mutations to accumulate that refine functional differentiation. Alternatively, a duplicate locus can instead first fix a null allele, becoming a pseudogene. Walsh 1995 (p. 426) Duplicated genes persist only if mutations create new and essential protein functions, an event that is predicted to occur rarely. Nadeau and Sankoff 1997 (p. 1259) Thus overall, with complex metazoans, the major mechanism for retention of ancient gene duplicates would appear to have been the acquisition of novel expression sites for developmental genes, with its accompanying opportunity for new gene roles underlying the progressive extension of development itself. Cooke et al. 1997 (p. 362)