Institute of Cytology and Genetics

An annotated reference sequence representing the hexaploid bread wheat genome in 21 pseudomolecules has been analyzed to identify the distribution and genomic context of coding and noncoding elements across the A, B, and D subgenomes. With an estimated coverage of 94% of the genome and containing 107,891 high-confidence gene models, this assembly enabled the discovery of tissue- and developmental stage-related coexpression networks by providing a transcriptome atlas representing major stages of wheat development. Dynamics of complex gene families involved in environmental adaptation and end-use quality were revealed at subgenome resolution and contextualized to known agronomic single-gene or quantitative trait loci. This community resource establishes the foundation for accelerating wheat research and application through improved understanding of wheat biology and genomics-assisted breeding.

Anopheles gambiae is the principal vector of malaria, a disease that afflicts more than 500 million people and causes more than 1 million deaths each year. Tenfold shotgun sequence coverage was obtained from the PEST strain of A. gambiae and assembled into scaffolds that span 278 million base pairs. A total of 91% of the genome was organized in 303 scaffolds; the largest scaffold was 23.1 million base pairs. There was substantial genetic variation within this strain, and the apparent existence of two haplotypes of approximately equal frequency ("dual haplotypes") in a substantial fraction of the genome likely reflects the outbred nature of the PEST strain. The sequence produced a conservative inference of more than 400,000 single-nucleotide polymorphisms that showed a markedly bimodal density distribution. Analysis of the genome sequence revealed strong evidence for about 14,000 protein-encoding transcripts. Prominent expansions in specific families of proteins likely involved in cell adhesion and immunity were noted. An expressed sequence tag analysis of genes regulated by blood feeding provided insights into the physiological adaptations of a hematophagous insect.

Deep sequencing of the gut microbiomes of 1135 participants from a Dutch population-based cohort shows relations between the microbiome and 126 exogenous and intrinsic host factors, including 31 intrinsic factors, 12 diseases, 19 drug groups, 4 smoking categories, and 60 dietary factors. These factors collectively explain 18.7% of the variation seen in the interindividual distance of microbial composition. We could associate 110 factors to 125 species and observed that fecal chromogranin A (CgA), a protein secreted by enteroendocrine cells, was exclusively associated with 61 microbial species whose abundance collectively accounted for 53% of microbial composition. Low CgA concentrations were seen in individuals with a more diverse microbiome. These results are an important step toward a better understanding of environment-diet-microbe-host interactions.

We report here the genome sequence of an ancient human. Obtained from ∼4,000-year-old permafrost-preserved hair, the genome represents a male individual from the first known culture to settle in Greenland. Sequenced to an average depth of 20×, we recover 79% of the diploid genome, an amount close to the practical limit of current sequencing technologies. We identify 353,151 high-confidence single-nucleotide polymorphisms (SNPs), of which 6.8% have not been reported previously. We estimate raw read contamination to be no higher than 0.8%. We use functional SNP assessment to assign possible phenotypic characteristics of the individual that belonged to a culture whose location has yielded only trace human remains. We compare the high-confidence SNPs to those of contemporary populations to find the populations most closely related to the individual. This provides evidence for a migration from Siberia into the New World some 5,500 years ago, independent of that giving rise to the modern Native Americans and Inuit. For the first time, the sequence of a near-complete nuclear genome has been obtained from the tissue of an ancient human. It comes from permafrost-preserved hair, about 4,000 years old, of a male palaeo-Eskimo of the Saqqaq culture, the earliest known settlers in Greenland. Functional single-nucleotide polymorphism (SNP) assessment was used to assign possible phenotypic characteristics. The analysis provides evidence for a migration from Siberia into the New World some 5,500 years ago, independent of the migration that gave rise to the modern Native Americans and Inuit. Elsewhere in the issue we profile the paper's last author Eske Willerslev, who headed the project and found the lock of hair in a Copenhagen museum basement — after a fruitless search among the archaeological sites of Peary Land. The first genome sequence of an ancient human is reported. It comes from an approximately 4,000-year-old permafrost-preserved hair from a male from the first known culture to settle in Greenland. Functional single-nucleotide polymorphism (SNP) assessment is used to assign possible phenotypic characteristics and high-confidence SNPs are compared to those of contemporary populations to find those most closely related to the individual.

The origins of neural systems remain unresolved. In contrast to other basal metazoans, ctenophores (comb jellies) have both complex nervous and mesoderm-derived muscular systems. These holoplanktonic predators also have sophisticated ciliated locomotion, behaviour and distinct development. Here we present the draft genome of Pleurobrachia bachei, Pacific sea gooseberry, together with ten other ctenophore transcriptomes, and show that they are remarkably distinct from other animal genomes in their content of neurogenic, immune and developmental genes. Our integrative analyses place Ctenophora as the earliest lineage within Metazoa. This hypothesis is supported by comparative analysis of multiple gene families, including the apparent absence of HOX genes, canonical microRNA machinery, and reduced immune complement in ctenophores. Although two distinct nervous systems are well recognized in ctenophores, many bilaterian neuron-specific genes and genes of ‘classical’ neurotransmitter pathways either are absent or, if present, are not expressed in neurons. Our metabolomic and physiological data are consistent with the hypothesis that ctenophore neural systems, and possibly muscle specification, evolved independently from those in other animals. The draft genome of the ctenophore Pleurobrachia bachei (Pacific sea gooseberry) is presented, together with ten other ctenophore transcriptomes — these genomes have a very different neurogenic, immune and developmental gene content when compared with other animal genomes, and it is proposed that ctenophore neural systems, and possibly muscle specification, evolved independently from those in other animals. Ctenophores (comb jellies) are enigmatic animals that combine two distinct nervous systems with an elementary brain-like centre and possess mesoderm-derived muscles appropriate to their predatory lifestyle. Leonid Moroz et al. present the draft genome of the ctenophore Pleurobrachia bachei (Pacific sea gooseberry), together with ten other ctenophore transcriptomes. These genomes have a neurogenic, immune and developmental gene content that differs markedly from other animal genomes: HOX genes and canonical microRNA machinery are absent, and the immune gene complement is reduced. Many bilaterian neuron-specific genes and genes of 'classical' neurotransmitter pathways either are absent or not expressed in neurons. The authors propose that ctenophore neural systems, and possibly muscle specification, evolved independently from those in other animals.

BACKGROUND: Statins increase the risk of new-onset type 2 diabetes mellitus. We aimed to assess whether this increase in risk is a consequence of inhibition of 3-hydroxy-3-methylglutaryl-CoA reductase (HMGCR), the intended drug target. METHODS: We used single nucleotide polymorphisms in the HMGCR gene, rs17238484 (for the main analysis) and rs12916 (for a subsidiary analysis) as proxies for HMGCR inhibition by statins. We examined associations of these variants with plasma lipid, glucose, and insulin concentrations; bodyweight; waist circumference; and prevalent and incident type 2 diabetes. Study-specific effect estimates per copy of each LDL-lowering allele were pooled by meta-analysis. These findings were compared with a meta-analysis of new-onset type 2 diabetes and bodyweight change data from randomised trials of statin drugs. The effects of statins in each randomised trial were assessed using meta-analysis. FINDINGS: Data were available for up to 223 463 individuals from 43 genetic studies. Each additional rs17238484-G allele was associated with a mean 0·06 mmol/L (95% CI 0·05-0·07) lower LDL cholesterol and higher body weight (0·30 kg, 0·18-0·43), waist circumference (0·32 cm, 0·16-0·47), plasma insulin concentration (1·62%, 0·53-2·72), and plasma glucose concentration (0·23%, 0·02-0·44). The rs12916 SNP had similar effects on LDL cholesterol, bodyweight, and waist circumference. The rs17238484-G allele seemed to be associated with higher risk of type 2 diabetes (odds ratio [OR] per allele 1·02, 95% CI 1·00-1·05); the rs12916-T allele association was consistent (1·06, 1·03-1·09). In 129 170 individuals in randomised trials, statins lowered LDL cholesterol by 0·92 mmol/L (95% CI 0·18-1·67) at 1-year of follow-up, increased bodyweight by 0·24 kg (95% CI 0·10-0·38 in all trials; 0·33 kg, 95% CI 0·24-0·42 in placebo or standard care controlled trials and -0·15 kg, 95% CI -0·39 to 0·08 in intensive-dose vs moderate-dose trials) at a mean of 4·2 years (range 1·9-6·7) of follow-up, and increased the odds of new-onset type 2 diabetes (OR 1·12, 95% CI 1·06-1·18 in all trials; 1·11, 95% CI 1·03-1·20 in placebo or standard care controlled trials and 1·12, 95% CI 1·04-1·22 in intensive-dose vs moderate dose trials). INTERPRETATION: The increased risk of type 2 diabetes noted with statins is at least partially explained by HMGCR inhibition. FUNDING: The funding sources are cited at the end of the paper.

Native Americans derive from a small number of Asian founders who likely arrived to the Americas via Beringia. However, additional details about the initial colonization of the Americas remain unclear. To investigate the pioneering phase in the Americas we analyzed a total of 623 complete mtDNAs from the Americas and Asia, including 20 new complete mtDNAs from the Americas and seven from Asia. This sequence data was used to direct high-resolution genotyping from 20 American and 26 Asian populations. Here we describe more genetic diversity within the founder population than was previously reported. The newly resolved phylogenetic structure suggests that ancestors of Native Americans paused when they reached Beringia, during which time New World founder lineages differentiated from their Asian sister-clades. This pause in movement was followed by a swift migration southward that distributed the founder types all the way to South America. The data also suggest more recent bi-directional gene flow between Siberia and the North American Arctic.

We review the evolution of domestic animals, emphasizing the effect of the earliest steps of domestication on its course. Using the first domesticated species, the dog (Canis familiaris), for illustration, we describe the evolutionary peculiarities during the historical domestication, such as the high level and wide range of diversity. We suggest that the process of earliest domestication via unconscious and later conscious selection of human-defined behavioral traits may accelerate phenotypic variations. The review is based on the results of a long-term experiment designed to reproduce early mammalian domestication in the silver fox (Vulpes vulpes) selected for tameability or amenability to domestication. We describe changes in behavior, morphology and physiology that appeared in the fox during its selection for tameability, which were similar to those observed in the domestic dog. Based on the data of the fox experiment and survey of relevant data, we discuss the developmental, genetic and possible molecular genetic mechanisms underlying these changes. We ascribe the causative role in evolutionary transformation of domestic animals to the selection for behavior and to the neurospecific regulatory genes it affects.

The sublethal effects of high predation risk on both prey behavior and physiology may have long-term consequences for prey population dynamics. We tested the hypothesis that snowshoe hares during the population decline are chronically stressed because of high predation risk whereas those during the population low are not, and that this has negative effects on both their physiology and demography. Snowshoe hares exhibit 10-yr population cycles; during declines, virtually every hare that dies is killed by a predator. We assessed the physiological responsiveness of the stress axis and of energy mobilization by subjecting hares during the population decline and low to a hormonal-challenge protocol. We monitored the population demography through live-trapping and assessed reproduction through a maternal-cage technique. During the 1990s' decline in the Yukon, Canada, hares were chronically stressed—as indicated by higher levels of free cortisol, reduced maximum corticosteroid-binding capacity, reduced testosterone response, reduced index of body condition, reduced leucocyte counts, increased overwinter body-mass loss, and increased glucose mobilization, relative to hares during the population low. This evidence is consistent with the explanation that predation risk, not high hare density or poor nutritional condition, accounted for the chronic stress and for the marked deterioration of reproduction during the decline. Reproduction and indices of stress physiology did not improve until predation risk declined. These findings may also account for the lag in recovery of hare reproduction after predator densities have declined and thus may implicate the long-term consequences of predation risk on prey populations beyond the immediate effects of predators on prey behavior and physiology.

How and when the Americas were populated remains contentious. Using ancient and modern genome-wide data, we found that the ancestors of all present-day Native Americans, including Athabascans and Amerindians, entered the Americas as a single migration wave from Siberia no earlier than 23 thousand years ago (ka) and after no more than an 8000-year isolation period in Beringia. After their arrival to the Americas, ancestral Native Americans diversified into two basal genetic branches around 13 ka, one that is now dispersed across North and South America and the other restricted to North America. Subsequent gene flow resulted in some Native Americans sharing ancestry with present-day East Asians (including Siberians) and, more distantly, Australo-Melanesians. Putative "Paleoamerican" relict populations, including the historical Mexican Pericúes and South American Fuego-Patagonians, are not directly related to modern Australo-Melanesians as suggested by the Paleoamerican Model.

Comparison of the genomes and proteomes of the two diptera Anopheles gambiae and Drosophila melanogaster, which diverged about 250 million years ago, reveals considerable similarities. However, numerous differences are also observed; some of these must reflect the selection and subsequent adaptation associated with different ecologies and life strategies. Almost half of the genes in both genomes are interpreted as orthologs and show an average sequence identity of about 56%, which is slightly lower than that observed between the orthologs of the pufferfish and human (diverged about 450 million years ago). This indicates that these two insects diverged considerably faster than vertebrates. Aligned sequences reveal that orthologous genes have retained only half of their intron/exon structure, indicating that intron gains or losses have occurred at a rate of about one per gene per 125 million years. Chromosomal arms exhibit significant remnants of homology between the two species, although only 34% of the genes colocalize in small "microsyntenic" clusters, and major interarm transfers as well as intra-arm shuffling of gene order are detected.

It is commonly thought that human genetic diversity in non-African populations was shaped primarily by an out-of-Africa dispersal 50-100 thousand yr ago (kya). Here, we present a study of 456 geographically diverse high-coverage Y chromosome sequences, including 299 newly reported samples. Applying ancient DNA calibration, we date the Y-chromosomal most recent common ancestor (MRCA) in Africa at 254 (95% CI 192-307) kya and detect a cluster of major non-African founder haplogroups in a narrow time interval at 47-52 kya, consistent with a rapid initial colonization model of Eurasia and Oceania after the out-of-Africa bottleneck. In contrast to demographic reconstructions based on mtDNA, we infer a second strong bottleneck in Y-chromosome lineages dating to the last 10 ky. We hypothesize that this bottleneck is caused by cultural changes affecting variance of reproductive success among males.

Importance: Blood pressure (BP) is a known risk factor for overall mortality and cardiovascular (CV)-specific fatal and nonfatal outcomes. It is uncertain which BP index is most strongly associated with these outcomes. Objective: To evaluate the association of BP indexes with death and a composite CV event. Design, Setting, and Participants: Longitudinal population-based cohort study of 11 135 adults from Europe, Asia, and South America with baseline observations collected from May 1988 to May 2010 (last follow-ups, August 2006-October 2016). Exposures: Blood pressure measured by an observer or an automated office machine; measured for 24 hours, during the day or the night; and the dipping ratio (nighttime divided by daytime readings). Main Outcomes and Measures: Multivariable-adjusted hazard ratios (HRs) expressed the risk of death or a CV event associated with BP increments of 20/10 mm Hg. Cardiovascular events included CV mortality combined with nonfatal coronary events, heart failure, and stroke. Improvement in model performance was assessed by the change in the area under the curve (AUC). Results: Among 11 135 participants (median age, 54.7 years, 49.3% women), 2836 participants died (18.5 per 1000 person-years) and 2049 (13.4 per 1000 person-years) experienced a CV event over a median of 13.8 years of follow-up. Both end points were significantly associated with all single systolic BP indexes (P < .001). For nighttime systolic BP level, the HR for total mortality was 1.23 (95% CI, 1.17-1.28) and for CV events, 1.36 (95% CI, 1.30-1.43). For the 24-hour systolic BP level, the HR for total mortality was 1.22 (95% CI, 1.16-1.28) and for CV events, 1.45 (95% CI, 1.37-1.54). With adjustment for any of the other systolic BP indexes, the associations of nighttime and 24-hour systolic BP with the primary outcomes remained statistically significant (HRs ranging from 1.17 [95% CI, 1.10-1.25] to 1.87 [95% CI, 1.62-2.16]). Base models that included single systolic BP indexes yielded an AUC of 0.83 for mortality and 0.84 for the CV outcomes. Adding 24-hour or nighttime systolic BP to base models that included other BP indexes resulted in incremental improvements in the AUC of 0.0013 to 0.0027 for mortality and 0.0031 to 0.0075 for the composite CV outcome. Adding any systolic BP index to models already including nighttime or 24-hour systolic BP did not significantly improve model performance. These findings were consistent for diastolic BP. Conclusions and Relevance: In this population-based cohort study, higher 24-hour and nighttime blood pressure measurements were significantly associated with greater risks of death and a composite CV outcome, even after adjusting for other office-based or ambulatory blood pressure measurements. Thus, 24-hour and nighttime blood pressure may be considered optimal measurements for estimating CV risk, although statistically, model improvement compared with other blood pressure indexes was small.

For pedigree-based quantitative trait loci (QTL) association analysis, a range of methods utilizing within-family variation such as transmission-disequilibrium test (TDT)-based methods have been developed. In scenarios where stratification is not a concern, methods exploiting between-family variation in addition to within-family variation, such as the measured genotype (MG) approach, have greater power. Application of MG methods can be computationally demanding (especially for large pedigrees), making genomewide scans practically infeasible. Here we suggest a novel approach for genomewide pedigree-based quantitative trait loci (QTL) association analysis: genomewide rapid association using mixed model and regression (GRAMMAR). The method first obtains residuals adjusted for family effects and subsequently analyzes the association between these residuals and genetic polymorphisms using rapid least-squares methods. At the final step, the selected polymorphisms may be followed up with the full measured genotype (MG) analysis. In a simulation study, we compared type 1 error, power, and operational characteristics of the proposed method with those of MG and TDT-based approaches. For moderately heritable (30%) traits in human pedigrees the power of the GRAMMAR and the MG approaches is similar and is much higher than that of TDT-based approaches. When using tabulated thresholds, the proposed method is less powerful than MG for very high heritabilities and pedigrees including large sibships like those observed in livestock pedigrees. However, there is little or no difference in empirical power of MG and the proposed method. In any scenario, GRAMMAR is much faster than MG and enables rapid analysis of hundreds of thousands of markers.

Fine structural details of glycans attached to the conserved N-glycosylation site significantly not only affect function of individual immunoglobulin G (IgG) molecules but also mediate inflammation at the systemic level. By analyzing IgG glycosylation in 5,117 individuals from four European populations, we have revealed very complex patterns of changes in IgG glycosylation with age. Several IgG glycans (including FA2B, FA2G2, and FA2BG2) changed considerably with age and the combination of these three glycans can explain up to 58% of variance in chronological age, significantly more than other markers of biological age like telomere lengths. The remaining variance in these glycans strongly correlated with physiological parameters associated with biological age. Thus, IgG glycosylation appears to be closely linked with both chronological and biological ages. Considering the important role of IgG glycans in inflammation, and because the observed changes with age promote inflammation, changes in IgG glycosylation also seem to represent a factor contributing to aging. SIGNIFICANCE STATEMENT: Glycosylation is the key posttranslational mechanism that regulates function of immunoglobulins, with multiple systemic repercussions to the immune system. Our study of IgG glycosylation in 5,117 individuals from four European populations has revealed very extensive and complex changes in IgG glycosylation with age. The combined index composed of only three glycans explained up to 58% of variance in age, considerably more than other biomarkers of age like telomere lengths. The remaining variance in these glycans strongly correlated with physiological parameters associated with biological age; thus, IgG glycosylation appears to be closely linked with both chronological and biological ages. The ability to measure human biological aging using molecular profiling has practical applications for diverse fields such as disease prevention and treatment, or forensics.

BACKGROUND: Over the last few years, genome-wide association (GWA) studies became a tool of choice for the identification of loci associated with complex traits. Currently, imputed single nucleotide polymorphisms (SNP) data are frequently used in GWA analyzes. Correct analysis of imputed data calls for the implementation of specific methods which take genotype imputation uncertainty into account. RESULTS: We developed the ProbABEL software package for the analysis of genome-wide imputed SNP data and quantitative, binary, and time-till-event outcomes under linear, logistic, and Cox proportional hazards models, respectively. For quantitative traits, the package also implements a fast two-step mixed model-based score test for association in samples with differential relationships, facilitating analysis in family-based studies, studies performed in human genetically isolated populations and outbred animal populations. CONCLUSIONS: ProbABEL package provides fast efficient way to analyze imputed data in genome-wide context and will facilitate future identification of complex trait loci.

In order to explore the diversity and selective signatures of duplication and deletion human copy-number variants (CNVs), we sequenced 236 individuals from 125 distinct human populations. We observed that duplications exhibit fundamentally different population genetic and selective signatures than deletions and are more likely to be stratified between human populations. Through reconstruction of the ancestral human genome, we identify megabases of DNA lost in different human lineages and pinpoint large duplications that introgressed from the extinct Denisova lineage now found at high frequency exclusively in Oceanic populations. We find that the proportion of CNV base pairs to single-nucleotide-variant base pairs is greater among non-Africans than it is among African populations, but we conclude that this difference is likely due to unique aspects of non-African population history as opposed to differences in CNV load.

Glycosylation of immunoglobulin G (IgG) influences IgG effector function by modulating binding to Fc receptors. To identify genetic loci associated with IgG glycosylation, we quantitated N-linked IgG glycans using two approaches. After isolating IgG from human plasma, we performed 77 quantitative measurements of N-glycosylation using ultra-performance liquid chromatography (UPLC) in 2,247 individuals from four European discovery populations. In parallel, we measured IgG N-glycans using MALDI-TOF mass spectrometry (MS) in a replication cohort of 1,848 Europeans. Meta-analysis of genome-wide association study (GWAS) results identified 9 genome-wide significant loci (P<2.27 × 10(-9)) in the discovery analysis and two of the same loci (B4GALT1 and MGAT3) in the replication cohort. Four loci contained genes encoding glycosyltransferases (ST6GAL1, B4GALT1, FUT8, and MGAT3), while the remaining 5 contained genes that have not been previously implicated in protein glycosylation (IKZF1, IL6ST-ANKRD55, ABCF2-SMARCD3, SUV420H1, and SMARCB1-DERL3). However, most of them have been strongly associated with autoimmune and inflammatory conditions (e.g., systemic lupus erythematosus, rheumatoid arthritis, ulcerative colitis, Crohn's disease, diabetes type 1, multiple sclerosis, Graves' disease, celiac disease, nodular sclerosis) and/or haematological cancers (acute lymphoblastic leukaemia, Hodgkin lymphoma, and multiple myeloma). Follow-up functional experiments in haplodeficient Ikzf1 knock-out mice showed the same general pattern of changes in IgG glycosylation as identified in the meta-analysis. As IKZF1 was associated with multiple IgG N-glycan traits, we explored biomarker potential of affected N-glycans in 101 cases with SLE and 183 matched controls and demonstrated substantial discriminative power in a ROC-curve analysis (area under the curve = 0.842). Our study shows that it is possible to identify new loci that control glycosylation of a single plasma protein using GWAS. The results may also provide an explanation for the reported pleiotropy and antagonistic effects of loci involved in autoimmune diseases and haematological cancer.

We surveyed nine diallelic polymorphic sites on the Y chromosomes of 1,544 individuals from Africa, Asia, Europe, Oceania, and the New World. Phylogenetic analyses of these nine sites resulted in a tree for 10 distinct Y haplotypes with a coalescence time of approximately 150,000 years. The 10 haplotypes were unevenly distributed among human populations: 5 were restricted to a particular continent, 2 were shared between Africa and Europe, 1 was present only in the Old World, and 2 were found in all geographic regions surveyed. The ancestral haplotype was limited to African populations. Random permutation procedures revealed statistically significant patterns of geographical structuring of this paternal genetic variation. The results of a nested cladistic analysis indicated that these geographical associations arose through a combination of processes, including restricted, recurrent gene flow (isolation by distance) and range expansions. We inferred that one of the oldest events in the nested cladistic analysis was a range expansion out of Africa which resulted in the complete replacement of Y chromosomes throughout the Old World, a finding consistent with many versions of the Out of Africa Replacement Model. A second and more recent range expansion brought Asian Y chromosomes back to Africa without replacing the indigenous African male gene pool. Thus, the previously observed high levels of Y chromosomal genetic diversity in Africa may be due in part to bidirectional population movements. Finally, a comparison of our results with those from nested cladistic analyses of human mtDNA and beta-globin data revealed different patterns of inferences for males and females concerning the relative roles of population history (range expansions) and population structure (recurrent gene flow), thereby adding a new sex-specific component to models of human evolution.

The three Drosophila EcR isoforms differ only at their N termini; thus, they share the conserved ligand-binding domain transcriptional activation function (AF2) and only differ in the unconserved A/B region, which contains a second, isoform-specific, activation function (AF1). We have developed a dominant-negative mutant EcR (EcR-DN), expressed it in flies with the GAL4/UAS system, and used it to block ecdysone signaling in eight tissues or groups of tissues. Localized EcR-DN arrests ecdysone-dependent development in the target cells and often--because of a molting checkpoint--arrests development globally. Simultaneously expressing individual wild-type EcR isoforms in the same target tissues suppresses the EcR-DN phenotype and identifies the rescuing isoform as sufficient to support the development of the target. Every isoform, and even an N-terminal truncated EcR that lacks any AF1, supports development in the fat body, eye discs, salivary glands, EH-secreting neurosecretory cells and in the dpp expression domain, implying that AF1 is dispensable in these tissues. By contrast, only EcR-A is able to support development in the margins of the wing discs, and only EcR-B2 can do so in the larval epidermis and the border cells of the developing egg chamber. In light of our results, the simplest explanations for the widespread spatial and temporal variations in EcR isoform titers appear untenable.