Epigenetic mechanisms of osteoarthritis risk in human skeletal development.

The epigenome, including the methylation of cytosine bases at CG dinucleotides, is intrinsically linked to transcriptional regulation. The tight regulation of gene expression during skeletal development is essential, with ~1/500 individuals born with skeletal abnormalities. Furthermore, increasing evidence is emerging to link age-associated complex genetic musculoskeletal diseases, including osteoarthritis (OA), to developmental factors including joint shape. Multiple studies have shown a functional role for DNA methylation in the genetic mechanisms of OA risk using articular cartilage samples taken from aged patients. Despite this, our knowledge of temporal changes to the methylome during human cartilage development has been limited. We quantified DNA methylation at ~700,000 individual CpGs across the epigenome of developing human articular cartilage in 72 samples ranging from 7-21 post-conception weeks, a time period that includes cavitation of the developing knee joint. We identified significant changes in 8% of all CpGs, and >9400 developmental differentially methylated regions (dDMRs). The largest hypermethylated dDMRs mapped to transcriptional regulators of early skeletal patterning including MEIS1 and IRX1. Conversely, the largest hypomethylated dDMRs mapped to genes encoding extracellular matrix proteins including SPON2 and TNXB and were enriched in chondrocyte enhancers. Significant correlations were identified between the expression of these genes and methylation within the hypomethylated dDMRs. We further identified 811 CpGs at which significant dimorphism was present between the male and female samples, with the majority (68%) being hypermethylated in female samples. Following imputation, we captured the genotype of these samples at >5 million variants and performed epigenome-wide methylation quantitative trait locus (mQTL) analysis. Colocalization analysis identified 26 loci at which genetic variants exhibited shared impacts upon methylation and OA genetic risk. This included loci which have been previously reported to harbour OA-mQTLs (including GDF5 and ALDH1A2), yet the majority (73%) were novel (including those mapping to CHST3, FGF1 and TEAD1). To our knowledge, this is the first extensive study of DNA methylation across human articular cartilage development. We identify considerable methylomic plasticity within the development of knee cartilage and report active epigenomic mediators of OA risk operating in prenatal joint tissues.

Epigenetic Contributions to Clinical Risk Prediction of Cardiovascular Disease.

BACKGROUND: Cardiovascular disease (CVD) is among the leading causes of death worldwide. The discovery of new omics biomarkers could help to improve risk stratification algorithms and expand our understanding of molecular pathways contributing to the disease. Here, ASSIGN-a cardiovascular risk prediction tool recommended for use in Scotland-was examined in tandem with epigenetic and proteomic features in risk prediction models in ≥12 657 participants from the Generation Scotland cohort. METHODS: Previously generated DNA methylation-derived epigenetic scores (EpiScores) for 109 protein levels were considered, in addition to both measured levels and an EpiScore for cTnI (cardiac troponin I). The associations between individual protein EpiScores and the CVD risk were examined using Cox regression (ncases≥1274; ncontrols≥11 383) and visualized in a tailored R application. Splitting the cohort into independent training (n=6880) and test (n=3659) subsets, a composite CVD EpiScore was then developed. RESULTS: Sixty-five protein EpiScores were associated with incident CVD independently of ASSIGN and the measured concentration of cTnI (P<0.05), over a follow-up of up to 16 years of electronic health record linkage. The most significant EpiScores were for proteins involved in metabolic, immune response, and tissue development/regeneration pathways. A composite CVD EpiScore (based on 45 protein EpiScores) was a significant predictor of CVD risk independent of ASSIGN and the concentration of cTnI (hazard ratio, 1.32; P=3.7×10-3; 0.3% increase in C-statistic). CONCLUSIONS: EpiScores for circulating protein levels are associated with CVD risk independent of traditional risk factors and may increase our understanding of the etiology of the disease.

Blood-based epigenome-wide analyses of 19 common disease states: A longitudinal, population-based linked cohort study of 18,413 Scottish individuals.

BACKGROUND: DNA methylation is a dynamic epigenetic mechanism that occurs at cytosine-phosphate-guanine dinucleotide (CpG) sites. Epigenome-wide association studies (EWAS) investigate the strength of association between methylation at individual CpG sites and health outcomes. Although blood methylation may act as a peripheral marker of common disease states, previous EWAS have typically focused only on individual conditions and have had limited power to discover disease-associated loci. This study examined the association of blood DNA methylation with the prevalence of 14 disease states and the incidence of 19 disease states in a single population of over 18,000 Scottish individuals. METHODS AND FINDINGS: DNA methylation was assayed at 752,722 CpG sites in whole-blood samples from 18,413 volunteers in the family-structured, population-based cohort study Generation Scotland (age range 18 to 99 years). EWAS tested for cross-sectional associations between baseline CpG methylation and 14 prevalent disease states, and for longitudinal associations between baseline CpG methylation and 19 incident disease states. Prevalent cases were self-reported on health questionnaires at the baseline. Incident cases were identified using linkage to Scottish primary (Read 2) and secondary (ICD-10) care records, and the censoring date was set to October 2020. The mean time-to-diagnosis ranged from 5.0 years (for chronic pain) to 11.7 years (for Coronavirus Disease 2019 (COVID-19) hospitalisation). The 19 disease states considered in this study were selected if they were present on the World Health Organisation's 10 leading causes of death and disease burden or included in baseline self-report questionnaires. EWAS models were adjusted for age at methylation typing, sex, estimated white blood cell composition, population structure, and 5 common lifestyle risk factors. A structured literature review was also conducted to identify existing EWAS for all 19 disease states tested. The MEDLINE, Embase, Web of Science, and preprint servers were searched to retrieve relevant articles indexed as of March 27, 2023. Fifty-four of approximately 2,000 indexed articles met our inclusion criteria: assayed blood-based DNA methylation, had >20 individuals in each comparison group, and examined one of the 19 conditions considered. First, we assessed whether the associations identified in our study were reported in previous studies. We identified 69 associations between CpGs and the prevalence of 4 conditions, of which 58 were newly described. The conditions were breast cancer, chronic kidney disease, ischemic heart disease, and type 2 diabetes mellitus. We also uncovered 64 CpGs that associated with the incidence of 2 disease states (COPD and type 2 diabetes), of which 56 were not reported in the surveyed literature. Second, we assessed replication across existing studies, which was defined as the reporting of at least 1 common site in >2 studies that examined the same condition. Only 6/19 disease states had evidence of such replication. The limitations of this study include the nonconsideration of medication data and a potential lack of generalizability to individuals that are not of Scottish and European ancestry. CONCLUSIONS: We discovered over 100 associations between blood methylation sites and common disease states, independently of major confounding risk factors, and a need for greater standardisation among EWAS on human disease.

Development and validation of DNA methylation scores in two European cohorts augment 10-year risk prediction of type 2 diabetes.

Type 2 diabetes mellitus (T2D) presents a major health and economic burden that could be alleviated with improved early prediction and intervention. While standard risk factors have shown good predictive performance, we show that the use of blood-based DNA methylation information leads to a significant improvement in the prediction of 10-year T2D incidence risk. Previous studies have been largely constrained by linear assumptions, the use of cytosine-guanine pairs one-at-a-time and binary outcomes. We present a flexible approach (via an R package, MethylPipeR) based on a range of linear and tree-ensemble models that incorporate time-to-event data for prediction. Using the Generation Scotland cohort (training set ncases = 374, ncontrols = 9,461; test set ncases = 252, ncontrols = 4,526) our best-performing model (area under the receiver operating characteristic curve (AUC) = 0.872, area under the precision-recall curve (PRAUC) = 0.302) showed notable improvement in 10-year onset prediction beyond standard risk factors (AUC = 0.839, precision-recall AUC = 0.227). Replication was observed in the German-based KORA study (n = 1,451, ncases = 142, P = 1.6 × 10-5).

Refining epigenetic prediction of chronological and biological age.

BACKGROUND: Epigenetic clocks can track both chronological age (cAge) and biological age (bAge). The latter is typically defined by physiological biomarkers and risk of adverse health outcomes, including all-cause mortality. As cohort sample sizes increase, estimates of cAge and bAge become more precise. Here, we aim to develop accurate epigenetic predictors of cAge and bAge, whilst improving our understanding of their epigenomic architecture. METHODS: First, we perform large-scale (N = 18,413) epigenome-wide association studies (EWAS) of chronological age and all-cause mortality. Next, to create a cAge predictor, we use methylation data from 24,674 participants from the Generation Scotland study, the Lothian Birth Cohorts (LBC) of 1921 and 1936, and 8 other cohorts with publicly available data. In addition, we train a predictor of time to all-cause mortality as a proxy for bAge using the Generation Scotland cohort (1214 observed deaths). For this purpose, we use epigenetic surrogates (EpiScores) for 109 plasma proteins and the 8 component parts of GrimAge, one of the current best epigenetic predictors of survival. We test this bAge predictor in four external cohorts (LBC1921, LBC1936, the Framingham Heart Study and the Women's Health Initiative study). RESULTS: Through the inclusion of linear and non-linear age-CpG associations from the EWAS, feature pre-selection in advance of elastic net regression, and a leave-one-cohort-out (LOCO) cross-validation framework, we obtain cAge prediction with a median absolute error equal to 2.3 years. Our bAge predictor was found to slightly outperform GrimAge in terms of the strength of its association to survival (HRGrimAge = 1.47 [1.40, 1.54] with p = 1.08 × 10-52, and HRbAge = 1.52 [1.44, 1.59] with p = 2.20 × 10-60). Finally, we introduce MethylBrowsR, an online tool to visualise epigenome-wide CpG-age associations. CONCLUSIONS: The integration of multiple large datasets, EpiScores, non-linear DNAm effects, and new approaches to feature selection has facilitated improvements to the blood-based epigenetic prediction of biological and chronological age.

Feasibility and ethics of using data from the Scottish newborn blood spot archive for research.

Background: Newborn heel prick blood spots are routinely used to screen for inborn errors of metabolism and life-limiting inherited disorders. The potential value of secondary data from newborn blood spot archives merits ethical consideration and assessment of feasibility for public benefit. Early life exposures and behaviours set health trajectories in childhood and later life. The newborn blood spot is potentially well placed to create an unbiased and cost-effective population-level retrospective birth cohort study. Scotland has retained newborn blood spots for all children born since 1965, around 3 million in total. However, a moratorium on research access is currently in place, pending public consultation. Methods: We conducted a Citizens' Jury as a first step to explore whether research use of newborn blood spots was in the public interest. We also assessed the feasibility and value of extracting research data from dried blood spots for predictive medicine. Results: Jurors delivered an agreed verdict that conditional research access to the newborn blood spots was in the public interest. The Chief Medical Officer for Scotland authorised restricted lifting of the current research moratorium to allow a feasibility study. Newborn blood spots from consented Generation Scotland volunteers were retrieved and their potential for both epidemiological and biological research demonstrated. Conclusions: Through the Citizens' Jury, we have begun to identify under what conditions, if any, should researchers in Scotland be granted access to the archive. Through the feasibility study, we have demonstrated the potential value of research access for health data science and predictive medicine.

A comparison of blood and brain-derived ageing and inflammation-related DNA methylation signatures and their association with microglial burdens.

Inflammation and ageing-related DNA methylation patterns in the blood have been linked to a variety of morbidities, including cognitive decline and neurodegenerative disease. However, it is unclear how these blood-based patterns relate to patterns within the brain and how each associates with central cellular profiles. In this study, we profiled DNA methylation in both the blood and in five post mortem brain regions (BA17, BA20/21, BA24, BA46 and hippocampus) in 14 individuals from the Lothian Birth Cohort 1936. Microglial burdens were additionally quantified in the same brain regions. DNA methylation signatures of five epigenetic ageing biomarkers ('epigenetic clocks'), and two inflammatory biomarkers (methylation proxies for C-reactive protein and interleukin-6) were compared across tissues and regions. Divergent associations between the inflammation and ageing signatures in the blood and brain were identified, depending on region assessed. Four out of the five assessed epigenetic age acceleration measures were found to be highest in the hippocampus (ß range = 0.83-1.14, p ≤ 0.02). The inflammation-related DNA methylation signatures showed no clear variation across brain regions. Reactive microglial burdens were found to be highest in the hippocampus (ß = 1.32, p = 5 × 10-4 ); however, the only association identified between the blood- and brain-based methylation signatures and microglia was a significant positive association with acceleration of one epigenetic clock (termed DNAm PhenoAge) averaged over all five brain regions (ß = 0.40, p = 0.002). This work highlights a potential vulnerability of the hippocampus to epigenetic ageing and provides preliminary evidence of a relationship between DNA methylation signatures in the brain and differences in microglial burdens.

DNA methylation in relation to gestational age and brain dysmaturation in preterm infants.

Preterm birth is associated with dysconnectivity of structural brain networks and is a leading cause of neurocognitive impairment in childhood. Variation in DNA methylation is associated with early exposure to extrauterine life but there has been little research exploring its relationship with brain development. Using genome-wide DNA methylation data from the saliva of 258 neonates, we investigated the impact of gestational age on the methylome and performed functional analysis to identify enriched gene sets from probes that contributed to differentially methylated probes or regions. We tested the hypothesis that variation in DNA methylation could underpin the association between low gestational age at birth and atypical brain development by linking differentially methylated probes with measures of white matter connectivity derived from diffusion MRI metrics: peak width skeletonized mean diffusivity, peak width skeletonized fractional anisotropy and peak width skeletonized neurite density index. Gestational age at birth was associated with widespread differential methylation at term equivalent age, with genome-wide significant associations observed for 8870 CpG probes (P < 3.6 × 10-8) and 1767 differentially methylated regions. Functional analysis identified 14 enriched gene ontology terms pertaining to cell-cell contacts and cell-extracellular matrix contacts. Principal component analysis of probes with genome-wide significance revealed a first principal component that explained 23.5% of the variance in DNA methylation, and this was negatively associated with gestational age at birth. The first principal component was associated with peak width of skeletonized mean diffusivity (ß = 0.349, P = 8.37 × 10-10) and peak width skeletonized neurite density index (ß = 0.364, P = 4.15 × 10-5), but not with peak width skeletonized fraction anisotropy (ß = -0.035, P = 0.510); these relationships mirrored the imaging metrics' associations with gestational age at birth. Low gestational age at birth has a profound and widely distributed effect on the neonatal saliva methylome that is apparent at term equivalent age. Enriched gene ontology terms related to cell-cell contacts reveal pathways that could mediate the effect of early life environmental exposures on development. Finally, associations between differential DNA methylation and image markers of white matter tract microstructure suggest that variation in DNA methylation may provide a link between preterm birth and the dysconnectivity of developing brain networks that characterizes atypical brain development in preterm infants.

Epigenetic predictors of lifestyle traits applied to the blood and brain.

Modifiable lifestyle factors influence the risk of developing many neurological diseases. These factors have been extensively linked with blood-based genome-wide DNA methylation, but it is unclear if the signatures from blood translate to the target tissue of interest-the brain. To investigate this, we apply blood-derived epigenetic predictors of four lifestyle traits to genome-wide DNA methylation from five post-mortem brain regions and the last blood sample prior to death in 14 individuals in the Lothian Birth Cohort 1936. Using these matched samples, we found that correlations between blood and brain DNA methylation scores for smoking, high-density lipoprotein cholesterol, alcohol and body mass index were highly variable across brain regions. Smoking scores in the dorsolateral prefrontal cortex had the strongest correlations with smoking scores in blood (r = 0.5, n = 14, P = 0.07) and smoking behaviour (r = 0.56, n = 9, P = 0.12). This was also the brain region which exhibited the largest correlations for DNA methylation at site cg05575921 - the single strongest correlate of smoking in blood-in relation to blood (r = 0.61, n = 14, P = 0.02) and smoking behaviour (r = -0.65, n = 9, P = 0.06). This suggested a particular vulnerability to smoking-related differential methylation in this region. Our work contributes to understanding how lifestyle factors affect the brain and suggest that lifestyle-related DNA methylation is likely to be both brain region dependent and in many cases poorly proxied for by blood. Though these pilot data provide a rarely-available opportunity for the comparison of methylation patterns across multiple brain regions and the blood, due to the limited sample size available our results must be considered as preliminary and should therefore be used as a basis for further investigation.

Genetic mechanisms of critical illness in COVID-19.

Host-mediated lung inflammation is present1, and drives mortality2, in the critical illness caused by coronavirus disease 2019 (COVID-19). Host genetic variants associated with critical illness may identify mechanistic targets for therapeutic development3. Here we report the results of the GenOMICC (Genetics Of Mortality In Critical Care) genome-wide association study in 2,244 critically ill patients with COVID-19 from 208 UK intensive care units. We have identified and replicated the following new genome-wide significant associations: on chromosome 12q24.13 (rs10735079, P = 1.65 × 10-8) in a gene cluster that encodes antiviral restriction enzyme activators (OAS1, OAS2 and OAS3); on chromosome 19p13.2 (rs74956615, P = 2.3 × 10-8) near the gene that encodes tyrosine kinase 2 (TYK2); on chromosome 19p13.3 (rs2109069, P = 3.98 ×  10-12) within the gene that encodes dipeptidyl peptidase 9 (DPP9); and on chromosome 21q22.1 (rs2236757, P = 4.99 × 10-8) in the interferon receptor gene IFNAR2. We identified potential targets for repurposing of licensed medications: using Mendelian randomization, we found evidence that low expression of IFNAR2, or high expression of TYK2, are associated with life-threatening disease; and transcriptome-wide association in lung tissue revealed that high expression of the monocyte-macrophage chemotactic receptor CCR2 is associated with severe COVID-19. Our results identify robust genetic signals relating to key host antiviral defence mechanisms and mediators of inflammatory organ damage in COVID-19. Both mechanisms may be amenable to targeted treatment with existing drugs. However, large-scale randomized clinical trials will be essential before any change to clinical practice.

Assessment of dried blood spots for DNA methylation profiling.

Background: DNA methylation reflects health-related environmental exposures and genetic risk, providing insights into aetiological mechanisms and potentially predicting disease onset, progression and treatment response. An increasingly recognised need for large-scale, longitudinally-profiled samples collected world-wide has made the development of efficient and straightforward sample collection and storage procedures a pressing issue. An alternative to the low-temperature storage of EDTA tubes of venous blood samples, which are frequently the source of the DNA used in such studies, is to collect and store at room temperature blood samples using purpose built filter paper, such as Whatman FTA® cards. Our goal was to determine whether DNA stored in this manner can be used to generate DNA methylation profiles comparable to those generated using blood samples frozen in EDTA tubes. Methods: DNA methylation profiles were obtained from matched EDTA tube and Whatman FTA® card whole-blood samples from 62 Generation Scotland: Scottish Family Health Study participants using the Infinium HumanMethylation450 BeadChip. Multiple quality control procedures were implemented, the relationship between the two sample types assessed, and epigenome-wide association studies (EWASs) performed for smoking status, age and the interaction between these variables and sample storage method. Results: Dried blood spot (DBS) DNA methylation profiles were of good quality and DNA methylation profiles from matched DBS and EDTA tube samples were highly correlated (mean r = 0.991) and could distinguish between participants. EWASs replicated established associations for smoking and age, with no evidence for moderation by storage method. Conclusions: Our results support the use of Whatman FTA® cards for collecting and storing blood samples for DNA methylation profiling. This approach is likely to be particularly beneficial for large-scale studies and those carried out in areas where freezer access is limited. Furthermore, our results will inform consideration of the use of newborn heel prick DBSs for research use.

Age-related gene expression changes, and transcriptome wide association study of physical and cognitive aging traits, in the Lothian Birth Cohort 1936.

Gene expression is influenced by both genetic variants and the environment. As individuals age, changes in gene expression may be associated with decline in physical and cognitive abilities. We measured transcriptome-wide expression levels in lymphoblastoid cell lines derived from members of the Lothian Birth Cohort 1936 at mean ages 70 and 76 years. Changes in gene expression levels were identified for 1,741 transcripts in 434 individuals. Gene Ontology enrichment analysis indicated an enrichment of biological processes involved in the immune system. Transcriptome-wide association analysis was performed for eleven cognitive, fitness, and biomedical aging-related traits at age 70 years (N=665 to 781) and with mortality. Transcripts for genes (F2RL3, EMILIN1 and CDC42BPA) previously identified as being differentially methylated or expressed in smoking or smoking-related cancers were overexpressed in smokers compared to non-smokers and the expression of transcripts for genes (HERPUD1, GAB2, FAM167A and GLS) previously associated with stress response, autoimmune disease and cancer were associated with telomere length. No associations between expression levels and other traits, or mortality were identified.

Transcription controls growth, cell kinetics and cholesterol supply to sustain ACTH responses.

Chronic ACTH exposure is associated with adrenal hypertrophy and steroidogenesis. The underlying molecular processes in mice have been analysed by microarray, histological and immunohistochemical techniques. Synacthen infused for 2 weeks markedly increased adrenal mass and plasma corticosterone levels. Microarray analysis found greater than 2-fold changes in expression of 928 genes (P < 0.001; 397 up, 531 down). These clustered in pathways involved in signalling, sterol/lipid metabolism, cell proliferation/hypertrophy and apoptosis. Signalling genes included some implicated in adrenal adenomas but also upregulated genes associated with cyclic AMP and downregulated genes associated with aldosterone synthesis. Sterol metabolism genes were those promoting cholesterol supply (Scarb1, Sqle, Apoa1) and disposal (Cyp27a1, Cyp7b1). Oil red O staining showed lipid depletion consistent with reduced expression of genes involved in lipid synthesis. Genes involved in steroidogenesis (Star, Cyp11a1, Cyp11b1) were modestly affected (P < 0.05; <1.3-fold). Increased Ki67, Ccna2, Ccnb2 and Tk1 expression complemented immunohistochemical evidence of a 3-fold change in cell proliferation. Growth arrest genes, Cdkn1a and Cdkn1c, which are known to be active in hypertrophied cells, were increased >4-fold and cross-sectional area of fasciculata cells was 2-fold greater. In contrast, genes associated with apoptosis (eg Casp12, Clu,) were downregulated and apoptotic cells (Tunel staining) were fewer (P < 0.001) and more widely distributed throughout the cortex. In summary, long-term steroidogenesis with ACTH excess is sustained by genes controlling cholesterol supply and adrenal mass. ACTH effects on adrenal morphology and genes controlling cell hypertrophy, proliferation and apoptosis suggest the involvement of different cell types and separate molecular pathways.

How many biological replicates are needed in an RNA-seq experiment and which differential expression tool should you use?

RNA-seq is now the technology of choice for genome-wide differential gene expression experiments, but it is not clear how many biological replicates are needed to ensure valid biological interpretation of the results or which statistical tools are best for analyzing the data. An RNA-seq experiment with 48 biological replicates in each of two conditions was performed to answer these questions and provide guidelines for experimental design. With three biological replicates, nine of the 11 tools evaluated found only 20%-40% of the significantly differentially expressed (SDE) genes identified with the full set of 42 clean replicates. This rises to >85% for the subset of SDE genes changing in expression by more than fourfold. To achieve >85% for all SDE genes regardless of fold change requires more than 20 biological replicates. The same nine tools successfully control their false discovery rate at ≲5% for all numbers of replicates, while the remaining two tools fail to control their FDR adequately, particularly for low numbers of replicates. For future RNA-seq experiments, these results suggest that at least six biological replicates should be used, rising to at least 12 when it is important to identify SDE genes for all fold changes. If fewer than 12 replicates are used, a superior combination of true positive and false positive performances makes edgeR and DESeq2 the leading tools. For higher replicate numbers, minimizing false positives is more important and DESeq marginally outperforms the other tools.

Statistical models for RNA-seq data derived from a two-condition 48-replicate experiment.

MOTIVATION: High-throughput RNA sequencing (RNA-seq) is now the standard method to determine differential gene expression. Identifying differentially expressed genes crucially depends on estimates of read-count variability. These estimates are typically based on statistical models such as the negative binomial distribution, which is employed by the tools edgeR, DESeq and cuffdiff. Until now, the validity of these models has usually been tested on either low-replicate RNA-seq data or simulations. RESULTS: A 48-replicate RNA-seq experiment in yeast was performed and data tested against theoretical models. The observed gene read counts were consistent with both log-normal and negative binomial distributions, while the mean-variance relation followed the line of constant dispersion parameter of ∼0.01. The high-replicate data also allowed for strict quality control and screening of 'bad' replicates, which can drastically affect the gene read-count distribution. AVAILABILITY AND IMPLEMENTATION: RNA-seq data have been submitted to ENA archive with project ID PRJEB5348. CONTACT: g.j.barton@dundee.ac.uk.

Linking genomics and ecology to investigate the complex evolution of an invasive Drosophila pest.

Drosophilid fruit flies have provided science with striking cases of behavioral adaptation and genetic innovation. A recent example is the invasive pest Drosophila suzukii, which, unlike most other Drosophila, lays eggs and feeds on undamaged, ripening fruits. This not only poses a serious threat for fruit cultivation but also offers an interesting model to study evolution of behavioral innovation. We developed genome and transcriptome resources for D. suzukii. Coupling analyses of these data with field observations, we propose a hypothesis of the origin of its peculiar ecology. Using nuclear and mitochondrial phylogenetic analyses, we confirm its Asian origin and reveal a surprising sister relationship between the eugracilis and the melanogaster subgroups. Although the D. suzukii genome is comparable in size and repeat content to other Drosophila species, it has the lowest nucleotide substitution rate among the species analyzed in this study. This finding is compatible with the overwintering diapause of D. suzukii, which results in a reduced number of generations per year compared with its sister species. Genome-scale relaxed clock analyses support a late Miocene origin of D. suzukii, concomitant with paleogeological and climatic conditions that suggest an adaptation to temperate montane forests, a hypothesis confirmed by field trapping. We propose a causal link between the ecological adaptations of D. suzukii in its native habitat and its invasive success in Europe and North America.

The genome of the heartworm, Dirofilaria immitis, reveals drug and vaccine targets.

The heartworm Dirofilaria immitis is an important parasite of dogs. Transmitted by mosquitoes in warmer climatic zones, it is spreading across southern Europe and the Americas at an alarming pace. There is no vaccine, and chemotherapy is prone to complications. To learn more about this parasite, we have sequenced the genomes of D. immitis and its endosymbiont Wolbachia. We predict 10,179 protein coding genes in the 84.2 Mb of the nuclear genome, and 823 genes in the 0.9-Mb Wolbachia genome. The D. immitis genome harbors neither DNA transposons nor active retrotransposons, and there is very little genetic variation between two sequenced isolates from Europe and the United States. The differential presence of anabolic pathways such as heme and nucleotide biosynthesis hints at the intricate metabolic interrelationship between the heartworm and Wolbachia. Comparing the proteome of D. immitis with other nematodes and with mammalian hosts, we identify families of potential drug targets, immune modulators, and vaccine candidates. This genome sequence will support the development of new tools against dirofilariasis and aid efforts to combat related human pathogens, the causative agents of lymphatic filariasis and river blindness.

Cyp11b1 null mouse, a model of congenital adrenal hyperplasia.

Patients with congenital adrenal hyperplasia arising from mutations of 11beta-hydroxylase, the final enzyme in the glucocorticoid biosynthetic pathway, exhibit glucocorticoid deficiency, adrenal hyperplasia driven by unsuppressed hypothalamo-pituitary-adrenal activity, and excess mineralocorticoid activity caused by the accumulation of deoxycorticosterone. A mouse model, in which exons 3-7 of Cyp11b1 (the gene encoding 11beta-hydroxylase) were replaced with cDNA encoding enhanced cyan fluorescent protein, was generated to investigate the underlying disease mechanisms. Enhanced cyan fluorescent protein was expressed appropriately in the zona fasciculata of the adrenal gland, and targeted knock-out was confirmed by urinary steroid profiles and, immunocytochemically, by the absence of 11beta-hydroxylase. The null mice exhibited glucocorticoid deficiency, mineralocorticoid excess, adrenal hyperplasia, mild hypertension, and hypokalemia. They also displayed glucose intolerance. Because rodents do not synthesize adrenal androgens, changes in reproductive function such as genital virilization of females were not anticipated. However, adult homozygote females were infertile, their ovaries showing an absence of corpora lutea and a central proliferation of disorganized steroidogenic tissue. Null females responded normally to superovulation, suggesting that raised systemic progesterone levels also contribute to infertility problems. The model reveals previously unrecognized phenotypic subtleties of congenital adrenal hyperplasia.

A switch in the mechanism of hypertension in the syndrome of apparent mineralocorticoid excess.

The syndrome of apparent mineralocorticoid excess arises from nonfunctional mutations in 11beta-hydroxysteroid dehydrogenase type 2 (11betaHSD2), an enzyme that inactivates cortisol and confers aldosterone specificity on the mineralocorticoid receptor. Loss of 11betaHSD2 permits glucocorticoids to activate the mineralocorticoid receptor, and the hypertension in the syndrome is presumed to arise from volume expansion secondary to renal sodium retention. An 11betaHSD2 null mouse was generated on an inbred C57BL/6J genetic background, allowing survival to adulthood. 11betaHSD2(-/-) mice had BP approximately 20 mmHg higher on average compared with wild-type mice but were volume contracted, not volume expanded as expected. Initially, impaired sodium excretion associated with increased activity of the epithelial sodium channel was observed. By 80 days of age, however, channel activity was abolished and 11betaHSD2(-/-) mice lost salt. Despite the natriuresis, hypertension remained but was not attributable to intrinsic vascular dysfunction. Instead, urinary catecholamine levels in 11betaHSD2(-/-) mice were double those in wild-type mice, and alpha1-adrenergic receptor blockade rescued the hypertensive phenotype, suggesting that vasoconstriction contributes to the sustained hypertension in this model. In summary, it is proposed that renal sodium retention remains a key event in apparent mineralocorticoid excess but that the accompanying hypertension changes from a renal to a vascular etiology over time.