Accurate long-read sequencing allows assembly of the duplicated RHD and RHCE genes harboring variants relevant to blood transfusion.

Next-generation sequencing (NGS) technologies have transformed medical genetics. However, short-read lengths pose a limitation on identification of structural variants, sequencing repetitive regions, phasing of distant nucleotide changes, and distinguishing highly homologous genomic regions. Long-read sequencing technologies may offer improvements in the characterization of genes that are currently difficult to assess. We used a combination of targeted DNA capture, long-read sequencing, and a customized bioinformatics pipeline to fully assemble the RH region, which harbors variation relevant to red cell donor-recipient mismatch, particularly among patients with sickle cell disease. RHD and RHCE are a pair of duplicated genes located within an ∼175 kb region on human chromosome 1 that have high sequence similarity and frequent structural variations. To achieve the assembly, we utilized palindrome repeats in PacBio SMRT reads to obtain consensus sequences of 2.1 to 2.9 kb average length with over 99% accuracy. We used these long consensus sequences to identify 771 assembly markers and to phase the RHD-RHCE region with high confidence. The dataset enabled direct linkage between coding and intronic variants, phasing of distant SNPs to determine RHD-RHCE haplotypes, and identification of known and novel structural variations along with the breakpoints. A limiting factor in phasing is the frequency of heterozygous assembly markers and therefore was most successful in samples from African Black individuals with increased heterogeneity at the RH locus. Overall, this approach allows RH genotyping and de novo assembly in an unbiased and comprehensive manner that is necessary to expand application of NGS technology to high-resolution RH typing.

Long-insert sequence capture detects high copy numbers in a defence-related beta-glucosidase gene ßglu-1 with large variations in white spruce but not Norway spruce.

Conifers are long-lived and slow-evolving, thus requiring effective defences against their fast-evolving insect natural enemies. The copy number variation (CNV) of two key acetophenone biosynthesis genes Ugt5/Ugt5b and ßglu-1 may provide a plausible mechanism underlying the constitutively variable defence in white spruce (Picea glauca) against its primary defoliator, spruce budworm. This study develops a long-insert sequence capture probe set (Picea_hung_p1.0) for quantifying copy number of ßglu-1-like, Ugt5-like genes and single-copy genes on 38 Norway spruce (Picea abies) and 40 P. glauca individuals from eight and nine provenances across Europe and North America respectively. We developed local assemblies (Piabi_c1.0 and Pigla_c.1.0), full-length transcriptomes (PIAB_v1 and PIGL_v1), and gene models to characterise the diversity of ßglu-1 and Ugt5 genes. We observed very large copy numbers of ßglu-1, with up to 381 copies in a single P. glauca individual. We observed among-provenance CNV of ßglu-1 in P. glauca but not P. abies. Ugt5b was predominantly single-copy in both species. This study generates critical hypotheses for testing the emergence and mechanism of extreme CNV, the dosage effect on phenotype, and the varying copy number of genes with the same pathway. We demonstrate new approaches to overcome experimental challenges in genomic research in conifer defences.

The effect of missing data on evolutionary analysis of sequence capture bycatch, with application to an agricultural pest.

Sequence capture is a genomic technique that selectively enriches target sequences before high throughput next-generation sequencing, to generate specific sequences of interest. Off-target or 'bycatch' data are often discarded from capture experiments, but can be leveraged to address evolutionary questions under some circumstances. Here, we investigated the effects of missing data on a variety of evolutionary analyses using bycatch from an exon capture experiment on the global pest moth, Helicoverpa armigera. We added > 200 new samples from across Australia in the form of mitogenomes obtained as bycatch from targeted sequence capture, and combined these into an additional larger dataset to total > 1000 mitochondrial cytochrome c oxidase subunit I (COI) sequences across the species' global distribution. Using discriminant analysis of principal components and Bayesian coalescent analyses, we showed that mitogenomes assembled from bycatch with up to 75% missing data were able to return evolutionary inferences consistent with higher coverage datasets and the broader literature surrounding H. armigera. For example, low-coverage sequences broadly supported the delineation of two H. armigera subspecies and also provided new insights into the potential for geographic turnover among these subspecies. However, we also identified key effects of dataset coverage and composition on our results. Thus, low-coverage bycatch data can offer valuable information for population genetic and phylodynamic analyses, but caution is required to ensure the reduced information does not introduce confounding factors, such as sampling biases, that drive inference. We encourage more researchers to consider maximizing the potential of the targeted sequence approach by examining evolutionary questions with their off-target bycatch where possible-especially in cases where no previous mitochondrial data exists-but recommend stratifying data at different genome coverage thresholds to separate sampling effects from genuine genomic signals, and to understand their implications for evolutionary research.

Four conserved amino acids on human papillomavirus E6 predict clinical high-risk types.

Human papillomavirus (HPV) types included in the genus alpha papillomavirus (alpha-HPVs) are subdivided into high- and low-risk HPVs associated with tumorigenicity. According to conventional risk classification, over 30 alpha-HPVs remain unclassified and HPV groups phylogenetically classified using the L1 gene do not exactly correspond to the conventional risk classification groups. Here, we propose a novel cervical lesion progression risk classification strategy. Using four E6 risk distinguishable amino acids (E6-RDAAs), we successfully expanded the conventional classification to encompass alpha-HPVs and resolve discrepancies. We validated our classification system using alpha-HPV-targeted sequence data of 325 cervical swab specimens from participants in Japan. Clinical outcomes significantly correlated with the E6-RDAA classification. Four of five HPV types in the data set that were not conventionally classified (HPV30, 34, 67, and 69) were high-risk according to our classification criteria. This report sheds light on the carcinogenicity of rare genital HPV types using a novel risk classification strategy.

Telomere integrated scoring system of myelodysplastic syndrome.

INTRODUCTION: Recently, multigene target sequencing is widely performed for the purpose of prognostic prediction and application of targeted therapy. Here, we proposed a new scoring system that encompasses gene variations, telomere length, and Revised International Prognostic Scoring System (IPSS-R) together in Asian myelodysplastic syndrome. METHODS: We developed a new scoring model of these variables: age ≥ 65 years + IPSS-R score + ASXL1 mutation + TP53 mutation + Telomere length (<5.37). According to this new scoring system, patients were divided into four groups: very good score cutoff (≤3.0), good (3.0-4.5), poor (4.5-7.0), and very poor (>7.0). RESULTS: The median OS was 170.1, 100.4, 46.0, and 12.0 months for very good, good, poor, and very poor, retrospectively (p < 0.001). Meanwhile, according to the conventional IPSS-R scoring system, the median OS was 141.3, 50.2, 93.0, 36.0, and 16.2 months for very low, low, intermediate, high, and very high, retrospectively (p < 0.001). CONCLUSIONS: The newly developed model incorporating molecular variations and TL yielded more clear separations of the survival curves. By adding the presence of gene mutation and telomere length to the existing IPSS-R, its predictive ability can be further improved in myelodysplastic syndrome.

The landscape of genetic aberrations in myxofibrosarcoma.

Myxofibrosarcoma (MFS) is a rare subtype of sarcoma, whose genetic basis is poorly understood. We analyzed 69 MFS cases using whole-genome (WGS), whole-exome (WES) and/or targeted-sequencing (TS). Newly sequenced genomic data were combined with additional deposited 116 MFS samples. WGS identified a high number of structural variations (SVs) per tumor most frequently affecting the TP53 and RB1 loci, 40% of tumors showed a BRCAness-associated mutation signature, and evidence of chromothripsis was found in all cases. Most frequently mutated/copy number altered genes affected known disease drivers such as TP53 (56.2%), CDKN2A/B (29.7%), RB1 (27.0%), ATRX (19.5%) and HDLBP (18.9%). Several previously unappreciated genetic aberrations including MUC17, FLG and ZNF780A were identified in more than 20% of patients. Longitudinal analysis of paired diagnosis and relapse time points revealed a 1.2-fold mutation number increase accompanied with substantial changes in clonal composition over time. Our study highlights the genetic complexity underlying sarcomagenesis of MFS.

Evolutionary origin of pathogenic GJB2 alleles in China.

The frequency of the pathogenic allele of the autosomal recessive deafness gene GJB2 varies among different populations in the world, and accumulates to a sufficiently high frequency in certain population. The purpose of this study is to investigate the origin and evolution of GJB2 pathogenic alleles in Chinese deaf patients. Children with non-syndromic hearing loss, and their parents, from 295 families were recruited. Customized capture probes targeted at 943 single nucleotide polymorphisms (SNPs) related to GJB2 gene were designed for sequencing of genomic DNA in blood samples. Haplotypes carrying pathogenic allele were analyzed through linkage disequilibrium block building, ancestry tracing, and extended haplotype heterozygosity calculation. Two pathogenic GJB2 alleles, c.235delC (18.41%) and c.109G > A (15.57%), were observed in 867 donors. For c.235delC allele, three different core haplotypes with one major haplotype (97.32%) were found, and their core SNPs were 100% conserved. For c.109G > A allele, six different haplotypes with one major haplotype (93.28%) were found and the major c.109G > A allele evolved from a specific ancestral haplotype. Geographical origins of donors carrying GJB2 c.109G > A and c.235delC core haplotypes centered between Qinghai and Neimenggu. GJB2 c.235delC has long-range linkage disequilibrium. No positive selection signature was found for GJB2 c.235delC or c.109G > A in the studied population. In conclusion, we discovered a single origin of GJB2 c.235delC allele and multiple independent origins of GJB2 c.109G > A allele. Alternative to positive selection or multiple independent recurrent mutation event, population bottleneck effect might account for the observed high population frequency of these pathogenic alleles.

Effects of sample age on data quality from targeted sequencing of museum specimens: what are we capturing in time?

BACKGROUND: Next generation sequencing (NGS) can recover DNA data from valuable extant and extinct museum specimens. However, archived or preserved DNA is difficult to sequence because of its fragmented, damaged nature, such that the most successful NGS methods for preserved specimens remain sub-optimal. Improving wet-lab protocols and comprehensively determining the effects of sample age on NGS library quality are therefore of vital importance. Here, I examine the relationship between sample age and several indicators of library quality following targeted NGS sequencing of ~ 1300 loci using 271 samples of pinned moth specimens (Helicoverpa armigera) ranging in age from 5 to 117 years. RESULTS: I find that older samples have lower DNA concentrations following extraction and thus require a higher number of indexing PCR cycles during library preparation. When sequenced reads are aligned to a reference genome or to only the targeted region, older samples have a lower number of sequenced and mapped reads, lower mean coverage, and lower estimated library sizes, while the percentage of adapters in sequenced reads increases significantly as samples become older. Older samples also show the poorest capture success, with lower enrichment and a higher improved coverage anticipated from further sequencing. CONCLUSIONS: Sample age has significant, measurable impacts on the quality of NGS data following targeted enrichment. However, incorporating a uracil-removing enzyme into the blunt end-repair step during library preparation could help to repair DNA damage, and using a method that prevents adapter-dimer formation may result in improved data yields.

Capturing the Resistome: a Targeted Capture Method To Reveal Antibiotic Resistance Determinants in Metagenomes.

Identification of the nucleotide sequences encoding antibiotic resistance elements and determination of their association with antibiotic resistance are critical to improve surveillance and monitor trends in antibiotic resistance. Current methods to study antibiotic resistance in various environments rely on extensive deep sequencing or laborious culturing of fastidious organisms, both of which are heavily time-consuming operations. An accurate and sensitive method to identify both rare and common resistance elements in complex metagenomic samples is needed. Referencing the sequences in the Comprehensive Antibiotic Resistance Database, we designed a set of 37,826 probes to specifically target over 2,000 nucleotide sequences associated with antibiotic resistance in clinically relevant bacteria. Testing of this probe set on DNA libraries generated from multidrug-resistant bacteria to selectively capture resistance genes reproducibly produced higher numbers of reads on target at a greater length of coverage than shotgun sequencing. We also identified additional resistance gene sequences from human gut microbiome samples that sequencing alone was not able to detect. Our method to capture the resistome enables a sensitive means of gene detection in diverse environments where genes encoding antibiotic resistance represent less than 0.1% of the metagenome.

Targeted capture-based NGS is superior to multiplex PCR-based NGS for hereditary BRCA1 and BRCA2 gene analysis in FFPE tumor samples.

BACKGROUND: With the introduction of Olaparib treatment for BRCA-deficient recurrent ovarian cancer, testing for somatic and/or germline mutations in BRCA1/2 genes in tumor tissues became essential for treatment decisions. In most cases only formalin-fixed paraffin-embedded (FFPE) samples, containing fragmented and chemically modified DNA of minor quality, are available. Thus, multiplex PCR-based sequencing is most commonly applied in routine molecular testing, which is predominantly focused on the identification of known hot spot mutations in oncogenes. METHODS: We compared the overall performance of an adjusted targeted capture-based enrichment protocol and a multiplex PCR-based approach for calling of pathogenic SNVs and InDels using DNA extracted from 13 FFPE tissue samples. We further applied both strategies to seven blood samples and five matched FFPE tumor tissues of patients with known germline exon-spanning deletions and gene-wide duplications in BRCA1/2 to evaluate CNV detection based solely on panel NGS data. Finally, we analyzed DNA from FFPE tissues of 11 index patients from families suspected of having hereditary breast and ovarian cancer, of whom no blood samples were available for testing, in order to identify underlying pathogenic germline BRCA1/2 mutations. RESULTS: The multiplex PCR-based protocol produced inhomogeneous coverage among targets of each sample and between samples as well as sporadic amplicon drop out, leading to insufficiently or non-covered nucleotides, which subsequently hindered variant detection. This protocol further led to detection of PCR-artifacts that could easily have been misinterpreted as pathogenic mutations. No such limitations were observed by application of an adjusted targeted capture-based protocol, which allowed for CNV calling with 86% sensitivity and 100% specificity. All pathogenic CNVs were confirmed in the five matched FFPE tumor samples from patients carrying known pathogenic germline mutations and we additionally identified somatic loss of the second allele in BRCA1/2. Furthermore we detected pathogenic BRCA1/2 variants in four the eleven FFPE samples from patients of whom no blood was available for analysis. CONCLUSIONS: We demonstrate that an adjusted targeted capture-based enrichment protocol is superior to commonly applied multiplex PCR-based protocols for reliable BRCA1/2 variant detection, including CNV-detection, using FFPE tumor samples.

The Microtus voles: Resolving the phylogeny of one of the most speciose mammalian genera using genomics.

Sequential rapid radiations pose some of the greatest difficulties in phylogenetics, especially when analysing only a small number of genetic markers. Given that most of the speciation events occur in quick succession at various points in time, this creates particular challenges in determining phylogenetic relationships, i.e. branching order and divergence times. With the development of high throughput sequencing, thousands of markers can now readily be used to tackle these issues. Microtus is a speciose genus currently composed of 65 species that evolved over the last 2 million years. Although it is a well-studied group, there is still phylogenetic uncertainty at various divergence levels. Building upon previous studies that generally used small numbers of mitochondrial and/or nuclear loci, in this genomic-scale study we used both mitochondrial and nuclear data to study the rapid radiation within Microtus, using partial mitogenomes and genotyping-by-sequencing (GBS) on seven species representing five Microtus subgenera and the main biogeographic ranges where this group occurs. Both types of genome (mitochondrial and nuclear) generated similar tree topologies, with a basal split of the Nearctic (M. ochrogaster) and Holarctic (M. oeconomus) species, and then a subdivision of the five Palearctic species into two subgroups. These data support the occurrence of two European radiations, one North American radiation, and a later expansion of M. oeconomus from Asia to both Europe and North America. We further resolved the positioning of M. cabrerae as sister group of M. agrestis and refute the claim that M. cabrerae should be elevated to its own genus (Iberomys). Finally, the data support ongoing speciation events, especially within M. agrestis, with high levels of genetic divergence between the three Evolutionarily Significant Units (ESUs) previously identified. Similar high levels of divergence were also found among ESUs within M. oeconomus and M. arvalis.

Comparative analysis of targeted long read sequencing approaches for characterization of a plant's immune receptor repertoire.

BACKGROUND: The Oxford Nanopore Technologies MinION™ sequencer is a small, portable, low cost device that is accessible to labs of all sizes and attractive for in-the-field sequencing experiments. Selective breeding of crops has led to a reduction in genetic diversity, and wild relatives are a key source of new genetic resistance to pathogens, usually via NLR immune receptor-encoding genes. Recent studies have demonstrated how crop NLR repertoires can be targeted for sequencing on Illumina or PacBio (RenSeq) and the specific gene conveying pathogen resistance identified. RESULTS: Sequence yields per MinION run are lower than Illumina, making targeted resequencing an efficient approach. While MinION generates long reads similar to PacBio it doesn't generate the highly accurate multipass consensus reads, which presents downstream bioinformatics challenges. Here we demonstrate how MinION data can be used for RenSeq achieving similar results to the PacBio and how novel NLR gene fusions can be identified via a Nanopore RenSeq pipeline. CONCLUSION: The described library preparation and bioinformatics methods should be applicable to other gene families or any targeted long DNA fragment nanopore sequencing project.

Microglandular adenosis associated with triple-negative breast cancer is a neoplastic lesion of triple-negative phenotype harbouring TP53 somatic mutations.

Microglandular adenosis (MGA) is a rare proliferative lesion of the breast composed of small glands lacking myoepithelial cells and lined by S100-positive, oestrogen receptor (ER)-negative, progesterone receptor (PR)-negative, and HER2-negative epithelial cells. There is evidence to suggest that MGA may constitute a non-obligate precursor of triple-negative breast cancer (TNBC). We sought to define the genomic landscape of pure MGA and of MGA, atypical MGA (AMGA) and associated TNBCs, and to determine whether synchronous MGA, AMGA, and TNBCs would be clonally related. Two pure MGAs and eight cases of MGA and/or AMGA associated with in situ or invasive TNBC were collected, microdissected, and subjected to massively parallel sequencing targeting all coding regions of 236 genes recurrently mutated in breast cancer or related to DNA repair. Pure MGAs lacked clonal non-synonymous somatic mutations and displayed limited copy number alterations (CNAs); conversely, all MGAs (n = 7) and AMGAs (n = 3) associated with TNBC harboured at least one somatic non-synonymous mutation (range 3-14 and 1-10, respectively). In all cases where TNBCs were analyzed, identical TP53 mutations and similar patterns of gene CNAs were found in the MGA and/or AMGA and in the associated TNBC. In the MGA/AMGA associated with TNBC lacking TP53 mutations, somatic mutations affecting PI3K pathway-related genes (eg PTEN, PIK3CA, and INPP4B) and tyrosine kinase receptor signalling-related genes (eg ERBB3 and FGFR2) were identified. At diagnosis, MGAs associated with TNBC were found to display subclonal populations, and clonal shifts in the progression from MGA to AMGA and/or to TNBC were observed. Our results demonstrate the heterogeneity of MGAs, and that MGAs associated with TNBC, but not necessarily pure MGAs, are genetically advanced, clonal, and neoplastic lesions harbouring recurrent mutations in TP53 and/or other cancer genes, supporting the notion that a subset of MGAs and AMGAs may constitute non-obligate precursors of TNBCs.

Reduction of bitter taste receptor gene family in folivorous colobine primates relative to omnivorous cercopithecine primates.

Bitter taste perception is important in preventing animals from ingesting potentially toxic compounds. Whole-genome assembly (WGA) data have revealed that bitter taste receptor genes (TAS2Rs) comprise a multigene family with dozens of intact and disrupted genes in primates. However, publicly available WGA data are often incomplete, especially for multigene families. In this study, we employed a targeted capture (TC) approach specifically probing TAS2Rs for ten species of cercopithecid primates with diverse diets, including eight omnivorous cercopithecine species and two folivorous colobine species. We designed RNA probes for all TAS2Rs that we modeled to be intact in the common ancestor of cercopithecids ("ancestral-cercopithecid TAS2R gene set"). The TC was followed by short-read and high-depth massive-parallel sequencing. TC retrieved more intact TAS2R genes than found in WGA databases. We confirmed a large number of gene "births" at the common ancestor of cercopithecids and found that the colobine common ancestor and the cercopithecine common ancestor had contrasting trajectories: four gene "deaths" and three gene births, respectively. The number of intact TAS2R genes was markedly reduced in colobines (25-28 detected via TC and 20-26 detected via WGA analysis) as compared with cercopithecines (27-36 via TC and 19-30 via WGA). Birth or death events occurred at almost every phylogenetic-tree branch, making the composition of intact genes variable among species. These results show that evolutionary change in intact TAS2R genes is a complex process, refute a simple general prediction that herbivory favors more TAS2R genes, and have implications for understanding dietary adaptations and the evolution of detoxification abilities.

Targeted enrichment of whole-genome SNPs from highly burned skeletal remains.

Genetic assessment of highly incinerated and/or degraded human skeletal material is a persistent challenge in forensic DNA analysis, including identifying victims of mass disasters. Few studies have investigated the impact of thermal degradation on whole-genome single-nucleotide polymorphism (SNP) quality and quantity using next-generation sequencing (NGS). We present whole-genome SNP data obtained from the bones and teeth of 27 fire victims using two DNA extraction techniques. Extracts were converted to double-stranded DNA libraries then enriched for whole-genome SNPs using unpublished biotinylated RNA baits and sequenced on an Illumina NextSeq 550 platform. Raw reads were processed using the EAGER (Efficient Ancient Genome Reconstruction) pipeline, and the SNPs filtered and called using FreeBayes and GATK (v. 3.8). Mixed-effects modeling of the data suggest that SNP variability and preservation is predominantly determined by skeletal element and burn category, and not by extraction type. Whole-genome SNP data suggest that selecting long bones, hand and foot bones, and teeth subjected to temperatures <350°C are the most likely sources for higher genomic DNA yields. Furthermore, we observed an inverse correlation between the number of captured SNPs and the extent to which samples were burned, as well as a significant decrease in the total number of SNPs measured for samples subjected to temperatures >350°C. Our data complement previous analyses of burned human remains that compare extraction methods for downstream forensic applications and support the idea of adopting a modified Dabney extraction technique when traditional forensic methods fail to produce DNA yields sufficient for genetic identification.

Genes and mutations causing retinitis pigmentosa.

Retinitis pigmentosa (RP) is a heterogeneous set of inherited retinopathies with many disease-causing genes, many known mutations, and highly varied clinical consequences. Progress in finding treatments is dependent on determining the genes and mutations causing these diseases, which includes both gene discovery and mutation screening in affected individuals and families. Despite the complexity, substantial progress has been made in finding RP genes and mutations. Depending on the type of RP, and the technology used, it is possible to detect mutations in 30-80% of cases. One of the most powerful approaches to genetic testing is high-throughput 'deep sequencing', that is, next-generation sequencing (NGS). NGS has identified several novel RP genes but a substantial fraction of previously unsolved cases have mutations in genes that are known causes of retinal disease but not necessarily RP. Apparent discrepancy between the molecular defect and clinical findings may warrant reevaluation of patients and families. In this review, we summarize the current approaches to gene discovery and mutation detection for RP, and indicate pitfalls and unsolved problems. Similar considerations apply to other forms of inherited retinal disease.

Description of a new species of black coral in the family Aphanipathidae (Anthozoa, Antipatharia) from Puerto Rico.

Black corals (Anthozoa: Antipatharia) are an anthozoan lineage in the class Hexacorallia that occur across a wide range of habitats from the tropics to the poles and from surface waters to depths deeper than 8000 m. A new species of black coral, Aphanipathespuertoricoensissp. nov., collected with a remotely operated vehicle 357 m deep off Puerto Rico is recognized in the family Aphanipathidae. The new species is characterized by very long and loosely coiled primary branches and up to 0.5 mm tall spines with as many as 40 or more small conical tubercles. A phylogeny composed of 13 taxa that are closely related to the new species was reconstructed from 793 nuclear loci to show their systematic relationships. Our study integrated morphological and genomic data to show that this new species is distinct from other species in the genus Aphanipathes. Furthermore, our results add to the growing knowledge of black coral diversity, while further demonstrating the need for exploration in deep waters of the Caribbean Sea.

Linking breadfruit cultivar names across the globe connects histories after 230 years of separation.

Every crop has a story. The story of breadfruit (Artocarpus altilis), an increasingly valued staple crop in tropical agroforestry systems, is filled with intrigue, oppression, and remains incomplete. The Caribbean is a major producer and consumer of breadfruit, yet most breadfruit there came from a single 1793 introduction aimed at providing a cheap food source for slaves forced to work on British plantations. St. Vincent was the first significant point of Caribbean introduction and played a vital role in subsequent breadfruit distribution throughout the region. Hundreds of cultivars are documented in breadfruit's native Oceania. It remains a mystery, however, which ones were introduced to the Caribbean 230 years ago-still comprising the vast diversity found there today. Integrating local knowledge, historical documents and specimens, morphological data, and DNA, we identify eight major global breadfruit lineages-five of which are found in the Caribbean and likely represent the original 1793 introduction. Genetic data were able to match two Caribbean cultivar names confidently to their Oceania counterparts. Genetics and morphology together enabled additional possible matches. Many other named cultivars within lineages are too genetically similar to differentiate, highlighting difficulties of defining and identifying variation among clonally propagated triploid crops. Breadfruit is important in resilient agroforestry in tropical islands predicted to be especially affected by climate change. Findings reveal global links, building upon collective knowledge that can be used to inform breadfruit management. Results are also summarized in a brochure about breadfruit history and diversity in St. Vincent, and the Caribbean more broadly.

Capturing the antibiotic resistome of preterm infants reveals new benefits of probiotic supplementation.

BACKGROUND: Probiotic use in preterm infants can mitigate the impact of antibiotic exposure and reduce rates of certain illnesses; however, the benefit on the gut resistome, the collection of antibiotic resistance genes, requires further investigation. We hypothesized that probiotic supplementation of early preterm infants (born < 32-week gestation) while in hospital reduces the prevalence of antibiotic resistance genes associated with pathogenic bacteria in the gut. We used a targeted capture approach to compare the resistome from stool samples collected at the term corrected age of 40 weeks for two groups of preterm infants (those that routinely received a multi-strain probiotic during hospitalization and those that did not) with samples from full-term infants at 10 days of age to identify if preterm birth or probiotic supplementation impacted the resistome. We also compared the two groups of preterm infants up to 5 months of age to identify persistent antibiotic resistance genes. RESULTS: At the term corrected age, or 10 days of age for the full-term infants, we found over 80 antibiotic resistance genes in the preterm infants that did not receive probiotics that were not identified in either the full-term or probiotic-supplemented preterm infants. More genes associated with antibiotic inactivation mechanisms were identified in preterm infants unexposed to probiotics at this collection time-point compared to the other infants. We further linked these genes to mobile genetic elements and Enterobacteriaceae, which were also abundant in their gut microbiomes. Various genes associated with aminoglycoside and beta-lactam resistance, commonly found in pathogenic bacteria, were retained for up to 5 months in the preterm infants that did not receive probiotics. CONCLUSIONS: This pilot survey of preterm infants shows that probiotics administered after preterm birth during hospitalization reduced the diversity and prevented persistence of antibiotic resistance genes in the gut microbiome. The benefits of probiotic use on the microbiome and the resistome should be further explored in larger groups of infants. Due to its high sensitivity and lower sequencing cost, our targeted capture approach can facilitate these surveys to further address the implications of resistance genes persisting into infancy without the need for large-scale metagenomic sequencing. Video Abstract.

Targeted screening of genetic associations with COVID-19 susceptibility and severity.

The COVID-19 pandemic has resulted in great morbidity and mortality worldwide and human genetic factors have been implicated in the susceptibility and severity of COVID-19. However, few replicate researches have been performed, and studies on associated genes mainly focused on genic regions while regulatory regions were a lack of in-depth dissection. Here, based on previously reported associated variants and genes, we designed a capture panel covering 1,238 candidate variants and 25 regulatory regions of 19 candidate genes and targeted-sequenced 96 mild and 145 severe COVID-19 patients. Genetic association analysis was conducted between mild and severe COVID-19 patients, between all COVID-19 patients and general population, or between severe COVID-19 patients and general population. A total of 49 variants were confirmed to be associated with susceptibility or severity of COVID-19 (p < 0.05), corresponding to 18 independent loci. Specifically, rs1799964 in the promoter of inflammation-related gene TNF, rs9975538 in the intron of interferon receptor gene IFNAR2, rs429358 in the exon of APOE, rs1886814 in the intron of FOXP4-AS1 and a list of variants in the widely reported 3p21.31 and ABO gene were confirmed. It is worth noting that, for the confirmed variants, the phenotypes of the cases and controls were highly consistent between our study and previous reports, and the confirmed variants identified between mild and severe patients were quite different from those identified between patients and general population, suggesting the genetic basis of susceptibility and severity of SARS-CoV-2 infection might be quite different. Moreover, we newly identified 67 significant associated variants in the 12 regulatory regions of 11 candidate genes (p < 0.05). Further annotation by RegulomeDB database and GTEx eQTL data filtered out two variants (rs11246060 and rs28655829) in the enhancer of broad-spectrum antiviral gene IFITM3 that might affect disease severity by regulating the gene expression. Collectively, we confirmed a list of previously reported variants and identified novel regulatory variants associated with susceptibility and severity of COVID-19, which might provide biological and clinical insights into COVID-19 pathogenesis and treatment.