ALBERT

Description: Frontonasal dysplasia (FND) refers to a class of midline facial malformations caused by abnormal development of the facial primordia. The term encompasses a spectrum of severities but characteristic features include combinations of ocular hypertelorism, malformations of the nose and forehead and clefting of the facial midline. Several recent studies have drawn attention to the importance of Alx homeobox transcription factors during craniofacial development. Most notably, loss of Alx1 has devastating consequences resulting in severe orofacial clefting and extreme microphthalmia. In contrast, mutations of Alx3 or Alx4 cause milder forms of FND. Whilst Alx1 , Alx3 and Alx4 are all known to be expressed in the facial mesenchyme of vertebrate embryos, little is known about the function of these proteins during development. Here, we report the establishment of a zebrafish model of Alx -related FND. Morpholino knock-down of zebrafish alx1 expression causes a profound craniofacial phenotype including loss of the facial cartilages and defective ocular development. We demonstrate for the first time that Alx1 plays a crucial role in regulating the migration of cranial neural crest (CNC) cells into the frontonasal primordia. Abnormal neural crest migration is coincident with aberrant expression of foxd3 and sox10 , two genes previously suggested to play key roles during neural crest development, including migration, differentiation and the maintenance of progenitor cells. This novel function is specific to Alx1, and likely explains the marked clinical severity of Alx1 mutation within the spectrum of Alx -related FND.

Description: Activating somatic and germline mutations of closely related RAS genes (H, K, N) have been found in various types of cancer and in patients with developmental disorders, respectively. The involvement of the RAS signalling pathways in developmental disorders has recently emerged as one of the most important drivers in RAS research. In the present study, we investigated the biochemical and cell biological properties of two novel missense KRAS mutations (Y71H and K147E). Both mutations affect residues that are highly conserved within the RAS family. KRAS Y71H showed no clear differences to KRAS wt , except for an increased binding affinity for its major effector, the RAF1 kinase. Consistent with this finding, even though we detected similar levels of active KRAS Y71H when compared with wild-type protein, we observed an increased activation of MEK1/2, irrespective of the stimulation conditions. In contrast, KRAS K147E exhibited a tremendous increase in nucleotide dissociation generating a self-activating RAS protein that can act independently of upstream signals. As a consequence, levels of active KRAS K147E were strongly increased regardless of serum stimulation and similar to the oncogenic KRAS G12V . In spite of this, KRAS K147E downstream signalling did not reach the level triggered by oncogenic KRAS G12V , especially because KRAS K147E was downregulated by RASGAP and moreover exhibited a 2-fold lower affinity for RAF kinase. Here, our findings clearly emphasize that individual RAS mutations, despite being associated with comparable phenotypes of developmental disorders in patients, can cause remarkably diverse biochemical effects with a common outcome, namely a rather moderate gain-of-function.

Description: Spinocerebellar ataxia type 2 (SCA2) is an autosomal dominantly inherited disorder, which is caused by a pathological expansion of a polyglutamine (polyQ) tract in the coding region of the ATXN2 gene. Like other ataxias, SCA2 most overtly affects Purkinje cells (PCs) in the cerebellum. Using a transgenic mouse model expressing a full-length ATXN2 Q127 -complementary DNA under control of the Pcp2 promoter (a PC-specific promoter), we examined the time course of behavioral, morphologic, biochemical and physiological changes with particular attention to PC firing in the cerebellar slice. Although motor performance began to deteriorate at 8 weeks of age, reductions in PC number were not seen until after 12 weeks. Decreases in the PC firing frequency first showed at 6 weeks and paralleled deterioration of motor performance with progression of disease. Transcription changes in several PC-specific genes such as Calb1 and Pcp2 mirrored the time course of changes in PC physiology with calbindin-28 K changes showing the first small, but significant decreases at 4 weeks. These results emphasize that in this model of SCA2, physiological and behavioral phenotypes precede morphological changes by several weeks and provide a rationale for future studies examining the effects of restoration of firing frequency on motor function and prevention of future loss of PCs.

Description: Balancing selection has maintained human leukocyte antigen (HLA) allele diversity, but it is unclear whether this selection is symmetric (all heterozygotes are comparable and all homozygotes are comparable in terms of fitness) or asymmetric (distinct heterozygote genotypes display greater fitness than others). We tested the hypothesis that HLA is under asymmetric balancing selection in populations by estimating allelic branch lengths from genetic sequence data encoding peptide-binding domains. Significant deviations indicated changes in the ratio of terminal to internal branch lengths. Such deviations could arise even if no individual alleles present a strikingly altered branch length (e.g. if there is an overall distortion, with all or many terminal branches being longer than expected). DQ and DP loci were also analyzed as haplotypes. Using allele frequencies for 419 distinct populations in 10 geographical regions, we examined population differentiation in alleles within and between regions, and the relationship between allelic branch length and frequency. The strongest evidence for asymmetrical balancing selection was observed for HLA-DRB1 , HLA-B and HLA-DPA1 , with significant deviation ( P ≤ 1.1 x 10 –4 ) in about half of the populations. There were significant results at all loci except HLA-DQB1 / DQA1 . We observed moderate genetic variation within and between geographic regions, similar to the rest of the genome. Branch length was not correlated with allele frequency. In conclusion, sequence data suggest that balancing selection in HLA is asymmetric (some heterozygotes enjoy greater fitness than others). Because HLA polymorphism is crucial for pathogen resistance, this may manifest as a frequency-dependent selection with fluctuation in the fitness of specific heterozygotes over time.

Description: Birt–Hogg–Dubé syndrome (BHD) is a human cancer disorder caused by mutations in the tumor suppressor gene Folliculin ( FLCN ) with unknown biological functions. Here, we show that the Drosophila homolog of FLCN, dFLCN (a.k.a. dBHD ) localizes to the nucleolus and physically interacts with the 19S proteasomal ATPase, Rpt4, a nucleolar resident and known regulator of rRNA transcription. Downregulation of dFLCN resulted in an increase in nucleolar volume and upregulation of rRNA synthesis, whereas dFLCN overexpression reduced rRNA transcription and counteracted the effects of Rpt4 on rRNA production by preventing the association of Rpt4 with the rDNA locus. We further show that human FLCN exhibited evolutionarily conserved function and that Rpt4 knockdown inhibits the growth of FLCN-deficient human renal cancer cells in mouse xenografts. Our study suggests that FLCN functions as a tumor suppressor by negatively regulating rRNA synthesis.

Description: KitL, via its receptor cKit, supports primordial germ cell (PGC) growth, survival, migration and reprogramming to pluripotent embryonic germ cells (EGCs). However, the signaling downstream of KitL and its regulation in PGCs remain unclear. A constitutively activating mutation, cKit V558 , causes gain-of-function phenotypes in mast cells and intestines, and gastrointestinal stromal tumors (GISTs) when heterozygous. Unexpectedly, we find that PGC growth is not significantly affected in cKit V558 heterozygotes, whereas in homozygotes, increased apoptosis and inefficient migration lead to the depletion of PGCs. Through genetic studies, we reveal that this oncogenic cKit allele exhibits loss-of-function behavior in PGCs distinct from that in GIST development. Examination of downstream signaling in GISTs from cKit V558/+ mice confirmed hyperphosphorylation of AKT and ERK, but both remain unperturbed in cKit V558/+ PGCs and EGCs. In contrast, we find reduced activation of ERK1/2 and JNK1 in cKit V558 homozygous PGCs and EGCs. Inhibiting JNK, though not ERK1/2, increased apoptosis of wild-type PGCs, but did not further affect the already elevated apoptosis of cKit V558 / V558 PGCs. These results demonstrate a cell-context-dependent response to the cKit V558 mutation. We propose that AKT overload protection and JNK-mediated survival comprise PGC-specific mechanisms for regulating cKit signaling.

Description: TDP-43 is an evolutionarily conserved RNA-binding protein currently under intense investigation for its involvement in the molecular pathogenesis of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). TDP-43 is normally localized in the nucleus, but translocated to the cytoplasm in diseased neurons. The endogenous functions of TDP-43 in the nervous system remain poorly understood. Here, we show that the loss of Drosophila TDP-43 (dTDP-43) results in an increased production of sensory bristles and sensory organ precursor (SOP) cells on the notum of some but not all flies. The location of ectopic SOPs varies among mutant flies. The penetrance of this novel phenotype is dependent on the gender and sensitive to environmental influences. A similar SOP phenotype was also observed on the wing and in the embryos. Overexpression of dTDP-43 causes both loss and ectopic production of SOPs. Ectopic expression of ALS-associated mutant human TDP-43 (hTDP-43 M337V and hTDP-43 Q331K ) produces a less severe SOP phenotype than hTDP-43 WT , indicating a partial loss of function of mutant hTDP-43. In dTDP-43 mutants, miR-9a expression is significantly reduced. Genetic interaction studies further support the notion that dTDP-43 acts through miR-9a to control the precision of SOP specification. These findings reveal a novel role for endogenous TDP-43 in neuronal specification and suggest that the FTD/ALS-associated RNA-binding protein TDP-43 functions to ensure the robustness of genetic control programs.

Description: Mutations in leucine-rich repeat kinase 2 (LRRK2) are the most frequent known cause of late-onset Parkinson's disease (PD). To explore the therapeutic potential of small molecules targeting the LRRK2 kinase domain, we characterized two LRRK2 kinase inhibitors, TTT-3002 and LRRK2-IN1, for their effects against LRRK2 activity in vitro and in Caenorhabditis elegans models of LRRK2-linked neurodegeneration. TTT-3002 and LRRK2-IN1 potently inhibited in vitro kinase activity of LRRK2 wild-type and mutant proteins, attenuated phosphorylation of cellular LRRK2 and rescued neurotoxicity of mutant LRRK2 in transfected cells. To establish whether LRRK2 kinase inhibitors can mitigate pathogenesis caused by different mutations including G2019S and R1441C located within and outside of the LRRK2 kinase domain, respectively, we evaluated effects of TTT-3002 and LRRK2-IN1 against R1441C- and G2019S-induced neurodegeneration in C. elegans models. TTT-3002 and LRRK2-IN1 rescued the behavioral deficit characteristic of dopaminergic impairment in transgenic C. elegans expressing human R1441C- and G2019S-LRRK2. The inhibitors displayed nanomolar to low micromolar rescue potency when administered either pre-symptomatically or post-symptomatically, indicating both prevention and reversal of the dopaminergic deficit. The same treatments also led to long-lasting prevention and rescue of neurodegeneration. In contrast, TTT-3002 and LRRK2-IN1 were ineffective against the neurodegenerative phenotype in transgenic worms carrying the inhibitor-resistant A2016T mutation of LRRK2, suggesting that they elicit neuroprotective effects in vivo by targeting LRRK2 specifically. Our findings indicate that the LRRK2 kinase activity is critical for neurodegeneration caused by R1441C and G2019S mutations, suggesting that kinase inhibition of LRRK2 may represent a promising therapeutic strategy for PD.

Description: Mitochondrial DNA (mtDNA) mutations leading to the disruption of respiratory complex I (CI) have been shown to exhibit anti-tumorigenic effects, at variance with those impairing only the function but not the assembly of the complex, which appear to contribute positively to cancer development. Owing to the challenges in the analysis of the multi-copy mitochondrial genome, it is yet to be determined whether tumour-associated mtDNA lesions occur as somatic modifying factors or as germ-line predisposing elements. Here we investigated the whole mitochondrial genome sequence of 20 pituitary adenomas with oncocytic phenotype and identified pathogenic and/or novel mtDNA mutations in 60% of the cases. Using highly sensitive techniques, namely fluorescent PCR and allele-specific locked nucleic acid quantitative PCR, we identified the most likely somatic nature of these mutations in our sample set, since none of the mutations was detected in the corresponding blood tissue of the patients analysed. Furthermore, we have subjected a series of 48 pituitary adenomas to a high-resolution array comparative genomic hybridization analysis, which revealed that CI disruptive mutations, and the oncocytic phenotype, significantly correlate with low number of chromosomal aberrations in the nuclear genome. We conclude that CI disruptive mutations in pituitary adenomas are somatic modifiers of tumorigenesis most likely contributing not only to the development of oncocytic change, but also to a less aggressive tumour phenotype, as indicated by a stable karyotype.

Description: Functional loss of SMN1 causes proximal spinal muscular atrophy (SMA), the most common genetic condition accounting for infant lethality. Hence, the hypomorphic copy gene SMN2 is the only resource of functional SMN protein in SMA patients and influences SMA severity in a dose-dependent manner. Consequently, current therapeutic approaches focus on SMN2 . Histone deacetylase inhibitors (HDACi), such as the short chain fatty acid VPA (valproic acid), ameliorate the SMA phenotype by activating the SMN2 expression. By analyzing blood SMN2 expression in 16 VPA-treated SMA patients, about one-third of individuals were identified as positive responders presenting increased SMN2 transcript levels. In 66% of enrolled patients, a concordant response was detected in the respective fibroblasts. Most importantly, by taking the detour of reprograming SMA patients' fibroblasts, we showed that the VPA response was maintained even in GABAergic neurons derived from induced pluripotent stem cells (iPS) cells. Differential expression microarray analysis revealed a complete lack of response to VPA in non-responders, which was associated with an increased expression of the fatty acid translocase CD36. The pivotal role of CD36 as the cause of non-responsiveness was proven in various in vitro approaches. Most importantly, knockdown of CD36 in SMA fibroblasts converted non- into pos-responders. In summary, the concordant response from blood to the central nervous system (CNS) to VPA may allow selection of pos-responders prior to therapy. Increased CD36 expression accounts for VPA non-responsiveness. These findings may be essential not only for SMA but also for other diseases such as epilepsy or migraine frequently treated with VPA.

Description: Rett syndrome (RTT) is a neurodevelopmental disorder caused primarily by mutations of the X-linked MECP2 gene. Although the loss of MeCP2 function affects many neural systems, impairments of catecholaminergic function have been hypothesized to underlie several of the cardinal behavioral deficits of RTT patients and Mecp2-deficient mice. Although recent Mecp2 reactivation studies indicate that RTT may be a reversible condition, it remains unclear whether specifically preserving Mecp2 function within a specific system will be sufficient to convey beneficial effects. Here, we test whether the selective preservation of Mecp2 within catecholaminergic cells will improve the phenotype of Mecp2-deficient mice. Our results show that this targeted preservation of Mecp2 significantly improves the lifespan, phenotypic severity and cortical epileptiform discharge activity of both male and female Mecp2-deficient mice. Further, we found that the catecholaminergic preservation of Mecp2 also improves the ambulatory rate, rearing activity, motor coordination, anxiety and nest-building performances of Mecp2-deficient mice of each gender. Interestingly, our results also revealed a gender-specific improvement, as specific cortical and hippocampal electroencephalographic abnormalities were significantly improved in male, but not female, rescue mice. Collectively, these results support the role of the catecholaminergic system in the pathogenesis of RTT and provide proof-of-principle that restoring MeCP2 function within this specific system could represent a treatment strategy for RTT.

Description: Mutations in COL4A1 have been identified in families with hereditary small vessel disease of the brain presumably due to a dominant-negative mechanism. Here, we report on two novel mutations in COL4A1 in two families with porencephaly, intracerebral hemorrhage and severe white matter disease caused by haploinsufficiency. Two families with various clinical presentations of cerebral microangiopathy and autosomal dominant inheritance were examined. Clinical, neuroradiological and genetic investigations were performed. Electron microscopy of the skin was also performed. In one of the families, sequence analysis revealed a one base deletion, c.2085del, leading to a frameshift and a premature stopcodon, p.(Gly696fs). In the other family, a splice site mutation was identified, c.2194-1G〉A, which most likely leads to skipping of an exon with a frameshift and premature termination as a result. In fibroblasts of affected individuals from both the families, nonsense-mediated decay (NMD) of the mutant COL4A1 messenger RNAs (mRNAs) and a clear reduction of COL4A1 protein expression were demonstrated, indicating haploinsufficiency of COL4A1. Moreover, thickening of the capillary basement membrane in the skin was documented, similar to reports in patients with COL4A1 missense mutations. These findings suggest haploinsufficiency, a different mechanism from the commonly assumed dominant-negative effect, for COL4A1 mutations as a cause of (antenatal) intracerebral hemorrhage and white matter disease.

Description: Angiogenesis, a formation of neovessels, is regulated by the local balance between angiogenesis stimulators and inhibitors. A number of such endogenous regulators of angiogenesis have been found in the body. Recently, vasohibin-1 (VASH1) was isolated as a negative feedback regulator of angiogenesis produced by endothelial cells (ECs) and subsequently vasohibin-2 (VASH2) as a homologue of VASH1. It was then explored that VASH1 is expressed in ECs to terminate angiogenesis, whereas VASH2 is expressed in cells other than ECs to promote angiogenesis in the mouse model of angiogenesis. This review will focus on the vasohibin family members, which are novel regulators of angiogenesis.

Description: Tertiary dentin is deposited inside teeth after various stimuli and serves as a major defensive wall to preserve pulp cells. However, the molecular mechanisms of the activation of quiescent odontoblasts, immature pulp cells and tertiary dentin formation are still unclear. Therefore, we performed a comprehensive gene expression analysis of pulp cells after cavity preparation of 9-week-old rat molars to clarify the critical molecules in tertiary dentinogenesis. As a result, mRNA expression of various molecules was up- or down-regulated. Notably, several members of the matrix metalloprotease family and their endogenous inhibitors were up-regulated after cavity preparation. In situ hybridization showed that tissue inhibitor of metalloprotease 1 ( Timp1 ) was widely and continuously distributed in the pulp beneath the cavity in vivo . We also observed accumulation of β-catenin in the pulp cells beneath the cavity by fluorescence immunohistochemistry. Furthermore, Timp1 transcription was repressed by a dominant-negative TCF4 in immature undifferentiated mesenchymal cells, but not altered in mature odontoblast-like cells. These results indicate that cavity preparation may activate the Wnt/β-catenin pathway and the Wnt/β-catenin pathway and Timp1 may be correlatively involved in pulp repair. Timp1 might play crucial roles in reactivation of immature pulp cells for tertiary dentinogenesis.

Description: The activity of biological molecules is often affected by their phosphorylation state. Regulatory phosphorylation operates as a binary switch and is usually controlled by counteracting kinases and phosphatases. However, phosphatidylinositol (PtdIns) has three phosphorylation sites on its inositol ring. The phosphorylation status of PtdIns is controlled by multiple kinases and phosphatases with distinct substrate specificities, serving as a ‘lipid code’ or ‘phosphoinositide code’. Class I phosphoinositide 3-kinase (PI3K) converts PtdIns(4,5)P 2 to PtdIns(3,4,5)P 3 , which plays a pivotal role in signals controlling glucose uptake, cytoskeletal reorganization, cell proliferation and apoptosis. PI3K is pro-oncogenic, whereas phosphoinositide phosphatases that degrade PtdIns(3,4,5)P 3 are not always anti-oncogenic. Recent studies have revealed the unique characteristics of phosphoinositide 5-phosphatases.

Description: Vascular endothelial growth factors (VEGFs) belong to the platelet-derived growth factor supergene family, and they play central roles in the regulation of angiogenesis and lymphangiogenesis. VEGF-A, the major factor for angiogenesis, binds to two tyrosine kinase (TK) receptors, VEGFR-1 (Flt-1) and VEGFR-2 (KDR/Flk-1), and regulates endothelial cell proliferation, migration, vascular permeability, secretion and other endothelial functions. VEGFR-2 exhibits a strong TK activity towards pro-angiogenic signals, whereas the soluble VEGFR-1 (sFlt-1) functions as an endogenous VEGF inhibitor. sFlt-1 is abnormally overexpressed in the placenta of preeclampsia patients, resulting in the major symptoms of the disease due to abnormal trapping of VEGFs. The VEGF-VEGFR system is crucial for tumour angiogenesis, and anti-VEGF-VEGFR molecules are now widely used in the clinical field to treat cancer patients. The efficacy of these molecules in prolonging the overall survival of patients has been established; however, some cancers do not respond well and reduced tumour sensitivity to anti-VEGF signals may occur after long-term treatment. The molecular basis of tumour refractoriness should be determined to improve anti-angiogenic therapy.

Description: The actual levels of steroid hormones in organs are vital for endocrine, reproductive and neuronal health and disorders. We developed an accurate method to determine the levels of steroid hormones and steroid conjugates in various organs by an efficient preparation using a solid-phase-extraction cartridge. Each steroid was identified by the precursor ion spectra using liquid chromatography–electrospray ionization time-of-flight mass spectrometry, and the respective steroids were quantitatively analysed in the selected reaction monitoring mode by liquid chromatograph-mass spectrometry/mass spectrometry (LC-MS/MS). The data showed that significant levels of testosterone, corticosterone and precursors of both hormones were detected in all organs except liver. The glucuronide conjugates of steroid hormones and the precursors were detected in all organs except liver, but sulfate conjugates of these steroids were observed only in the target organs of the hormones and kidney. Interestingly, these steroids and the conjugates were not observed in the liver except pregnenolone. In conclusion, an accurate determination of tissue steroids was developed using LC-MS analysis. Biosynthesis of steroid hormones from the precursors was estimated even in the target organs, and the delivery of these steroid conjugates was also suggested via the circulation without any significant hepatic participation.

Description: Papain-like cysteine protease activity that shows a unique transient expression profile in cotyledons of daikon radish during germination was detected. The enzyme showed a distinct elution pattern on DEAE-cellulose compared with cathepsin B-like and Responsive to dessication-21 cysteine protease. Although this activity was not detected in seed prior to imbibition, the activity increased markedly and reached a maximum at 2 days after imbibition and then decreased rapidly and completely disappeared after 5 days. Using cystatin-Sepharose, the 26 kDa cysteine protease (DRCP26) was isolated from cotyledons at 2 days after imbibition. The deduced amino acid sequence from the cDNA nucleotide sequence indicated that DRCP26 is an orthologue of Arabidopsis unidentified protein, germination-specific cysteine protease-1, belonging to the C1 family of cysteine protease predicted from genetic information. In an effort to characterize the enzymatic properties of DRCP26, the enzyme was purified to homogeneity from cotyledons at 48 h after imbibition. The best synthetic substrate for the enzyme was carbobenzoxy-Phe-Arg-4-methylcoumaryl-7-amide. All model peptides were digested to small peptides by the enzyme, suggesting that DRCP26 possesses broad cleavage specificity. These results indicated that DRCP26 plays a role in the mobilization of storage proteins in the early phase of seed germination.

Description: Proteins in the superfamily of voltage-gated ion channels mediate behavior across the tree of life. These proteins regulate the movement of ions across cell membranes by opening and closing a central pore that controls ion flow. The best-known members of this superfamily are the voltage-gated potassium, calcium (Ca v ), and sodium (Na v ) channels, which underlie impulse conduction in nerve and muscle. Not all members of this family are opened by changes in voltage, however. NALCN (NA + leak channel nonselective) channels, which encode a voltage-insensitive "sodium leak" channel, have garnered a growing interest. This study examines the phylogenetic relationship among Na v /Ca v voltage-gated and voltage-insensitive channels in the eukaryotic group Opisthokonta, which includes animals, fungi, and their unicellular relatives. We show that NALCN channels diverged from voltage-gated channels before the divergence of fungi and animals and that the closest relatives of NALCN channels are fungal calcium channels, which they functionally resemble.

Description: Newer parts of sex chromosomes, neo-sex chromosomes, offer unique possibilities for studying gene degeneration and sequence evolution in response to loss of recombination and population size decrease. We have recently described a neo-sex chromosome system in Sylvioidea passerines that has resulted from a fusion between the first half (10 Mb) of chromosome 4a and the ancestral sex chromosomes. In this study, we report the results of molecular analyses of neo-Z and neo-W gametologs and intronic parts of neo-Z and autosomal genes on the second half of chromosome 4a in three species within different Sylvioidea lineages (Acrocephalidea, Timaliidae, and Alaudidae). In line with hypotheses of neo-sex chromosome evolution, we observe 1) lower genetic diversity of neo-Z genes compared with autosomal genes, 2) moderate synonymous and weak nonsynonymous sequence divergence between neo-Z and neo-W gametologs, and 3) lower GC content on neo-W than neo-Z gametologs. Phylogenetic reconstruction of eight neo-Z and neo-W gametologs suggests that recombination continued after the split of Alaudidae from the rest of the Sylvioidea lineages (i.e., after ~42.2 Ma) and with some exceptions also after the split of Acrocephalidea and Timaliidae (i.e., after ~39.4 Ma). The Sylvioidea neo-sex chromosome shares classical evolutionary features with the ancestral sex chromosomes but, as expected from its more recent origin, shows weaker divergence between gametologs.

Description: In bacteria, physiological change may be effected by a single gene acquisition, producing ecological differentiation without genetic isolation. Natural selection acting on such differences can reduce the frequency of genotypes that arise from recombination at these loci. However, gene acquisition can only account for recombination interference in the fraction of the genome that is tightly linked to the integration site. To identify additional loci that contribute to adaptive differences, we examined orthologous genes in species of Enterobacteriaceae to identify significant differences in the degree of codon selection. Significance was assessed using the Adaptive Codon Enrichment metric, which accounts for the variation in codon usage bias that is expected to arise from mutation and drift; large differences in codon usage bias were identified in more genes than would be expected to arise from stochastic processes alone. Genes in the same operon showed parallel differences in codon usage bias, suggesting that changes in the overall levels of gene expression led to changes in the degree of adaptive codon usage. Most significant differences between orthologous operons were found among those involved with specific environmental adaptations, whereas "housekeeping" genes rarely showed significant changes. When considered together, the loci experiencing significant changes in codon selection outnumber potentially adaptive gene acquisition events. The identity of genes under strong codon selection seems to be influenced by the habitat from which the bacteria were isolated. We propose a two-stage model for how adaptation to different selective regimes can drive bacterial speciation. Initially, gene acquisitions catalyze rapid ecological differentiation, which modifies the utilization of genes, thereby changing the strength of codon selection on them. Alleles develop fitness variation by substitution, producing recombination interference at these loci in addition to those flanking acquired genes, allowing sequences to diverge across the entire genome and establishing genetic isolation (i.e., protection from frequent homologous recombination).

Description: Gene expression levels correlate with multiple aspects of gene sequence and gene structure in phylogenetically diverse taxa, suggesting an important role of gene expression levels in the evolution of protein-coding genes. Here we present results of a genome-wide study of the influence of gene expression on synonymous codon usage, amino acid composition, and gene structure in the red flour beetle, Tribolium castaneum . Consistent with the action of translational selection, we find that synonymous codon usage bias increases with gene expression. However, the correspondence between tRNA gene copy number and optimal codons is weak. At the amino acid level, translational selection is suggested by the positive correlation between tRNA gene numbers and amino acid usage, which is stronger for highly expressed genes. In addition, there is a clear trend for increased use of metabolically cheaper, less complex amino acids as gene expression increases. tRNA gene numbers also correlate negatively with amino acid size/complexity (S/C) score indicating the coupling between translational selection and selection to minimize the use of large/complex amino acids. Interestingly, the analysis of 10 additional genomes suggests that the correlation between tRNA gene numbers and amino acid S/C score is widespread and might be explained by selection against negative consequences of protein misfolding. At the level of gene structure, three major trends are detected: 1) complete coding region length increases across low and intermediate expression levels but decreases in highly expressed genes; 2) the average intron size shows the opposite trend, first decreasing with expression, followed by a slight increase in highly expressed genes; and 3) intron density remains nearly constant across all expression levels. These changes in gene architecture are only in partial agreement with selection favoring reduced cost of biosynthesis.

Description: New protein-coding genes can originate either through modification of existing genes or de novo. Recently, the importance of de novo origination has been recognized in eukaryotes, although eukaryotic genes originated de novo are relatively rare and difficult to identify. In contrast, viruses contain many de novo genes , namely those in which an existing gene has been "overprinted" by a new open reading frame, a process that generates a new protein-coding gene overlapping the ancestral gene. We analyzed the evolution of 12 experimentally validated viral genes that originated de novo and estimated their relative ages. We found that young de novo genes have a different codon usage from the rest of the genome. They evolve rapidly and are under positive or weak purifying selection. Thus, young de novo genes might have strain-specific functions, or no function, and would be difficult to detect using current genome annotation methods that rely on the sequence signature of purifying selection. In contrast to young de novo genes, older de novo genes have a codon usage that is similar to the rest of the genome. They evolve slowly and are under stronger purifying selection. Some of the oldest de novo genes evolve under stronger selection pressure than the ancestral gene they overlap, suggesting an evolutionary tug of war between the ancestral and the de novo gene.

Description: Here, we present a study of the molecular evolution of the pheromone receptor genes ( pre-1 and pre-2 ) in Neurospora taxa with different mating systems. We focus on comparisons between heterothallic and homothallic taxa, reproducing sexually by outcrossing and by intrahaploid selfing, respectively. Our general aim was to use a phylogenetic framework to investigate whether the evolutionary trajectory of the pheromone and receptor genes in Neurospora differs between heterothallic and homothallic taxa, and among the homothallic lineages/clades previously indicated to represent independent switches from heterothallism to homothallism in the evolutionary history of the genus. We complemented molecular evolution analyses with an expression study of the pre genes and their upstream regulators, the mating-type ( mat ) genes, in homothallic taxa. Our analyses suggest that the pheromone receptor gene pre-1 is functionally conserved in both heterothallic and homothallic taxa. Moreover, we found evidence of positive selection for a small fraction of codons in the cytoplasmic signal-transducing C-terminal region of the protein PRE-1. Distribution of positively selected codons differs between heterothallic and homothallic groups, suggesting functional divergence associated with mating system. The gene pre-2 was shown to evolve under high selective constraints, with no strong evidence for positive selection. Although our data suggest that both pre-1 and pre-2 are overall functional in homothallic taxa, individual taxa display frame-shift mutations causing premature stop codons, which might indicate loss of function. Transcriptional patterns of pre and mat genes in six homothallic taxa, selected to represent six different switches from heterothallism to homothallism, do not support a universal pattern of regulation of these genes during reproductive tissue development. Taken together, our analyses suggest that the pheromone receptor genes pre-1 and pre-2 are in general functional in homothallic Neurospora taxa, in contrast with the situation for the mat genes that are generally degenerate in these taxa.

Description: Comparative genome biology has unveiled the polyploid origin of all angiosperms and the role of recurrent polyploidization in the amplification of gene families and the structuring of genomes. Which species share certain ancient polyploidy events, and which do not, is ill defined because of the limited number of sequenced genomes and transcriptomes and their uneven phylogenetic distribution. Previously, it has been suggested that most, but probably not all, of the eudicots have shared an ancient hexaploidy event, referred to as the gamma triplication. In this study, detailed phylogenies of subfamilies of MADS-box genes suggest that the gamma triplication has occurred before the divergence of Gunnerales but after the divergence of Buxales and Trochodendrales. Large-scale phylogenetic and K S - based approaches on the inflorescence transcriptomes of Gunnera manicata (Gunnerales) and Pachysandra terminalis (Buxales) provide further support for this placement, enabling us to position the gamma triplication in the stem lineage of the core eudicots. This triplication likely initiated the functional diversification of key regulators of reproductive development in the core eudicots, comprising 75% of flowering plants. Although it is possible that the gamma event triggered early core eudicot diversification, our dating estimates suggest that the event occurred early in the stem lineage, well before the rapid speciation of the earliest core eudicot lineages. The evolutionary significance of this paleopolyploidy event may thus rather lie in establishing a species lineage that was resilient to extinction, but with the genomic potential for later diversification. We consider that the traits generated from this potential characterize extant core eudicots both chemically and morphologically.

Description: Genome reduction in obligately intracellular bacteria is one of the most well-established patterns in the field of molecular evolution. In the extreme, many sap-feeding insects harbor nutritional symbionts with genomes that are so reduced that it is not clear how they perform basic cellular functions. For example, the primary symbiont of psyllids ( Carsonella ) maintains one of the smallest and most AT-rich bacterial genomes ever identified and has surprisingly lost many genes that are thought to be essential for its role in provisioning its host with amino acids. However, our understanding of this extreme case of genome reduction is limited, as genomic data for Carsonella are available from only a single host species, and little is known about the functional role of "secondary" bacterial symbionts in psyllids. To address these limitations, we analyzed complete Carsonella genomes from pairs of congeneric hosts in three divergent genera within the Psyllidae ( Ctenarytaina , Heteropsylla , and Pachypsylla ) as well as complete secondary symbiont genomes from two of these host species ( Ctenarytaina eucalypti and Heteropsylla cubana ). Although the Carsonella genomes are generally conserved in size, structure, and GC content and exhibit genome-wide signatures of purifying selection, we found that gene loss has remained active since the divergence of the host species and had a particularly large impact on the amino acid biosynthesis pathways that define the symbiotic role of Carsonella . In some cases, the presence of additional bacterial symbionts may compensate for gene loss in Carsonella , as functional gene content indicates a high degree of metabolic complementarity between co-occurring symbionts. The genomes of the secondary symbionts also show signatures of long-term evolution as vertically transmitted, intracellular bacteria, including more extensive genome reduction than typically observed in facultative symbionts. Therefore, a history of co-evolution with secondary bacterial symbionts can partially explain the ongoing genome reduction in Carsonella . However, the absence of these secondary symbionts in other host lineages indicates that the relationships are dynamic and that other mechanisms, such as changes in host diet or functional coordination with the host genome, must also be at play.

Description: Gene duplication is a major driver of organismal adaptation and evolution and plays an important role in multiple human diseases. Whole-genome analyses have shown similar and high rates of gene duplication across a variety of eukaryotic species. Most of these studies, however, did not address the possible impact of interlocus gene conversion (IGC) on the evolution of gene duplicates. Because IGC homogenizes pairs of duplicates, widespread conversion would cause gene duplication events that happened long ago to appear more recent, resulting in artificially high estimates of duplication rates. Although the majority of genome-wide studies (including in the budding yeast Saccharomyces cerevisiae [Scer]) point to levels of IGC between paralogs ranging from 2% to 18%, Gao and Innan (Gao LZ, Innan H. 2004. Very low gene duplication rate in the yeast genome. Science 306:1367–1370.) found that gene conversion in yeast affected 〉80% of paralog pairs. If conversion rates really are this high, it would imply that the rate of gene duplication in eukaryotes is much lower than previously reported. In this work, we apply four different methodologies—including one approach that closely mirrors Gao and Innan’s method—to estimate the level of IGC in Scer. Our analyses point to a maximum conversion level of 13% between paralogs in this species, in close agreement with most estimates of IGC in eukaryotes. We also show that the exceedingly high levels of conversion found previously derive from application of an accurate method to an inappropriate data set. In conclusion, our work provides the most striking evidence to date supporting the reduced incidence of IGC among Scer paralogs and sets up a framework for future analyses in other eukaryotes.

Description: Cyclin B3 evolution has the unique peculiarity of an abrupt 3-fold increase of the protein size in the mammalian lineage due to the extension of a single exon. We have analyzed the evolution of the gene to define the modalities of this event and the possible consequences on the function of the protein. Database searches can trace the appearance of the gene to the origin of metazoans. Most introns were already present in early metazoans, and the intron–exon structure as well as the protein size were fairly conserved in invertebrates and nonmammalian vertebrates. Although intron gains are considered as rare events, we identified two cases, one at the prochordate–chordate transition and one in murids, resulting from different mechanisms. At the emergence of mammals, the gene was relocated from chromosome 6 of platypus to the X chromosome in marsupials, but the exon extension occurred only in placental mammals. A repetitive structure of 18 amino acids, of uncertain origin, is detectable in the 3,000-nt mammalian exon-encoded sequence, suggesting an extension by multiple internal duplications, some of which are still detectable in the primate lineage. Structure prediction programs suggest that the repetitive structure has no associated three-dimensional structure but rather a tendency for disorder. Splice variant isoforms were detected in several mammalian species but without conserved pattern, notably excluding the constant coexistence of premammalian-like transcripts, without the extension. The yeast two-hybrid method revealed that, in human, the extension allowed new interactions with ten unrelated proteins, most of them with specific three-dimensional structures involved in protein–protein interactions, and some highly expressed in testis, as is cyclin B3. The interactions with activator of cAMP-responsive element modulator in testis (ACT), germ cell-less homolog 1, and chromosome 1 open reading frame 14 remain to be verified in vivo since they may not be expressed at the same stages of spermatogenesis as cyclin B3.

Description: The X chromosome has a large effect on hybrid dysfunction, particularly on hybrid male sterility. Although the evidence for this so-called large-X effect is clear, its molecular causes are not yet fully understood. One possibility is that, under certain conditions, evolution proceeds faster in X-linked than in autosomal loci (i.e., faster-X effect) due to both natural selection and their hemizygosity in males, an effect that is expected to be greatest in genes with male-biased expression. Here, I study genome-wide variation in transcript abundance between Drosophila yakuba and D. santomea , within these species and in their hybrid males to evaluate both the faster-X and large-X effects at the level of expression. I find that in X-linked male-biased genes (MBGs) expression evolves faster than in their autosomal counterparts, an effect that is accompanied by a unique reduction in expression polymorphism. This suggests that Darwinian selection is driving expression differences between species, likely enhanced by the hemizygosity of the X chromosome in males. Despite the recent split of the two sister species under study, abundant changes in both cis - and trans -regulatory elements underlie expression divergence in the majority of the genes analyzed, with significant differences in allelic ratios of transcript abundance between the two reciprocal F 1 hybrid males. Cis–trans coevolution at molecular level, evolved shortly after populations become isolated, may therefore contribute to explain the breakdown of the regulation of gene expression in hybrid males. Additionally, the X chromosome plays a large role in this hybrid male misexpression, which affects not only MBG but also, to a lesser degree, nonsex-biased genes. Interestingly, hybrid male misexpression is concentrated mostly in autosomal genes, likely facilitated by the rapid evolution of sex-linked trans -acting factors. I suggest that the faster evolution of X-linked MBGs, at both protein and expression levels, contributes to explain the large effect of the X chromosome on hybrid male sterility, likely mediating widespread autosomal misexpression through the preferential recognition of cis -regulatory elements by conspecific trans -acting factors (i.e., cis–trans conspecific recognition).

Description: The sex chromosomes of the tropical crop papaya ( Carica papaya ) are evolutionarily young and consequently allow for the examination of evolutionary mechanisms that drive early sex chromosome divergence. We conducted a molecular population genetic analysis of four X/Y gene pairs from a collection of 45 wild papaya accessions. These population genetic analyses reveal striking differences in the patterns of polymorphism between the X and Y chromosomes that distinguish them from other sex chromosome systems. In most sex chromosome systems, the Y chromosome displays significantly reduced polymorphism levels, whereas the X chromosome maintains a level of polymorphism that is comparable to autosomal loci. However, the four papaya sex-linked loci that we examined display diversity patterns that are opposite this trend: the papaya X alleles exhibit significantly reduced polymorphism levels, whereas the papaya Y alleles maintain greater than expected levels of diversity. Our analyses suggest that selective sweeps in the regions of the X have contributed to this pattern while also revealing geographically restricted haplogroups on the Y. We discuss the possible role sexual selection and/or genomic conflict have played in shaping the contrasting patterns of polymorphism found for the papaya X and Y chromosomes.

Description: The HLA region shows diversity concerning the number and content of DRB genes present per haplotype. Similar observations are made for the equivalent regions in other primate species. To elucidate the evolutionary history of the various HLA-DRB genes, a large panel of intron sequences obtained from humans, chimpanzees, rhesus macaques, and common marmosets has been subjected to phylogenetic analyses. Special attention was paid to the presence and absence of particular transposable elements and/or to their segments. The sharing of different parts of the same long interspersed nuclear element-2 (LINE2, L2) and various Alu insertions by the species studied demonstrates that one precursor gene must have been duplicated several times before the Old World monkey (OWM) and hominid (HOM) divergence. At least four ancestral DRB gene families appear to have been present before the radiation of OWM and HOM, and one of these even predates the speciation of Old and New World primates. Two of these families represent the pseudogenes DRB6/DRB2 and DRB7 , which have been locked in the genomes of various primate species over long evolutionary time spans. Furthermore, all phylogenies of different intron segments show consistently that, apart from the pseudogenes, only DRB5 genes are shared by OWM and HOM, and they demonstrate the common history of certain DRB genes/lineages of humans and chimpanzees. In contrast, the evolutionary history of some other DRB loci is difficult to decipher, thus illustrating the complex history of the evolution of DRB genes due to a combination of mutations and recombination-like events. The selected approach allowed us to shed light on the ancestral DRB gene pool in primates and on the evolutionary relationship of the various HLA-DRB genes.

Description: Chemosensory genes are frequently the target of positive selection and are often present in large gene families, but little is known about heterogeneity of selection in these cases and its relation to function. Here, we use the vomeronasal-1 receptor ( V1R ) repertoire of mouse lemurs ( Microcebus spp.) as a model system to study patterns of selection of chemosensory genes at several different levels. Mouse lemurs are small nocturnal strepsirrhine primates and have a large (~200 loci) repertoire of V1R loci that are likely important for intraspecific pheromonal communication and interspecific interactions, for example, recognition of predator cues. We investigated signals and patterns of positive selection among the 105 identified full length V1R loci in the gray mouse lemur and within 7 V1R loci amplified across multiple mouse lemur species. Phylogenetic reconstructions of published sequences revealed at least nine monophyletic clusters of V1R s in gray mouse lemurs that have diversified since the split between lemurs and lorisoid primates. A large majority of clusters evolved under significant positive selection. Similar results were found in V1R s of closely related greater galagos. Comparison with function of related V1R clusters in mice suggested a potential relationship between receptor function and strength of selection. Interestingly, most codons identified as being under positive selection are located in the extracellular domains of the receptors and hence likely indicate the position of residues involved in ligand binding. Positive selection was also detected within five V1R loci (=71% of analyzed loci) sequenced from 6 to 10 mouse lemur species, indicating ongoing selection within the genus, which may be related to sexual selection and, potentially, speciation processes. Variation in strength of positive selection on V1R s showed no simple relationship to cluster size. The diversity of V1R loci in mouse lemurs reflects their adaptive evolution and is most likely related to the fundamental relevance of olfactory communication and predator recognition in these primates. Overall, adaptive evolution is the predominant mode of evolution of V1R loci at all levels, and the substantial heterogeneity in the strength of selection may be related to receptor function.

Description: The large number of sexually transmitted diseases and ocular trachoma cases that are caused globally each year by Chlamydia trachomatis has made this organism a World Health Organization priority for vaccine development. However, there is no gene transfer system for Chlamydia to help identify potential vaccine targets. To accelerate discoveries toward this goal, here we analyzed the broadest diversity of C. trachomatis genomes to date, including 25 geographically dispersed clinical and seven reference strains representing 14 of the 19 known serotypes. Strikingly, all 32 genomes were found to have evidence of DNA acquisition by homologous recombination in their history. Four distinct clades were identified, which correspond to all C. trachomatis disease phenotypes: lymphogranuloma venereum (LGV; Clade 1); noninvasive urogenital infections (Clade 2); ocular trachoma (Clade 3); and protocolitis (Clade 4; also includes some noninvasive urogenital infections). Although the ancestral relationship between clades varied, most strains acted as donor and recipient of recombination with no evidence for barriers to genetic exchange. The niche-specific LGV and trachoma clades have undergone less recombination, although the opportunity for mixing with strains from other clades that infect the rectal and ocular mucosa, respectively, is evident. Furthermore, there are numerous occasions for gene conversion events through sequential infections at the same anatomic sites. The size of recombinant segments is relatively small (~357 bp) compared with in vitro experiments of various C. trachomatis strains but is consistent with in vitro estimates for other bacterial species including Escherichia coli and Helicobacter pylori . Selection has also played a crucial role during the diversification of the organism. Clade 2 had the lowest nonsynonymous to synonymous ratio (d N /d S ) but the highest effect of recombination, which is consistent with the widespread occurrence of synonymous substitutions in recombined genomic segments. The trachoma Clade 3 had the highest d N /d S estimates, which may be caused by an increased effect of genetic drift from niche specialization and a reduced effective population size. The degree of drift, selection, and recombination in C. trachomatis suggests that the challenge will remain to identify genomic regions that are stable and cross protective for the development of an efficacious vaccine.

Description: We surveyed genetic variation in alr2 , an allodeterminant of the colonial hydroid Hydractinia symbiolongicarpus . We generated cDNA from a sample of 239 Hydractinia colonies collected at Lighthouse Point, Connecticut, and identified 473 alr2 alleles, 198 of which were unique. Rarefaction analysis suggested that the sample was near saturation. Most alleles were rare, with 86% occurring at frequencies of 1% or less. Alleles were highly variable, diverging on average by 18% of the amino acids in a predicted extracellular domain of the molecule. Analysis of 152 full-length alleles confirmed the existence of two structural types, defined by exons 4–8 of the gene. Several residues of the predicted immunoglobulin superfamily-like domains display signatures of positive selection. We also identified 77 unique alr2 pseudogene sequences from 85 colonies. Twenty-seven of these sequences matched expressed alr2 sequences from other colonies. This observation is consistent with pseudogenes contributing to alr2 diversification through sequence donation. A more limited collection of animals was made from a distant, relict population of H. symbiolongicarpus . Sixty percent of the unique sequences identified in this sample were found to match sequences from the Lighthouse Point population. The large number of alr2 alleles, their degree of divergence, the predominance of rare alleles in the population, their persistence over broad spatial and temporal scales, and the signatures of positive selection in multiple residues of the putative recognition domain paint a consistent picture of negative-frequency-dependent selection operating in this system. The genetic diversity observed at alr2 is comparable to that of the most highly polymorphic genetic systems known to date.

Description: The mechanosensory lateral line, found only in fishes and amphibians, is an important sense organ associated with aquatic life. Lateral line patterns differ among teleost, the most diverse vertebrate taxa, hypothetically in response to selective pressures from different aquatic habitats. In this article, we conduct evolutionary genomic analyses of 34 genes associated with lateral line system development in teleosts to elucidate the significance of contrasting evolutionary rates and changes in the protein coding sequences. We find that duplicated copies of these genes are preferentially retained in the teleost genomes and that episodic events of positive selection have occurred in 22 of the 30 postduplication branches. In general, teleost genes evolved at a faster rate relative to their tetrapod counterparts, and the mutation rates of 26 of the 34 genes differed among teleosts and tetrapods. We conclude that following whole genome duplication, evolutionary rates and episodic events of positive selection on the lateral line system development genes might have been one of the factors favoring the subsequent adaptive radiation of teleosts into diverse habitats. These results provide the foundation for further detailed explorations into lateral line system genes and the evolution of diverse phenotypes and adaptations.

Description: In any comparative studies striving to understand the similarities and differences of the living organisms at the molecular genetic level, the crucial first step is to establish the homology (orthology and paralogy) of genes between different organisms. Determination of the homology of genes becomes complicated when the genes have undergone a rapid divergence in sequence or when the involved genes are members of a gene family that has experienced a differential gain or loss of its constituents in different taxonomic groups. Organisms with duplicated genomes such as teleost fishes might have been especially prone to these problems because the functional redundancies provided by the duplicate copies of genes would have allowed a rapid divergence or loss of genes during evolution. In this study, we will demonstrate that much of the ambiguities in the determination of the homology between fish and tetrapod genes resulting from the problems like these can be eliminated by complementing the sequence-based phylogenies with nonsequence information, such as the exon–intron structure of a gene or the composition of a gene’s genomic neighbors. We will use the Tbx6/16 subfamily genes of zebrafish ( tbx6 , tbx16 , tbx24 , and mga genes), which have been well known for the ambiguity of their evolutionary relationships to the Tbx6/16 subfamily genes of tetrapods, as an illustrative example. We will show that, despite the similarity of sequence and expression to the tetrapod Tbx6 genes, zebrafish tbx6 gene is actually a novel T-box gene more closely related to the tetrapod Tbx16 genes, whereas the zebrafish tbx24 gene, hitherto considered to be a novel gene due to the high level of sequence divergence, is actually an ortholog of tetrapod Tbx6 genes. We will also show that, after their initial appearance by the multiplication of a common ancestral gene at the beginning of vertebrate evolution, the Tbx6/16 subfamily of vertebrate T-box genes might have experienced differential losses of member genes in different vertebrate groups and gradual pooling of member gene’s functions in surviving members, which might have prevented the revelation of the true identity of member genes by way of the comparison of sequence and function.

Description: Cyanobacteria are among the most ancient organisms known to have circadian rhythms. The cpmA gene is involved in controlling the circadian output signal. We studied polymorphism and divergence of this gene in six populations of a stress-tolerant cyanobacterium, Chroococcidiopsis sp., sampled in extreme habitats across the globe. Despite high haplotype diversity (0.774), nucleotide diversity of cpmA is very low ( = 0.0034): the gene appears to be even more conserved than housekeeping genes. Even though the populations were sampled thousands kilometers apart, they manifested virtually no genetic differentiation at this locus ( F ST = 0.0228). Using various tests for neutrality, we determined that evolution of cpmA significantly departures from the neutral model and is governed by episodic positive selection.

Description: Environmental drivers of biodiversity can be identified by relating patterns of community similarity to ecological factors. Community variation has traditionally been assessed by considering changes in species composition and more recently by incorporating phylogenetic information to account for the relative similarity of taxa. Here, we describe how an important class of measures including Bray–Curtis, Canberra, and UniFrac can be extended to allow community variation to be computed on a phylogenetic network. We focus on phylogenetic split systems, networks that are produced by the widely used median network and neighbor-net methods, which can represent incongruence in the evolutionary history of a set of taxa. Calculating β diversity over a split system provides a measure of community similarity averaged over uncertainty or conflict in the available phylogenetic signal. Our freely available software, Network Diversity, provides 11 qualitative (presence–absence, unweighted) and 14 quantitative (weighted) network-based measures of community similarity that model different aspects of community richness and evenness. We demonstrate the broad applicability of network-based diversity approaches by applying them to three distinct data sets: pneumococcal isolates from distinct geographic regions, human mitochondrial DNA data from the Indonesian island of Nias, and proteorhodopsin sequences from the Sargasso and Mediterranean Seas. Our results show that major expected patterns of variation for these data sets are recovered using network-based measures, which indicates that these patterns are robust to phylogenetic uncertainty and conflict. Nonetheless, network-based measures of community similarity can differ substantially from measures ignoring phylogenetic relationships or from tree-based measures when incongruent signals are present in the underlying data. Network-based measures provide a methodology for assessing the robustness of β-diversity results in light of incongruent phylogenetic signal and allow β diversity to be calculated over widely used network structures such as median networks.

Description: Although biallelic mutations in non-collagen genes account for 〈10% of individuals with osteogenesis imperfecta, the characterization of these genes has identified new pathways and potential interventions that could benefit even those with mutations in type I collagen genes. We identified mutations in FKBP10 , which encodes the 65 kDa prolyl cis–trans isomerase, FKBP65, in 38 members of 21 families with OI. These include 10 families from the Samoan Islands who share a founder mutation. Of the mutations, three are missense; the remainder either introduce premature termination codons or create frameshifts both of which result in mRNA instability. In four families missense mutations result in loss of most of the protein. The clinical effects of these mutations are short stature, a high incidence of joint contractures at birth and progressive scoliosis and fractures, but there is remarkable variability in phenotype even within families. The loss of the activity of FKBP65 has several effects: type I procollagen secretion is slightly delayed, the stabilization of the intact trimer is incomplete and there is diminished hydroxylation of the telopeptide lysyl residues involved in intermolecular cross-link formation in bone. The phenotype overlaps with that seen with mutations in PLOD2 (Bruck syndrome II), which encodes LH2, the enzyme that hydroxylates the telopeptide lysyl residues. These findings define a set of genes, FKBP10 , PLOD2 and SERPINH1 , that act during procollagen maturation to contribute to molecular stability and post-translational modification of type I procollagen, without which bone mass and quality are abnormal and fractures and contractures result.

Description: The GUCY2D gene encodes retinal membrane guanylyl cyclase (RetGC1), a key component of the phototransduction machinery in photoreceptors. Mutations in GUCY2D cause Leber congenital amaurosis type 1 (LCA1), an autosomal recessive human retinal blinding disease. The effects of RetGC1 deficiency on human rod and cone photoreceptor structure and function are currently unknown. To move LCA1 closer to clinical trials, we characterized a cohort of patients (ages 6 months—37 years) with GUCY2D mutations. In vivo analyses of retinal architecture indicated intact rod photoreceptors in all patients but abnormalities in foveal cones. By functional phenotype, there were patients with and those without detectable cone vision. Rod vision could be retained and did not correlate with the extent of cone vision or age. In patients without cone vision, rod vision functioned unsaturated under bright ambient illumination. In vitro analyses of the mutant alleles showed that in addition to the major truncation of the essential catalytic domain in RetGC1, some missense mutations in LCA1 patients result in a severe loss of function by inactivating its catalytic activity and/or ability to interact with the activator proteins, GCAPs. The differences in rod sensitivities among patients were not explained by the biochemical properties of the mutants. However, the RetGC1 mutant alleles with remaining biochemical activity in vitro were associated with retained cone vision in vivo . We postulate a relationship between the level of RetGC1 activity and the degree of cone vision abnormality, and argue for cone function being the efficacy outcome in clinical trials of gene augmentation therapy in LCA1.

Description: We implement an isolation with migration model for three species, with migration occurring between two closely related species while an out-group species is used to provide further information concerning gene trees and model parameters. The model is implemented in the likelihood framework for analyzing multilocus genomic sequence alignments, with one sequence sampled from each of the three species. The prior distribution of gene tree topology and branch lengths at every locus is calculated using a Markov chain characterization of the genealogical process of coalescent and migration, which integrates over the histories of migration events analytically. The likelihood function is calculated by integrating over branch lengths in the gene trees (coalescent times) numerically. We analyze the model to study the gene tree-species tree mismatch probability and the time to the most recent common ancestor at a locus. The model is used to construct a likelihood ratio test (LRT) of speciation with gene flow. We conduct computer simulations to evaluate the LRT and found that the test is in general conservative, with the false positive rate well below the significance level. For the test to have substantial power, hundreds of loci are needed. Application of the test to a human–chimpanzee–gorilla genomic data set suggests gene flow around the time of speciation of the human and the chimpanzee.

Description: Although both genotypes with elevated mutation rate (mutators) and mobilization of insertion sequence (IS) elements have substantial impact on genome diversification, their potential interactions are unknown. Moreover, the evolutionary forces driving gradual accumulation of these elements are unclear: Do these elements spread in an initially transposon-free bacterial genome as they enable rapid adaptive evolution? To address these issues, we inserted an active IS 1 element into a reduced Escherichia coli genome devoid of all other mobile DNA. Evolutionary laboratory experiments revealed that IS elements increase mutational supply and occasionally generate variants with especially large phenotypic effects. However, their impact on adaptive evolution is small compared with mismatch repair mutator alleles, and hence, the latter impede the spread of IS-carrying strains. Given their ubiquity in natural populations, such mutator alleles could limit early phase of IS element evolution in a new bacterial host. More generally, our work demonstrates the existence of an evolutionary conflict between mutation-promoting mechanisms.

Description: Muscles are composed of multinucleated muscle fibers with different contractile and physiological properties, which result from specific slow or fast gene expression programs in the differentiated muscle cells. In the zebra fish embryo, the slow program is under the control of Hedgehog signaling from the notochord and floor plate. This pathway activates the expression of the conserved transcriptional repressor, Prdm1 (Blimp1), which in turn represses the fast program and promotes the slow program in adaxial cells of the somite and their descendants. In the mouse embryo, myogenesis is also initiated in the myotomal compartment of the somite, but the slow muscle program is not confined to a specific subset of cells. We now show that Prdm1 is expressed in the first differentiated myocytes of the early myotome from embryonic day (E)9.5–E11.5. During this period, muscle formation depends on the myogenic regulatory factors, Myf5 and Mrf4. In their absence, Prdm1 is not activated, in apparent contrast to zebra fish where Prdm1 is expressed in the absence of Myf5 and MyoD that drive myogenesis in adaxial cells. However, as in zebra fish, Prdm1 expression in the mouse myotome does not occur in the absence of Hedgehog signaling. Analysis of the muscle phenotype of Prdm1 mutant embryos shows that myogenesis appears to proceed normally. Notably, there is no requirement for Prdm1 activation of the slow muscle program in the mouse myotome. Furthermore, the gene for the transcriptional repressor, Sox6, which is repressed by Prdm1 to permit slow muscle differentiation in zebra fish, is not expressed in the mouse myotome. We propose that the lack of functional conservation for mouse Prdm1, that can nevertheless partially rescue the adaxial cells of zebra fish Prdm1 mutants, reflects differences in the evolution of the role of key regulators such as Prdm1 or Sox6, in initiating the onset of the slow muscle program, between teleosts and mammals.

Description: The specific recognition of antigen by T cells is critical to the generation of adaptive immune responses in vertebrates. T cells recognize antigen using a somatically diversified T-cell receptor (TCR). All jawed vertebrates use four TCR chains called α, β, , and , which are expressed as either a αβ or heterodimer. Nonplacental mammals (monotremes and marsupials) are unusual in that their genomes encode a fifth TCR chain, called TCRµ, whose function is not known but is also somatically diversified like the conventional chains. The origins of TCRµ are also unclear, although it appears distantly related to TCR. Recent analysis of avian and amphibian genomes has provided insight into a model for understanding the evolution of the TCR genes in tetrapods that was not evident from humans, mice, or other commonly studied placental (eutherian) mammals. An analysis of the genes encoding the TCR chains in the duckbill platypus revealed the presence of a highly divergent variable (V) gene, indistinguishable from immunoglobulin heavy (IgH) chain V genes (VH) and related to V genes used in TCRµ. They are expressed as part of TCR repertoire (VH) and similar to what has been found in frogs and birds. This, however, is the first time a VH has been found in a mammal and provides a critical link in reconstructing the evolutionary history of TCRµ. The current structure of TCR and TCRµ genes in tetrapods suggests ancient and possibly recurring translocations of gene segments between the IgH and TCR genes, as well as translocations of TCR genes out of the TCRα/ locus early in mammals, creating the TCRµ locus.

Description: Starch is the most widespread and abundant storage carbohydrate in crops and its production is critical to both crop yield and quality. In regard to the starch content in the seeds of crop plants, there is a distinct difference between grasses (Poaceae) and dicots. However, few studies have described the evolutionary pattern of genes in the starch biosynthetic pathway in these two groups of plants. In this study, therefore, an attempt was made to compare evolutionary rate, gene duplication, and selective pattern of the key genes involved in this pathway between the two groups, using five grasses and five dicots as materials. The results showed 1) distinct differences in patterns of gene duplication and loss between grasses and dicots; duplication in grasses mainly occurred before the divergence of grasses, whereas duplication mostly occurred in individual species within the dicots; there is less gene loss in grasses than in dicots, 2) a considerably higher evolutionary rate in grasses than in dicots in most gene families analyzed, and 3) evidence of a different selective pattern between grasses and dicots; positive selection may have occurred asymmetrically in grasses in some gene families, for example, ADP-glucose pyrophosphorylase small subunit. Therefore, we deduced that gene duplication contributes to, and a higher evolutionary rate is associated with, the higher starch content in grasses. In addition, two novel aspects of the evolution of the starch biosynthetic pathway were observed.

Description: From studies investigating the differences in evolutionary rates between genes, gene compactness and gene expression level have been identified as important determinants of gene-level protein evolutionary rate, as represented by nonsynonymous to synonymous substitution rate ( d N / d S ) ratio. However, the causes of exon-level variances in d N / d S are less understood. Here, we use principal component regression to examine to what extent 13 exon features explain the variance in d N , d S , and the d N / d S ratio of human–rhesus macaque or human–mouse orthologous exons. The exon features were grouped into six functional categories: expression features, mRNA splicing features, structural–functional features, compactness features, exon duplicability, and other features, including G + C content and exon length. Although expression features are important for determining d N and d N / d S between exons of different genes, structural–functional features and splicing features explained more of the variance for exons of the same genes. Furthermore, we show that compactness features can explain only a relatively small percentage of variance in exon-level d N or d N / d S in either between-gene or within-gene comparison. By contrast, d S yielded inconsistent results in the human–mouse comparison and the human–rhesus macaque comparison. This inconsistency may suggest rapid evolutionary changes of the mutation landscape in mammals. Our results suggest that between-gene and within-gene variation in d N / d S (and d N ) are driven by different evolutionary forces and that the role of mRNA splicing in causing the variation in evolutionary rates of coding sequences may be underappreciated.

Description: Extensive synonymous codon modification of viral genomes appears to be an effective way of attenuating strains for use as live vaccines. An assumption of this method is that codon changes have individually small effects, such that codon-attenuated viruses will be slow to evolve back to high fitness (and thus to high virulence). The major capsid gene of the bacterial virus T7 was modified to have varying levels of suboptimal synonymous codons in different constructs, and fitnesses declined linearly with the number of changes. Adaptation of the most extreme design, with 182 codon changes, resulted in a slow fitness recovery by standards of previous experimental evolution with this virus, although fitness effects of substitutions were higher than expected from the average effect of an engineered codon modification. Molecular evolution during recovery was modest, and changes evolved both within the modified gene and outside it. Some changes within the modified gene evolved in parallel across replicates, but with no obvious explanation. Overall, the study supports the premise that codon-modified viruses recover fitness slowly, although the evolution is substantially more rapid than expected from the design principle.

Description: In the fungal kingdom, the evolution of mating systems is highly dynamic, varying even among closely related species. Rearrangements in the mating-type ( mat ) locus, which contains the major regulators of sexual development, are expected to underlie the transitions between self-sterility (heterothallism) and self-fertility (homothallism). However, both the genetic mechanisms and the direction of evolutionary transitions in fungal mating systems are under debate. Here, we present new sequences of the mat locus of four homothallic and one heterothallic species of the model genus Neurospora (Ascomycota). By examining the patterns of synteny among these sequences and previously published data, we show that the locus is conserved among heterothallic species belonging to distinct phylogenetic clades, while different gene arrangements characterize the four homothallic species. These results allowed us to ascertain a heterothallic ancestor for the genus, confirming the prediction of the dead-end theory on unidirectional transitions toward selfing. We show that at least four shifts from heterothallism to homothallism have occurred in Neurospora , three of which involve the acquisition of sequences of both mating types into the same haploid genome. We present evidence for two genetic mechanisms allowing these shifts: translocation and unequal crossover. Finally, we identified two novel retrotransposons and suggest that these have played a major role in mating-system transitions, by facilitating multiple rearrangements of the mat locus.

Description: Male patients with Peutz–Jeghers syndrome (PJS) have defective spermatogenesis and are at increased risk of developing Sertoli cell tumors. Mutations in the Liver Kinase B1 ( LKB1/STK11) gene are associated with the pathogenesis of PJS and have been identified in non-PJS patients with sporadic testicular cancers. The mechanisms controlled by LKB1 signaling in Sertoli cell functions and testicular biology have not been described. We have conditionally deleted the Lkb1 gene ( Lkb1 cko ) in somatic testicular cells to define the molecular mechanisms involved in the development of the testicular phenotype observed in PJS patients. Focal vacuolization in some of the seminiferous tubules was observed in 4-week-old mutant testes but germ cell development appeared to be normal. However, similar to PJS patients, we observed progressive germ cell loss and Sertoli cell only tubules in Lkb1 cko testes from mice older than 10 weeks, accompanied by defects in Sertoli cell polarity and testicular junctional complexes and decreased activation of the MAP/microtubule affinity regulating and focal adhesion kinases. Suppression of AMP kinase and activation of mammalian target of rapamycin (mTOR) signaling were also observed in Lkb1 cko testes. Loss of Tsc1 or Tsc2 copies the progressive Lkb1 cko phenotype, suggesting that dysregulated activation of mTOR contributes to the pathogenesis of the Lkb1 cko testicular phenotype. Pten cko mice had a normal testicular phenotype, which could be explained by the comparative lack of mTOR activation detected. These studies describe the importance of LKB1 signaling in testicular biology and the possible molecular mechanisms driving the pathogenesis of the testicular defects observed in PJS patients.

Description: The human X-linked macrosatellite DXZ4 is a large tandem repeat located at Xq23 that is packaged into heterochromatin on the male X chromosome and female active X chromosome and, in response to X chromosome, inactivation is organized into euchromatin bound by the insulator protein CCCTC-binding factor (CTCF) on the inactive X chromosome (Xi). The purpose served by this unusual epigenetic regulation is unclear, but suggests a Xi-specific gain of function for DXZ4. Other less extensive bands of euchromatin can be observed on the Xi, but the identity of the underlying DNA sequences is unknown. Here, we report the identification of two novel human X-linked tandem repeats, located 58 Mb proximal and 16 Mb distal to the macrosatellite DXZ4. Both tandem repeats are entirely contained within the transcriptional unit of novel spliced transcripts. Like DXZ4, the tandem repeats are packaged into Xi-specific CTCF-bound euchromatin. These sequences undergo frequent CTCF-dependent interactions with DXZ4 on the Xi, implicating DXZ4 as an epigenetically regulated Xi-specific structural element and providing the first putative functional attribute of a macrosatellite in the human genome.

Description: The apolipoprotein E ( APOE ) genotype is the major genetic risk factor for Alzheimer's disease (AD). We have access to cerebrospinal fluid (CSF) and plasma APOE protein levels from 641 individuals and genome-wide genotyped data from 570 of these samples. The aim of this study was to test whether CSF or plasma APOE levels could be a useful endophenotype for AD and to identify genetic variants associated with APOE levels. We found that CSF ( P = 8.15 x 10 –4 ) but not plasma ( P = 0.071) APOE protein levels are significantly associated with CSF Aβ 42 levels. We used Mendelian randomization and genetic variants as instrumental variables to confirm that the association of CSF APOE with CSF Aβ 42 levels and clinical dementia rating (CDR) is not because of a reverse causation or confounding effect. In addition the association of CSF APOE with Aβ 42 levels was independent of the APOE 4 genotype, suggesting that APOE levels in CSF may be a useful endophenotype for AD. We performed a genome-wide association study to identify genetic variants associated with CSF APOE levels: the APOE 4 genotype was the strongest single-genetic factor associated with CSF APOE protein levels ( P = 6.9 x 10 –13 ). In aggregate, the Illumina chip single nucleotide polymorphisms explain 72% of the variability in CSF APOE protein levels, whereas the APOE 4 genotype alone explains 8% of the variability. No other genetic variant reached the genome-wide significance threshold, but nine additional variants exhibited a P -value 〈10 –6 . Pathway mining analysis indicated that these nine additional loci are involved in lipid metabolism ( P = 4.49 x 10 –9 ).

Description: Single-molecule imaging is a powerful technique to visualize molecular interactions and movements. Translation is one of the most interesting targets for researchers with the molecular-imaging skills, since mRNA, tRNA and translation factors interact with or move inside or on the ribosome in an ordered manner. Trans -translation is a bacterial quality control system to rescue the ribosomes stalled at the 3' end of the mRNA, and this phenomenon is recapitulated in vitro with defined factors including two trans -translation-specific entities tmRNA and SmpB. Zhou et al. (Single molecule imaging of the trans -translation entry process via anchoring of the tagged ribosome. J Biochem 2011;149:609-618.) successfully visualized the interaction of the tmRNA–SmpB complex with the ribosome by immobilizing the ribosome on the quartz surface with the HaloTag technology. This ribosome-anchoring system may be useful for the imaging analysis of other processes of translation.

Description: The lipid mediator sphingosine-1-phosphate (S1P) is generated within cells from sphingosine by two sphingosine kinases (SPHK1 and SPHK2). Intracellularly synthesized S1P is released into the extracellular fluid by S1P transporters, including SPNS2. Released S1P binds specifically to the G protein-coupled S1P receptors (S1PR1/S1P 1 –S1PR5/S1P 5 ), which activate a diverse range of downstream signalling pathways. Recent studies have proposed that one of the central physiological functions of intercellular S1P signalling is in lymphocyte trafficking in vivo because genetic disruption of SPHK1/2, SPNS2 or S1PR1/S1P 1 in mice induces a lymphopenia phenotype. In this review, we discuss the current understanding of intercellular S1P signalling in the context of immunity.

Description: Cdc6 is the AAA+ ATPase that assembles prereplicative complexes on replication origins in eukaryotic chromosomes. Recently, the same Cdc6 protein was found to exert two more functions in mammalian cells to promote cell proliferation and survival: ATP-dependent activation of p21 CIP1 - or p27 KIP1 -bound Cdk2-cyclin A/E complexes and obstruction of apoptosome assembly and consequent cell death by forming stable complexes with activated Apaf-1 molecules. These findings not only redefined the biological role of mammalian Cdc6 but also led the discovery of entirely new mechanisms controlling Cdk2 activity and apoptosis. This review focuses on this multi-functional AAA+ ATPase and the newly discovered mechanisms by which it controls the G 1 –S transition and cell survival during proliferation.

Description: The pyrimidine reductive catabolic pathway is important for the utilization of uracil and thymine as sources of nitrogen and carbon. The pathway is controlled by three enzymes: dihydropyrimidine dehydrogenase (DPD), dihydropyrimidinase and β-alanine synthase. The putative DPD genes, pydX and pydA , are tandemly arranged in the Pseudomonas putida genome. Intriguingly, a putative transcriptional regulator, PydR, homologous to Escherichia coli RutR, a repressor of the Rut-dependent pyrimidine degradation pathway, is located downstream of pydX and pydA . In this study, we show that a pydA strain of P. putida fails to grow on a minimal media containing uracil or thymine as a sole nitrogen source, demonstrating the physiological importance of DPD in the reductive pathway. The expression of pydA and DPD activity in the absence of uracil were significantly higher in a pydR strain than in the wild-type strain, indicating that PydR acts as a repressor of the pyrimidine reductive pathway in P. putida . Phylogenetic analysis of RutR and PydR suggests that these homologous repressors may have evolved from a common ancestral protein that regulates pyrimidine degradation.

Description: Acute inflammation is an indispensable host response to foreign challenges or tissue injury. In healthy conditions, inflammatory processes are self-limiting and self-resolving, suggesting the existence of endogenous mechanisms for the control of inflammation and resolution. A comprehensive understanding of the cellular and molecular events of a well-orchestrated inflammatory response is required, and recent studies have uncovered the roles of endogenous lipid mediators derived from polyunsaturated fatty acids (i.e. lipoxins, resolvins, protectins) in controlling the resolution of inflammation. This review presents recent advances in understanding the formation and action of these mediators, especially focusing on the LC-MS/MS-based lipidomics approach and the emerging roles of eosinophils and eosinophil-derived lipid mediators in controlling acute inflammation and resolution.

Description: The cytokine transforming growth factor-beta (TGF-β) has multiple effects in both physiological and pathological conditions. TGF-β is secreted as part of a tripartite complex from which it must be released in order to bind to its receptor. Sequestration of latent TGF-β in the extracellular matrix (ECM) is crucial for proper mobilization of the latent cytokine and its activation. However, contrary to expectation, loss-of-function mutations in genes encoding certain matrix proteins that bind TGF-β yield elevated, rather than decreased, TGF-β levels, posing a ‘TGF-β paradox.’ In this review, we discuss recent findings concerning the relationship of TGF-β, ECM molecules, and latent TGF-β activation and propose a model to resolve the ‘TGF-β paradox.’

Description: Colicin E5 cleaves tRNAs for Tyr, His, Asn and Asp in their anticodons to abolish protein synthesis in Escherichia coli . We previously showed how its C-terminal RNase domain, E5-CRD, recognizes the anticodon bases but the catalytic mechanism remained to be elucidated. Although the reaction products with 5'-OH and 2',3'-cyclic phosphate ends suggested a similar mechanism to those of RNases A and T1, E5-CRD does not have the His residues necessary as a catalyst in usual RNases. To identify residues important for the catalytic reaction, mutants as to all residues within 5 Å from the central phosphorus of the scissile phosphodiester bond were prepared. Evaluation of the killing activities of the mutant colicins and the RNase activities of the mutant E5-CRDs suggested direct involvement of Arg33, Lys25, Gln29 and Lys60 in the reaction. Particularly, Arg33 plays a critical role and Ile94 provides a structural support of Arg33. Crystal structure of the complex of E5-CRD(R33Q)/dGpdUp showed structural and binding functional integrity of this mutant protein, suggesting involvement of Arg33 in the catalytic reaction. The structure of the free E5–CRD, we also determined, showed great flexibility of a flap region, which facilitates the access of Lys60 to the substrate in an induced-fit manner.

Description: Sulfatide (HSO 3 -3-galactosylceramide), which enriched in lipid rafts of plasma membranes in various epithelial cell lines, is a critical component of host cells for effective production of influenza A virus. However, the function of sulfatide in other virus infections targeting epithelial cells remains unknown. In this study, the effect of sulfatide on infection of human parainfluenza virus type 3 (hPIV3) was demonstrated by using genetically produced sulfatide-enriched cells and by treatment of hPIV3-infected cells with anti-sulfatide monoclonal antibody (GS-5) as well as by addition of sulfatide to the cells. hPIV3 was found to bind to sulfatide in a virus overlay assay and a solid-phase binding assay. Genetic expression of sulfatide in COS-7 cells defective in sulfatide suppressed initial hPIV3 infection and formation of multinucleate virus-infected cells. Treatment of virus-infected LLC-MK2 cells with GS-5 promoted formation of multinucleate cells. In contrast, exogenous addition of sulfatide to hPIV3-infected COS-7 cells and cells expressing the hPIV3- hemagglutinin-neuraminidase ( HN ) gene and fusion ( F ) gene conspicuously reduced the formation of multinucleate cells. The results suggest that sulfatide negatively regulates the fusion process of hPIV3, possibly through interaction with HN or F glycoprotein on the cell surface.

Description: Annexin A3 is a protein belonging to the annexin family, and it is mainly present in cellular membranes as a phospholipid-binding protein that binds via the calcium ion. However, its physiological function remains to be clarified. We examined the expression of annexin A3 in mouse tissues and found for the first time that annexin A3 mRNA and its protein were expressed more strongly in adipose tissues than in other tissues. In adipose tissues, annexin A3-expressing cells were present in the stromal vascular fraction, and precisely identical to Pref-1-positive preadipocytes, Pref-1 being an epidermal growth factor repeat-containing transmembrane protein that inhibits adipogenesis. In 3T3-L1 cells, used as a model of adipogenesis, annexin A3 was down-regulated at an early phase of adipocyte differentiation, and this pattern paralleled that of Pref-1. Suppression of annexin A3 in these cells with siRNA caused elevation of the PPAR2 mRNA level and lipid droplet accumulation. In conclusion, our data suggest that annexin A3 is a negative regulator of adipocyte differentiation.

Description: Protein phosphorylation by protein tyrosine (Tyr) kinases plays important roles in a variety of signalling pathways in cell growth, differentiation and oncogenesis in animals. Despite the absence of classical Tyr kinases in plants, a similar ratio of phosphotyrosine residues in phosphorylated proteins was found in Arabidopsis thaliana as in human. However, protein kinases responsible for tyrosine phosphorylation in plants except some dedicated dual-specificity kinases still remain unclear. In this study, we found that PKL01, a nuclear Dbf2-related (Ndr) kinase homologue in Lotus japonicus , was autophosphorylated at a tyrosine residue when it was expressed in Escherichia coli , but kinase-dead mutant of PKL01 was not. Tyrosine phophorylation site in PKL01 was identified as Tyr-56 by LC-MS/MS analysis. Recombinant PKL01, which had been dephosphorylated by an alkaline phosphatase, could be phosphorylated again at the Tyr residue when it was incubated with ATP. Furthermore, other Ndr kinases in plants and PKL01 phosphorylated on Tyr residues in the exogenous substrates such as poly(Glu, Tyr) 4:1 and casein. Therefore, the Ndr kinases in plants, which had been assumed as protein serine (Ser)/threonine (Thr) kinases, turned out to be dual-specificity kinases responsible for phosphorylation of Tyr residues and Ser/Thr residues in plant proteins.

Description: Mouse UDP-glucuronosyltransferase 1a6 (Ugt1a6) contains two functional copies of 1a6a and 1a6b that share high sequence homology (98%). Only 10 amino acids located around the substrate recognition region are different out of 531 total residues. Although Ugt1a6 plays important roles in conjugating phenolic compounds, the functional characteristics of these isozymes are unclear. We performed functional analyses of mouse Ugt1a6a and Ugt1a6b using two isomeric polyphenols ( trans - and cis -resveratrol). The cDNAs of mouse Ugt1a6a and Ugt1a6b were cloned and constructed as recombinant proteins using a yeast expression system, and kinetic parameters were evaluated. The wild-type Ugt1a6a and Ugt1a6b proteins catalysed trans - and cis -resveratrol 3- O -glucuronidation. Although the K m value for trans -resveratrol was significantly lower for Ugt1a6a compared with Ugt1a6b, the K m values for cis -resveratrol were comparable for the isozymes. Despite high sequence homology, significant kinetic differences were observed between the isozymes. To identify the critical residues for resveratrol glucuronidation, we constructed 10 variants of Ugt1a6a (T81P, N96R, H98Q, L100V, S104P, N115S, I117L, V118T, V119L and D120E). The I117L variant had Ugt1a6b-like enzymatic properties of K m in trans -resveratrol, and V max and K si in cis -form, suggesting that the residues located at position 117 of Ugt1a6a and Ugt1a6b play an important role in resveratrol glucuronidation.

Description: Estimates of the proportion of amino acid substitutions that have been fixed by selection ( α ) vary widely among taxa, ranging from zero in humans to over 50% in Drosophila . This wide range may reflect differences in the efficacy of selection due to differences in the effective population size ( N e ). However, most comparisons have been made among distantly related organisms that differ not only in N e but also in many other aspects of their biology. Here, we estimate α in three closely related lineages of house mice that have a similar ecology but differ widely in N e : Mus musculus musculus ( N e ~ 25,000–120,000), M. m. domesticus ( N e ~ 58,000–200,000), and M. m. castaneus ( N e ~ 200,000–733,000). Mice were genotyped using a high-density single nucleotide polymorphism array, and the proportions of replacement and silent mutations within subspecies were compared with those fixed between each subspecies and an outgroup, Mus spretus . There was significant evidence of positive selection in M. m. castaneus , the lineage with the largest N e , with α estimated to be approximately 40%. In contrast, estimates of α for M. m. domesticus ( α = 13%) and for M. m. musculus ( α = 12 %) were much smaller. Interestingly, the higher estimate of α for M. m. castaneus appears to reflect not only more adaptive fixations but also more effective purifying selection. These results support the hypothesis that differences in N e contribute to differences among species in the efficacy of selection.

Description: Neanderthals have been shown to share more genetic variants with present-day non-Africans than Africans. Recent admixture between Neanderthals and modern humans outside of Africa was proposed as the most parsimonious explanation for this observation. However, the hypothesis of ancient population structure within Africa could not be ruled out as an alternative explanation. We use simulations to test whether the site frequency spectrum, conditioned on a derived Neanderthal and an ancestral Yoruba (African) nucleotide (the doubly conditioned site frequency spectrum [ dcfs ]), can distinguish between models that assume recent admixture or ancient population structure. We compare the simulations to the dcfs calculated from data taken from populations of European, Chinese, and Japanese descent in the Complete Genomics Diversity Panel. Simulations under a variety of plausible demographic parameters were used to examine the shape of the dcfs for both models. The observed shape of the dcfs cannot be explained by any set of parameter values used in the simulations of the ancient structure model. The dcfs simulations for the recent admixture model provide a good fit to the observed dcfs for non-Africans, thereby supporting the hypothesis that recent admixture with Neanderthals accounts for the greater similarity of Neanderthals to non-Africans than Africans.

Description: The Trinidadian pike cichlid ( Crenicichla frenata ) is a major predator of the guppy ( Poecilia reticulata ), a model system for visual ecology research, and visual predation by the pike cichlid is known to select for male guppies with reduced short-wavelength reflectance. However, an early study of the pike cichlid’s visual system suggested a lack of short-wavelength–sensitive cone photoreceptors, a surprising finding as many African cichlids have highly developed short-wavelength vision. In this study, we found evidence for only four expressed cone opsins (LWS, RH2a, SWS2a, and SWS2b), plus one pseudogene (RH2b). Taken together with our microspectrophotometry data, which revealed the presence of three types of cone photoreceptor, including one sensitive to short-wavelength light, this would indicate a broader spectral capacity than previously believed from earlier visual studies of this fish. Relative to the highly diverse African cichlids, however, this Neotropical cichlid appears to have a greatly reduced opsin complement, reflecting both gene loss along the Neotropical lineage (lacking functional RH2b and, possibly, SWS1 opsins) and gene duplication within the African clade (which possesses paralogous RH2aα and RH2aβ opsins). Molecular evolutionary analyses show that positive selection has shaped the SWS2b and RH1 opsins along the Neotropical lineage, which may be indicative of adaptive evolution to alter nonspectral aspects of opsin biology. These results represent the first molecular evolutionary study of visual pigments in a Neotropical cichlid and thus provide a foundation for further study of a morphologically and ecologically diverse clade that has been understudied with respect to the link between visual ecology and diversification.

Description: On an evolutionary time scale, polymorphic alleles are believed to have a short life, persisting at most tens of millions of years even under long-term balancing selection. Here, we report highly diverged trans -species dimorphism of the proteasome subunit beta type 8 ( PSMB8 ) gene, which encodes a catalytic subunit of the immunoproteasome responsible for the generation of peptides presented by major histocompatibility complex (MHC) class I molecules, in lower teleosts including Cypriniformes (zebrafish and loach) and Salmoniformes (trout and salmon), whose last common ancestor dates to 300 Ma. Moreover, phylogenetic analyses indicated that these dimorphic alleles share lineages with two shark paralogous genes, suggesting that these two lineages have been maintained for more than 500 My either as alleles or as paralogs, and that conversion between alleles and paralogs has occurred at least once during vertebrate evolution. Two lineages termed PSMB8A and PSMB8F show an A 31 F substitution that would probably affect their cleaving specificity, and whereas the PSMB8A lineage has been retained by all analyzed jawed vertebrates, the PSMB8F lineage has been lost by most jawed vertebrates except for cartilaginous fish and basal teleosts. However, a possible functional equivalent of the PSMB8F lineage has been revived as alleles within the PSMB8A lineage at least twice during vertebrate evolution in the amphibian Xenopus and teleostean Oryzias species. Dynamic evolution of the PSMB8 polymorphism through long-term persistence, loss, and regaining of dimorphism and conversion between alleles and paralogs implies the presence of strong selective pressure for functional polymorphism of this gene.

Description: Adaptive immune systems are present only in vertebrates. How do all the remaining animals withstand continuous attacks of permanently evolving pathogens? Even in the absence of adaptive immunity, every organism must be able to unambiguously distinguish "self" cells from any imaginable "nonself." Here, we analyzed the function of highly polymorphic gene vCRL1 , which is expressed in follicle and blood cells of Ciona intestinalis , pointing to possible recognition roles either during fertilization or in immune reactions. By using segregation analysis, we demonstrate that vCRL1 locus is not involved in the control of self-sterility. Interestingly, genetic knockdown of vCRL1 in all tissues or specifically in hemocytes results in a drastic developmental arrest during metamorphosis exactly when blood system formation in Ciona normally occurs. Our data demonstrate that vCRL1 gene might be essential for the establishment of a functional blood system in Ciona . Presumably, presence of the vCRL1 receptor on the surface of blood cells renders them as self, whereas any cell lacking it is referred to as nonself and will be consequently destroyed. We propose that individual-specific receptor vCRL1 might be utilized to facilitate somatic self/nonself discrimination.

Description: In the age of whole-genome population genetics, so-called genomic scan studies often conclude with a long list of putatively selected loci. These lists are then further scrutinized to annotate these regions by gene function, corresponding biological processes, expression levels, or gene networks. Such annotations are often used to assess and/or verify the validity of the genome scan and the statistical methods that have been used to perform the analyses. Furthermore, these results are frequently considered to validate "true-positives" if the identified regions make biological sense a posteriori. Here, we show that this approach can be potentially misleading. By simulating neutral evolutionary histories, we demonstrate that it is possible not only to obtain an extremely high false-positive rate but also to make biological sense out of the false-positives and construct a sensible biological narrative. Results are compared with a recent polymorphism data set from Drosophila melanogaster .

Description: Synonymous codons are widely selected for various biological mechanisms in both prokaryotes and eukaryotes. Recent evidence suggests that microRNA (miRNA) function may affect synonymous codon choices near miRNA target sites. To better understand this, we perform genome-wide analysis on synonymous codon usage around miRNA target sites in four plant genomes. We observed a general trend of increased site accessibility around miRNA target sites in plants. Guanine-cytosine (GC)–poor codons are preferred in the flank region of miRNA target sites. Within-genome analyses show significant variation among miRNA targets in species. GC content of the target gene can partly explain the variation of site accessibility among miRNA targets. miRNA targets in GC-rich genes show stronger selection signals than those in GC-poor genes. Gene’s codon usage bias and the conservation level of miRNA and its target also have some effects on site accessibility, but the expression level of miRNA or its target and the mechanism of miRNA activity do not contribute to site accessibility differences among miRNA targets. We suggest that synonymous codons near miRNA targets are selected for efficient miRNA binding and proper miRNA function. Our results present a new dimension of natural selection on synonymous codons near miRNA target sites in plants, which will have important implications of coding sequence evolution.

Description: Chloroplasts have evolved from a cyanobacterial endosymbiont and their continuity has been maintained by chloroplast division, which is performed by the constriction of a ring-like division complex at the division site. It is believed that the synchronization of the endosymbiotic and host cell division events was a critical step in establishing a permanent endosymbiotic relationship, such as is commonly seen in existing algae. In the majority of algal species, chloroplasts divide once per specific period of the host cell division cycle. In order to understand both the regulation of the timing of chloroplast division in algal cells and how the system evolved, we examined the expression of chloroplast division genes and proteins in the cell cycle of algae containing chloroplasts of cyanobacterial primary endosymbiotic origin (glaucophyte, red, green, and streptophyte algae). The results show that the nucleus-encoded chloroplast division genes and proteins of both cyanobacterial and eukaryotic host origin are expressed specifically during the S phase, except for FtsZ in one graucophyte alga. In this glaucophyte alga, FtsZ is persistently expressed throughout the cell cycle, whereas the expression of the nucleus-encoded MinD and MinE as well as FtsZ ring formation are regulated by the phases of the cell cycle. In contrast to the nucleus-encoded division genes, it has been shown that the expression of chloroplast-encoded division genes is not regulated by the host cell cycle. The endosymbiotic gene transfer of minE and minD from the chloroplast to the nuclear genome occurred independently on multiple occasions in distinct lineages, whereas the expression of nucleus-encoded MIND and MINE is regulated by the cell cycle in all lineages examined in this study. These results suggest that the timing of chloroplast division in algal cell cycle is restricted by the cell cycle-regulated expression of some but not all of the chloroplast division genes. In addition, it is suggested that the regulation of each division-related gene was established shortly after the endosymbiotic gene transfer, and this event occurred multiple times independently in distinct genes and in distinct lineages.

Description: During allopolyploid speciation, two divergent nuclear genomes merge, yet only one (usually the maternal) of the two sets of progenitor organellar genomes is maintained. Rubisco (1,5-bisphosphate carboxylase/oxygenase) is composed of nuclear-encoded small subunits (SSUs) and plastome-encoded large subunits (LSUs), providing an ideal system to explore the evolutionary process of cytonuclear accommodation. Here, we take initial steps in this direction, using Gossypium allopolyploids as our model. SSU copies from divergent (5–10 My) progenitor diploids ("A" and "D" genomes) were combined at the time of polyploid formation 1–2 Ma, with the LSU encoded by the maternal A-genome parent. LSU genes from A- and D-genome diploids and AD-genome allopolyploids were sequenced, revealing several nonsynonymous substitutions and suggesting the possibility of differential selection on the nuclear-encoded rbcS partner following allopolyploid formation. Sequence data for the rbcS gene family revealed nonreciprocal homoeologous recombination between A- and D- rbcS homoeologs in all polyploid species but not in a synthetic intergenomic F1 hybrid, demonstrating "gene conversion" during allopolyploid evolution. All progenitor rbcS genes are retained and expressed in the five extant allopolyploid species, but analysis of the leaf transcriptome showed that A-homoeologs are preferentially expressed in both the allopolyploid and hybrid, consistent with the maternal origin of rbcL . Although rbcS genes from both progenitor genomes are expressed, some appear to have experienced mutations that may represent cytonuclear coevolution.

Description: In the genome of Artiodactyla (cow, sheep, pigs, camels, and whales), a major retroposon group originated from a presumable horizontal transfer of BovB, a retrotransposon-like element retroposon, between 52 and 70 million years ago. Since then, BovB retroposons have proliferated and today occupy a quarter of the cow’s genome sequence. The BovB-related short interspersed elements (SINEs) were used for resolving the phylogeny of Bovinae (cows, spiral-horned antelopes, and nilgais) and their relatives. In silico screening of 55,000 intronic retroposon insertions in the cow genome and experimental validation of 126 introns resulted in 29 informative retroposon markers for resolving bovine evolutionary relationships. A transposition-in-transposition analysis identifies three different phases of SINE activity and show how BovB elements have expanded in the cattle genome.

Description: The use of codon substitution models to compare synonymous and nonsynonymous substitution rates is a widely used approach to detecting positive Darwinian selection affecting protein evolution. However, in several recent papers, Hughes and colleagues claim that codon-based likelihood-ratio tests (LRTs) are logically flawed as they lack prior hypotheses and fail to accommodate random fluctuations in synonymous and nonsynonymous substitutions Friedman and Hughes (2007) also used site-based LRTs to analyze 605 gene families consisting of human and mouse paralogues. They found that the outcome of the tests was largely determined by irrelevant factors such as the GC content at the third codon positions and the synonymous rate d S , but not by the nonsynonymous rate d N or the d N / d S ratio, factors that should be related to selection. Here, we reanalyze those data. Contra Friedman and Hughes, we found that the test results are related to sequence length and the average d N / d S ratio. We examine the criticisms of Hughes and suggest that they are based on misunderstandings of the codon models and on statistical errors. Our analyses suggest that codon-based tests are useful tools for comparative analysis of genomic data sets.

Description: Giardia intestinalis is a major cause of waterborne enteric disease in humans. The species is divided into eight assemblages suggested to represent separate Giardia species based on host specificities and the genetic divergence of marker genes. We have investigated whether genome-wide recombination occurs between assemblages using the three available G. intestinalis genomes. First, the relative nonsynonymous substitution rates of the homologs were compared for 4,009 positional homologs. The vast majority of these comparisons indicate genetic isolation without interassemblage recombinations. Only a region of 6 kbp suggests genetic exchange between assemblages A and E, followed by gene conversion events. Second, recombination-detecting software fails to identify within-gene recombination between the different assemblages for most of the homologs. Our results indicate very low frequency of recombination between the syntenic core genes, suggesting that G. intestinalis assemblages are genetically isolated lineages and thus should be viewed as separated Giardia species.

Description: During the early stages of speciation, interspecific gene flow may be impeded by deleterious epistatic interactions in hybrids, which maintain parental allelic combinations at the speciation genes. The resulting semipermeable nature of the barrier to interspecific gene flow provides a valuable framework to identify the genes involved in hybrid mortality or sterility, as well as the evolutionary mechanisms that initially caused their divergence. The two Atlantic eels Anguilla anguilla and A. rostrata are partially isolated sister species that naturally hybridize, but whose genetic basis of postzygotic isolation remains unknown. We collected high-throughput sequencing data from the transcriptomes of 58 individuals and discovered 94 genes showing differentially fixed mutations between species. Evidence for positive selection at nuclear diagnostic genes was obtained using multilocus extensions of the McDonald–Kreitman test with polymorphism data from each species. In contrast, mitochondrial protein-coding genes experienced strong purifying selection and mostly diverged at synonymous sites, except for the mt-atp6 gene, which showed an atypically high nonsynonymous to synonymous rate ratio. Nuclear-encoded protein interactors of the mt-atp6 gene in the ATP synthase complex were significantly overrepresented in the list of nuclear diagnostic genes. Further analysis of resequencing data showed that positive selection has operated at both the mt-atp6 gene and its nuclear interactor atp5c1 . These findings suggest that a cytonuclear incompatibility caused by a disruption of normal ATP synthase function in hybrids contributes to partial reproductive isolation between European and American eels.

Description: A number of mouse models for spinal muscular atrophy (SMA) have been genetically engineered to recapitulate the severity of human SMA by using a targeted null mutation at the mouse Smn1 locus coupled with the transgenic addition of varying copy numbers of human SMN2 genes. Although this approach has been useful in modeling severe SMA and very mild SMA, a mouse model of the intermediate form of the disease would provide an additional research tool amenable for drug discovery. In addition, many of the previously engineered SMA strains are multi-allelic by design, containing a combination of transgenes and targeted mutations in the homozygous state, making further genetic manipulation difficult. A new genetic engineering approach was developed whereby variable numbers of SMN2 sequences were incorporated directly into the murine Smn1 locus. Using combinations of these alleles, we generated an allelic series of SMA mouse strains harboring no, one, two, three, four, five, six or eight copies of SMN2. We report here the characterization of SMA mutants in this series that displayed a range in disease severity from embryonic lethal to viable with mild neuromuscular deficits.

Description: Insulin resistance (IR) is a key determinant of type 2 diabetes (T2D) and other metabolic disorders. This genome-wide association study (GWAS) was designed to shed light on the genetic basis of fasting insulin (FI) and IR in 927 non-diabetic African Americans. 5 396 838 single-nucleotide polymorphisms (SNPs) were tested for associations with FI or IR with adjustments for age, sex, body mass index, hypertension status and first two principal components. Genotyped SNPs ( n = 12) with P 〈 5 x 10 –6 in African Americans were carried forward for de novo genotyping in 570 non-diabetic West Africans. We replicated SNPs in or near SC4MOL and TCERG1L in West Africans. The meta-analysis of 1497 African Americans and West Africans yielded genome-wide significant associations for SNPs in the SC4MOL gene: rs17046216 ( P = 1.7 x 10 –8 and 2.9 x 10 –8 for FI and IR, respectively); and near the TCERG1L gene with rs7077836 as the top scoring ( P = 7.5 x 10 –9 and 4.9 x 10 –10 for FI and IR, respectively). In silico replication in the MAGIC study ( n = 37 037) showed weak but significant association (adjusted P -value of 0.0097) for rs34602777 in the MYO5A gene. In addition, we replicated previous GWAS findings for IR and FI in Europeans for GCKR , and for variants in four T2D loci ( FTO , IRS1 , KLF14 and PPARG ) which exert their action via IR. In summary, variants in/near SC4MOL , and TCERG1L were associated with FI and IR in this cohort of African Americans and were replicated in West Africans. SC4MOL is under-expressed in an animal model of T2D and plays a key role in lipid biosynthesis, with implications for the regulation of energy metabolism, obesity and dyslipidemia. TCERG1L is associated with plasma adiponectin, a key modulator of obesity, inflammation, IR and diabetes.

Description: Loss of dystrophin protein due to mutations in the DMD gene causes Duchenne muscular dystrophy. Dystrophin loss also leads to the loss of the dystrophin glycoprotein complex (DGC) from the sarcolemma which contributes to the dystrophic phenotype. Tyrosine phosphorylation of dystroglycan has been identified as a possible signal to promote the proteasomal degradation of the DGC. In order to test the role of tyrosine phosphorylation of dystroglycan in the aetiology of DMD, we generated a knock-in mouse with a phenylalanine substitution at a key tyrosine phosphorylation site in dystroglycan, Y890. Dystroglycan knock-in mice ( Dag1 Y890F/Y890F ) had no overt phenotype. In order to examine the consequence of blocking dystroglycan phosphorylation on the aetiology of dystrophin-deficient muscular dystrophy, the Y890F mice were crossed with mdx mice an established model of muscular dystrophy. Dag1 Y890F/Y890F / mdx mice showed a significant improvement in several parameters of muscle pathophysiology associated with muscular dystrophy, including a reduction in centrally nucleated fibres, less Evans blue dye infiltration and lower serum creatine kinase levels. With the exception of dystrophin, other DGC components were restored to the sarcolemma including α-sarcoglycan, α-/β-dystroglycan and sarcospan. Furthermore, Dag1 Y890F/Y890F / mdx showed a significant resistance to muscle damage and force loss following repeated eccentric contractions when compared with mdx mice. While the Y890F substitution may prevent dystroglycan from proteasomal degradation, an increase in sarcolemmal plectin appeared to confer protection on Dag1 Y890F/Y890F / mdx mouse muscle. This new model confirms dystroglycan phosphorylation as an important pathway in the aetiology of DMD and provides novel targets for therapeutic intervention.

Description: Facioscapulohumeral muscular dystrophy (FSHD), the most prevalent myopathy afflicting both children and adults, is predominantly associated with contractions in the 4q35-localized macrosatellite D4Z4 repeat array. Recent studies have proposed that FSHD pathology is caused by the misexpression of the DUX4 (double homeobox 4) gene resulting in production of a pathogenic protein, DUX4-FL, which has been detected in FSHD, but not in unaffected control myogenic cells and muscle tissue. Here, we report the analysis of DUX4 mRNA and protein expression in a much larger collection of myogenic cells and muscle biopsies derived from biceps and deltoid muscles of FSHD affected subjects and their unaffected first-degree relatives. We confirmed that stable DUX4-fl mRNA and protein were expressed in myogenic cells and muscle tissues derived from FSHD affected subjects, including several genetically diagnosed adult FSHD subjects yet to show clinical manifestations of the disease in the assayed muscles. In addition, we report DUX4-fl mRNA and protein expression in muscle biopsies and myogenic cells from genetically unaffected relatives of the FSHD subjects, although at a significantly lower frequency. These results establish that DUX4-fl expression per se is not sufficient for FSHD muscle pathology and indicate that quantitative modifiers of DUX4-fl expression and/or function and family genetic background are determinants of FSHD muscle disease progression.

Description: Abnormal presence of autophagic vacuoles is evident in brains of patients with Parkinson's disease (PD), in contrast to the rare detection of autophagosomes in a normal brain. However, the actual cause and pathological significance of these observations remain unknown. Here, we demonstrate a role for mitochondrial metabolism in the regulation of the autophagy-lysosomal pathway in ex vivo and in vitro models of PD. We show that transferring mitochondria from PD patients into cells previously depleted of mitochondrial DNA is sufficient to reproduce the alterations in the autophagic system observed in PD patient brains. Although the initial steps of this pathway are not compromised, there is an increased accumulation of autophagosomes associated with a defective autophagic activity. We prove that this functional decline was originated from a deficient mobilization of autophagosomes from their site of formation toward lysosomes due to disruption in microtubule-dependent trafficking. This contributed directly to a decreased proteolytic flux of α-synuclein and other autophagic substrates. Our results lend strong support for a direct impact of mitochondria in autophagy as defective autophagic clearance ability secondary to impaired microtubule trafficking is driven by dysfunctional mitochondria. We uncover mitochondria and mitochondria-dependent intracellular traffic as main players in the regulation of autophagy in PD.

Description: Resistin is a polypeptide hormone that was reported to be associated with insulin resistance, inflammation and risk of type 2 diabetes and cardiovascular disease. We conducted a genome-wide association (GWA) study on circulating resistin levels in individuals of European ancestry drawn from the two independent studies: the Nurses' Health Study ( n = 1590) and the Health, Aging and Body Composition Study ( n = 1658). Single-nucleotide polymorphisms (SNPs) identified in the GWA analysis were replicated in an independent cohort of Europeans: the Gargano Family Study ( n = 659). We confirmed the association with a previously known locus, the RETN gene (19p13.2), and identified two novel loci near the TYW3/CRYZ gene (1p31) and the NDST4 gene (4q25), associated with resistin levels at a genome-wide significant level, best represented by SNP rs3931020 ( P = 6.37 x 10 –12 ) and SNP rs13144478 ( P = 6.19 x 10 –18 ), respectively. Gene expression quantitative trait loci analyses showed a significant cis association between the SNP rs3931020 and CRYZ gene expression levels ( P = 3.68 x 10 –7 ). We also found that both of these two SNPs were significantly associated with resistin gene ( RETN ) mRNA levels in white blood cells from 68 subjects with type 2 diabetes (both P = 0.02). In addition, the resistin-rising allele of the TYW3/CRYZ SNP rs3931020, but not the NDST4 SNP rs13144478, showed a consistent association with increased coronary heart disease risk [odds ratio = 1.18 (95% CI, 1.03–1.34); P = 0.01]. Our results suggest that genetic variants in TYW3/CRYZ and NDST4 loci may be involved in the regulation of circulating resistin levels. More studies are needed to verify the associations of the SNP rs13144478 with NDST4 gene expression and resistin-related disease.

Description: Gibbons (Hylobatidae) are small, arboreal apes indigenous to Southeast Asia that diverged from other apes ~15–18 Ma. Extant lineages radiated rapidly 6–10 Ma and are organized into four genera ( Hylobates , Hoolock , Symphalangus , and Nomascus ) consisting of 12–19 species. The use of short interspersed elements (SINEs) as phylogenetic markers has seen recent popularity due to several desirable characteristics: the ancestral state of a locus is known to be the absence of an element, rare potentially homoplasious events are relatively easy to resolve, and samples can be quickly and inexpensively genotyped. During radiation of primates, one particular family of SINEs, the Alu family, has proliferated in primate genomes. Nomascus leucogenys (northern white-cheeked gibbon) sequences were analyzed for repetitive content with RepeatMasker using a custom library. The sequences containing Alu elements identified as members of a gibbon-specific subfamily were then compared with orthologous positions in other primate genomes. A primate phylogenetic panel consisting of 18 primate species, including 13 gibbon species representing all four extant genera, was assayed for all loci, and a total of 125 gibbon-specific Alu insertions were identified. The resulting amplification patterns were used to generate a phylogenetic tree. We demonstrate significant support for Symphalangus as the most basal lineage within the family. Our findings also place Nomascus as a derived lineage, sister to Hoolock , with the Nomascus–Hoolock clade sister to Hylobates . Further, our analysis groups N. leucogenys and Nomascus siki as sister taxa to the exclusion of the other Nomascus species assayed. This study represents the first use of SINEs to determine the genus level phylogenetic relationships within the family Hylobatidae. These relationships have been resolved with robust support at most internal nodes, demonstrating the utility of SINE-based phylogenetic analysis. We postulate that hybridization and rapid radiation may have contributed to the complex and contradictory findings of the previous studies. Our findings will aid in the conservation of these threatened primates and inform future studies of the biogeographical history and distribution of modern gibbon species.

Description: Translation termination is accomplished by proteins of the Class I release factor family (RF) that recognize stop codons and catalyze the ribosomal release of the newly synthesized peptide. Bacteria have two canonical RFs: RF1 recognizes UAA and UAG, RF2 recognizes UAA and UGA. Despite that these two release factor proteins are sufficient for de facto translation termination, the eukaryotic organellar RF protein family, which has evolved from bacterial release factors, has expanded considerably, comprising multiple subfamilies, most of which have not been functionally characterized or formally classified. Here, we integrate multiple sources of information to analyze the remarkable differentiation of the RF family among organelles. We document the origin, phylogenetic distribution and sequence structure features of the mitochondrial and plastidial release factors: mtRF1a, mtRF1, mtRF2a, mtRF2b, mtRF2c, ICT1, C12orf65, pRF1, and pRF2, and review published relevant experimental data. The canonical release factors (mtRF1a, mtRF2a, pRF1, and pRF2) and ICT1 are derived from bacterial ancestors, whereas the others have resulted from gene duplications of another release factor. These new RF family members have all lost one or more specific motifs relevant for bona fide release factor function but are mostly targeted to the same organelle as their ancestor. We also characterize the subset of canonical release factor proteins that bear nonclassical PxT/SPF tripeptide motifs and provide a molecular-model-based rationale for their retained ability to recognize stop codons. Finally, we analyze the coevolution of canonical RFs with the organellar genetic code. Although the RF presence in an organelle and its stop codon usage tend to coevolve, we find three taxa that encode an RF2 without using UGA stop codons, and one reverse scenario, where mamiellales green algae use UGA stop codons in their mitochondria without having a mitochondrial type RF2. For the latter, we put forward a "stop-codon reinvention" hypothesis that involves the retargeting of the plastid release factor to the mitochondrion.

Description: The population history of the indigenous populations in island Southeast Asia is generally accepted to have been shaped by two major migrations: the ancient "Out of Africa" migration ~50,000 years before present (YBP) and the relatively recent "Out of Taiwan" expansion of Austronesian agriculturalists approximately 5,000 YBP. The Negritos are believed to have originated from the ancient migration, whereas the majority of island Southeast Asians are associated with the Austronesian expansion. We determined 86 mitochondrial DNA (mtDNA) complete genome sequences in four indigenous Malaysian populations, together with a reanalysis of published autosomal single-nucleotide polymorphism (SNP) data of Southeast Asians to test the plausibility and impact of those migration models. The three Austronesian groups (Bidayuh, Selatar, and Temuan) showed high frequencies of mtDNA haplogroups, which originated from the Asian mainland ~30,000–10,000 YBP, but low frequencies of "Out of Taiwan" markers. Principal component analysis and phylogenetic analysis using autosomal SNP data indicate a dichotomy between continental and island Austronesian groups. We argue that both the mtDNA and autosomal data suggest an "Early Train" migration originating from Indochina or South China around the late-Pleistocene to early-Holocene period, which predates, but may not necessarily exclude, the Austronesian expansion.

Description: Haploid genomes greater than 25,000 Mb are rare, within the animals only the lungfish and some of the salamanders and crustaceans are known to have genomes this large. There is very little data on the structure of genomes this size. It is known, however, that for animal genomes up to 3,000 Mb, there is in general a good correlation between genome size and the percent of the genome composed of repetitive sequence and that this repetitive component is highly dynamic. In this study, we sampled the Australian lungfish genome using three mini-genomic libraries and found that with very little sequence, the results converged on an estimate of 40% of the genome being composed of recognizable transposable elements (TEs), chiefly from the CR1 and L2 long interspersed nuclear element clades. We further characterized the CR1 and L2 elements in the lungfish genome and show that although most CR1 elements probably represent recent amplifications, the L2 elements are more diverse and are more likely the result of a series of amplifications. We suggest that our sampling method has probably underestimated the recognizable TE content. However, on the basis of the most likely sources of error, we suggest that this very large genome is not largely composed of recently amplified, undetected TEs but may instead include a large component of older degenerate TEs. Based on these estimates, and on Thomson’s (Thomson K. 1972. An attempt to reconstruct evolutionary changes in the cellular DNA content of lungfish. J Exp Zool. 180:363–372) inference that in the lineage leading to the extant Australian lungfish, there was massive increase in genome size between 350 and 200 mya, after which the size of the genome changed little, we speculate that the very large Australian lungfish genome may be the result of a massive amplification of TEs followed by a long period with a very low rate of sequence removal and some ongoing TE activity.

Description: L1 elements are mammalian non–long terminal repeat retrotransposons, or long interspersed elements (LINEs), that significantly influence the dynamics and fluidity of the genome. A series of observations suggest that plant L1-clade LINEs, just as mammalian L1s, mobilize both short interspersed elements (SINEs) and certain messenger RNA by recognizing the 3'-poly(A) tail of RNA. However, one L1 lineage in monocots was shown to possess a conserved 3'-end sequence with a solid RNA structure also observed in maize and sorghum SINEs. This strongly suggests that plant LINEs require a particular 3'-end sequence during initiation of reverse transcription. As one L1-clade LINE was also found to share the 3'-end sequence with a SINE in a green algal genome, I propose that the ancestral L1-clade LINE in the common ancestor of green plants may have recognized the specific RNA template, with stringent recognition then becoming relaxed during the course of plant evolution.

Description: Plant artificial micro-RNAs (amiRs) have been engineered to target viral genomes and induce their degradation. However, the exceptional evolutionary plasticity of RNA viruses threatens the durability of the resistance conferred by these amiRs. It has recently been shown that viral populations not experiencing strong selective pressure from an antiviral amiR may already contain enough genetic variability in the target sequence to escape plant resistance in an almost deterministic manner. Furthermore, it has also been shown that viral populations exposed to subinhibitory concentrations of the antiviral amiR speed up this process. In this article, we have characterized the molecular evolutionary dynamics of an amiR target sequence in a viral genome under both conditions. The use of Illumina ultradeep sequencing has allowed us to identify virus sequence variants at frequencies as low as 2 x 10 –6 and to track their variation in time before and after the viral population was able of successfully infecting plants fully resistant to the ancestral virus. We found that every site in the amiR-target sequence of the viral genome presented variation and that the variant that eventually broke resistance was sampled among the many coexisting ones. In this system, viral evolution in fully susceptible plants results from an equilibrium between mutation and genetic drift, whereas evolution in partially resistant plants originates from more complex dynamics involving mutation, selection, and drift.

Description: Molecular evolutionary rate estimates have been shown to depend on the time period over which they are estimated. Factors such as demographic processes, calibration errors, purifying selection, and the heterogeneity of substitution rates among sites (RHAS) are known to affect the accuracy with which rates of evolution are estimated. We use mathematical modeling and Bayesian analyses of simulated sequence alignments to explore how mutational hotspots can lead to time-dependent rate estimates. Mathematical modeling shows that underestimation of molecular rates over increasing time scales is inevitable when RHAS is ignored. Although a gamma distribution is commonly used to model RHAS, we show that when the actual RHAS deviates from a gamma-like distribution, rates can either be under- or overestimated in a time-dependent manner. Simulations performed under different scenarios of RHAS confirm the mathematical modeling and demonstrate the impacts of time-dependent rates on estimates of divergence times. Most notably, erroneous rate estimates can have narrow credibility intervals, leading to false confidence in biased estimates of rates, and node ages. Surprisingly, large errors in estimates of overall molecular rate do not necessarily generate large errors in divergence time estimates. Finally, we illustrate the correlation between time-dependent rate patterns and differential saturation between quickly and slowly evolving sites. Our results suggest that data partitioning or simple nonparametric mixture models of RHAS significantly improve the accuracy with which node ages and substitution rates can be estimated.

Description: Multiple visual pigments, prerequisites for color vision, are found in arthropods, but the evolutionary origin of their diversity remains obscure. In this study, we explore the opsin genes in five distantly related species of Onychophora, using deep transcriptome sequencing and screening approaches. Surprisingly, our data reveal the presence of only one opsin gene (onychopsin) in each onychophoran species, and our behavioral experiments indicate a maximum sensitivity of onychopsin to blue–green light. In our phylogenetic analyses, the onychopsins represent the sister group to the monophyletic clade of visual r-opsins of arthropods. These results concur with phylogenomic support for the sister-group status of the Onychophora and Arthropoda and provide evidence for monochromatic vision in velvet worms and in the last common ancestor of Onychophora and Arthropoda. We conclude that the diversification of visual pigments and color vision evolved in arthropods, along with the evolution of compound eyes—one of the most sophisticated visual systems known.

Description: Cis- regulatory DNA has been suspected to play a preeminent role in adaptive evolution, but understanding the role of cis -regulatory mutations in gene expression divergence first requires an accurate analysis of the functional differences associated with these regions. We analyzed allele-specific expression (ASE) in leaf and floral tissues of F1 interspecific hybrids generated between the two closely related species Arabidopsis thaliana and A. lyrata with a whole-genome SNP (single nucleotide polymorphism) tiling array. We observed 2,205 genes showing ASE pattern in at least one tissue. Nearly 90% of genes displaying ASE preferentially expressed the allele of A. lyrata . Genome-wide comparison of sequence divergence revealed that genes displaying ASE had a higher ratio of nonsynonymous to synonymous substitutions in coding regions. We further observe that the epigenetic landscape of histone methylation in A. thaliana genome associate with ASE. The asymmetry in the direction of allele-specific expression suggests interspecific differences in the efficiency of gene silencing in F1 hybrids.

Description: The olfactory receptor (OR) gene family is the largest gene family found in mammalian genomes. It is known to evolve through a birth-and-death process. Here, we characterized the sequences of 16 segregating OR pseudogenes in the samples of the wolf and the Chinese village dog (CVD) and compared them with the sequences from dogs of different breeds. Our results show that the segregating OR pseudogenes in breed dogs are under strong purifying selection, while evolving neutrally in the CVD, and show a more complicated pattern in the wolf. In the wolf, we found a trend to remove deleterious polymorphisms and accumulate nondeleterious polymorphisms. On the basis of protein structure of the ORs, we found that the distribution of different types of polymorphisms (synonymous, nonsynonymous, tolerated, and untolerated) varied greatly between the wolf and the breed dogs. In summary, our results suggest that different forms of selection have acted on the segregating OR pseudogenes in the CVD since domestication, breed dogs after breed formation, and ancestral wolf population, which has driven the evolution of these genes in different directions.

Description: Although gene duplications occur at a higher rate, only a small fraction of these are retained. The position of a gene’s encoded product in the protein–protein interaction network has recently emerged as a determining factor of gene duplicability. However, the direction of the relationship between network centrality and duplicability is not universal: In Escherichia coli , yeast, fly, and worm, duplicated genes more often act at the periphery of the network, whereas in humans, such genes tend to occupy the most central positions. Herein, we have inferred duplication events that took place in the different branches of the primate phylogeny. In agreement with previous observations, we found that duplications generally affected the most central network genes, which is presumably the process that has most influenced the trend in humans. However, the opposite trend—that is, duplication being more common in genes whose encoded products are peripheral in the network—is observed for three recent branches, including, quite counterintuitively, the external branch leading to humans. This indicates a shift in the relationship between centrality and duplicability during primate evolution. Furthermore, we found that genes encoding interacting proteins exhibit phylogenetic tree topologies that are more similar than expected for random pairs and that genes duplicated in a given branch of the phylogeny tend to interact with those that duplicated in the same lineage. These results indicate that duplication of a gene increases the likelihood of duplication of its interacting partners. Our observations indicate that the structure of the primate protein–protein interaction network affects gene duplicability in previously unrecognized ways.