Admix-kit: an integrated toolkit and pipeline for genetic analyses of admixed populations.

SUMMARY: Admixed populations, with their unique and diverse genetic backgrounds, are often underrepresented in genetic studies. This oversight not only limits our understanding but also exacerbates existing health disparities. One major barrier has been the lack of efficient tools tailored for the special challenges of genetic studies of admixed populations. Here, we present admix-kit, an integrated toolkit and pipeline for genetic analyses of admixed populations. Admix-kit implements a suite of methods to facilitate genotype and phenotype simulation, association testing, genetic architecture inference, and polygenic scoring in admixed populations. AVAILABILITY AND IMPLEMENTATION: Admix-kit package is open-source and available at https://github.com/KangchengHou/admix-kit. Additionally, users can use the pipeline designed for admixed genotype simulation available at https://github.com/UW-GAC/admix-kit_workflow.

Novel ancestry-specific primary open-angle glaucoma loci and shared biology with vascular mechanisms and cell proliferation.

Primary open-angle glaucoma (POAG), a leading cause of irreversible blindness globally, shows disparity in prevalence and manifestations across ancestries. We perform meta-analysis across 15 biobanks (of the Global Biobank Meta-analysis Initiative) (n = 1,487,441: cases = 26,848) and merge with previous multi-ancestry studies, with the combined dataset representing the largest and most diverse POAG study to date (n = 1,478,037: cases = 46,325) and identify 17 novel significant loci, 5 of which were ancestry specific. Gene-enrichment and transcriptome-wide association analyses implicate vascular and cancer genes, a fifth of which are primary ciliary related. We perform an extensive statistical analysis of SIX6 and CDKN2B-AS1 loci in human GTEx data and across large electronic health records showing interaction between SIX6 gene and causal variants in the chr9p21.3 locus, with expression effect on CDKN2A/B. Our results suggest that some POAG risk variants may be ancestry specific, sex specific, or both, and support the contribution of genes involved in programmed cell death in POAG pathogenesis.

Building a vertically integrated genomic learning health system: The biobank at the Colorado Center for Personalized Medicine.

Precision medicine initiatives across the globe have led to a revolution of repositories linking large-scale genomic data with electronic health records, enabling genomic analyses across the entire phenome. Many of these initiatives focus solely on research insights, leading to limited direct benefit to patients. We describe the biobank at the Colorado Center for Personalized Medicine (CCPM Biobank) that was jointly developed by the University of Colorado Anschutz Medical Campus and UCHealth to serve as a unique, dual-purpose research and clinical resource accelerating personalized medicine. This living resource currently has more than 200,000 participants with ongoing recruitment. We highlight the clinical, laboratory, regulatory, and HIPAA-compliant informatics infrastructure along with our stakeholder engagement, consent, recontact, and participant engagement strategies. We characterize aspects of genetic and geographic diversity unique to the Rocky Mountain region, the primary catchment area for CCPM Biobank participants. We leverage linked health and demographic information of the CCPM Biobank participant population to demonstrate the utility of the CCPM Biobank to replicate complex trait associations in the first 33,674 genotyped individuals across multiple disease domains. Finally, we describe our current efforts toward return of clinical genetic test results, including high-impact pathogenic variants and pharmacogenetic information, and our broader goals as the CCPM Biobank continues to grow. Bringing clinical and research interests together fosters unique clinical and translational questions that can be addressed from the large EHR-linked CCPM Biobank resource within a HIPAA- and CLIA-certified environment.

Thyroid Cancer Polygenic Risk Score Improves Classification of Thyroid Nodules as Benign or Malignant.

CONTEXT: Thyroid nodule ultrasound-based risk stratification schemas rely on the presence of high-risk sonographic features. However, some malignant thyroid nodules have benign appearance on thyroid ultrasound. New methods for thyroid nodule risk assessment are needed. OBJECTIVE: We investigated polygenic risk score (PRS) accounting for inherited thyroid cancer risk combined with ultrasound-based analysis for improved thyroid nodule risk assessment. METHODS: The convolutional neural network classifier was trained on thyroid ultrasound still images and cine clips from 621 thyroid nodules. Phenome-wide association study (PheWAS) and PRS PheWAS were used to optimize PRS for distinguishing benign and malignant nodules. PRS was evaluated in 73 346 participants in the Colorado Center for Personalized Medicine Biobank. RESULTS: When the deep learning model output was combined with thyroid cancer PRS and genetic ancestry estimates, the area under the receiver operating characteristic curve (AUROC) of the benign vs malignant thyroid nodule classifier increased from 0.83 to 0.89 (DeLong, P value = .007). The combined deep learning and genetic classifier achieved a clinically relevant sensitivity of 0.95, 95% CI [0.88-0.99], specificity of 0.63 [0.55-0.70], and positive and negative predictive values of 0.47 [0.41-0.58] and 0.97 [0.92-0.99], respectively. AUROC improvement was consistent in European ancestry-stratified analysis (0.83 and 0.87 for deep learning and deep learning combined with PRS classifiers, respectively). Elevated PRS was associated with a greater risk of thyroid cancer structural disease recurrence (ordinal logistic regression, P value = .002). CONCLUSION: Augmenting ultrasound-based risk assessment with PRS improves diagnostic accuracy.

Genomic Selection Signals in Andean Highlanders Reveal Adaptive Placental Metabolic Phenotypes That Are Disrupted in Preeclampsia.

BACKGROUND: The chronic hypoxia of high-altitude residence poses challenges for tissue oxygen supply and metabolism. Exposure to high altitude during pregnancy increases the incidence of hypertensive disorders of pregnancy and fetal growth restriction and alters placental metabolism. High-altitude ancestry protects against altitude-associated fetal growth restriction, indicating hypoxia tolerance that is genetic in nature. Yet, not all babies are protected and placental pathologies associated with fetal growth restriction occur in some Andean highlanders. METHODS: We examined placental metabolic function in 79 Andeans (18-45 years; 39 preeclamptic and 40 normotensive) living in La Paz, Bolivia (3600-4100 m) delivered by unlabored Cesarean section. Using a selection-nominated approach, we examined links between putatively adaptive genetic variation and phenotypes related to oxygen delivery or placental metabolism. RESULTS: Mitochondrial oxidative capacity was associated with fetal oxygen delivery in normotensive but not preeclamptic placenta and was also suppressed in term preeclamptic pregnancy. Maternal haplotypes in or within 200 kb of selection-nominated genes were associated with lower placental mitochondrial respiratory capacity (PTPRD [protein tyrosine phosphatase receptor-δ]), lower maternal plasma erythropoietin (CPT2 [carnitine palmitoyl transferase 2], proopiomelanocortin, and DNMT3 [DNA methyltransferase 3]), and lower VEGF (vascular endothelial growth factor) in umbilical venous plasma (TBX5 [T-box transcription factor 5]). A fetal haplotype within 200 kb of CPT2 was associated with increased placental mitochondrial complex II capacity, placental nitrotyrosine, and GLUT4 (glucose transporter type 4) protein expression. CONCLUSIONS: Our findings reveal novel associations between putatively adaptive gene regions and phenotypes linked to oxygen delivery and placental metabolic function in highland Andeans, suggesting that such effects may be of genetic origin. Our findings also demonstrate maladaptive metabolic mechanisms in the context of preeclampsia, including dysregulation of placental oxygen consumption.

Admix-kit: An Integrated Toolkit and Pipeline for Genetic Analyses of Admixed Populations.

Admixed populations, with their unique and diverse genetic backgrounds, are often underrepresented in genetic studies. This oversight not only limits our understanding but also exacerbates existing health disparities. One major barrier has been the lack of efficient tools tailored for the special challenges of genetic study of admixed populations. Here, we present admix-kit, an integrated toolkit and pipeline for genetic analyses of admixed populations. Admix-kit implements a suite of methods to facilitate genotype and phenotype simulation, association testing, genetic architecture inference, and polygenic scoring in admixed populations.

Altered placental ion channel gene expression in preeclamptic high-altitude pregnancies.

High-altitude (>2,500 m) residence increases the risk of pregnancy vascular disorders such as fetal growth restriction and preeclampsia, each characterized by impaired placental function. Genetic attributes of highland ancestry confer relative protection against vascular disorders of pregnancy at high altitudes. Although ion channels have been implicated in placental function regulation, neither their expression in high-altitude placentas nor their relationship to high-altitude preeclampsia has been determined. Here, we measured the expression of 26 ion-channel genes in placentas from preeclampsia cases and normotensive controls in La Paz, Bolivia (3,850 m). In addition, we correlated gene transcription to maternal and infant ancestry proportions. Gene expression was assessed by PCR, genetic ancestry evaluated by ADMIXTURE, and ion channel proteins localized by immunofluorescence. In preeclamptic placentas, 11 genes were downregulated (ABCC9, ATP2A2, CACNA1C, KCNE1, KCNJ8, KCNK3, KCNMA1, KCNQ1, KCNQ4, PKD2, and TRPV6) and two were upregulated (KCNQ3 and SCNN1G). KCNE1 expression was positively correlated with high-altitude Amerindian ancestry and negatively correlated with non-high altitude. SCNN1G was negatively correlated with African ancestry, despite minimal African admixture. Most ion channels were localized in syncytiotrophoblasts (Cav1.2, TRPP2, TRPV6, and Kv7.1), whereas expression of Kv7.4 was primarily in microvillous membranes, Kir6.1 in chorionic plate and fetal vessels, and MinK in stromal cells. Our findings suggest a role for differential placental ion channel expression in the development of preeclampsia. Functional studies are needed to determine processes affected by these ion channels in the placenta and whether therapies directed at modulating their activity could influence the onset or severity of preeclampsia.

Disease risk and healthcare utilization among ancestrally diverse groups in the Los Angeles region.

An individual's disease risk is affected by the populations that they belong to, due to shared genetics and environmental factors. The study of fine-scale populations in clinical care is important for identifying and reducing health disparities and for developing personalized interventions. To assess patterns of clinical diagnoses and healthcare utilization by fine-scale populations, we leveraged genetic data and electronic medical records from 35,968 patients as part of the UCLA ATLAS Community Health Initiative. We defined clusters of individuals using identity by descent, a form of genetic relatedness that utilizes shared genomic segments arising due to a common ancestor. In total, we identified 376 clusters, including clusters with patients of Afro-Caribbean, Puerto Rican, Lebanese Christian, Iranian Jewish and Gujarati ancestry. Our analysis uncovered 1,218 significant associations between disease diagnoses and clusters and 124 significant associations with specialty visits. We also examined the distribution of pathogenic alleles and found 189 significant alleles at elevated frequency in particular clusters, including many that are not regularly included in population screening efforts. Overall, this work progresses the understanding of health in understudied communities and can provide the foundation for further study into health inequities.

COVID-19 Mortality in the Colorado Center for Personalized Medicine Biobank.

Over 6.37 million people have died from COVID-19 worldwide, but factors influencing COVID-19-related mortality remain understudied. We aimed to describe and identify risk factors for COVID-19 mortality in the Colorado Center for Personalized Medicine (CCPM) Biobank using integrated data sources, including Electronic Health Records (EHRs). We calculated cause-specific mortality and case-fatality rates for COVID-19 and common pre-existing health conditions defined by diagnostic phecodes and encounters in EHRs. We performed multivariable logistic regression analyses of the association between each pre-existing condition and COVID-19 mortality. Of the 155,859 Biobank participants enrolled as of July 2022, 20,797 had been diagnosed with COVID-19. Of 5334 Biobank participants who had died, 190 were attributed to COVID-19. The case-fatality rate was 0.91% and the COVID-19 mortality rate was 122 per 100,000 persons. The odds of dying from COVID-19 were significantly increased among older men, and those with 14 of the 61 pre-existing conditions tested, including hypertensive chronic kidney disease (OR: 10.14, 95% CI: 5.48, 19.16) and type 2 diabetes with renal manifestations (OR: 5.59, 95% CI: 3.42, 8.97). Male patients who are older and have pre-existing kidney diseases may be at higher risk for death from COVID-19 and may require special care.

COVID-19 Surveillance in the Biobank at the Colorado Center for Personalized Medicine: Observational Study.

BACKGROUND: Characterizing the experience and impact of the COVID-19 pandemic among various populations remains challenging due to the limitations inherent in common data sources, such as electronic health records (EHRs) or cross-sectional surveys. OBJECTIVE: This study aims to describe testing behaviors, symptoms, impact, vaccination status, and case ascertainment during the COVID-19 pandemic using integrated data sources. METHODS: In summer 2020 and 2021, we surveyed participants enrolled in the Biobank at the Colorado Center for Personalized Medicine (CCPM; N=180,599) about their experience with COVID-19. The prevalence of testing, symptoms, and impacts of COVID-19 on employment, family life, and physical and mental health were calculated overall and by demographic categories. Survey respondents who reported receiving a positive COVID-19 test result were considered a "confirmed case" of COVID-19. Using EHRs, we compared COVID-19 case ascertainment and characteristics in EHRs versus the survey. Positive cases were identified in EHRs using the International Statistical Classification of Diseases, 10th revision (ICD-10) diagnosis codes, health care encounter types, and encounter primary diagnoses. RESULTS: Of the 25,063 (13.9%) survey respondents, 10,661 (42.5%) had been tested for COVID-19, and of those, 1366 (12.8%) tested positive. Nearly half of those tested had symptoms or had been exposed to someone who was infected. Young adults (18-29 years) and Hispanics were more likely to have positive tests compared to older adults and persons of other racial/ethnic groups. Mental health (n=13,688, 54.6%) and family life (n=12,233, 48.8%) were most negatively affected by the pandemic and more so among younger groups and women; negative impacts on employment were more commonly reported among Black respondents. Of the 10,249 individuals who responded to vaccination questions from version 2 of the survey (summer 2021), 9770 (95.3%) had received the vaccine. After integration with EHR data up to the time of the survey completion, 1006 (4%) of the survey respondents had a discordant COVID-19 case status between EHRs and the survey. Using all longitudinal EHR and survey data, we identified 11,472 (6.4%) COVID-19-positive cases among Biobank participants. In comparison to COVID-19 cases identified through the survey, EHR-identified cases were younger and more likely to be Hispanic. CONCLUSIONS: We found that the COVID-19 pandemic has had far-reaching and varying effects among our Biobank participants. Integrated data assets, such as the Biobank at the CCPM, are key resources for population health monitoring in response to public health emergencies, such as the COVID-19 pandemic.

Allele imputation for the killer cell immunoglobulin-like receptor KIR3DL1/S1.

Highly polymorphic interaction of KIR3DL1 and KIR3DS1 with HLA class I ligands modulates the effector functions of natural killer (NK) cells and some T cells. This genetically determined diversity affects severity of infections, immune-mediated diseases, and some cancers, and impacts the course of immunotherapies, including transplantation. KIR3DL1 is an inhibitory receptor, and KIR3DS1 is an activating receptor encoded by the KIR3DL1/S1 gene that has more than 200 diverse and divergent alleles. Determination of KIR3DL1/S1 genotypes for medical application is hampered by complex sequence and structural variation, requiring targeted approaches to generate and analyze high-resolution allele data. To overcome these obstacles, we developed and optimized a model for imputing KIR3DL1/S1 alleles at high-resolution from whole-genome SNP data. We designed the model to represent a substantial component of human genetic diversity. Our Global imputation model is effective at genotyping KIR3DL1/S1 alleles with an accuracy ranging from 88% in Africans to 97% in East Asians, with mean specificity of 99% and sensitivity of 95% for alleles >1% frequency. We used the established algorithm of the HIBAG program, in a modification named Pulling Out Natural killer cell Genomics (PONG). Because HIBAG was designed to impute HLA alleles also from whole-genome SNP data, PONG allows combinatorial diversity of KIR3DL1/S1 with HLA-A and -B to be analyzed using complementary techniques on a single data source. The use of PONG thus negates the need for targeted sequencing data in very large-scale association studies where such methods might not be tractable.

Ancestral patterns of recessive dystrophic epidermolysis bullosa mutations in Hispanic populations suggest sephardic ancestry.

Recessive dystrophic epidermolysis bullosa (RDEB) is a rare genodermatosis caused by mutations in the gene coding for type VII collagen (COL7A1). More than 800 different pathogenic mutations in COL7A1 have been described to date; however, the ancestral origins of many of these mutations have not been precisely identified. In this study, 32 RDEB patient samples from the Southwestern United States, Mexico, Chile, and Colombia carrying common mutations in the COL7A1 gene were investigated to determine the origins of these mutations and the extent to which shared ancestry contributes to disease prevalence. The results demonstrate both shared European and American origins of RDEB mutations in distinct populations in the Americas and suggest the influence of Sephardic ancestry in at least some RDEB mutations of European origins. Knowledge of ancestry and relatedness among RDEB patient populations will be crucial for the development of future clinical trials and the advancement of novel therapeutics.

Soluble CD14 Levels in the Jackson Heart Study: Associations With Cardiovascular Disease Risk and Genetic Variants.

[Figure: see text].

Population genomic analyses of schistosome parasites highlight critical challenges facing endgame elimination efforts.

Schistosomiasis persists in Asian regions despite aggressive elimination measures. To identify factors enabling continued parasite transmission, we performed reduced representation genome sequencing on Schistosoma japonicum miracidia collected across multiple years from transmission hotspots in Sichuan, China. We discovered strong geographic structure, suggesting that local, rather than imported, reservoirs are key sources of persistent infections in the region. At the village level, parasites collected after referral for praziquantel treatment are closely related to local pre-treatment populations. Schistosomes within villages are also highly related, suggesting that only a few parasites from a limited number of hosts drive re-infection. The close familial relationships among miracidia from different human hosts also implicate short transmission routes among humans. At the individual host level, genetic evidence indicates that multiple humans retained infections following referral for treatment. Our findings suggest that end-game schistosomiasis control measures should focus on completely extirpating local parasite reservoirs and confirming successful treatment of infected human hosts.

Applicability of ancestral genotyping in pharmacogenomic research with hormonal contraception.

To compare etonogestrel pharmacokinetic and pharmacodynamic outcomes by both self-reported race/ethnicity and genetically determined ancestry among contraceptive implant users. We conducted a secondary analysis of our parent pharmacogenomic study of 350 implant users. We genotyped these reproductive-aged (18-45 years) women for 88 ancestry-informative single nucleotide polymorphisms. We then assigned each participant a proportion value for African (AFR), European (EUR), and Indigenous American (AMR) ancestry based on reference population data. We correlated genetic ancestry with self-reported race/ethnicity and utilized genetic ancestry proportion values as variables for previously performed association analyses with serum etonogestrel concentrations and progestin-related side effects (e.g., bothersome bleeding and subjective weight gain). We successfully estimated genetically determined ancestry for 332 participants. EUR, AFR, and AMR ancestry were each highly correlated with self-reported White/non-Hispanic race (r = 0.64, p = 4.14 × 10-40 ), Black/African American race (r = 0.88, p = 1.36 × 10-107 ), and Hispanic/Latina ethnicity (r = 0.68, p = 4.03 × 10-47 ), respectively. Neither genetically determined ancestry nor self-reported race/ethnicity were significantly associated with serum etonogestrel concentrations. AFR ancestry and self-reported Black race had similar associations with reporting monthly periods (odds ratio [OR] 2.18, p = 0.09 vs. OR 2.22, p = 0.02) and having received treatment for bothersome bleeding (OR 5.19, p = 0.005 vs. OR 4.73, p = 2.0 × 10-4 ). In multivariable logistic regression for subjective weight gain, AMR ancestry dropped out of the model in preference for self-reported Hispanic/Latina ethnicity. We found no new associations between genetically determined ancestry and contraceptive implant pharmacodynamics/pharmacokinetics. Self-reported race/ethnicity were strong surrogates for genetically determined ancestry among this population of contraceptive implant users. Our data suggest that self-reported race/ethnicity, capturing societal and cultural aspects, remain important to the investigation of progestin-related side effects.

Adaptive Admixture of HLA Class I Allotypes Enhanced Genetically Determined Strength of Natural Killer Cells in East Asians.

Human natural killer (NK) cells are essential for controlling infection, cancer, and fetal development. NK cell functions are modulated by interactions between polymorphic inhibitory killer cell immunoglobulin-like receptors (KIR) and polymorphic HLA-A, -B, and -C ligands expressed on tissue cells. All HLA-C alleles encode a KIR ligand and contribute to reproduction and immunity. In contrast, only some HLA-A and -B alleles encode KIR ligands and they focus on immunity. By high-resolution analysis of KIR and HLA-A, -B, and -C genes, we show that the Chinese Southern Han (CHS) are significantly enriched for interactions between inhibitory KIR and HLA-A and -B. This enrichment has had substantial input through population admixture with neighboring populations, who contributed HLA class I haplotypes expressing the KIR ligands B*46:01 and B*58:01, which subsequently rose to high frequency by natural selection. Consequently, over 80% of Southern Han HLA haplotypes encode more than one KIR ligand. Complementing the high number of KIR ligands, the CHS KIR locus combines a high frequency of genes expressing potent inhibitory KIR, with a low frequency of those expressing activating KIR. The Southern Han centromeric KIR region encodes strong, conserved, inhibitory HLA-C-specific receptors, and the telomeric region provides a high number and diversity of inhibitory HLA-A and -B-specific receptors. In all these characteristics, the CHS represent other East Asians, whose NK cell repertoires are thus enhanced in quantity, diversity, and effector strength, likely augmenting resistance to endemic viral infections.

Finding and extending ancient simple sequence repeat-derived regions in the human genome.

BACKGROUND: Previously, 3% of the human genome has been annotated as simple sequence repeats (SSRs), similar to the proportion annotated as protein coding. The origin of much of the genome is not well annotated, however, and some of the unidentified regions are likely to be ancient SSR-derived regions not identified by current methods. The identification of these regions is complicated because SSRs appear to evolve through complex cycles of expansion and contraction, often interrupted by mutations that alter both the repeated motif and mutation rate. We applied an empirical, kmer-based, approach to identify genome regions that are likely derived from SSRs. RESULTS: The sequences flanking annotated SSRs are enriched for similar sequences and for SSRs with similar motifs, suggesting that the evolutionary remains of SSR activity abound in regions near obvious SSRs. Using our previously described P-clouds approach, we identified 'SSR-clouds', groups of similar kmers (or 'oligos') that are enriched near a training set of unbroken SSR loci, and then used the SSR-clouds to detect likely SSR-derived regions throughout the genome. CONCLUSIONS: Our analysis indicates that the amount of likely SSR-derived sequence in the human genome is 6.77%, over twice as much as previous estimates, including millions of newly identified ancient SSR-derived loci. SSR-clouds identified poly-A sequences adjacent to transposable element termini in over 74% of the oldest class of Alu (roughly, AluJ), validating the sensitivity of the approach. Poly-A's annotated by SSR-clouds also had a length distribution that was more consistent with their poly-A origins, with mean about 35 bp even in older Alus. This work demonstrates that the high sensitivity provided by SSR-Clouds improves the detection of SSR-derived regions and will enable deeper analysis of how decaying repeats contribute to genome structure.

Conservation, Extensive Heterozygosity, and Convergence of Signaling Potential All Indicate a Critical Role for KIR3DL3 in Higher Primates.

Natural killer (NK) cell functions are modulated by polymorphic killer cell immunoglobulin-like receptors (KIR). Among 13 human KIR genes, which vary by presence and copy number, KIR3DL3 is ubiquitously present in every individual across diverse populations. No ligand or function is known for KIR3DL3, but limited knowledge of expression suggests involvement in reproduction, likely during placentation. With 157 human alleles, KIR3DL3 is also highly polymorphic and we show heterozygosity exceeds that of HLA-B in many populations. The external domains of catarrhine primate KIR3DL3 evolved as a conserved lineage distinct from other KIR. Accordingly, and in contrast to other KIR, we show the focus of natural selection does not correspond exclusively to known ligand binding sites. Instead, a strong signal for diversifying selection occurs in the D1 Ig domain at a site involved in receptor aggregation, which we show is polymorphic in humans worldwide, suggesting differential ability for receptor aggregation. Meanwhile in the cytoplasmic tail, the first of two inhibitory tyrosine motifs (ITIM) is conserved, whereas independent genomic events have mutated the second ITIM of KIR3DL3 alleles in all great apes. Together, these findings suggest that KIR3DL3 binds a conserved ligand, and a function requiring both receptor aggregation and inhibitory signal attenuation. In this model KIR3DL3 resembles other NK cell inhibitory receptors having only one ITIM, which interact with bivalent downstream signaling proteins through dimerization. Due to the extensive conservation across species, selection, and other unusual properties, we consider elucidating the ligand and function of KIR3DL3 to be a pressing question.

Squamate reptiles challenge paradigms of genomic repeat element evolution set by birds and mammals.

Broad paradigms of vertebrate genomic repeat element evolution have been largely shaped by analyses of mammalian and avian genomes. Here, based on analyses of genomes sequenced from over 60 squamate reptiles (lizards and snakes), we show that patterns of genomic repeat landscape evolution in squamates challenge such paradigms. Despite low variance in genome size, squamate genomes exhibit surprisingly high variation among species in abundance (ca. 25-73% of the genome) and composition of identifiable repeat elements. We also demonstrate that snake genomes have experienced microsatellite seeding by transposable elements at a scale unparalleled among eukaryotes, leading to some snake genomes containing the highest microsatellite content of any known eukaryote. Our analyses of transposable element evolution across squamates also suggest that lineage-specific variation in mechanisms of transposable element activity and silencing, rather than variation in species-specific demography, may play a dominant role in driving variation in repeat element landscapes across squamate phylogeny.

Whole Genome Amplification and Reduced-Representation Genome Sequencing of Schistosoma japonicum Miracidia.

BACKGROUND: In areas where schistosomiasis control programs have been implemented, morbidity and prevalence have been greatly reduced. However, to sustain these reductions and move towards interruption of transmission, new tools for disease surveillance are needed. Genomic methods have the potential to help trace the sources of new infections, and allow us to monitor drug resistance. Large-scale genotyping efforts for schistosome species have been hindered by cost, limited numbers of established target loci, and the small amount of DNA obtained from miracidia, the life stage most readily acquired from humans. Here, we present a method using next generation sequencing to provide high-resolution genomic data from S. japonicum for population-based studies. METHODOLOGY/PRINCIPAL FINDINGS: We applied whole genome amplification followed by double digest restriction site associated DNA sequencing (ddRADseq) to individual S. japonicum miracidia preserved on Whatman FTA cards. We found that we could effectively and consistently survey hundreds of thousands of variants from 10,000 to 30,000 loci from archived miracidia as old as six years. An analysis of variation from eight miracidia obtained from three hosts in two villages in Sichuan showed clear population structuring by village and host even within this limited sample. CONCLUSIONS/SIGNIFICANCE: This high-resolution sequencing approach yields three orders of magnitude more information than microsatellite genotyping methods that have been employed over the last decade, creating the potential to answer detailed questions about the sources of human infections and to monitor drug resistance. Costs per sample range from $50-$200, depending on the amount of sequence information desired, and we expect these costs can be reduced further given continued reductions in sequencing costs, improvement of protocols, and parallelization. This approach provides new promise for using modern genome-scale sampling to S. japonicum surveillance, and could be applied to other schistosome species and other parasitic helminthes.