Genomic divergence and adaptive convergence in Drosophila simulans from Evolution Canyon, Israel.

Biodiversity refugia formed by unique features of the Mediterranean arid landscape, such as the dramatic ecological contrast of "Evolution Canyon," provide a natural laboratory in which local adaptations to divergent microclimate conditions can be investigated. Significant insights have been provided by studies of Drosophila melanogaster diversifying along the thermal gradient in Evolution Canyon, but a comparative framework to survey adaptive convergence across sister species at the site has been lacking. To fill this void, we present an analysis of genomic polymorphism and evolutionary divergence of Drosophila simulans, a close relative of Drosophila melanogaster with which it co-occurs on both slopes of the canyon. Our results show even deeper interslope divergence in D. simulans than in D. melanogaster, with extensive signatures of selective sweeps present in flies from both slopes but enhanced in the population from the hotter and drier south-facing slope. Interslope divergence was enriched for genes related to electrochemical balance and transmembrane transport, likely in response to increased selection for dehydration resistance on the hotter slope. Both species shared genomic regions that underwent major selective sweeps, but the overall level of adaptive convergence was low, demonstrating no shortage of alternative genomic solutions to cope with the challenges of the microclimate contrast. Mobile elements were a major source of genetic polymorphism and divergence, affecting all parts of the genome, including coding sequences of mating behavior-related genes.

CAGm: a repository of germline microsatellite variations in the 1000 genomes project.

The human genome harbors an abundance of repetitive DNA; however, its function continues to be debated. Microsatellites-a class of short tandem repeat-are established as an important source of genetic variation. Array length variants are common among microsatellites and affect gene expression; but, efforts to understand the role and diversity of microsatellite variation has been hampered by several challenges. Without adequate depth, both long-read and short-read sequencing may not detect the variants present in a sample; additionally, large sample sizes are needed to reveal the degree of population-level polymorphism. To address these challenges we present the Comparative Analysis of Germline Microsatellites (CAGm): a database of germline microsatellites from 2529 individuals in the 1000 genomes project. A key novelty of CAGm is the ability to aggregate microsatellite variation by population, ethnicity (super population) and gender. The database provides advanced searching for microsatellites embedded in genes and functional elements. All data can be downloaded as Microsoft Excel spreadsheets. Two use-case scenarios are presented to demonstrate its utility: a mononucleotide (A) microsatellite at the BAT-26 locus and a dinucleotide (CA) microsatellite in the coding region of FGFRL1. CAGm is freely available at http://www.cagmdb.org/.

Estimating the number of genetic mutations (hits) required for carcinogenesis based on the distribution of somatic mutations.

Individual instances of cancer are primarily a result of a combination of a small number of genetic mutations (hits). Knowing the number of such mutations is a prerequisite for identifying specific combinations of carcinogenic mutations and understanding the etiology of cancer. We present a mathematical model for estimating the number of hits based on the distribution of somatic mutations. The model is fundamentally different from previous approaches, which are based on cancer incidence by age. Our somatic mutation based model is likely to be more robust than age-based models since it does not require knowing or accounting for the highly variable mutation rate, which can vary by over three orders of magnitude. In fact, we find that the number of somatic mutations at diagnosis is weakly correlated with age at cancer diagnosis, most likely due to the extreme variability in mutation rates between individuals. Comparing the distribution of somatic mutations predicted by our model to the actual distribution from 6904 tumor samples we estimate the number of hits required for carcinogenesis for 17 cancer types. We find that different cancer types exhibit distinct somatic mutational profiles corresponding to different numbers of hits. Why might different cancer types require different numbers of hits for carcinogenesis? The answer may provide insight into the unique etiology of different cancer types.

DNA sequencing of anatomy lab cadavers to provide hands-on precision medicine introduction to medical students.

BACKGROUND: Medical treatment informed by Precision Medicine is becoming a standard practice for many diseases, and patients are curious about the consequences of genomic variants in their genome. However, most medical students' understanding of Precision Medicine derives from classroom lectures. This format does little to foster an understanding for the potential and limitations of Precision Medicine. To close this gap, we implemented a hands-on Precision Medicine training program utilizing exome sequencing to prepare a clinical genetic report of cadavers studied in the anatomy lab. The program reinforces Precision Medicine related learning objectives for the Genetics curriculum. METHODS: Pre-embalmed blood samples and embalmed tissue were obtained from cadavers (donors) used in the anatomy lab. DNA was isolated and sequenced and illustrative genetic reports provided to the students. The reports were used to facilitate discussion with students on the implications of pathogenic genomic variants and the potential correlation of these variants in each "donor" with any anatomical anomalies identified during cadaver dissection. RESULTS: In 75% of cases, analysis of whole exome sequencing data identified a variant associated with increased risk for a disease/abnormal condition noted in the donor's cause of death or in the students' anatomical findings. This provided students with real-world examples of the potential relationship between genomic variants and disease risk. Our students also noted that diseases associated with 92% of the pathogenic variants identified were not related to the anatomical findings, demonstrating the limitations of Precision Medicine. CONCLUSION: With this study, we have established protocols and classroom procedures incorporating hands-on Precision Medicine training in the medical student curriculum and a template for other medical educators interested in enhancing their Precision Medicine training program. The program engaged students in discovering variants that were associated with the pathophysiology of the cadaver they were studying, which led to more exposure and understanding of the potential risks and benefits of genomic medicine.

Genomic signatures of experimental adaptive radiation in Drosophila.

Abiotic environmental factors play a fundamental role in determining the distribution, abundance and adaptive diversification of species. Empowered by new technologies enabling rapid and increasingly accurate examination of genomic variation in populations, researchers may gain new insights into the genomic background of adaptive radiation and stress resistance. We investigated genomic variation across generations of large-scale experimental selection regimes originating from a single founder population of Drosophila melanogaster, diverging in response to ecologically relevant environmental stressors: heat shock, heat knock down, cold shock, desiccation and starvation. When compared to the founder population, and to parallel unselected controls, there were more than 100,000 single nucleotide polymorphisms (SNPs) displaying consistent allelic changes in response to selective pressures across generations. These SNPs were found in both coding and noncoding sequences, with the highest density in promoter regions, and involved a broad range of functionalities, including molecular chaperoning by heat-shock proteins. The SNP patterns were highly stressor-specific despite considerable variation among line replicates within each selection regime, as reflected by a principal component analysis, and co-occurred with selective sweep regions. Only ~15% of SNPs with putatively adaptive changes were shared by at least two selective regimes, while less than 1% of SNPs diverged in opposite directions. Divergent stressors driving evolution in the experimental system of adaptive radiation left distinct genomic signatures, most pronounced in starvation and heat-shock selection regimes.

The Kub5-Hera/RPRD1B interactome: a novel role in preserving genetic stability by regulating DNA mismatch repair.

Ku70-binding protein 5 (Kub5)-Hera (K-H)/RPRD1B maintains genetic integrity by concomitantly minimizing persistent R-loops and promoting repair of DNA double strand breaks (DSBs). We used tandem affinity purification-mass spectrometry, co-immunoprecipitation and gel-filtration chromatography to define higher-order protein complexes containing K-H scaffolding protein to gain insight into its cellular functions. We confirmed known protein partners (Ku70, RNA Pol II, p15RS) and discovered several novel associated proteins that function in RNA metabolism (Topoisomerase 1 and RNA helicases), DNA repair/replication processes (PARP1, MSH2, Ku, DNA-PKcs, MCM proteins, PCNA and DNA Pol δ) and in protein metabolic processes, including translation. Notably, this approach directed us to investigate an unpredicted involvement of K-H in DNA mismatch repair (MMR) where K-H depletion led to concomitant MMR deficiency and compromised global microsatellite stability. Mechanistically, MMR deficiency in K-H-depleted cells was a consequence of reduced stability of the core MMR proteins (MLH1 and PMS2) caused by elevated basal caspase-dependent proteolysis. Pan-caspase inhibitor treatment restored MMR protein loss. These findings represent a novel mechanism to acquire MMR deficiency/microsatellite alterations. A significant proportion of colon, endometrial and ovarian cancers exhibit k-h expression/copy number loss and may have severe mutator phenotypes with enhanced malignancies that are currently overlooked based on sporadic MSI+ screening.

A Test for Gene Flow among Sympatric and Allopatric Hawaiian Picture-Winged Drosophila.

The Hawaiian Drosophila are one of the most species-rich endemic groups in Hawaii and a spectacular example of adaptive radiation. Drosophila silvestris and D. heteroneura are two closely related picture-winged Drosophila species that occur sympatrically on Hawaii Island and are known to hybridize in nature, yet exhibit highly divergent behavioral and morphological traits driven largely through sexual selection. Their closest-related allopatric species, D. planitibia from Maui, exhibits hybrid male sterility and reduced behavioral reproductive isolation when crossed experimentally with D. silvestris or D. heteroneura. A modified four-taxon test for gene flow was applied to recently obtained genomes of the three Hawaiian Drosophila species. The analysis indicates recent gene flow in sympatry, but also, although less extensive, between allopatric species. This study underscores the prevalence of gene flow, even in taxonomic groups considered classic examples of allopatric speciation on islands. The potential confounding effects of gene flow in phylogenetic and population genetics inference are discussed, as well as the implications for conservation.

Divergence of Drosophila melanogaster repeatomes in response to a sharp microclimate contrast in Evolution Canyon, Israel.

Repeat sequences, especially mobile elements, make up large portions of most eukaryotic genomes and provide enormous, albeit commonly underappreciated, evolutionary potential. We analyzed repeatomes of Drosophila melanogaster that have been diverging in response to a microclimate contrast in Evolution Canyon (Mount Carmel, Israel), a natural evolutionary laboratory with two abutting slopes at an average distance of only 200 m, which pose a constant ecological challenge to their local biotas. Flies inhabiting the colder and more humid north-facing slope carried about 6% more transposable elements than those from the hot and dry south-facing slope, in parallel to a suite of other genetic and phenotypic differences between the two populations. Nearly 50% of all mobile element insertions were slope unique, with many of them disrupting coding sequences of genes critical for cognition, olfaction, and thermotolerance, consistent with the observed patterns of thermotolerance differences and assortative mating.

EvoCor: a platform for predicting functionally related genes using phylogenetic and expression profiles.

The wealth of publicly available gene expression and genomic data provides unique opportunities for computational inference to discover groups of genes that function to control specific cellular processes. Such genes are likely to have co-evolved and be expressed in the same tissues and cells. Unfortunately, the expertise and computational resources required to compare tens of genomes and gene expression data sets make this type of analysis difficult for the average end-user. Here, we describe the implementation of a web server that predicts genes involved in affecting specific cellular processes together with a gene of interest. We termed the server 'EvoCor', to denote that it detects functional relationships among genes through evolutionary analysis and gene expression correlation. This web server integrates profiles of sequence divergence derived by a Hidden Markov Model (HMM) and tissue-wide gene expression patterns to determine putative functional linkages between pairs of genes. This server is easy to use and freely available at http://pilot-hmm.vbi.vt.edu/.

Nucleolar dominance and maternal control of 45S rDNA expression.

Using a system of interspecies hybrids, trihybrids, and recombinants with varying proportions of genomes from three distinct Xenopus species, we provide evidence for de novo epigenetic silencing of paternal 45 S ribosomal ribonucleic acid (rRNA) genes and their species-dependent expression dominance that escapes transcriptional inactivation after homologous recombination. The same pattern of imprinting is maintained in the offspring from mothers being genetic males (ZZ) sex-reversed to females, indicating that maternal control of ribosomal deoxyribonucleic acid (rDNA) expression is not sex-chromosome linked. Nucleolar dominance (nucleolus underdevelopment) in Xenopus hybrids appears to be associated with a major non-Mendelian reduction in the number of 45 S rDNA gene copies rather than a specific pattern of their expression. The loss of rRNA gene copies in F1 hybrids was non-random with respect to the parental species, with the transcriptionally dominant variant preferentially removed from hybrid zygotes. This dramatic disruption in the structure and function of 45 S rDNA impacts transcriptome patterns of small nucleolar RNAs and messenger RNAs, with genes from the ribosome and oxidative stress pathways being among the most affected. Unorthodoxies of rDNA inheritance and expression may be interpreted as hallmarks of genetic conflicts between parental genomes, as well as defensive epigenetic mechanisms employed to restore genome integrity.

Discretized Gaussian mixture for genotyping of microsatellite loci containing homopolymer runs.

MOTIVATION: Inferring lengths of inherited microsatellite alleles with single base pair resolution from short sequence reads is challenging due to several sources of noise caused by the repetitive nature of microsatellites and the technologies used to generate raw sequence data. RESULTS: We have developed a program, GenoTan, using a discretized Gaussian mixture model combined with a rules-based approach to identify inherited variation of microsatellite loci from short sequence reads without paired-end information. It effectively distinguishes length variants from noise including insertion/deletion errors in homopolymer runs by addressing the bidirectional aspect of insertion and deletion errors in sequence reads. Here we first introduce a homopolymer decomposition method which estimates error bias toward insertion or deletion in homopolymer sequence runs. Combining these approaches, GenoTan was able to genotype 94.9% of microsatellite loci accurately from simulated data with 40x sequence coverage quickly while the other programs showed <90% correct calls for the same data and required 5∼30× more computational time than GenoTan. It also showed the highest true-positive rate for real data using mixed sequence data of two Drosophila inbred lines, which was a novel validation approach for genotyping. AVAILABILITY: GenoTan is open-source software available at http://genotan.sourceforge.net.

Large scale comparison of non-human sequences in human sequencing data.

Several studies have demonstrated that unmapped reads in next generation sequencing data could be used to identify infectious agents or structural variants, but there has been no intensive effort to analyze and classify all non-human sequences found in individual large data sets. To identify commonality in non-human sequences by infectious agents and putative contamination events, we analyzed non-human sequences in 150 genomic sequencing data files from the 1000 Genomes Project and observed that 0.13% of reads on average showed similarities to non-human genomes. We compared results among different sample groups divided based on ethnicities, sequencing centers and enrichment methods (whole genome sequencing vs. exome sequencing) and found that sequencing centers had specific signatures of contaminating genomes as 'time stamps'. We also observed many unmapped reads that falsely indicated contamination because of the high similarity of human sequences to sequences in non-human genome assemblies such as mouse and Nicotiana.

ReviSTER: an automated pipeline to revise misaligned reads to simple tandem repeats.

MOTIVATION: Simple tandem repeats are highly variable genetic elements and widespread in genomes of many organisms. Next-generation sequencing technologies have enabled a robust comparison of large numbers of simple tandem repeat loci; however, analysis of their variation using traditional sequence analysis approaches still remains limiting and problematic due to variants occurring in repeat sequences confusing alignment programs into mapping sequence reads to incorrect loci when the sequence reads are significantly different from the reference sequence. RESULTS: We have developed a program, ReviSTER, which is an automated pipeline using a 'local mapping reference reconstruction method' to revise mismapped or partially misaligned reads at simple tandem repeat loci. RevisSTER estimates alleles of repeat loci using a local alignment method and creates temporary local mapping reference sequences, and finally remaps reads to the local mapping references. Using this approach, ReviSTER was able to successfully revise reads misaligned to repeat loci from both simulated data and real data. AVAILABILITY: ReviSTER is open-source software available at http://revister.sourceforge.net. CONTACT: garner@vbi.vt.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Common outpatient diagnoses and associated treatments logged by osteopathic medical students within a geriatric population.

CONTEXT: Clinical clerkships provide osteopathic medical students the opportunity to participate in the diagnosis and treatment of commonly encountered medical conditions. Appropriate management of these conditions may include pharmacotherapy and/or nonpharmacologic interventions, such as osteopathic manipulative treatment (OMT). Opportunities may exist to expand the utilization of OMT in the management of common conditions, particularly for geriatric patients, who are at increased risk for adverse outcomes from pharmacologic treatments. OBJECTIVES: This study aimed to assess the most common diagnoses and corresponding treatments logged by osteopathic medical students within an ambulatory geriatric population. METHODS: Patient encounters logged electronically by osteopathic medical students were retrospectively reviewed to determine the most commonly reported diagnostic codes and their treatments. Logged interventions were filtered to include patients over the age of 65 years who were seen on family medicine rotations within an ambulatory setting. The top 10 diagnoses were sorted and assessed to determine the associated treatments, including medications, procedures, and OMT. RESULTS: Between January 2018 and June 2020, a total of 11,185 primary diagnoses were logged pertaining to the defined patient population. The most frequently documented diagnoses were essential hypertension (n=1,420; 12.7 %), encounter for well examination (n=1,144; 10.2 %), type 2 diabetes mellitus (n=837; 7.5 %), hyperlipidemia (n=346; 3.1 %), chronic obstructive pulmonary disease (COPD; n=278; 2.5 %), osteoarthritis (OA; n=221; 2.0 %), low back pain (LBP; n=202; 1.8 %), pain in joint (n=187; 1.7 %), hypothyroidism (n=164; 1.5 %), and urinary tract infections (n=160; 1.4 %). Three of the top 10 logged diagnoses were musculoskeletal in nature (OA, LBP, and pain in joint). Pharmacotherapy was reported as the predominant treatment for musculoskeletal conditions, with OMT being logged as a treatment for 10.9 % (n=50) of those cases. The most commonly logged medication class in the management of patients with those musculoskeletal conditions was nonsteroidal anti-inflammatory drugs (NSAIDs; n=128; 27.9 %), while opioids were the second most frequently documented class of medications (n=65; 14.2 %). CONCLUSIONS: Musculoskeletal complaints were commonly logged by osteopathic medical students within the studied population. Opioids were documented as a treatment for musculoskeletal conditions more frequently than OMT. As such, opportunities exist to expand the utilization of OMT during clinical clerkships and to decrease the frequency of prescribed medications for pain management.

Improved variation calling via an iterative backbone remapping and local assembly method for bacterial genomes.

Sequencing data analysis remains limiting and problematic, especially for low complexity repeat sequences and transposon elements due to inherent sequencing errors and short sequence read lengths. We have developed a program, ReviSeq, which uses a hybrid method composed of iterative remapping and local assembly upon a bacterial sequence backbone. Application of this method to six Brucella suis field isolates compared to the newly revised B. suis 1330 reference genome identified on average 13, 15, 19 and 9 more variants per sample than STAMPY/SAMtools, BWA/SAMtools, iCORN and BWA/PINDEL pipelines, and excluded on average 4, 2, 3 and 19 variants per sample, respectively. In total, using this iterative approach, we identified on average 87 variants including SNVs, short INDELs and long INDELs per strain when compared to the reference. Our program outperforms other methods especially for long INDEL calling. The program is available at http://reviseq.sourceforge.net.

Complete genome sequence of Brucella suis VBI22, isolated from bovine milk.

Brucella suis is the causative agent of swine brucellosis and is known to be able to infect several different hosts, including cattle, dogs, and horses, without causing disease symptoms. Here we report the complete genome sequence of Brucella suis VBI22, which was isolated from raw milk from an infected cow.

An isogenetic myoblast expression screen identifies DUX4-mediated FSHD-associated molecular pathologies.

Facioscapulohumeral muscular dystrophy (FSHD) is caused by an unusual deletion with neomorphic activity. This deletion derepresses genes in cis; however which candidate gene causes the FSHD phenotype, and through what mechanism, is unknown. We describe a novel genetic tool, inducible cassette exchange, enabling rapid generation of isogenetically modified cells with conditional and variable transgene expression. We compare the effects of expressing variable levels of each FSHD candidate gene on myoblasts. This screen identified only one gene with overt toxicity: DUX4 (double homeobox, chromosome 4), a protein with two homeodomains, each similar in sequence to Pax3 and Pax7. DUX4 expression recapitulates key features of the FSHD molecular phenotype, including repression of MyoD and its target genes, diminished myogenic differentiation, repression of glutathione redox pathway components, and sensitivity to oxidative stress. We further demonstrate competition between DUX4 and Pax3/Pax7: when either Pax3 or Pax7 is expressed at high levels, DUX4 is no longer toxic. We propose a hypothesis for FSHD in which DUX4 expression interferes with Pax7 in satellite cells, and inappropriately regulates Pax targets, including myogenic regulatory factors, during regeneration.

Sporadic breast cancer patients' germline DNA exhibit an AT-rich microsatellite signature.

Using a custom CGH-like oligonucleotide array to measure the global microsatellite content in the genomes of 72 cancer, cancer-free, and high risk patient and cell line samples (56 germline DNA and 16 in tumor or tumor cell line DNA) we found a unique, reproducible, and statistically significant pattern of 18 motif-specific microsatellite families (out of 962 possible 1-6 mer repeats) in breast cancer patient germline and tumor DNA, but not in germline DNA of cancer-free volunteer controls or in breast cancer patients with BRCA1/2 mutations. These high-similarity A/T rich repetitive motifs were also more pronounced in the germlines and tumors of colon cancer tumor patients (3/6 samples) and microsatellite unstable colon cancer cell lines; however, germline DNA of sporadic breast cancer patients exhibited the largest global content shift for those motifs with extreme AT/GC ratios. These results indicate that global microsatellite variability is complex, suggest the existence of a previously unknown genomic destabilization mechanism in breast cancer patients' germline DNA, and warrant further testing of such microsatellite variability as a predictor of future breast cancer development.

Tracing the Genetic Evolution of Canine Parvovirus Type 2 (CPV-2) in Thailand.

Canine parvovirus type 2 (CPV-2) is responsible for hemorrhagic gastroenteritis in dogs worldwide. High genomic substitution rates in CPV-2 contribute to the progressive emergence of novel variants with increased ability to evade the host immune response. Three studies have analyzed the genomic mutations of CPV-2 variants in Thailand. These investigations were independently conducted at different timepoints. Thus, a retrospective integrated analysis of CPV-2 genomic mutations has not been fully performed. Our study aimed at evaluating the evolutionary changes in CPV-2 in Thailand from 2003 to 2019. Two hundred and sixty-eight Thai CPV-2 nucleotide sequences were used for multiple amino acid sequence alignment and phylogenetic analyses. From 2003 to 2010, CPV-2a and -2b were the only variants detected. CPV-2c, emerged in 2014, replacing CPV-2a and -2b, and has become a major variant in 2019. Phylogenetic analysis revealed that the proposed mutation pattern of VP2 amino acid residues could help distinguish Thai CPV-2 variants. This comprehensive examination provides insight into the genomic evolution of CPV-2 in Thailand since its first reporting in 2003, which may facilitate the surveillance of the potential genetic alteration of emergent CPV-2 variants.