Genome Sequence Variation Research Articles

Genetic architecture of quantitative traits in beef cattle revealed by genome wide association studies of imputed whole genome sequence variants: II: carcass merit traits

BackgroundGenome wide association studies (GWAS) were conducted on 7,853,211 imputed whole genome sequence variants in a population of 3354 to 3984 animals from multiple beef cattle breeds for five carcass merit traits including hot carcass weight (HCW), average backfat thickness (AFAT), rib eye area (REA), lean meat yield (LMY) and carcass marbling score (CMAR). Based on the GWAS results, genetic architectures of the carcass merit traits in beef cattle were elucidated.ResultsThe distributions of DNA variant allele substitution effects approximated a bell-shaped distribution for all the traits while the distribution of additive genetic variances explained by single DNA variants conformed to a scaled inverse chi-squared distribution to a greater extent. At a threshold of P-value < 10–5, 51, 33, 46, 40, and 38 lead DNA variants on multiple chromosomes were significantly associated with HCW, AFAT, REA, LMY, and CMAR, respectively. In addition, lead DNA variants with potentially large pleiotropic effects on HCW, AFAT, REA, and LMY were found on chromosome 6. On average, missense variants, 3’UTR variants, 5’UTR variants, and other regulatory region variants exhibited larger allele substitution effects on the traits in comparison to other functional classes. The amounts of additive genetic variance explained per DNA variant were smaller for intergenic and intron variants on all the traits whereas synonymous variants, missense variants, 3’UTR variants, 5’UTR variants, downstream and upstream gene variants, and other regulatory region variants captured a greater amount of additive genetic variance per sequence variant for one or more carcass merit traits investigated. In total, 26 enriched cellular and molecular functions were identified with lipid metabolisms, small molecular biochemistry, and carbohydrate metabolism being the most significant for the carcass merit traits.ConclusionsThe GWAS results have shown that the carcass merit traits are controlled by a few DNA variants with large effects and many DNA variants with small effects. Nucleotide polymorphisms in regulatory, synonymous, and missense functional classes have relatively larger impacts per sequence variant on the variation of carcass merit traits. The genetic architecture as revealed by the GWAS will improve our understanding on genetic controls of carcass merit traits in beef cattle.

Genetic architecture of quantitative traits in beef cattle revealed by genome wide association studies of imputed whole genome sequence variants: I: feed efficiency and component traits

BackgroundGenome wide association studies (GWAS) on residual feed intake (RFI) and its component traits including daily dry matter intake (DMI), average daily gain (ADG), and metabolic body weight (MWT) were conducted in a population of 7573 animals from multiple beef cattle breeds based on 7,853,211 imputed whole genome sequence variants. The GWAS results were used to elucidate genetic architectures of the feed efficiency related traits in beef cattle.ResultsThe DNA variant allele substitution effects approximated a bell-shaped distribution for all the traits while the distribution of additive genetic variances explained by single DNA variants followed a scaled inverse chi-squared distribution to a greater extent. With a threshold of P-value < 1.00E-05, 16, 72, 88, and 116 lead DNA variants on multiple chromosomes were significantly associated with RFI, DMI, ADG, and MWT, respectively. In addition, lead DNA variants with potentially large pleiotropic effects on DMI, ADG, and MWT were found on chromosomes 6, 14 and 20. On average, missense, 3’UTR, 5’UTR, and other regulatory region variants exhibited larger allele substitution effects in comparison to other functional classes. Intergenic and intron variants captured smaller proportions of additive genetic variance per DNA variant. Instead 3’UTR and synonymous variants explained a greater amount of genetic variance per DNA variant for all the traits examined while missense, 5’UTR and other regulatory region variants accounted for relatively more additive genetic variance per sequence variant for RFI and ADG, respectively. In total, 25 to 27 enriched cellular and molecular functions were identified with lipid metabolism and carbohydrate metabolism being the most significant for the feed efficiency traits.ConclusionsRFI is controlled by many DNA variants with relatively small effects whereas DMI, ADG, and MWT are influenced by a few DNA variants with large effects and many DNA variants with small effects. Nucleotide polymorphisms in regulatory region and synonymous functional classes play a more important role per sequence variant in determining variation of the feed efficiency traits. The genetic architecture as revealed by the GWAS of the imputed 7,853,211 DNA variants will improve our understanding on the genetic control of feed efficiency traits in beef cattle.

Assessing genomic diversity and signatures of selection in Original Braunvieh cattle using whole-genome sequencing data

BackgroundAutochthonous cattle breeds are an important source of genetic variation because they might carry alleles that enable them to adapt to local environment and food conditions. Original Braunvieh (OB) is a local cattle breed of Switzerland used for beef and milk production in alpine areas. Using whole-genome sequencing (WGS) data of 49 key ancestors, we characterize genomic diversity, genomic inbreeding, and signatures of selection in Swiss OB cattle at nucleotide resolution.ResultsWe annotated 15,722,811 SNPs and 1,580,878 Indels including 10,738 and 2763 missense deleterious and high impact variants, respectively, that were discovered in 49 OB key ancestors. Six Mendelian trait-associated variants that were previously detected in breeds other than OB, segregated in the sequenced key ancestors including variants causal for recessive xanthinuria and albinism. The average nucleotide diversity (1.6  × 10− 3) was higher in OB than many mainstream European cattle breeds. Accordingly, the average genomic inbreeding derived from runs of homozygosity (ROH) was relatively low (FROH = 0.14) in the 49 OB key ancestor animals. However, genomic inbreeding was higher in OB cattle of more recent generations (FROH = 0.16) due to a higher number of long (> 1 Mb) runs of homozygosity. Using two complementary approaches, composite likelihood ratio test and integrated haplotype score, we identified 95 and 162 genomic regions encompassing 136 and 157 protein-coding genes, respectively, that showed evidence (P < 0.005) of past and ongoing selection. These selection signals were enriched for quantitative trait loci related to beef traits including meat quality, feed efficiency and body weight and pathways related to blood coagulation, nervous and sensory stimulus.ConclusionsWe provide a comprehensive overview of sequence variation in Swiss OB cattle genomes. With WGS data, we observe higher genomic diversity and less inbreeding in OB than many European mainstream cattle breeds. Footprints of selection were detected in genomic regions that are possibly relevant for meat quality and adaptation to local environmental conditions. Considering that the population size is low and genomic inbreeding increased in the past generations, the implementation of optimal mating strategies seems warranted to maintain genetic diversity in the Swiss OB cattle population.

Analysis of preferred codon usage in the coronavirus N genes and their implications for genome evolution and vaccine design

The nucleocapsid (N) protein of a coronavirus plays a crucial role in virus assembly and in its RNA transcription. It is important to characterize a virus at the nucleotide level to discover the virus’s genomic sequence variations and similarities relative to other viruses that could have an impact on the functions of its genes and proteins. This entails a comprehensive and comparative analysis of the viral genomes of interest for preferred nucleotides, codon bias, nucleotide changes at the 3rd position (NT3s), synonymous codon usage and relative synonymous codon usage. In this study, the variations in the N proteins among 13 different coronaviruses (CoVs) were analysed at the nucleotide and amino acid levels in an attempt to reveal how these viruses adapt to their hosts relative to their preferred codon usage in the N genes. The results revealed that, overall, eighteen amino acids had different preferred codons and eight of these were over-biased. The N genes had a higher AT% over GC% and the values of their effective number of codons ranged from 40.43 to 53.85, indicating a slight codon bias. Neutrality plots and correlation analyses showed a very high level of GC3s/GC correlation in porcine epidemic diarrhea CoV (pedCoV), followed by Middle East respiratory syndrome-CoV (MERS CoV), porcine delta CoV (dCoV), bat CoV (bCoV) and feline CoV (fCoV) with r values 0.81, 0.68, -0.47, 0.98 and 0.58, respectively. These data implied a high rate of evolution of the CoV genomes and a strong influence of mutation on evolutionary selection in the CoV N genes. This type of genetic analysis would be useful for evaluating a virus’s host adaptation, evolution and is thus of value to vaccine design strategies.

Genome reconstruction and haplotype phasing using chromosome conformation capture methodologies.

Genomic analysis of individuals or organisms is predicated on the availability of high-quality reference and genotype information. With the rapidly dropping costs of high-throughput DNA sequencing, this is becoming readily available for diverse organisms and for increasingly large populations of individuals. Despite these advances, there are still aspects of genome sequencing that remain challenging for existing sequencing methods. This includes the generation of long-range contiguity during genome assembly, identification of structural variants in both germline and somatic tissues, the phasing of haplotypes in diploid organisms and the resolution of genome sequence for organisms derived from complex samples. These types of information are valuable for understanding the role of genome sequence and genetic variation on genome function, and numerous approaches have been developed to address them. Recently, chromosome conformation capture (3C) experiments, such as the Hi-C assay, have emerged as powerful tools to aid in these challenges for genome reconstruction. We will review the current use of Hi-C as a tool for aiding in genome sequencing, addressing the applications, strengths, limitations and potential future directions for the use of 3C data in genome analysis. We argue that unique features of Hi-C experiments make this data type a powerful tool to address challenges in genome sequencing, and that future integration of Hi-C data with alternative sequencing assays will facilitate the continuing revolution in genomic analysis and genome sequencing.

Acute exposure and chronic retention of aluminum in three vaccine schedules and effects of genetic and environmental variation

Like the mechanisms of action as adjuvants, the pharmacodynamics of injected forms of aluminum commonly used in vaccines are not well-characterized, particularly with respect to how differences in schedules impact accumulation and how factors such as genetics and environmental influences on detoxification influence clearance. Previous modeling efforts are based on very little empirical data, with the model by Priest based on whole-body clearance rates estimated from a study involving a single human subject. In this analysis, we explore the expected acute exposures and longer-term whole-body accumulation/clearance across three vaccination schedules: the current US Centers for Disease Control and Prevention (CDC) schedule, the current CDC schedule using low aluminum or no aluminum vaccines, and Dr. Paul Thomas’ “Vaccine Friendly Plan” schedule. We then study the effects of an implicit assumption of the Priest model on whether clearance dynamics from successive doses are influenced by the current level of aluminum or modeled by the assumption that a new dose has its own whole-body dynamics “reset” on the day of injection. We model two additional factors: variation (deficiency) in aluminum detoxification, and a factor added to the Priest equation to model the potential impact of aluminum itself on cellular and whole-body detoxification. These explorations are compared to a previously estimated pediatric dose limit (PDL) of whole-body aluminum exposure and provide a new statistic: %alumTox, the (expected) percentage of days (or weeks) an infant is in aluminum toxicity, reflecting chronic toxicity. We show that among three schedules, the CDC schedule results in the highest %alumTox regardless of model assumptions, and the Vaccine Friendly Plan schedule, which avoids >1 ACV per office visit results in the lowest (expected) %alumTox. These results are conservative, as the MSL is derived from data used by FDA to estimate safety of aluminum in adult humans. These results demonstrate high potential utility of modeling variation in patient responses to aluminum. More empirical data from individuals who are suspected of being intolerant of aluminum from vaccines, evidenced by high aluminum retention, neurodevelopmental disorders and/or a myriad of chronic illnesses would help answer questions on whether the model predictions can be used to estimate parameter values tied to genetic factors including genomic sequence variation and family history of chronic illnesses tied to aluminum exposure.

Host genome depletion to determine the evolution, genetic diversity and transmission patterns of full genome sequences of African swine fever genotype IX from Uganda

African swine fever (ASF) is a highly contagious viral haemorrhagic disease of domestic pigs, wild boar and feral swine. It is caused by the only DNA arbovirus, African swine fever virus (ASFV), which is the only member of the family Asfarviridae and is transmitted by the soft ticks in the genus Ornithodoros. ASF results in mortalities of up to 100% and exhibits different epidemiological patterns in Africa, Europe and Asia. This ASFV has recently spread to other countries in Europe and is also now present in Asia where it has spread like wildfire, resulting into high economic losses. The lack of a treatment or vaccine against this disease makes its control and prevention extremely difficult. This is partly due to the large size of the ASFV genome (170-190 kbp) and the many genes (187–190) that it encodes. There are currently 24 genotypes, with genotype IX predominating in Uganda and Kenya. Understanding the genetic characteristics of ASFV have focused mainly on partial fragments of selected genes, thus missing the opportunity to investigate the evolution, genetic diversity and patterns of transmission at the whole genomic level. Using the host-genome depletion approach we show an increase in the amount of ASFV DNA recovered directly from clinical material and generate genotype IX complete genomes on an NGS platform. We further show the inter-epidemic relationships of ASFV outbreaks in Uganda. The data contribute towards understanding the patterns and dynamics of transmission and maintenance of ASFV in a country in which ASF is endemic.

Genepanel.iobio - an easy to use web tool for generating disease- and phenotype-associated gene lists

When ordering genetic testing or triaging candidate variants in exome and genome sequencing studies, it is critical to generate and test a comprehensive list of candidate genes that succinctly describe the complete and objective phenotypic features of disease. Significant efforts have been made to curate gene:disease associations both in academic research and commercial genetic testing laboratory settings. However, many of these valuable resources exist as islands and must be used independently, generating static, single-resource gene:disease association lists. Here we describe genepanel.iobio (https://genepanel.iobio.io) an easy to use, free and open-source web tool for generating disease- and phenotype-associated gene lists from multiple gene:disease association resources, including the NCBI Genetic Testing Registry (GTR), Phenolyzer, and the Human Phenotype Ontology (HPO). We demonstrate the utility of genepanel.iobio by applying it to complex, rare and undiagnosed disease cases that had reached a diagnostic conclusion. We find that genepanel.iobio is able to correctly prioritize the gene containing the diagnostic variant in roughly half of these challenging cases. Importantly, each component resource contributed diagnostic value, showing the benefits of this aggregate approach. We expect genepanel.iobio will improve the ease and diagnostic value of generating gene:disease association lists for genetic test ordering and whole genome or exome sequencing variant prioritization.

Species Identification of Oaks (Quercus L., Fagaceae) from Gene to Genome

Species identification of oaks (Quercus) is always a challenge because many species exhibit variable phenotypes that overlap with other species. Oaks are notorious for interspecific hybridization and introgression, and complex speciation patterns involving incomplete lineage sorting. Therefore, accurately identifying Quercus species barcodes has been unsuccessful. In this study, we used chloroplast genome sequence data to identify molecular markers for oak species identification. Using next generation sequencing methods, we sequenced 14 chloroplast genomes of Quercus species in this study and added 10 additional chloroplast genome sequences from GenBank to develop a DNA barcode for oaks. Chloroplast genome sequence divergence was low. We identified four mutation hotspots as candidate Quercus DNA barcodes; two intergenic regions (matK-trnK-rps16 and trnR-atpA) were located in the large single copy region, and two coding regions (ndhF and ycf1b) were located in the small single copy region. The standard plant DNA barcode (rbcL and matK) had lower variability than that of the newly identified markers. Our data provide complete chloroplast genome sequences that improve the phylogenetic resolution and species level discrimination of Quercus. This study demonstrates that the complete chloroplast genome can substantially increase species discriminatory power and resolve phylogenetic relationships in plants.

Gene gain and loss push prokaryotes beyond the homologous recombination barrier and accelerate genome sequence divergence

Bacterial and archaeal evolution involve extensive gene gain and loss. Thus, phylogenetic trees of prokaryotes can be constructed both by traditional sequence-based methods (gene trees) and by comparison of gene compositions (genome trees). Comparing the branch lengths in gene and genome trees with identical topologies for 34 clusters of closely related bacterial and archaeal genomes, we show here that terminal branches of gene trees are systematically compressed compared to those of genome trees. Thus, sequence evolution is delayed compared to genome evolution by gene gain and loss. The extent of this delay differs widely among bacteria and archaea. Mathematical modeling shows that the divergence delay can result from sequence homogenization by homologous recombination. The model explains how homologous recombination maintains the cohesiveness of the core genome of a species while allowing extensive gene gain and loss within the accessory genome. Once evolving genomes become isolated by barriers impeding homologous recombination, gene and genome evolution processes settle into parallel trajectories, and genomes diverge, resulting in speciation.

Prediction of impacts of mutations on protein structure and interactions: SDM, a statistical approach, and mCSM, using machine learning.

Next-generation sequencing methods have not only allowed an understanding of genome sequence variation during the evolution of organisms but have also provided invaluable information about genetic variants in inherited disease and the emergence of resistance to drugs in cancers and infectious disease. A challenge is to distinguish mutations that are drivers of disease or drug resistance, from passengers that are neutral or even selectively advantageous to the organism. This requires an understanding of impacts of missense mutations in gene expression and regulation, and on the disruption of protein function by modulating protein stability or disturbing interactions with proteins, nucleic acids, small molecule ligands, and other biological molecules. Experimental approaches to understanding differences between wild-type and mutant proteins are most accurate but are also time-consuming and costly. Computational tools used to predict the impacts of mutations can provide useful information more quickly. Here, we focus on two widely used structure-based approaches, originally developed in the Blundell lab: site-directed mutator (SDM), a statistical approach to analyze amino acid substitutions, and mutation cutoff scanning matrix (mCSM), which uses graph-based signatures to represent the wild-type structural environment and machine learning to predict the effect of mutations on protein stability. Here, we describe DUET that uses machine learning to combine the two approaches. We discuss briefly the development of mCSM for understanding the impacts of mutations on interfaces with other proteins, nucleic acids, and ligands, and we exemplify the wide application of these approaches to understand human genetic disorders and drug resistance mutations relevant to cancer and mycobacterial infections. STATEMENT FOR A BROADER AUDIENCE: Genetic or somatic changes in genes can lead to mutations in human proteins, which give rise to genetic disorders or cancer, or to genes of pathogens leading to drug resistance. Computer software described here, using statistical approaches or machine learning, uses the information from genome sequencing of humans and pathogens, together with experimental or modeled 3D structures of gene products, the proteins, to predict impacts of mutations in genetic disease, cancer and drug resistance.

Target-Triggered Polymerization of Branched DNA Enables Enzyme-free and Fast Discrimination of Single-Base Changes.

SummarySingle-base changes lead to important biological and biomedical implications; however, the discrimination of single-base changes from normal DNA always remains a grand challenge. Herein we developed a DNA recognition and amplification system based on artificial branched DNA, namely, target-triggered polymerization (TTP), to realize enzyme-free and fast discrimination of single-base changes. Branched DNA as monomers rapidly polymerized into DNA nanospheres only with the trigger of specific DNA. Our TTP system worked reliably over a wide range of conditions. Remarkably, our TTP system was capable of discriminating base-changing DNA from normal DNA, including distinguishing 1–4 nucleotide changes and positions of single base, which was attributed to the significant amplification of small differences in hybridization thermodynamics and kinetics. We further proposed a theoretical method for calculating the hybridization probability of nucleic acids, which performed highly consistent with experimental results.

Identification of human genetic variants controlling circular RNA expression.

Circular RNAs (circRNAs) are abundant in eukaryotic transcriptomes and have been linked to various human disorders. However, understanding genetic control of circular RNA expression is in the early stages. Here we present the first integrated analysis of circRNAs and genome sequence variation from lymphoblastoid cell lines of the 1000 Genomes Project. We identified thousands of circRNAs in the RNA-seq data and show their association with local single-nucleotide polymorphic sites, referred to as circQTLs, which influence the circRNA transcript abundance. Strikingly, we found that circQTLs exist independently of eQTLs with most circQTLs having no effect on mRNA expression. Only a fraction of the polymorphic sites are shared and linked to both circRNA and mRNA expression with mostly similar effects on circular and linear RNA. A shared intronic QTL, rs55928920, of HMSD gene drives the circular and linear expression in opposite directions, potentially modulating circRNA levels at the expense of mRNA. Finally, circQTLs and eQTLs are largely independent and exist in separate linkage disequilibrium (LD) blocks with circQTLs highly enriched for functional genomic elements and regulatory regions. This study reveals a previously uncharacterized role of DNA sequence variation in human circular RNA regulation.

Diversity and Phylogenetic Relationships of Full Genome Sequences of Cassava Brown Streak Viruses in Kenya

Cassava brown streak disease is caused by cassava brown streak virus (CBSV) and Uganda cassava brown streak virus (UCBSV). Many of the CBSV and UCBSV diversity studies utilize partial coat protein sequences due to the unavailability of representative full genome sequences. Hence, there is little information on the diversity of cassava brown streak viruses in the rest of the genomes of the two species that are present in the farmers’ fields. The aim of this study was to determine Kenyan full CBSV and UCBSV genomes, and their sequence diversity and phylogenetic relationships within various genome and genome segments. Twenty four CBSVs positive samples tested by RT PCR from major cassava producing regions in Kenya were sequenced using Illumina MiSeq. Quality assessment of the output reads was done using the CLC Genomics 5.5.1 software programs. Genome assembly was done by de novo and reference guided assembly. Nucleotide sequence similarity of CBSV and UCBSV was determined. Phylogenetic relationships between CBSV and UCBSV were determined by performing the neighbour-joining analysis using MEGA 6.0 software. Six CBSV and 9 UCBSV genomes were generated from this study. The coat protein of the CBSV sequences had nucleotide sequence similarity of 92-100% while P1 and P3 gene had 75-100% and 76-100%, respectively. The coat protein of the UCBSV sequences had nucleotide sequence similarity of 91-99%. P1 and P3 gene had 83-100% and 86-99%, respectively. The phylogenetic analysis of full genomes revealed two distinct clusters one for UCBSV and another cluster for CBSV. Individual gene segments phylogenetic tree resembled that of the whole genome by clustering the nucleotide sequences into two clusters, one belonging to UCBSV and the other CBSV. The study revealed an average genome nucleotide diversity of 21.5% and 15.8% in CBSV and UCBSV respectively. A vast genetic diversity observed in CBSV and UCBSV in this study portends a lot of challenges in developing molecular diagnostic techniques as well as control strategies against the viruses.

Safety evaluation of myostatin-edited Meishan pigs by whole genome resequencing analyses

Genome editing technology can make specifically target genomic modifications, resulting in site specific DNA insertion, deletion or replacement in the genome of an organism. We have recently produced genetically engineered (GE) Meishan pigs containing a ZFN-edited myostatin (MSTN) loss-of-function mutation that leads to a clear “double muscle” phenotype as observed for Belgian cattle. In this study, whole genome resequencing was used as an approach to evaluate the safety risk, if any, associated with the introduction of a ZFN-edited myostatin (MSTN) loss-of-function mutation in a local pig breed, the Meishan pigs. The results of resequencing analyses show that the effective data from pigs of wild-type group and MSTN-edited GE group is greater than 99%. The 1× coverage rate is &amp;gt; 98%, and the 4× coverage rate is &amp;gt; 96%. The genetic variation on each chromosome is close to 1. From this whole genome resequencing study, our results demonstrated that 99.7% of single nucleotide polymorphisms (SNPs) are the same in the same genetic variation from both wild-type group and MSTN-edited GE group, implying genomic sequence variations are highly similar between the two groups of pigs.

Diversifying Evolution of the Ubiquitin-26S Proteasome System in Brassicaceae and Poaceae.

Genome amplification and sequence divergence provides raw materials to allow organismal adaptation. This is exemplified by the large expansion of the ubiquitin-26S proteasome system (UPS) in land plants, which primarily rely on intracellular signaling and biochemical metabolism to combat biotic and abiotic stresses. While a handful of functional genomic studies have demonstrated the adaptive role of the UPS in plant growth and development, many UPS members remain unknown. In this work, we applied a comparative genomic study to address the functional divergence of the UPS at a systematic level. We first used a closing-target-trimming annotation approach to identify most, if not all, UPS members in six species from each of two evolutionarily distant plant families, Brassicaceae and Poaceae. To reduce age-related errors, the two groups of species were selected based on their similar chronological order of speciation. Through size comparison, chronological expansion inference, evolutionary selection analyses, duplication mechanism prediction, and functional domain enrichment assays, we discovered significant diversities within the UPS, particularly between members from its three largest ubiquitin ligase gene families, the F-box (FBX), the Really Interesting New Gene (RING), and the Bric-a-Brac/Tramtrack/Broad Complex (BTB) families, between Brassicaceae and Poaceae. Uncovering independent Arabidopsis and Oryza genus–specific subclades of the 26S proteasome subunits from a comprehensive phylogenetic analysis further supported a diversifying evolutionary model of the UPS in these two genera, confirming its role in plant adaptation.

Genomic erosion and extensive horizontal gene transfer in gut-associated Acetobacteraceae

BackgroundSymbiotic relationships between animals and bacteria have profound impacts on the evolutionary trajectories of each partner. Animals and gut bacteria engage in a variety of relationships, occasionally persisting over evolutionary timescales. Ants are a diverse group of animals that engage in many types of associations with taxonomically distinct groups of bacterial associates. Here, we bring into culture and characterize two closely-related strains of gut associated Acetobacteraceae (AAB) of the red carpenter ant, Camponotus chromaiodes.ResultsGenome sequencing, assembly, and annotation of both strains delineate stark patterns of genomic erosion and sequence divergence in gut associated AAB. We found widespread horizontal gene transfer (HGT) in these bacterial associates and report elevated gene acquisition associated with energy production and conversion, amino acid and coenzyme transport and metabolism, defense mechanisms, and lysine export. Both strains have acquired the complete NADH-quinone oxidoreductase complex, plausibly from an Enterobacteriaceae origin, likely facilitating energy production under diverse conditions. Conservation of several lysine biosynthetic and salvage pathways and accumulation of lysine export genes via HGT implicate L-lysine supplementation by both strains as a potential functional benefit for the host. These trends are contrasted by genome-wide erosion of several amino acid biosynthetic pathways and pathways in central metabolism. We perform phylogenomic analyses on both strains as well as several free living and host associated AAB. Based on their monophyly and deep divergence from other AAB, these C. chromaiodes gut associates may represent a novel genus. Together, our results demonstrate how extensive horizontal transfer between gut associates along with genome-wide deletions leads to mosaic metabolic pathways. More broadly, these patterns demonstrate that HGT and genomic erosion shape metabolic capabilities of persistent gut associates and influence their genomic evolution.ConclusionsUsing comparative genomics, our study reveals substantial changes in genomic content in persistent associates of the insect gastrointestinal tract and provides evidence for the evolutionary pressures inherent to this environment. We describe patterns of genomic erosion and horizontal acquisition that result in mosaic metabolic pathways. Accordingly, the phylogenetic position of both strains of these associates form a divergent, monophyletic clade sister to gut associates of honey bees and more distantly to Gluconobacter.

Analysis of the genetic diversity of genome sequences of variants of apple hammerhead viroid

A dataset of approximately 85 apple hammerhead viroid (AHVd) genome sequences was analysed to determine genetic diversity. Sequences analysed included 18 new AHVd variants detected in a 2017/2018 survey, as well as suitable sequences deposited in GenBank. Two phylogenetically distinct clades were observed. A range of indel events was identified with most AHVd sequences ranging in size from 433 to 435 nt. None of these events appear to affect the predicted secondary structure of AHVd, which was well supported by a prevalence of unambiguous covariation events and or substitutions resulting in the conversion of canonical to wobble sequences, or vice versa. All AHVd variants analysed possess the 13 nt that are conserved in most hammerheads of viroids and viroid-like satellite RNA, except at positions 146 and 149. No significant recombination event was identified.

Essential role of inverted repeat in Epstein-Barr virus IR-1 in B cell transformation; geographical variation of the viral genome.

Many regions of the Epstein–Barr virus (EBV) genome, repeated and unique sequences, contribute to the geographical variation observed between strains. Here we use a large alignment of curated EBV genome sequences to identify major sites of variation in the genome of type 1 EBV strains; the CAO deletion in latent membrane protein 1 (LMP1) is the most frequent major indel present in the unique regions of EBV strains from various parts of the world. Principal component analysis was used to identify patterns of sequence variation and nucleotide positions in the sequences that can distinguish EBV from some different geographical regions. Viral genome sequence variation also affects interpretation of genetic content; known genes, origins of replication and gene expression control regions explain most of the viral genome but there are still a few sections of unknown function. One of these EBV genome regions contains a large inverted repeat sequence (invR) within the IR-1 major internal repeat array. We deleted this invR sequence and showed that this abolished the ability of the virus to transform human B cells into lymphoblastoid cell lines.This article is part of the theme issue ‘Silent cancer agents: multi-disciplinary modelling of human DNA oncoviruses’.

Copy-number variation in goat genome sequence: A comparative analysis of the different litter size trait groups

Copy number variation (CNV), as an important component of genomic structural variation (SV), plays essential roles in phenotypic variability, disease susceptibility and species evolution. To investigate whether critical CNVs exist in dairy goats with differing fecundity, we performed genome-wide sequencing of two populations of Laoshan dairy goats with large differences in litter size. After reference genome aligning, CNV calling, and annotation, we screened identified CNVs in the high-fecundity (HF) and low-fecundity (LF) groups to identify discrepant CNVs and their distribution within the genome. Prolactin-related protein 1 and 6 (PRP1 and PRP6), important factors regulating reproductive processes, were demonstrated to be duplicated in the HF group. In summary, based on the differences in CNVs between goats with differing litter sizes, it suggests CNVs may contribute to litter size in Laoshan dairy goats.