Local ecological divergence of two closely related stag beetles based on genetic,morphological, and environmental analyses

The process of phenotypic adaptation to the environments is widely recognized. However, comprehensive studies integrating phylogenetic, phenotypic, and ecological approaches to assess this process are scarce. Our study aims to assess whether local adaptation may explain intraspecific differentiation by quantifying multidimensional differences among populations in closely related lucanid species, Platycerus delicatulus and Platycerus kawadai, which are endemic saproxylic beetles in Japan. First, we determined intraspecific analysis units based on nuclear and mitochondrial gene analyses of Platycerus delicatulus and Platycerus kawadai under sympatric and allopatric conditions. Then, we compared differences in morphology and environmental niche between populations (analysis units) within species. We examined the relationship between morphology and environmental niche via geographic distance. P. kawadai was subdivided into the “No introgression” and “Introgression” populations based on mitochondrial COI gene – nuclear ITS region discordance. P. delicatulus was subdivided into “Allopatric” and “Sympatric” populations. Body length differed significantly among the populations of each species. For P. delicatulus, character displacement was suggested. For P. kawadai, the morphological difference was likely caused by geographic distance or genetic divergence rather than environmental differences. The finding showed that the observed mitochondrial–nuclear discordance is likely due to historical mitochondrial introgression following a range of expansion. Our results show that morphological variation among populations of P. delicatulus and P. kawadai reflects an ecological adaptation process based on interspecific interactions, geographic distance, or genetic divergence. Our results will deepen understanding of ecological specialization processes across the distribution and adaptation of species in natural systems.     

Bidirectional-Genetics Platform,a Dual-Purpose Mutagenesis Strategy for Filamentous Fungi

Rapidly increasing fungal genome sequences call for efficient ways of generating mutants to translate quickly gene sequences into their functions. A reverse genetic strategy via targeted gene replacement (TGR) has been inefficient for many filamentous fungi due to dominant production of undesirable ectopic transformants. Although large-scale random insertional mutagenesis via transformation (i.e., forward genetics) facilitates high-throughput uncovering of novel genes of interest, generating a huge number of transformants, which is necessary to ensure the likelihood of mutagenizing most genes, is time-consuming. We propose a new strategy, entitled the Bidirectional-Genetics (BiG) platform, which combines both forward and reverse genetic strategies by recycling ectopic transformants derived from TGR as a source for random insertional mutants. The BiG platform was evaluated using the rice blast fungus Magnaporthe oryzae as a model. Over 10% of >1,000 M. oryzae ectopic transformants, generated during disruption of specific genes, displayed abnormality in vegetative growth, pigmentation, and/or asexual reproduction. In this pool of putative mutants, we isolated insertional mutants with mutations in three genes involved in histidine biosynthesis (MoHIS5), vegetative growth (MoVPS74), or conidiophore formation (MoFRQ) (where “Mo” indicates “M. oryzae”), supporting the utility of this platform for systematic gene function studies.     

Genome Wide Association for Addiction: Replicated Results and Comparisons of Two Analytic Approaches

Background

Vulnerabilities to dependence on addictive substances are substantially heritable complex disorders whose underlying genetic architecture is likely to be polygenic, with modest contributions from variants in many individual genes. “Nontemplate” genome wide association (GWA) approaches can identity groups of chromosomal regions and genes that, taken together, are much more likely to contain allelic variants that alter vulnerability to substance dependence than expected by chance.

Methodology/Principal Findings

We report pooled “nontemplate” genome-wide association studies of two independent samples of substance dependent vs control research volunteers (n = 1620), one European-American and the other African-American using 1 million SNP (single nucleotide polymorphism) Affymetrix genotyping arrays. We assess convergence between results from these two samples using two related methods that seek clustering of nominally-positive results and assess significance levels with Monte Carlo and permutation approaches. Both “converge then cluster” and “cluster then converge” analyses document convergence between the results obtained from these two independent datasets in ways that are virtually never found by chance. The genes identified in this fashion are also identified by individually-genotyped dbGAP data that compare allele frequencies in cocaine dependent vs control individuals.

Conclusions/Significance

These overlapping results identify small chromosomal regions that are also identified by genome wide data from studies of other relevant samples to extents much greater than chance. These chromosomal regions contain more genes related to “cell adhesion” processes than expected by chance. They also contain a number of genes that encode potential targets for anti-addiction pharmacotherapeutics. “Nontemplate” GWA approaches that seek chromosomal regions in which nominally-positive associations are found in multiple independent samples are likely to complement classical, “template” GWA approaches in which “genome wide” levels of significance are sought for SNP data from single case vs control comparisons.     

KIAA1462, A Coronary Artery Disease Associated Gene,Is a Candidate Gene for Late Onset Alzheimer Disease in APOE Carriers

Alzheimer disease (AD) is a devastating neurodegenerative disease affecting more than five million Americans. In this study, we have used updated genetic linkage data from chromosome 10 in combination with expression data from serial analysis of gene expression to choose a new set of thirteen candidate genes for genetic analysis in late onset Alzheimer disease (LOAD). Results in this study identify the KIAA1462 locus as a candidate locus for LOAD in APOE4 carriers. Two genes exist at this locus, KIAA1462, a gene associated with coronary artery disease, and “rokimi”, encoding an untranslated spliced RNA The genetic architecture at this locus suggests that the gene product important in this association is either “rokimi”, or a different isoform of KIAA1462 than the isoform that is important in cardiovascular disease. Expression data suggests that isoform f of KIAA1462 is a more attractive candidate for association with LOAD in APOE4 carriers than “rokimi” which had no detectable expression in brain.     

Efficient Agrobacterium-mediated genetic transformation method using hypocotyl explants of radish (Raphanus sativus L.)

To investigate the gene function of radish (Raphanus sativus L.), several attempts have been made to generate genetically transformed radish. However, no efficient and relatively simple method for the genetic transformation of radish has been developed to date. In this study, we established an Agrobacterium-mediated genetic transformation method using the hypocotyl-derived explants of radish cultivar “Pirabikku”. Primarily based on the Brassica transformation procedure, we optimized it for radish transformation. Using this system, the transformation efficiency of radish hypocotyl explants by Agrobacterium tumefaciens strain GV3101 harboring pIG121-Hm was 13.3%. The copy number of transfer DNA integrated into the genome was either one or two in the four independent transgenic plants. Two of the four plants exhibited male sterility and did not produce self-pollinated seeds. Examination of the expression of the β-glucuronidase (GUS) gene in T₁ plants from fertile T₀ plants showed that the GUS genes were inherited. The improvement in the genetic transformation in this study might pave the way for accelerated molecular breeding and genetic analysis of radish.     

Discriminating between the Effects of Founding Events and Reproductive Mode on the Genetic Structure of Triops Populations (Branchiopoda: Notostraca)

Crustaceans that initially colonize a freshwater temporary pond can strongly bias the subsequent genetic composition of the population, causing nearby populations to be genetically distinct. In addition, these crustaceans have various reproductive modes that can influence genetic differentiation and diversity within and between populations. We report on two species of tadpole shrimp, Triops newberryi and Triops longicaudatus “short”, with different reproductive modes. Reproduction in the tadpole shrimp can occur clonally (parthenogenesis), with self fertilization (hermaphroditism), or through outcrossing of hermaphrodites with males (androdioecy). For all these reproductive modes, population genetic theory predicts decreased genetic diversity and increased population differentiation. Here we use mitochondrial control region (mtCR) sequences and nuclear microsatellite loci to determine if the difference in reproductive mode affects the high genetic structure typical of persistent founder effects. Previous authors indicated that T. newberryi is androdioecious because populations are composed of hermaphrodites and males, and T. longicaudatus “short” is hermaphroditic or parthenogenetic because males are absent. In our data, T. newberryi and T. longicaudatus “short” populations were highly structured genetically over short geographic distances for mtCR sequences and microsatellite loci (T. newberryi: Φ_ST = 0.644, F _ST = 0.252, respectively; T. l. “short”: invariant mtCR sequences, F _ST = 0.600). Differences between the two Triops species in a number of diversity measures were generally consistent with expectations from population genetic theory regarding reproductive mode; however, three of four comparisons were not statistically significant. We conclude the high genetic differentiation between populations is likely due to founder effects and results suggest both species are composed of selfing hermaphrodites with some level of outcrossing; the presence of males in T. newberryi does not appreciably reduce inbreeding. We cannot exclude the possibility that males in T. newberryi are non-reproductive individuals and the two species have the same mating system.     

Differences in Behavior and Activity Associated with a Poly(A) Expansion in the Dopamine Transporter in Belgian Malinois

In Belgian Malinois dogs, a 38-base pair variable number tandem repeat in the dopamine transporter gene (SLC6A3) is associated with behavior changes in Malinois. By additional sequencing in SLC6A3, we identified an intronic 12-nucleotide poly(A) insertion (“PolyA(22)”) before the terminal exon that was associated with seizure, “glazing over” behaviors, and episodic biting behaviors in a sample of 138 Malinois. We next investigated whether PolyA(22) was associated with 1) increased locomotor activity and 2) response to novelty. Using a sample of 22 Malinois and 25 dogs of other breeds, dogs’ activity was monitored in a novel and non-novel environment while wearing activity monitoring collars. All dogs were more active in novel compared with non-novel environments, and Malinois were more active overall than other breeds. There was an effect of PolyA(22) genotype on activity levels, and this effect appeared to underlie the difference detected between Malinois and other breeds. There was no effect of PolyA(22) genotype on the relative decrease in activity between novel and non-novel environments for either group or all dogs considered together. In addition to an association between PolyA(22) and owner reports of seizure, “glazing over” behaviors, and episodic biting behaviors, these findings support an effect of PolyA(22) on dopamine transporter function related to activity. Further investigation is required to confirm mechanistic effects of PolyA(22) on SLC6A3. The complex polygenic nature of behavior and the range of behaviors associated with this insertion predict that effects are likely also modified by additional genetic and environmental factors.     

“Candidatus Accumulibacter” Population Structure in Enhanced Biological Phosphorus Removal Sludges as Revealed by Polyphosphate Kinase Genes

We investigated the fine-scale population structure of the “Candidatus Accumulibacter” lineage in enhanced biological phosphorus removal (EBPR) systems using the polyphosphate kinase 1 gene (ppk1) as a genetic marker. We retrieved fragments of “Candidatus Accumulibacter” 16S rRNA and ppk1 genes from one laboratory-scale and several full-scale EBPR systems. Phylogenies reconstructed using 16S rRNA genes and ppk1 were largely congruent, with ppk1 granting higher phylogenetic resolution and clearer tree topology and thus serving as a better genetic marker than 16S rRNA for revealing population structure within the “Candidatus Accumulibacter” lineage. Sequences from at least five clades of “Candidatus Accumulibacter” were recovered by ppk1-targeted PCR, and subsequently, specific primer sets were designed to target the ppk1 gene for each clade. Quantitative real-time PCR (qPCR) assays using “Candidatus Accumulibacter”-specific 16S rRNA and “Candidatus Accumulibacter” clade-specific ppk1 primers were developed and conducted on three laboratory-scale and nine full-scale EBPR samples and two full-scale non-EBPR samples to determine the abundance of the total “Candidatus Accumulibacter” lineage and the relative distributions and abundances of the five “Candidatus Accumulibacter” clades. The qPCR-based estimation of the total “Candidatus Accumulibacter” fraction as a proportion of the bacterial community as measured using 16S rRNA genes was not significantly different from the estimation measured using ppk1, demonstrating the power of ppk1 as a genetic marker for detection of all currently defined “Candidatus Accumulibacter” clades. The relative distributions of “Candidatus Accumulibacter” clades varied among different EBPR systems and also temporally within a system. Our results suggest that the “Candidatus Accumulibacter” lineage is more diverse than previously realized and that different clades within the lineage are ecologically distinct.     

Evolution and Functional Characterisation of Melanopsins in a Deep-Sea Chimaera (Elephant Shark,Callorhinchus milii)

Non-visual photoreception in mammals is primarily mediated by two splice variants that derive from a single melanopsin (OPN4M) gene, whose expression is restricted to a subset of retinal ganglion cells. Physiologically, this sensory system regulates the photoentrainment of many biological rhythms, such as sleep via the melatonin endocrine system and pupil constriction. By contrast, melanopsin exists as two distinct lineages in non-mammals, opn4m and opn4x, and is broadly expressed in a wide range of tissue types, including the eye, brain, pineal gland and skin. Despite these findings, the evolution and function of melanopsin in early vertebrates are largely unknown. We, therefore, investigated the complement of opn4 classes present in the genome of a model deep-sea cartilaginous species, the elephant shark (Callorhinchus milii), as a representative vertebrate that resides at the base of the gnathostome (jawed vertebrate) lineage. We reveal that three melanopsin genes, opn4m1, opn4m2 and opn4x, are expressed in multiple tissues of the elephant shark. The two opn4m genes are likely to have arisen as a result of a lineage-specific duplication, whereas “long” and “short” splice variants are generated from a single opn4x gene. By using a heterologous expression system, we suggest that these genes encode functional photopigments that exhibit both “invertebrate-like” bistable and classical “vertebrate-like” monostable biochemical characteristics. We discuss the evolution and function of these melanopsin pigments within the context of the diverse photic and ecological environments inhabited by this chimaerid holocephalan, as well as the origin of non-visual sensory systems in early vertebrates.     

Learning a Prior on Regulatory Potential from eQTL Data

Genome-wide RNA expression data provide a detailed view of an organism's biological state; hence, a dataset measuring expression variation between genetically diverse individuals (eQTL data) may provide important insights into the genetics of complex traits. However, with data from a relatively small number of individuals, it is difficult to distinguish true causal polymorphisms from the large number of possibilities. The problem is particularly challenging in populations with significant linkage disequilibrium, where traits are often linked to large chromosomal regions containing many genes. Here, we present a novel method, Lirnet, that automatically learns a regulatory potential for each sequence polymorphism, estimating how likely it is to have a significant effect on gene expression. This regulatory potential is defined in terms of “regulatory features”—including the function of the gene and the conservation, type, and position of genetic polymorphisms—that are available for any organism. The extent to which the different features influence the regulatory potential is learned automatically, making Lirnet readily applicable to different datasets, organisms, and feature sets. We apply Lirnet both to the human HapMap eQTL dataset and to a yeast eQTL dataset and provide statistical and biological results demonstrating that Lirnet produces significantly better regulatory programs than other recent approaches. We demonstrate in the yeast data that Lirnet can correctly suggest a specific causal sequence variation within a large, linked chromosomal region. In one example, Lirnet uncovered a novel, experimentally validated connection between Puf3—a sequence-specific RNA binding protein—and P-bodies—cytoplasmic structures that regulate translation and RNA stability—as well as the particular causative polymorphism, a SNP in Mkt1, that induces the variation in the pathway.     

MET Gene Copy Number Predicts Worse Overall Survival in Patients with Non-Small Cell Lung Cancer (NSCLC); A Systematic Review and Meta-Analysis

Objectives

MET is a receptor present in the membrane of NSCLC cells and is known to promote cell proliferation, survival and migration. MET gene copy number is a common genetic alteration and inhibition o MET emerges as a promising targeted therapy in NSCLC. Here we aim to combine in a meta-analysis, data on the effect of high MET gene copy number on the overall survival of patients with resected NSCLC.

Methods

Two independent investigators applied parallel search strategies with the terms “MET AND lung cancer”, “MET AND NSCLC”, “MET gene copy number AND prognosis” in PubMed through January 2014. We selected the studies that investigated the association of MET gene copy number with survival, in patients who received surgery.

Results

Among 1096 titles that were identified in the initial search, we retrieved 9 studies on retrospective cohorts with adequate retrievable data regarding the prognostic impact of MET gene copy number on the survival of patients with NSCLC. Out of those, 6 used FISH and the remaining 3 used RT PCR to assess the MET gene copy number in the primary tumor. We calculated the I² statistic to assess heterogeneity (I² = 72%). MET gene copy number predicted worse overall survival when all studies were combined in a random effects model (HR = 1.78, 95% CI 1.22–2.60). When only the studies that had at least 50% of adenocarcinoma patients in their populations were included, the effect was significant (five studies, HR 1.55, 95% CI 1.23–1.94). This was not true when we included only the studies with no more than 50% of the patients having adenocarcinoma histology (four studies HR 2.18, 95% CI 0.97–4.90).

Conclusions

Higher MET gene copy number in the primary tumor at the time of diagnosis predicts worse outcome in patients with NSCLC. This prognostic impact may be adenocarcinoma histology specific.     

Higher Vulnerability and Stress Sensitivity of Neuronal Precursor Cells Carrying an Alpha-Synuclein Gene Triplication

Parkinson disease (PD) is a multi-factorial neurodegenerative disorder with loss of dopaminergic neurons in the substantia nigra and characteristic intracellular inclusions, called Lewy bodies. Genetic predisposition, such as point mutations and copy number variants of the SNCA gene locus can cause very similar PD-like neurodegeneration. The impact of altered α-synuclein protein expression on integrity and developmental potential of neuronal stem cells is largely unexplored, but may have wide ranging implications for PD manifestation and disease progression. Here, we investigated if induced pluripotent stem cell-derived neuronal precursor cells (NPCs) from a patient with Parkinson’s disease carrying a genomic triplication of the SNCA gene (SNCA-Tri). Our goal was to determine if these cells these neuronal precursor cells already display pathological changes and impaired cellular function that would likely predispose them when differentiated to neurodegeneration. To achieve this aim, we assessed viability and cellular physiology in human SNCA-Tri NPCs both under normal and environmentally stressed conditions to model in vitro gene-environment interactions which may play a role in the initiation and progression of PD. Human SNCA-Tri NPCs displayed overall normal cellular and mitochondrial morphology, but showed substantial changes in growth, viability, cellular energy metabolism and stress resistance especially when challenged by starvation or toxicant challenge. Knockdown of α-synuclein in the SNCA-Tri NPCs by stably expressed short hairpin RNA (shRNA) resulted in reversal of the observed phenotypic changes. These data show for the first time that genetic alterations such as the SNCA gene triplication set the stage for decreased developmental fitness, accelerated aging, and increased neuronal cell loss. The observation of this “stem cell pathology” could have a great impact on both quality and quantity of neuronal networks and could provide a powerful new tool for development of neuroprotective strategies for PD.     

The Effectiveness of Three Regions in Mitochondrial Genome for Aphid DNA Barcoding: A Case in Lachininae

Background

The mitochondrial gene COI has been widely used by taxonomists as a standard DNA barcode sequence for the identification of many animal species. However, the COI region is of limited use for identifying certain species and is not efficiently amplified by PCR in all animal taxa. To evaluate the utility of COI as a DNA barcode and to identify other barcode genes, we chose the aphid subfamily Lachninae (Hemiptera: Aphididae) as the focus of our study. We compared the results obtained using COI with two other mitochondrial genes, COII and Cytb. In addition, we propose a new method to improve the efficiency of species identification using DNA barcoding.

Methodology/Principal Findings

Three mitochondrial genes (COI, COII and Cytb) were sequenced and were used in the identification of over 80 species of Lachninae. The COI and COII genes demonstrated a greater PCR amplification efficiency than Cytb. Species identification using COII sequences had a higher frequency of success (96.9% in “best match” and 90.8% in “best close match”) and yielded lower intra- and higher interspecific genetic divergence values than the other two markers. The use of “tag barcodes” is a new approach that involves attaching a species-specific tag to the standard DNA barcode. With this method, the “barcoding overlap” can be nearly eliminated. As a result, we were able to increase the identification success rate from 83.9% to 95.2% by using COI and the “best close match” technique.

Conclusions/Significance

A COII-based identification system should be more effective in identifying lachnine species than COI or Cytb. However, the Cytb gene is an effective marker for the study of aphid population genetics due to its high sequence diversity. Furthermore, the use of “tag barcodes” can improve the accuracy of DNA barcoding identification by reducing or removing the overlap between intra- and inter-specific genetic divergence values.     

Genetic Dissection of the Drosophila melanogaster Female Head Transcriptome Reveals Widespread Allelic Heterogeneity

Modern genetic mapping is plagued by the “missing heritability” problem, which refers to the discordance between the estimated heritabilities of quantitative traits and the variance accounted for by mapped causative variants. One major potential explanation for the missing heritability is allelic heterogeneity, in which there are multiple causative variants at each causative gene with only a fraction having been identified. The majority of genome-wide association studies (GWAS) implicitly assume that a single SNP can explain all the variance for a causative locus. However, if allelic heterogeneity is prevalent, a substantial amount of genetic variance will remain unexplained. In this paper, we take a haplotype-based mapping approach and quantify the number of alleles segregating at each locus using a large set of 7922 eQTL contributing to regulatory variation in the Drosophila melanogaster female head. Not only does this study provide a comprehensive eQTL map for a major community genetic resource, the Drosophila Synthetic Population Resource, but it also provides a direct test of the allelic heterogeneity hypothesis. We find that 95% of cis-eQTLs and 78% of trans-eQTLs are due to multiple alleles, demonstrating that allelic heterogeneity is widespread in Drosophila eQTL. Allelic heterogeneity likely contributes significantly to the missing heritability problem common in GWAS studies.     

Differential Distribution of Genes Encoding the Virulence Factor Trans-Sialidase along Trypanosoma cruzi Discrete Typing Units

Trypanosoma cruzi the agent of Chagas disease is a monophyletic but heterogeneous group conformed by several Discrete Typing Units (DTUs) named TcI to TcVI characterized by genetic markers. The trans-sialidase (TS) is a virulence factor involved in cell invasion and pathogenesis that is differentially expressed in aggressive and less virulent parasite stocks. Genes encoding TS-related proteins are included in a large family divided in several groups but only one of them contains TS genes. Two closely related genes differing in a T/C transition encode the enzymatically active TS (aTS) and a lectin-like TS (iTS). We quantified the aTS/iTS genes from TcII and TcVI aggressive and TcI low virulent strains and found variable aTS number (1–32) per haploid genome. In spite of being low TS enzyme-expressers, TcI strains carry 28–32 aTS gene copies. The intriguing absence of iTS genes in TcI strains together with the presence of aTS/iTS in TcII and TcVI strains (virulent) were observed. Moreover, after sequencing aTS/iTS from 38 isolates collected along the Americas encompassing all DTUs, the persistent absence of the iTS gene in TcI, TcIII and TcIV was found. In addition, the sequence clustering together with T/C transition analysis correlated to DTUs of T. cruzi. The consistence of TS results with both evolutionary genome models proposed for T. cruzi, namely the “Two Hybridization” and the “Three Ancestor” was discussed and reviewed to fit present findings. Parasite stocks to attempt genetic KO or to assay the involvement of iTS in parasite biology and virulence are finally available.     

Genome-wide association scan shows genetic variants in the FTO gene are associated with obesity-related traits

The obesity epidemic is responsible for a substantial economic burden in developed countries and is a major risk factor for type 2 diabetes and cardiovascular disease. The disease is the result not only of several environmental risk factors, but also of genetic predisposition. To take advantage of recent advances in gene-mapping technology, we executed a genome-wide association scan to identify genetic variants associated with obesity-related quantitative traits in the genetically isolated population of Sardinia. Initial analysis suggested that several SNPs in the FTO and PFKP genes were associated with increased BMI, hip circumference, and weight. Within the FTO gene, rs9930506 showed the strongest association with BMI (p = 8.6 ×10⁻⁷), hip circumference (p = 3.4 × 10⁻⁸), and weight (p = 9.1 × 10⁻⁷). In Sardinia, homozygotes for the rare “G” allele of this SNP (minor allele frequency = 0.46) were 1.3 BMI units heavier than homozygotes for the common “A” allele. Within the PFKP gene, rs6602024 showed very strong association with BMI (p = 4.9 × 10⁻⁶). Homozygotes for the rare “A” allele of this SNP (minor allele frequency = 0.12) were 1.8 BMI units heavier than homozygotes for the common “G” allele. To replicate our findings, we genotyped these two SNPs in the GenNet study. In European Americans (N = 1,496) and in Hispanic Americans (N = 839), we replicated significant association between rs9930506 in the FTO gene and BMI (p-value for meta-analysis of European American and Hispanic American follow-up samples, p = 0.001), weight (p = 0.001), and hip circumference (p = 0.0005). We did not replicate association between rs6602024 and obesity-related traits in the GenNet sample, although we found that in European Americans, Hispanic Americans, and African Americans, homozygotes for the rare “A” allele were, on average, 1.0–3.0 BMI units heavier than homozygotes for the more common “G” allele. In summary, we have completed a whole genome–association scan for three obesity-related quantitative traits and report that common genetic variants in the FTO gene are associated with substantial changes in BMI, hip circumference, and body weight. These changes could have a significant impact on the risk of obesity-related morbidity in the general population.