The status of the human gene catalogue.

Scientists have been trying to identify every gene in the human genome since the initial draft was published in 2001. In the years since, much progress has been made in identifying protein-coding genes, currently estimated to number fewer than 20,000, with an ever-expanding number of distinct protein-coding isoforms. Here we review the status of the human gene catalogue and the efforts to complete it in recent years. Beside the ongoing annotation of protein-coding genes, their isoforms and pseudogenes, the invention of high-throughput RNA sequencing and other technological breakthroughs have led to a rapid growth in the number of reported non-coding RNA genes. For most of these non-coding RNAs, the functional relevance is currently unclear; we look at recent advances that offer paths forward to identifying their functions and towards eventually completing the human gene catalogue. Finally, we examine the need for a universal annotation standard that includes all medically significant genes and maintains their relationships with different reference genomes for the use of the human gene catalogue in clinical settings.

Prey Shifts Drive Venom Evolution in Cone Snails.

Venom systems are complex traits that have independently emerged multiple times in diverse plant and animal phyla. Within each venomous lineage there typically exists interspecific variation in venom composition where several factors have been proposed as drivers of variation, including phylogeny and diet. Understanding these factors is of broad biological interest and has implications for the development of anti-venom therapies and venom-based drug discovery. Because of their high species richness and the presence of several major evolutionary prey shifts, venomous marine cone snails (genus Conus) provide an ideal system to investigate drivers of interspecific venom variation. Here, by analyzing the venom gland expression profiles of ∼3,000 toxin genes from 42 species of cone snail, we elucidate the role of prey-specific selection pressures in shaping venom variation. By analyzing overall venom composition and individual toxin structures, we demonstrate that the shifts from vermivory to piscivory in Conus are complemented by distinct changes in venom composition independent of phylogeny. In vivo injections of venom from piscivorous cone snails in fish further showed a higher potency compared to venom of non-piscivores demonstrating a selective advantage. Together, our findings provide compelling evidence for the role of prey shifts in directing the venom composition of cone snails and expand our understanding of the mechanisms of venom variation and diversification.

The genomic basis of evolutionary differentiation among honey bees.

In contrast to the western honey bee, Apis mellifera, other honey bee species have been largely neglected despite their importance and diversity. The genetic basis of the evolutionary diversification of honey bees remains largely unknown. Here, we provide a genome-wide comparison of three honey bee species, each representing one of the three subgenera of honey bees, namely the dwarf (Apis florea), giant (A. dorsata), and cavity-nesting (A. mellifera) honey bees with bumblebees as an outgroup. Our analyses resolve the phylogeny of honey bees with the dwarf honey bees diverging first. We find that evolution of increased eusocial complexity in Apis proceeds via increases in the complexity of gene regulation, which is in agreement with previous studies. However, this process seems to be related to pathways other than transcriptional control. Positive selection patterns across Apis reveal a trade-off between maintaining genome stability and generating genetic diversity, with a rapidly evolving piRNA pathway leading to genomes depleted of transposable elements, and a rapidly evolving DNA repair pathway associated with high recombination rates in all Apis species. Diversification within Apis is accompanied by positive selection in several genes whose putative functions present candidate mechanisms for lineage-specific adaptations, such as migration, immunity, and nesting behavior.

Neurite outgrowth deficits caused by rare PLXNB1 mutation in pediatric bipolar disorder.

Pediatric bipolar disorder (PBD) is a severe mood dysregulation condition that affects 0.5-1% of children and teens in the United States. It is associated with recurrent episodes of mania and depression and an increased risk of suicidality. However, the genetics and neuropathology of PBD are largely unknown. Here, we used a combinatorial family-based approach to characterize cellular, molecular, genetic, and network-level deficits associated with PBD. We recruited a PBD patient and three unaffected family members from a family with a history of psychiatric illnesses. Using resting-state functional magnetic resonance imaging (rs-fMRI), we detected altered resting-state functional connectivity in the patient as compared to an unaffected sibling. Using transcriptomic profiling of patient and control induced pluripotent stem cell (iPSC)-derived telencephalic organoids, we found aberrant signaling in the molecular pathways related to neurite outgrowth. We corroborated the presence of neurite outgrowth deficits in patient iPSC-derived cortical neurons and identified a rare homozygous loss-of-function PLXNB1 variant (c.1360C>C; p.Ser454Arg) responsible for the deficits in the patient. Expression of wild-type PLXNB1, but not the variant, rescued neurite outgrowth in patient neurons, and expression of the variant caused the neurite outgrowth deficits in cortical neurons from PlxnB1 knockout mice. These results indicate that dysregulated PLXNB1 signaling may contribute to an increased risk of PBD and other mood dysregulation-related disorders by disrupting neurite outgrowth and functional brain connectivity. Overall, this study established and validated a novel family-based combinatorial approach for studying cellular and molecular deficits in psychiatric disorders and identified dysfunctional PLXNB1 signaling and neurite outgrowth as potential risk factors for PBD.

Identification of sodium channel toxins from marine cone snails of the subgenera Textilia and Afonsoconus.

Voltage-gated sodium (NaV) channels are transmembrane proteins that play a critical role in electrical signaling in the nervous system and other excitable tissues. µ-Conotoxins are peptide toxins from the venoms of marine cone snails (genus Conus) that block NaV channels with nanomolar potency. Most species of the subgenera Textilia and Afonsoconus are difficult to acquire; therefore, their venoms have yet to be comprehensively interrogated for µ-conotoxins. The goal of this study was to find new µ-conotoxins from species of the subgenera Textilia and Afonsoconus and investigate their selectivity at human NaV channels. Using RNA-seq of the venom gland of Conus (Textilia) bullatus, we identified 12 µ-conotoxin (or µ-conotoxin-like) sequences. Based on these sequences we designed primers which we used to identify additional µ-conotoxin sequences from DNA extracted from historical specimens of species from Textilia and Afonsoconus. We synthesized six of these µ-conotoxins and tested their activity on human NaV1.1-NaV1.8. Five of the six synthetic peptides were potent blockers of human NaV channels. Of these, two peptides (BuIIIB and BuIIIE) were potent blockers of hNaV1.3. Three of the peptides (BuIIIB, BuIIIE and AdIIIA) had submicromolar activity at hNaV1.7. This study serves as an example of the identification of new peptide toxins from historical DNA and provides new insights into structure-activity relationships of µ-conotoxins with activity at hNaV1.3 and hNaV1.7.

Settling the score: variant prioritization and Mendelian disease.

When investigating Mendelian disease using exome or genome sequencing, distinguishing disease-causing genetic variants from the multitude of candidate variants is a complex, multidimensional task. Many prioritization tools and online interpretation resources exist, and professional organizations have offered clinical guidelines for review and return of prioritization results. In this Review, we describe the strengths and weaknesses of widely used computational approaches, explain their roles in the diagnostic and discovery process and discuss how they can inform (and misinform) expert reviewers. We place variant prioritization in the wider context of gene prioritization, burden testing and genotype-phenotype association, and we discuss opportunities and challenges introduced by whole-genome sequencing.

A copy number variant is associated with a spectrum of pigmentation patterns in the rock pigeon (Columba livia).

Rock pigeons (Columba livia) display an extraordinary array of pigment pattern variation. One such pattern, Almond, is characterized by a variegated patchwork of plumage colors that are distributed in an apparently random manner. Almond is a sex-linked, semi-dominant trait controlled by the classical Stipper (St) locus. Heterozygous males (ZStZ+ sex chromosomes) and hemizygous Almond females (ZStW) are favored by breeders for their attractive plumage. In contrast, homozygous Almond males (ZStZSt) develop severe eye defects and often lack plumage pigmentation, suggesting that higher dosage of the mutant allele is deleterious. To determine the molecular basis of Almond, we compared the genomes of Almond pigeons to non-Almond pigeons and identified a candidate St locus on the Z chromosome. We found a copy number variant (CNV) within the differentiated region that captures complete or partial coding sequences of four genes, including the melanosome maturation gene Mlana. We did not find fixed coding changes in genes within the CNV, but all genes are misexpressed in regenerating feather bud collar cells of Almond birds. Notably, six other alleles at the St locus are associated with depigmentation phenotypes, and all exhibit expansion of the same CNV. Structural variation at St is linked to diversity in plumage pigmentation and gene expression, and thus provides a potential mode of rapid phenotypic evolution in pigeons.

Two Genomic Loci Control Three Eye Colors in the Domestic Pigeon (Columba livia).

The iris of the eye shows striking color variation across vertebrate species, and may play important roles in crypsis and communication. The domestic pigeon (Columba livia) has three common iris colors, orange, pearl (white), and bull (dark brown), segregating in a single species, thereby providing a unique opportunity to identify the genetic basis of iris coloration. We used comparative genomics and genetic mapping in laboratory crosses to identify two candidate genes that control variation in iris color in domestic pigeons. We identified a nonsense mutation in the solute carrier SLC2A11B that is shared among all pigeons with pearl eye color, and a locus associated with bull eye color that includes EDNRB2, a gene involved in neural crest migration and pigment development. However, bull eye is likely controlled by a heterogeneous collection of alleles across pigeon breeds. We also found that the EDNRB2 region is associated with regionalized plumage depigmentation (piebalding). Our study identifies two candidate genes for eye colors variation, and establishes a genetic link between iris and plumage color, two traits that vary widely in the evolution of birds and other vertebrates.

Elephant Genomes Reveal Accelerated Evolution in Mechanisms Underlying Disease Defenses.

Disease susceptibility and resistance are important factors for the conservation of endangered species, including elephants. We analyzed pathology data from 26 zoos and report that Asian elephants have increased neoplasia and malignancy prevalence compared with African bush elephants. This is consistent with observed higher susceptibility to tuberculosis and elephant endotheliotropic herpesvirus (EEHV) in Asian elephants. To investigate genetic mechanisms underlying disease resistance, including differential responses between species, among other elephant traits, we sequenced multiple elephant genomes. We report a draft assembly for an Asian elephant, and defined 862 and 1,017 conserved potential regulatory elements in Asian and African bush elephants, respectively. In the genomes of both elephant species, conserved elements were significantly enriched with genes differentially expressed between the species. In Asian elephants, these putative regulatory regions were involved in immunity pathways including tumor-necrosis factor, which plays an important role in EEHV response. Genomic sequences of African bush, forest, and Asian elephant genomes revealed extensive sequence conservation at TP53 retrogene loci across three species, which may be related to TP53 functionality in elephant cancer resistance. Positive selection scans revealed outlier genes related to additional elephant traits. Our study suggests that gene regulation plays an important role in the differential inflammatory response of Asian and African elephants, leading to increased infectious disease and cancer susceptibility in Asian elephants. These genomic discoveries can inform future functional and translational studies aimed at identifying effective treatment approaches for ill elephants, which may improve conservation.

Comprehensive Genomic Profiling of Cell-Free DNA in Men With Advanced Prostate Cancer: Differences in Genomic Landscape Based on Race.

Advanced prostate cancer (aPC) in Black men was reported to present with aggressive features and to be associated with poor prognosis. Herein, we compared the cell-free DNA (cfDNA) genomic landscape of aPC in Black vs White men. Patients (pts) with aPC from 6 academic institutions and available cfDNA comprehensive genomic profiling (CGP) were included. Association between mutated genes and race was evaluated using Barnard's test and a Probabilistic Graphical Model (PGM) machine learning approach. Analysis included 743 aPC pts (217 Black, 526 White) with available cfDNA CGP. The frequency of alterations in the androgen receptor gene was significantly higher in Black vs White men (55.3% vs 35% respectively, P < .001). Additionally, alterations in EGFR, MYC, FGFR1, and CTNNB1 were present at higher frequencies in Black men. PGM analysis and Barnard's test were concordant. Findings from the largest cohort of Black men with aPC undergoing cfDNA CGP may guide further drug development in these men.

Noncoding sequence variants define a novel regulatory element in the first intron of the N-acetylglutamate synthase gene.

N-acetylglutamate synthase deficiency is an autosomal recessive urea cycle disorder caused either by decreased expression of the NAGS gene or defective NAGS enzyme resulting in decreased production of N-acetylglutamate (NAG), an allosteric activator of carbamylphosphate synthetase 1 (CPS1). NAGSD is the only urea cycle disorder that can be effectively treated with a single drug, N-carbamylglutamate (NCG), a stable NAG analog, which activates CPS1 to restore ureagenesis. We describe three patients with NAGSD due to four novel noncoding sequence variants in the NAGS regulatory regions. All three patients had hyperammonemia that resolved upon treatment with NCG. Sequence variants NM_153006.2:c.427-222G>A and NM_153006.2:c.427-218A>C reside in the 547 bp-long first intron of NAGS and define a novel NAGS regulatory element that binds retinoic X receptor α. Sequence variants NC_000017.10:g.42078967A>T (NM_153006.2:c.-3065A>T) and NC_000017.10:g.42078934C>T (NM_153006.2:c.-3098C>T) reside in the NAGS enhancer, within known HNF1 and predicted glucocorticoid receptor binding sites, respectively. Reporter gene assays in HepG2 and HuH-7 cells demonstrated that all four substitutions could result in reduced expression of NAGS. These findings show that analyzing noncoding regions of NAGS and other urea cycle genes can reveal molecular causes of disease and identify novel regulators of ureagenesis.

Identification of Somatic Gene Signatures in Circulating Cell-Free DNA Associated with Disease Progression in Metastatic Prostate Cancer by a Novel Machine Learning Platform.

PURPOSE: Progression from metastatic castration-sensitive prostate cancer (mCSPC) to a castration-resistant (mCRPC) state heralds the lethal phenotype of prostate cancer. Identifying genomic alterations associated with mCRPC may help find new targets for drug development. In the majority of patients, obtaining a tumor biopsy is challenging because of the predominance of bone-only metastasis. In this study, we hypothesize that machine learning (ML) algorithms can identify clinically relevant patterns of genomic alterations (GAs) that distinguish mCRPC from mCSPC, as assessed by next-generation sequencing (NGS) of circulating cell-free DNA (cfDNA). EXPERIMENTAL DESIGN: Retrospective clinical data from men with metastatic prostate cancer were collected. Men with NGS of cfDNA performed at a Clinical Laboratory Improvement Amendments (CLIA)-certified laboratory at time of diagnosis of mCSPC or mCRPC were included. A combination of supervised and unsupervised ML algorithms was used to obtain biologically interpretable, potentially actionable insights into genomic signatures that distinguish mCRPC from mCSPC. RESULTS: GAs that distinguish patients with mCRPC (n = 187) from patients with mCSPC (n = 154) (positive predictive value = 94%, specificity = 91%) were identified using supervised ML algorithms. These GAs, primarily amplifications, corresponded to androgen receptor, Mitogen-activated protein kinase (MAPK) signaling, Phosphoinositide 3-kinase (PI3K) signaling, G1/S cell cycle, and receptor tyrosine kinases. We also identified recurrent patterns of gene- and pathway-level alterations associated with mCRPC by using Bayesian networks, an unsupervised machine learning algorithm. CONCLUSION: These results provide clinical evidence that progression from mCSPC to mCRPC is associated with stereotyped concomitant gain-of-function aberrations in these pathways. Furthermore, detection of these aberrations in cfDNA may overcome the challenges associated with obtaining tumor bone biopsies and allow contemporary investigation of combinatorial therapies that target these aberrations. IMPLICATIONS FOR PRACTICE: The progression from castration-sensitive to castration-resistant prostate cancer is characterized by worse prognosis and there is a pressing need for targeted drugs to prevent or delay this transition. This study used machine learning algorithms to examine the cell-free DNA of patients to identify alterations to specific pathways and genes associated with progression. Detection of these alterations in cell-free DNA may overcome the challenges associated with obtaining tumor bone biopsies and allow contemporary investigation of combinatorial therapies that target these aberrations.

Pulmonary Eosinophilic Granulomatosis with Polyangiitis Has IgG4 Plasma Cells and Immunoregulatory Features.

The immunologic mechanisms promoting eosinophilic granulomatosis with polyangiitis (EGPA) are unclear. To characterize the mechanisms underlying pulmonary EGPA, we examined and compared EGPA paraffin-embedded lung biopsies with normal lung biopsies, using immunostaining, RNA sequencing, and RT-PCR. The results revealed novel type 2 as well as immuneregulatory features. These features included basophils and increased mast cell contents; increased immunostaining for tumor necrosis factor ligand superfamily member 14; sparse mast cell degranulation; numerous forkhead box protein P3 (FoxP3)+ regulatory T cells and IgG4 plasma cells; and abundant arachidonate 15-lipoxygenase and 25-hydroxyvitamin D-1 α hydroxylase, mitochondrial. Significantly decreased 15-hydroxyprostaglandin dehydrogenase [NAD(+)], which degrades eicosanoids, was observed in EGPA samples. In addition, there was significantly increased mRNA for chemokine (C-C motif) ligands 18 and 13 and major collagen genes, IgG4-rich immune complexes coating alveolar macrophages, and increased immunostaining for phosphorylated mothers against decapentaplegic homolog 2/SMAD2, suggesting transforming growth factor-ß activation. These findings suggest a novel self-promoting mechanism of activation of alveolar macrophages by arachidonate 15-lipoxygenase-derived eicosanoids to express chemokines that recruit a combined type 2/immunoregulatory immune response, which produces these eicosanoids. These results suggest that the pulmonary EGPA immune response resembles the immune response to a tissue-invasive parasite infection.

Regulatory Divergence in Wound-Responsive Gene Expression between Domesticated and Wild Tomato.

The evolution of transcriptional regulatory mechanisms is central to how stress response and tolerance differ between species. However, it remains largely unknown how divergence in cis-regulatory sites and, subsequently, transcription factor (TF) binding specificity contribute to stress-responsive expression divergence, particularly between wild and domesticated species. By profiling wound-responsive gene transcriptomes in wild Solanum pennellii and domesticated S. lycopersicum, we found extensive wound response divergence and identified 493 S. lycopersicum and 278 S. pennellii putative cis-regulatory elements (pCREs) that were predictive of wound-responsive gene expression. Only 24-52% of these wound response pCREs (depending on wound response patterns) were consistently enriched in the putative promoter regions of wound-responsive genes across species. In addition, between these two species, their differences in pCRE site sequences were significantly and positively correlated with differences in wound-responsive gene expression. Furthermore, ∼11-39% of pCREs were specific to only one of the species and likely bound by TFs from different families. These findings indicate substantial regulatory divergence in these two plant species that diverged ∼3-7 million years ago. Our study provides insights into the mechanistic basis of how the transcriptional response to wounding is regulated and, importantly, the contribution of cis-regulatory components to variation in wound-responsive gene expression between a wild and a domesticated plant species.

Comprehensive Analysis of AR Alterations in Circulating Tumor DNA from Patients with Advanced Prostate Cancer.

BACKGROUND: Somatic alterations in circulating tumor DNA (ctDNA) may be associated with treatment response or prognosis in prostate cancer (PCa). The goal was to characterize androgen receptor gene (AR) amplifications and mutations detected in ctDNA from patients with PCa and to further understand the somatic genetic heterogeneity of advanced prostate cancer. PATIENTS AND METHODS: This study included a heterogeneous group of 892 patients with advanced PCa (predominantly castrate-resistant prostate cancer) with AR alterations detected in ctDNA that underwent next-generation sequencing of 54 to 73 genes via Guardant360 testing (Guardant Health, Inc., Redwood City, CA). Distribution and summary of AR alterations detected, the association of AR alterations with other genes, and a pathway analysis are reported. RESULTS: The median absolute plasma copy number of AR amplifications was 3.3 (range, 1.2-165.2). Many patients had multiple AR mutations; a total of 112 unique mutations were identified in AR, including L702H (25%), T878A (14%), H875Y (11%), W742C (8%), W742L (4%), F877L (2%), and T878S (2%). Other ctDNA gene alterations in the Guardant assays included TP53 (50%), MYC (34%), BRAF (32%), PIK3CA (29%), MET (25%), CDK6 (26%), EGFR (24%), FGFR1 (21%), and APC (12%). Many of these non-AR alterations are not tissue verified in other studies. AR amplification cosegregated with alterations in MYC (p < .001), BRAF (p < .001), PIK3CA (p < .001), MET (p < .001), CDK6 (p < .001), EGFR (p < .001), FGFR1 (p = .391), and more. Alterations in APC were significantly associated with mutations in AR (p < .001). CONCLUSION: Several AR alterations and concomitant non-AR alterations that associate with drug resistance were detected. These findings provide additional insights into the heterogeneity of advanced prostate cancer. IMPLICATIONS FOR PRACTICE: The goal was to characterize androgen receptor gene (AR) amplifications and mutations detected in circulating tumor DNA (ctDNA) from patients with prostate cancer in relation to non-AR gene alterations detected in the ctDNA landscape. The study included 892 patients with prostate cancer with AR alterations in ctDNA. AR alterations were significantly associated with other gene alterations detected in ctDNA. The common AR mutations found are linked to resistance to abiraterone, enzalutamide, or bicalutamide. Characterization of the circulating AR landscape and gene alterations provides potential additional insight into the somatic genetic heterogeneity of advanced prostate cancer.

Identification of Polycystic Kidney Disease 1 Like 1 Gene Variants in Children With Biliary Atresia Splenic Malformation Syndrome.

Biliary atresia (BA) is the most common cause of end-stage liver disease in children and the primary indication for pediatric liver transplantation, yet underlying etiologies remain unknown. Approximately 10% of infants affected by BA exhibit various laterality defects (heterotaxy) including splenic abnormalities and complex cardiac malformations-a distinctive subgroup commonly referred to as the biliary atresia splenic malformation (BASM) syndrome. We hypothesized that genetic factors linking laterality features with the etiopathogenesis of BA in BASM patients could be identified through whole-exome sequencing (WES) of an affected cohort. DNA specimens from 67 BASM subjects, including 58 patient-parent trios, from the National Institute of Diabetes and Digestive and Kidney Diseases-supported Childhood Liver Disease Research Network (ChiLDReN) underwent WES. Candidate gene variants derived from a prespecified set of 2,016 genes associated with ciliary dysgenesis and/or dysfunction or cholestasis were prioritized according to pathogenicity, population frequency, and mode of inheritance. Five BASM subjects harbored rare and potentially deleterious biallelic variants in polycystic kidney disease 1 like 1 (PKD1L1), a gene associated with ciliary calcium signaling and embryonic laterality determination in fish, mice, and humans. Heterozygous PKD1L1 variants were found in 3 additional subjects. Immunohistochemical analysis of liver from the one BASM subject available revealed decreased PKD1L1 expression in bile duct epithelium when compared to normal livers and livers affected by other noncholestatic diseases. Conclusion: WES identified biallelic and heterozygous PKD1L1 variants of interest in 8 BASM subjects from the ChiLDReN data set; the dual roles for PKD1L1 in laterality determination and ciliary function suggest that PKD1L1 is a biologically plausible, cholangiocyte-expressed candidate gene for the BASM syndrome.

XPAT: a toolkit to conduct cross-platform association studies with heterogeneous sequencing datasets.

High-throughput sequencing data are increasingly being made available to the research community for secondary analyses, providing new opportunities for large-scale association studies. However, heterogeneity in target capture and sequencing technologies often introduce strong technological stratification biases that overwhelm subtle signals of association in studies of complex traits. Here, we introduce the Cross-Platform Association Toolkit, XPAT, which provides a suite of tools designed to support and conduct large-scale association studies with heterogeneous sequencing datasets. XPAT includes tools to support cross-platform aware variant calling, quality control filtering, gene-based association testing and rare variant effect size estimation. To evaluate the performance of XPAT, we conducted case-control association studies for three diseases, including 783 breast cancer cases, 272 ovarian cancer cases, 205 Crohn disease cases and 3507 shared controls (including 1722 females) using sequencing data from multiple sources. XPAT greatly reduced Type I error inflation in the case-control analyses, while replicating many previously identified disease-gene associations. We also show that association tests conducted with XPAT using cross-platform data have comparable performance to tests using matched platform data. XPAT enables new association studies that combine existing sequencing datasets to identify genetic loci associated with common diseases and other complex traits.

Correction: Genome, Functional Gene Annotation, and Nuclear Transformation of the Heterokont Oleaginous Alga Nannochloropsis oceanica CCMP1779.

[This corrects the article DOI: 10.1371/journal.pgen.1003064.].

Predicting gene structure changes resulting from genetic variants via exon definition features.

Motivation: Genetic variation that disrupts gene function by altering gene splicing between individuals can substantially influence traits and disease. In those cases, accurately predicting the effects of genetic variation on splicing can be highly valuable for investigating the mechanisms underlying those traits and diseases. While methods have been developed to generate high quality computational predictions of gene structures in reference genomes, the same methods perform poorly when used to predict the potentially deleterious effects of genetic changes that alter gene splicing between individuals. Underlying that discrepancy in predictive ability are the common assumptions by reference gene finding algorithms that genes are conserved, well-formed and produce functional proteins. Results: We describe a probabilistic approach for predicting recent changes to gene structure that may or may not conserve function. The model is applicable to both coding and non-coding genes, and can be trained on existing gene annotations without requiring curated examples of aberrant splicing. We apply this model to the problem of predicting altered splicing patterns in the genomes of individual humans, and we demonstrate that performing gene-structure prediction without relying on conserved coding features is feasible. The model predicts an unexpected abundance of variants that create de novo splice sites, an observation supported by both simulations and empirical data from RNA-seq experiments. While these de novo splice variants are commonly misinterpreted by other tools as coding or non-coding variants of little or no effect, we find that in some cases they can have large effects on splicing activity and protein products and we propose that they may commonly act as cryptic factors in disease. Availability and implementation: The software is available from geneprediction.org/SGRF. Supplementary information: Supplementary information is available at Bioinformatics online.

Antimicrobial Functions of Lactoferrin Promote Genetic Conflicts in Ancient Primates and Modern Humans.

Lactoferrin is a multifunctional mammalian immunity protein that limits microbial growth through sequestration of nutrient iron. Additionally, lactoferrin possesses cationic protein domains that directly bind and inhibit diverse microbes. The implications for these dual functions on lactoferrin evolution and genetic conflicts with microbes remain unclear. Here we show that lactoferrin has been subject to recurrent episodes of positive selection during primate divergence predominately at antimicrobial peptide surfaces consistent with long-term antagonism by bacteria. An abundant lactoferrin polymorphism in human populations and Neanderthals also exhibits signatures of positive selection across primates, linking ancient host-microbe conflicts to modern human genetic variation. Rapidly evolving sites in lactoferrin further correspond to molecular interfaces with opportunistic bacterial pathogens causing meningitis, pneumonia, and sepsis. Because microbes actively target lactoferrin to acquire iron, we propose that the emergence of antimicrobial activity provided a pivotal mechanism of adaptation sparking evolutionary conflicts via acquisition of new protein functions.