The genome of the leaf-cutting ant Acromyrmex echinatior suggests key adaptations to advanced social life and fungus farming.

We present a high-quality (>100× depth) Illumina genome sequence of the leaf-cutting ant Acromyrmex echinatior, a model species for symbiosis and reproductive conflict studies. We compare this genome with three previously sequenced genomes of ants from different subfamilies and focus our analyses on aspects of the genome likely to be associated with known evolutionary changes. The first is the specialized fungal diet of A. echinatior, where we find gene loss in the ant's arginine synthesis pathway, loss of detoxification genes, and expansion of a group of peptidase proteins. One of these is a unique ant-derived contribution to the fecal fluid, which otherwise consists of "garden manuring" fungal enzymes that are unaffected by ant digestion. The second is multiple mating of queens and ejaculate competition, which may be associated with a greatly expanded nardilysin-like peptidase gene family. The third is sex determination, where we could identify only a single homolog of the feminizer gene. As other ants and the honeybee have duplications of this gene, we hypothesize that this may partly explain the frequent production of diploid male larvae in A. echinatior. The fourth is the evolution of eusociality, where we find a highly conserved ant-specific profile of neuropeptide genes that may be related to caste determination. These first analyses of the A. echinatior genome indicate that considerable genetic changes are likely to have accompanied the transition from hunter-gathering to agricultural food production 50 million years ago, and the transition from single to multiple queen mating 10 million years ago.

A burst of segmental duplications in the genome of the African great ape ancestor

It is generally accepted that the extent of phenotypic change between human and great apes is dissonant with the rate of molecular change. Between these two groups, proteins are virtually identical, cytogenetically there are few rearrangements that distinguish ape-human chromosomes, and rates of single-base-pair change and retrotransposon activity have slowed particularly within hominid lineages when compared to rodents or monkeys. Studies of gene family evolution indicate that gene loss and gain are enriched within the primate lineage. Here, we perform a systematic analysis of duplication content of four primate genomes (macaque, orang-utan, chimpanzee and human) in an effort to understand the pattern and rates of genomic duplication during hominid evolution. We find that the ancestral branch leading to human and African great apes shows the most significant increase in duplication activity both in terms of base pairs and in terms of events. This duplication acceleration within the ancestral species is significant when compared to lineage-specific rate estimates even after accounting for copy-number polymorphism and homoplasy. We discover striking examples of recurrent and independent gene-containing duplications within the gorilla and chimpanzee that are absent in the human lineage. Our results suggest that the evolutionary properties of copy-number mutation differ significantly from other forms of genetic mutation and, in contrast to the hominid slowdown of single-base-pair mutations, there has been a genomic burst of duplication activity at this period during human evolution.

The genome sequencing of an albino Western lowland gorilla reveals inbreeding in the wild

BACKGROUND: The only known albino gorilla, named Snowflake, was a male wild born individual from Equatorial Guinea who lived at the Barcelona Zoo for almost 40 years. He was diagnosed with non-syndromic oculocutaneous albinism, i.e. white hair, light eyes, pink skin, photophobia and reduced visual acuity. Despite previous efforts to explain the genetic cause, this is still unknown. Here, we study the genetic cause of his albinism and making use of whole genome sequencing data we find a higher inbreeding coefficient compared to other gorillas.RESULTS: We successfully identified the causal genetic variant for Snowflake's albinism, a non-synonymous single nucleotide variant located in a transmembrane region of SLC45A2. This transporter is known to be involved in oculocutaneous albinism type 4 (OCA4) in humans. We provide experimental evidence that shows that this amino acid replacement alters the membrane spanning capability of this transmembrane region. Finally, we provide a comprehensive study of genome-wide patterns of autozygogosity revealing that Snowflake's parents were related, being this the first report of inbreeding in a wild born Western lowland gorilla.CONCLUSIONS: In this study we demonstrate how the use of whole genome sequencing can be extended to link genotype and phenotype in non-model organisms and it can be a powerful tool in conservation genetics (e.g., inbreeding and genetic diversity) with the expected decrease in sequencing cost.

Global Genomic Diversity of Human Papillomavirus 6 Based on 724 Isolates and 190 Complete Genome Sequences.

Human papillomavirus type 6 (HPV6) is the major etiological agent of anogenital warts and laryngeal papillomas and has been included in both the quadrivalent and nonavalent prophylactic HPV vaccines. This study investigated the global genomic diversity of HPV6, using 724 isolates and 190 complete genomes from six continents, and the association of HPV6 genomic variants with geographical location, anatomical site of infection/disease, and gender. Initially, a 2,800-bp E5a-E5b-L1-LCR fragment was sequenced from 492/530 (92.8%) HPV6-positive samples collected for this study. Among them, 130 exhibited at least one single nucleotide polymorphism (SNP), indel, or amino acid change in the E5a-E5b-L1-LCR fragment and were sequenced in full. A global alignment and maximum likelihood tree of 190 complete HPV6 genomes (130 fully sequenced in this study and 60 obtained from sequence repositories) revealed two variant lineages, A and B, and five B sublineages: B1, B2, B3, B4, and B5. HPV6 (sub)lineage-specific SNPs and a 960-bp representative region for whole-genome-based phylogenetic clustering within the L2 open reading frame were identified. Multivariate logistic regression analysis revealed that lineage B predominated globally. Sublineage B3 was more common in Africa and North and South America, and lineage A was more common in Asia. Sublineages B1 and B3 were associated with anogenital infections, indicating a potential lesion-specific predilection of some HPV6 sublineages. Females had higher odds for infection with sublineage B3 than males. In conclusion, a global HPV6 phylogenetic analysis revealed the existence of two variant lineages and five sublineages, showing some degree of ethnogeographic, gender, and/or disease predilection in their distribution. IMPORTANCE: This study established the largest database of globally circulating HPV6 genomic variants and contributed a total of 130 new, complete HPV6 genome sequences to available sequence repositories. Two HPV6 variant lineages and five sublineages were identified and showed some degree of association with geographical location, anatomical site of infection/disease, and/or gender. We additionally identified several HPV6 lineage- and sublineage-specific SNPs to facilitate the identification of HPV6 variants and determined a representative region within the L2 gene that is suitable for HPV6 whole-genome-based phylogenetic analysis. This study complements and significantly expands the current knowledge of HPV6 genetic diversity and forms a comprehensive basis for future epidemiological, evolutionary, functional, pathogenicity, vaccination, and molecular assay development studies.

Mitochondrial uncoupling as a regulator of life-history trajectories in birds: an experimental study in the zebra finch.

Mitochondria have a fundamental role in the transduction of energy from food into ATP. The coupling between food oxidation and ATP production is never perfect, but may nevertheless be of evolutionary significance. The 'uncoupling to survive' hypothesis suggests that 'mild' mitochondrial uncoupling evolved as a protective mechanism against the excessive production of damaging reactive oxygen species (ROS). Because resource allocation and ROS production are thought to shape animal life histories, alternative life-history trajectories might be driven by individual variation in the degree of mitochondrial uncoupling. We tested this hypothesis in a small bird species, the zebra finch (Taeniopygia guttata), by treating adults with the artificial mitochondrial uncoupler 2,4-dinitrophenol (DNP) over a 32-month period. In agreement with our expectations, the uncoupling treatment increased metabolic rate. However, we found no evidence that treated birds enjoyed lower oxidative stress levels or greater survival rates, in contrast to previous results in other taxa. In vitro experiments revealed lower sensitivity of ROS production to DNP in mitochondria isolated from skeletal muscles of zebra finch than mouse. In addition, we found significant reductions in the number of eggs laid and in the inflammatory immune response in treated birds. Altogether, our data suggest that the 'uncoupling to survive' hypothesis may not be applicable for zebra finches, presumably because of lower effects of mitochondrial uncoupling on mitochondrial ROS production in birds than in mammals. Nevertheless, mitochondrial uncoupling appeared to be a potential life-history regulator of traits such as fecundity and immunity at adulthood, even with food supplied ad libitum.

Mfn2 downregulation in excitotoxicity causes mitochondrial dysfunction and delayed neuronal death.

Mitochondrial fusion and fission is a dynamic process critical for the maintenance of mitochondrial function and cell viability. During excitotoxicity neuronal mitochondria are fragmented, but the mechanism underlying this process is poorly understood. Here, we show that Mfn2 is the only member of the mitochondrial fusion/fission machinery whose expression is reduced in in vitro and in vivo models of excitotoxicity. Whereas in cortical primary cultures, Drp1 recruitment to mitochondria plays a primordial role in mitochondrial fragmentation in an early phase that can be reversed once the insult has ceased, Mfn2 downregulation intervenes in a delayed mitochondrial fragmentation phase that progresses even when the insult has ceased. Downregulation of Mfn2 causes mitochondrial dysfunction, altered calcium homeostasis, and enhanced Bax translocation to mitochondria, resulting in delayed neuronal death. We found that transcription factor MEF2 regulates basal Mfn2 expression in neurons and that excitotoxicity-dependent degradation of MEF2 causes Mfn2 downregulation. Thus, Mfn2 reduction is a late event in excitotoxicity and its targeting may help to reduce excitotoxic damage and increase the currently short therapeutic window in stroke.

Integration of genome-wide association studies with biological knowledge identifies six novel genes related to kidney function.

In conducting genome-wide association studies (GWAS), analytical approaches leveraging biological information may further understanding of the pathophysiology of clinical traits. To discover novel associations with estimated glomerular filtration rate (eGFR), a measure of kidney function, we developed a strategy for integrating prior biological knowledge into the existing GWAS data for eGFR from the CKDGen Consortium. Our strategy focuses on single nucleotide polymorphism (SNPs) in genes that are connected by functional evidence, determined by literature mining and gene ontology (GO) hierarchies, to genes near previously validated eGFR associations. It then requires association thresholds consistent with multiple testing, and finally evaluates novel candidates by independent replication. Among the samples of European ancestry, we identified a genome-wide significant SNP in FBXL20 (P = 5.6 × 10(-9)) in meta-analysis of all available data, and additional SNPs at the INHBC, LRP2, PLEKHA1, SLC3A2 and SLC7A6 genes meeting multiple-testing corrected significance for replication and overall P-values of 4.5 × 10(-4)-2.2 × 10(-7). Neither the novel PLEKHA1 nor FBXL20 associations, both further supported by association with eGFR among African Americans and with transcript abundance, would have been implicated by eGFR candidate gene approaches. LRP2, encoding the megalin receptor, was identified through connection with the previously known eGFR gene DAB2 and extends understanding of the megalin system in kidney function. These findings highlight integration of existing genome-wide association data with independent biological knowledge to uncover novel candidate eGFR associations, including candidates lacking known connections to kidney-specific pathways. The strategy may also be applicable to other clinical phenotypes, although more testing will be needed to assess its potential for discovery in general.

De Novo Assembly of the Pneumocystis jirovecii Genome from a Single Bronchoalveolar Lavage Fluid Specimen from a Patient.

ABSTRACT Pneumocystis jirovecii is a fungus that causes severe pneumonia in immunocompromised patients. However, its study is hindered by the lack of an in vitro culture method. We report here the genome of P. jirovecii that was obtained from a single bronchoalveolar lavage fluid specimen from a patient. The major challenge was the in silico sorting of the reads from a mixture representing the different organisms of the lung microbiome. This genome lacks virulence factors and most amino acid biosynthesis enzymes and presents reduced GC content and size. Together with epidemiological observations, these features suggest that P. jirovecii is an obligate parasite specialized in the colonization of human lungs, which causes disease only in immune-deficient individuals. This genome sequence will boost research on this deadly pathogen. IMPORTANCE Pneumocystis pneumonia is a major cause of mortality in patients with impaired immune systems. The availability of the P. jirovecii genome sequence allows new analyses to be performed which open avenues to solve critical issues for this deadly human disease. The most important ones are (i) identification of nutritional supplements for development of culture in vitro, which is still lacking 100 years after discovery of the pathogen; (ii) identification of new targets for development of new drugs, given the paucity of present treatments and emerging resistance; and (iii) identification of targets for development of vaccines.

Genome-wide association study of blood pressure extremes identifies variant near UMOD associated with hypertension.

Hypertension is a heritable and major contributor to the global burden of disease. The sum of rare and common genetic variants robustly identified so far explain only 1%-2% of the population variation in BP and hypertension. This suggests the existence of more undiscovered common variants. We conducted a genome-wide association study in 1,621 hypertensive cases and 1,699 controls and follow-up validation analyses in 19,845 cases and 16,541 controls using an extreme case-control design. We identified a locus on chromosome 16 in the 5' region of Uromodulin (UMOD; rs13333226, combined P value of 3.6×10(-11)). The minor G allele is associated with a lower risk of hypertension (OR [95%CI]: 0.87 [0.84-0.91]), reduced urinary uromodulin excretion, better renal function; and each copy of the G allele is associated with a 7.7% reduction in risk of CVD events after adjusting for age, sex, BMI, and smoking status (H.R. = 0.923, 95% CI 0.860-0.991; p = 0.027). In a subset of 13,446 individuals with estimated glomerular filtration rate (eGFR) measurements, we show that rs13333226 is independently associated with hypertension (unadjusted for eGFR: 0.89 [0.83-0.96], p = 0.004; after eGFR adjustment: 0.89 [0.83-0.96], p = 0.003). In clinical functional studies, we also consistently show the minor G allele is associated with lower urinary uromodulin excretion. The exclusive expression of uromodulin in the thick portion of the ascending limb of Henle suggests a putative role of this variant in hypertension through an effect on sodium homeostasis. The newly discovered UMOD locus for hypertension has the potential to give new insights into the role of uromodulin in BP regulation and to identify novel drugable targets for reducing cardiovascular risk.

Control of mitochondrial pH by uncoupling protein 4 in astrocytes promotes neuronal survival.

Brain activity is energetically costly and requires a steady and highly regulated flow of energy equivalents between neural cells. It is believed that a substantial share of cerebral glucose, the major source of energy of the brain, will preferentially be metabolized in astrocytes via aerobic glycolysis. The aim of this study was to evaluate whether uncoupling proteins (UCPs), located in the inner membrane of mitochondria, play a role in setting up the metabolic response pattern of astrocytes. UCPs are believed to mediate the transmembrane transfer of protons, resulting in the uncoupling of oxidative phosphorylation from ATP production. UCPs are therefore potentially important regulators of energy fluxes. The main UCP isoforms expressed in the brain are UCP2, UCP4, and UCP5. We examined in particular the role of UCP4 in neuron-astrocyte metabolic coupling and measured a range of functional metabolic parameters including mitochondrial electrical potential and pH, reactive oxygen species production, NAD/NADH ratio, ATP/ADP ratio, CO2 and lactate production, and oxygen consumption rate. In brief, we found that UCP4 regulates the intramitochondrial pH of astrocytes, which acidifies as a consequence of glutamate uptake, with the main consequence of reducing efficiency of mitochondrial ATP production. The diminished ATP production is effectively compensated by enhancement of glycolysis. This nonoxidative production of energy is not associated with deleterious H2O2 production. We show that astrocytes expressing more UCP4 produced more lactate, which is used as an energy source by neurons, and had the ability to enhance neuronal survival.

Mitochondrial DNA (mtDNA) A 3243G mutation associated with an annular perimacular retinal atrophy

BACKGROUND: A point mutation at the locus 3243 of the mitonchondrial DNA (mtDNA) is associated with either the MIDD syndrome (maternally inherited diabetes, deafness), the MELAS syndrome (myopathy, encephalitis, lactic acidosis, stroke) or cardiac, digestive, endocrine or exocrine dysfunctions. We report a peculiar maculopathy in two patients with an mtDNA 3243 mutation. HISTORY AND SIGNS: Case 1: A visually asymptomatic 40-year-old woman was examined for screening of diabetic retinopathy. Visual acuity was 10 / 10 in both eyes. Case 2: A 54-year-old woman with deafness and diabetes complained of visual loss. Visual acuity was 6 / 10 for the right eye and 0.5 / 10 for the left eye. Both patients exhibited a chorioretinal areolar atrophy. Case 1 was followed over 15 years and exhibited a slow progression of the maculopathy with moderate loss of visual acuity to 6 / 10 in both eyes, but marked handicap from the annular scotoma. THERAPY AND OUTCOME: None. CONCLUSION: Both patients presented a perimacular annular retinal atrophy. Patients harbouring mtDNA 3243 mutation should be examined for the presence of a maculopathy, even if they are asymptomatic. Conversely, the finding of such a geographic maculopathy should suggest the possibility of a point mutation at the locus 3243 of the mitochondrial DNA, especially in the presences of diabetes mellitus and/or deafness

Genome-wide association study to identify common variants associated with brachial circumference: a meta-analysis of 14 cohorts.

Brachial circumference (BC), also known as upper arm or mid arm circumference, can be used as an indicator of muscle mass and fat tissue, which are distributed differently in men and women. Analysis of anthropometric measures of peripheral fat distribution such as BC could help in understanding the complex pathophysiology behind overweight and obesity. The purpose of this study is to identify genetic variants associated with BC through a large-scale genome-wide association scan (GWAS) meta-analysis. We used fixed-effects meta-analysis to synthesise summary results across 14 GWAS discovery and 4 replication cohorts comprising overall 22,376 individuals (12,031 women and 10,345 men) of European ancestry. Individual analyses were carried out for men, women, and combined across sexes using linear regression and an additive genetic model: adjusted for age and adjusted for age and BMI. We prioritised signals for follow-up in two-stages. We did not detect any signals reaching genome-wide significance. The FTO rs9939609 SNP showed nominal evidence for association (p<0.05) in the age-adjusted strata for men and across both sexes. In this first GWAS meta-analysis for BC to date, we have not identified any genome-wide significant signals and do not observe robust association of previously established obesity loci with BC. Large-scale collaborations will be necessary to achieve higher power to detect loci underlying BC.

BIRO1, a cell-permeable BH3 peptide, promotes mitochondrial fragmentation and death of retinoblastoma cells.

Retinoblastoma is the most common pediatric intraocular neoplasm. While retinoblastoma development requires the inactivation of both alleles of the retinoblastoma tumor suppressor gene (RB1) in the developing retina, additional genomic changes are involved in tumor progression, which progressively lead to resistance of tumor cells to death. Therapeutics acting at very downstream levels of death signaling pathways should therefore be interesting in killing retinoblastoma cells. The BH3-only proteins promote apoptosis by modulating the interaction between the pro- and antiapoptotic members of the BCL2 protein family, and this effect can be recapitulated by the BH3 domains. This report analyzes the effect of various BH3 peptides, corresponding to different BH3-only proteins, on two retinoblastoma cell lines, Y79 and WERI-Rb, as well as on the photoreceptor cell line 661W. The BH3 peptide BIRO1, derived from the BCL2L11 death domain, was very effective in promoting Y79 and WERI-Rb cell death without affecting the 661W photoreceptor cells. This cell death was efficient even in absence of BAX and was shown to be caspase independent. While ROS production or AIF release was not detected from mitochondria of treated cells, BIRO1 initiated mitochondria fragmentation in a short period of time following treatment. IMPLICATIONS: The BIRO1 peptide is highly effective at killing retinoblastoma cells and has potential as a peptidomimetic.

Genome-wide association analysis of blood-pressure traits in African-ancestry individuals reveals common associated genes in African and non-African populations.

High blood pressure (BP) is more prevalent and contributes to more severe manifestations of cardiovascular disease (CVD) in African Americans than in any other United States ethnic group. Several small African-ancestry (AA) BP genome-wide association studies (GWASs) have been published, but their findings have failed to replicate to date. We report on a large AA BP GWAS meta-analysis that includes 29,378 individuals from 19 discovery cohorts and subsequent replication in additional samples of AA (n = 10,386), European ancestry (EA) (n = 69,395), and East Asian ancestry (n = 19,601). Five loci (EVX1-HOXA, ULK4, RSPO3, PLEKHG1, and SOX6) reached genome-wide significance (p < 1.0 × 10(-8)) for either systolic or diastolic BP in a transethnic meta-analysis after correction for multiple testing. Three of these BP loci (EVX1-HOXA, RSPO3, and PLEKHG1) lack previous associations with BP. We also identified one independent signal in a known BP locus (SOX6) and provide evidence for fine mapping in four additional validated BP loci. We also demonstrate that validated EA BP GWAS loci, considered jointly, show significant effects in AA samples. Consequently, these findings suggest that BP loci might have universal effects across studied populations, demonstrating that multiethnic samples are an essential component in identifying, fine mapping, and understanding their trait variability.