Sanger Sequencing Data Research Articles

Scifer: An R/Bioconductor package for large-scale integration of Sanger sequencing and flow cytometry data of index-sorted single cells

Sanger sequencing remains widely used in various experimental contexts, often in combination with flow cytometry for indexing specific cell populations. However, existing software lacks the capability to automate quality control (QC) of raw Sanger sequencing data and integrate it with flow cytometry information on a large scale. Here, we introduce scifer, an R package now available in the latest release of Bioconductor (3.20) showcasing its effectiveness in seamlessly integrating these types of data as demonstrated by analyses of B cell and T cell receptor sequences. Scifer preprocesses raw data from index sorts and immune receptor Sanger sequencing. It identifies high-quality sequences based on selected parameters, such as length, Phred scores, and heavy-chain complementarity-determining region 3 (HCDR3) quality. As a result, the quality of germline assignments is significantly increased and spurious variable gene mutations are reduced. Scifer is automated and can process thousands of sequences in less than an hour. Its output provides quality control reports, FASTA files, summarized tables, and electropherograms for manual inspection. In summary, scifer is a user-friendly software that speeds up the analysis of immune receptor repertoire sequences, offering wide applicability.

HerpesDRG: a comprehensive resource for human herpesvirus antiviral drug resistance genotyping

The prevention and treatment of many herpesvirus associated diseases is based on the utilization of antiviral therapies, however therapeutic success is limited by the development of drug resistance. Currently no single database cataloguing resistance mutations exists, which hampers the use of sequence data for patient management. We therefore developed HerpesDRG, a drug resistance mutation database that incorporates all the known resistance genes and current treatment options, built from a systematic review of available genotype to phenotype literature. The database is released along with an R package that provides a simple approach to resistance variant annotation and clinical implication analysis from common sanger and next generation sequencing data. This represents the first openly available and community maintainable database of drug resistance mutations for the human herpesviruses (HHV), developed for the community of researchers and clinicians tackling HHV drug resistance.

Accurate and Automated Genotyping of the CFTR Poly-T/TG Tract with CFTR-TIPS.

Cystic fibrosis is caused by biallelic pathogenic variants in the CFTR gene, which contains a polymorphic (TG)mTn sequence (the "poly-T/TG tract") in intron 9. While T9 and T7 alleles are benign, T5 alleles with longer TG repeats, e.g., (TG)12T5 and (TG)13T5, are clinically significant. Thus, professional medical societies currently recommend reporting the TG repeat size when T5 is detected. Sanger sequencing is a cost-effective method of genotyping the (TG)mTn tract; however, its polymorphic length substantially complicates data analysis. We developed CFTR-TIPS, a freely available web-based software tool that infers the (TG)mTn genotype from Sanger sequencing data. This tool detects the (TG)mTn tract in the chromatograms, quantifies goodness of fit with expected patterns, and visualizes the results in a graphical user interface. It is broadly compatible with any Sanger chromatogram that contains the (TG)mTn tract ± 15 bp. We evaluated CFTR-TIPS using 835 clinical samples previously analyzed in a CLIA-certified, CAP-accredited laboratory. When operated fully automatically, CFTR-TIPS achieved 99.8% concordance with our clinically validated manual workflow, while generally taking less than 10 s per sample. There were two discordant samples: one due to a co-occurring heterozygous duplication that confounded the tool and the other due to incomplete (TG)mTn tract detection in the reverse chromatogram. No clinically significant misclassifications were observed. CFTR-TIPS is a free, accurate, and rapid tool for CFTR (TG)mTn tract genotyping using cost-effective Sanger sequencing. This tool is suitable both for automated use and as an aid to manual review to enhance accuracy and reduce analysis time.

UCE-based phylogenomics of the lepidopteran endoparasitoid wasp subfamily Rogadinae (Hymenoptera: Braconidae) unveils a new Neotropical tribe.

During the past two decades, the phylogenetic relationships and higher-level classification of the subfamily Rogadinae have received relevant contributions based on Sanger, mitogenome and genome-wide nuclear DNA sequence data. These studies have helped to update the circumscription and tribal classification of this subfamily, with six tribes currently recognised (Aleiodini, Betylobraconini, Clinocentrini, Rogadini, Stiropiini and Yeliconini). The tribal relationships within Rogadinae, however, are yet to be fully resolved, including the status of tribe Facitorini, previously regarded as betylobraconine, with respect to the members of Yeliconini. We conducted a phylogenomic analysis among the tribes of Rogadinae based on genomic ultraconserved element (UCE) data and extensive taxon sampling including three undescribed genera of uncertain tribal placement. Our almost fully supported estimate of phylogeny confirmed the basal position of Rogadini within the subfamily and a Facitorini clade (Yeliconini+Aleiodini) that led us to propose the former group as a valid rogadine tribe (Facitorini stat. res.). Stiropiini, however, was recovered for the first time as sister to the remaining rogadine tribes except Rogadini, and Clinocentrini as sister to a clade with Betylobraconini+the three undescribed genera. The relationships recovered and morphological examination of the material included led us to place the latter three new genera and recently described genus Gondwanocentrus within a new rogadine tribe, Gondwanocentrini Shimbori & Zaldívar-Riverón trib. nov. We described these genera (Ghibli Shimbori & Zaldívar-Riverón gen. nov., Racionais Shimbori & Zaldívar-Riverón gen. nov. and Soraya Shimbori gen. nov.) with two or three new species each (G. miyazakii Shimbori & Zaldívar-Riverón sp. nov., G. totoro Shimbori & Zaldívar-Riverón sp. nov., R. brunus Shimbori & Zaldívar-Riverón sp. nov., R. kaelejay Shimbori & Zaldívar-Riverón sp. nov., R. superstes Shimbori & Zaldívar-Riverón sp. nov., S. alencarae Shimbori sp. nov. and S. venus Shimbori & Zaldívar-Riverón sp. nov.). A new species of Facitorini, Jannya pasargadae Gadelha & Shimbori sp. nov., is also described. Our newly proposed classification expands the number of tribes and genera within Rogadinae to 8 and 66 respectively. ZooBank: urn:lsid:zoobank.org:pub:51951C78-069A-4D8B-B5F0-7EBD4D9D21CE.

Population-based nanopore sequencing of the HIV-1 pangenome to identify drug resistance mutations.

HIV-1 drug resistance genotypic tests have primarily been performed by Sanger sequencing of gene segments encoding different drug target proteins. Since the number of targets has increased with the addition of a new class of antiretroviral drugs, a simple high-throughput system for assessing nucleotide sequences throughout the HIV-1 genome is required. Here, we developed a new solution using nanopore sequencing of viral pangenomes amplified by PCR. Benchmark tests using HIV-1 molecular clones demonstrated an accuracy of up to 99.9%. In addition, validation tests of our protocol in 106 clinical samples demonstrated high concordance of drug resistance and tropism genotypes (92.5% and 98.1%, respectively) between the nanopore sequencing-based results and archived clinical determinations made based on Sanger sequencing data. These results suggest that our new approach will be a powerful solution for the comprehensive survey of HIV-1 drug resistance mutations in clinical settings.

Phylogenomics-driven host test list selection for weed biological control

A concern in weed biological control research is the potential for candidate biocontrol agents to impact not only the target weed but also native or economically important flora. The degree of evolutionary relatedness between the target weed species and a non-target species is a key predictor of the susceptibility of the non-target to the biocontrol agent. To manage this risk, biocontrol practitioners need to understand the phylogenetic position of the host weed relative to non-target plant species. However, comprehensively sampled phylogenetic trees are often unavailable, with incomplete information scattered across multiple publications. Further, older published phylogenies based on Sanger sequence data often lack branch resolution and support, which increases uncertainty in biocontrol decision making. Decreasing sequencing cost and technological advances have led to phylogenomic approaches being more widely used to understand evolutionary relationships between species. For example, target capture sequencing methods using bait kits such as Angiosperms353 enable cost-effective and timely phylogenomic-level analysis of flowering plant groups at different scales. Here, we introduce a workflow to embed a comprehensive understanding of evolutionary relationships into the efficient development of host test lists in weed biological control. We demonstrate the effectiveness of the workflow through a case study on the major crop weed flaxleaf fleabane (Erigeron bonariensis). Phylogenomic analysis was conducted on 280 species of the tribe Astereae (family Asteraceae) occurring in Australia and New Zealand, clarifying relationships between the target species and related clades of native and non-native Astereae. We consider the phylogenetic tree in the context of a previously proposed host test list and discuss taxonomic implications, highlighting avenues of future molecular-based work to uncover the origin of Australian fleabanes. This study provides a workflow and demonstrates the practical application of target sequence capture for phylogenomic inference to support risk analysis and decision making in classical weed biological control.

Mitogenome architecture supports the non-monophyly of the cosmopolitan parasitoid wasp subfamily Doryctinae (Hymenoptera: Braconidae) recovered by nuclear and mitochondrial phylogenomics.

Mitochondrial DNA gene organisation is an important source of phylogenetic information for various metazoan taxa at different evolutionary timescales, though this has not been broadly tested for all insect groups nor within a phylogenetic context. The cosmopolitan subfamily Doryctinae is a highly diverse group of braconid wasps mainly represented by ectoparasitoids of xylophagous beetle larvae. Previous molecular studies based on Sanger and genome-wide (ultraconserved elements, UCE; and mitochondrial genomes) sequence data have recovered a non-monophyletic Doryctinae, though the relationships involved have always been weakly supported. We characterised doryctine mitogenomes and conducted separate phylogenetic analyses based on mitogenome and UCE sequence data of ~100 representative doryctine genera to assess the monophyly and higher-level classification of the subfamily. We identified rearrangements of mitochondrial transfer RNAs (tRNAs) that support a non-monophyletic Doryctinae consisting of two separate non-related clades with strong geographic structure ('New World' and 'Old World' clades). This geographic structure was also consistently supported by the phylogenetic analyses preformed with mitogenome and UCE sequence data. These results highlight the utility of the mitogenome gene rearrangements as a potential source of phylogenetic information at different evolutionary timescales.

Evolutionary Characterization of tubulin Gene Family in the Desert Biomass Willow (Salix psammophila) and Expression of the β-tubulin Gene SpsTUB10 during Different Stresses

Microtubules, polymerized from α-tubulin (TUA) and β-tubulin (TUB) monomers, play a pivotal role in shaping plant morphogenesis according to developmental and environmental cues. Salix psammophila C. Wang & C. Y. Yang is an important shrub plant in sand-fixing afforestation in arid regions, with three significantly distinct plant types shaped under various environments, namely, upright, intermediate, and scattered types. However, how tubulin genes respond to the developmental and environmental signs in S. psammophila has been far less studied. Here, based on RNA-seq, Sanger sequencing, and real-time PCR (RT-PCR) data, we analyzed the phylogeny of tubulins and their expression profiles in S. psammophila among the three plant types. Furthermore, we analyzed the genetic structure and expression pattern of SpsTUB10 in S. psammophila under various abiotic stress treatments. In total, we identified 26 SpsTubulin genes in S. psammophila. The homologous alignment and phylogenetic analysis revealed that these SpsTubulin genes can be classified into two groups, corresponding to the TUA and TUB genes. The expression profiles of these SpsTubulin genes in various organs showed that most SpsTubulin genes were mainly expressed in the root. SpsTUB10 is a member of the TUB IIa group, consisting of two intros and three exons. The SpsTUB10 protein contains a typical GTPase domain and a C-terminal domain, with α-helix and random coil dominant in the secondary and tertiary structures. The RT-PCR results of SpsTUB10 showed an extremely significant difference in expression levels among the root and stem-developing organs between the upright and scattered types, and the transcript level of SpsTUB10 had a significantly negative correlation with the crown-height ratio. Under different treatments, we found that cold, osmotic stress, and short daylight could significantly increase SpsTUB10 expression levels compared to those in the controls, thereby supporting the positive role of SpsTUB10 in stress-induced responses. These results will provide evidence for the SpsTubulin genes’ response to the developmental and environmental cues in S. psammophila.

VirusPredictor: XGBoost-based software to predict virus-related sequences in human data.

Discovering disease causative pathogens, particularly viruses without reference genomes, poses a technical challenge as they are often unidentifiable through sequence alignment. Machine learning prediction of patient high-throughput sequences unmappable to human and pathogen genomes may reveal sequences originating from uncharacterized viruses. Currently, there is a lack of software specifically designed for accurately predicting such viral sequences in human data. We developed a fast XGBoost method and software VirusPredictor leveraging an in-house viral genome database. Our two-step XGBoost models first classify each query sequence into one of three groups: infectious virus, endogenous retrovirus (ERV) or non-ERV human. The prediction accuracies increased as the sequences became longer, i.e. 0.76, 0.93, and 0.98 for 150-350 (Illumina short reads), 850-950 (Sanger sequencing data), and 2000-5000 bp sequences, respectively. Then, sequences predicted to be from infectious viruses are further classified into one of six virus taxonomic subgroups, and the accuracies increased from 0.92 to >0.98 when query sequences increased from 150-350 to >850 bp. The results suggest that Illumina short reads should be de novo assembled into contigs (e.g. ∼1000 bp or longer) before prediction whenever possible. We applied VirusPredictor to multiple real genomic and metagenomic datasets and obtained high accuracies. VirusPredictor, a user-friendly open-source Python software, is useful for predicting the origins of patients' unmappable sequences. This study is the first to classify ERVs in infectious viral sequence prediction. This is also the first study combining virus sub-group predictions. www.dllab.org/software/VirusPredictor.html.

The origin and speciation of orchids.

Orchids constitute one of the most spectacular radiations of flowering plants. However, their origin, spread across the globe, and hotspots of speciation remain uncertain due to the lack of an up-to-date phylogeographic analysis. We present a new Orchidaceae phylogeny based on combined high-throughput and Sanger sequencing data, covering all five subfamilies, 17/22 tribes, 40/49 subtribes, 285/736 genera, and c. 7% (1921) of the 29 524 accepted species, and use it to infer geographic range evolution, diversity, and speciation patterns by adding curated geographical distributions from the World Checklist of Vascular Plants. The orchids' most recent common ancestor is inferred to have lived in Late Cretaceous Laurasia. The modern range of Apostasioideae, which comprises two genera with 16 species from India to northern Australia, is interpreted as relictual, similar to that of numerous other groups that went extinct at higher latitudes following the global climate cooling during the Oligocene. Despite their ancient origin, modern orchid species diversity mainly originated over the last 5 Ma, with the highest speciation rates in Panama and Costa Rica. These results alter our understanding of the geographic origin of orchids, previously proposed as Australian, and pinpoint Central America as a region of recent, explosive speciation.

Two complete 1918 influenza A/H1N1 pandemic virus genomes characterized by next-generation sequencing using RNA isolated from formalin-fixed, paraffin-embedded autopsy lung tissue samples along with evidence of secondary bacterial co-infection.

This study applied modern molecular biotechnology and high-throughput sequencing to formalin-fixed, paraffin-embedded autopsy lung samples from two fatal cases during the fall wave of the 1918 influenza A (H1N1) pandemic in the United States. Complete influenza genomes were obtained from both cases, which increases the total number of available complete or near-complete influenza genomes of the 1918 pandemic virus from four to six. Sequence analysis confirms that the 1918 pandemic virus was highly conserved during the main wave of the pandemic with geographic separation in North America and Europe. Metagenomic analyses revealed bacterial co-infections in both cases, including the first reported evidence of Rhodococcus-influenza co-infection. Overall, this study offers a detailed view at the molecular level of the very limited samples from the most devastating influenza pandemic in modern human history.

Phylogenetic affinities of the non-cyclostome subfamilies Amicrocentrinae and Dirrhopinae (Hymenoptera, Braconidae) confirmed by ultraconserved element data

The subfamilies Amicrocentrinae and Dirrhopinae (Hymenoptera, Braconidae) are two small, monogeneric braconid subfamilies whose species exclusively attack lepidopteran larvae. The phylogenetic placement of Amicrocentrinae as a member of the helconoid complex of subfamilies has been supported by morphological and nuclear Sanger sequence data. The subfamilial status of Dirrhopinae on the other hand has been subject to debate, although it has been suggested as closely related to the microgastroid complex based on morphology only. Here we generated for the first time genomic ultraconserved element data for members of the above subfamilies (Amicrocentrum seyrigi van Achterberg and Dirrhope americana Muesebeck) to assess their phylogenetic affinities using exhaustive taxon sampling that includes all but one of the currently valid braconid subfamilies. Our results strongly confirm the placement of both taxa within the non-cyclostome helconoid and microgastroid complexes, respectively.

Comparison of methylation estimates obtained via MinION nanopore sequencing and sanger bisulfite sequencing in the TRPA1 promoter region

BackgroundBisulfite sequencing has long been considered the gold standard for measuring DNA methylation at single CpG resolution. However, in recent years several new approaches like nanopore sequencing have been developed due to hints for a partial error-proneness of bisulfite sequencing. Since these errors were shown to be sequence-specific, we aimed to verify the methylation data of a particular region of the TRPA1 promoter from our previous studies obtained by bisulfite sequencing.MethodsWe compared methylation rates determined by direct bisulfite sequencing and nanopore sequencing following Cas9-mediated PCR-free enrichment.ResultsWe could show that CpG methylation levels above 20% corroborate well with our previous data. Within the range between 0 and 20% methylation, however, Sanger sequencing data have to be interpreted cautiously, at least in the investigated region of interest (TRPA1 promotor region).ConclusionBased on the investigation of the TRPA1- region as an example, the present work can help in choosing the right method out of the two current main approaches for methylation analysis for different individual settings regarding many factors like cohort size, costs and prerequisites that should be fulfilled for each method. All in all, both methods have their raison d’être. Furthermore, the present paper contains and illustrates some important basic information and explanation of how guide RNAs should be located for an optimal outcome in Cas9 mediated PCR free target enrichment.

PHFinder: assisted detection of point heteroplasmy in Sanger sequencing chromatograms.

Heteroplasmy is the presence of two or more organellar genomes (mitochondrial or plastid DNA) in an organism, tissue, cell or organelle. Heteroplasmy can be detected by visual inspection of Sanger sequencing chromatograms, where it appears as multiple peaks of fluorescence at a single nucleotide position. Visual inspection of chromatograms is both consuming and highly subjective, as heteroplasmy is difficult to differentiate from background noise. Few software solutions are available to automate the detection of point heteroplasmies, and those that are available are typically proprietary, lack customization or are unsuitable for automated heteroplasmy assessment in large datasets. Here, we present PHFinder, a Python-based, open-source tool to assist in the detection of point heteroplasmies in large numbers of Sanger chromatograms. PHFinder automatically identifies point heteroplasmies directly from the chromatogram trace data. The program was tested with Sanger sequencing data from 100 humpback whales (Megaptera novaeangliae) tissue samples with known heteroplasmies. PHFinder detected most (90%) of the known heteroplasmies thereby greatly reducing the amount of visual inspection required. PHFinder is flexible and enables explicit specification of key parameters to infer double peaks (i.e., heteroplasmies).

First Report of Tomato Spotted Wilt Virus infecting Lettuce in Yuma, Arizona.

Tomato spotted wilt virus (TSWV, family Tospoviridae, genus Orthotospovirus) is a thrips-vectored pathogen that infects lettuce (Lactuca sativa) and many vegetable crops (Kuo et al. 2014, Hasegawa et al. 2022). Another thrips-borne pathogen of lettuce, impatiens necrotic spot virus (INSV, Tospoviridae, Orthotospovirus), was first reported in 2021 in Yuma, Arizona (Hasegawa et al. 2022). Symptoms of both viruses in lettuce are similar and include necrotic spotting, leaf chlorosis and plant stunting (Kuo et al. 2014). Beginning February through April of 2022, lettuce displaying symptoms of orthotospovirus infection was collected from romaine lettuce (var. longifolia) fields in three regions of Yuma County. A total of 96 plants were collected (5 from Tacna on 2/21, 5 from Wellton on 2/21, 15 from Wellton on 3/23, 30 from Tacna on 4/4, 20 from Wellton on 4/4, and 21 from Yuma Valley on 4/4). The area of the fields ranged from 10 to 18 acres, and the percent disease incidence ranged from 0.8% (Tacna on 4/4) to 2.75% (Tacna on 2/21). Thrips vector were present in all fields were symptomatic plants were observed. One leaf disk per plant (8 mm in diameter) was sampled with a cork borer and grounded individually with a micro pestle in a 1.7 ml microcentrifuge tube with 150 ul of Tri-reagent (Molecular Research Center). Total RNA was extracted from each sample using the Zymo Direct-zol-96 kit (Zymo Research). Samples were diluted with water to a ratio of 1:10 after RNA extraction. RT-qPCR was performed in 20 ul reactions with 5 ul of input RNA using the PCR Biosystems qPCRBIO Probe 1-Step Go No-ROX for the cDNA/qPCR master mix. RT-qPCR assays were carried out in multiplex reactions using primers specific for TSWV and INSV, in addition to a lettuce internal control gene (LOC111918243), along with negative controls. Primer and probe sequence details are reported in supplemental Table 1. We used a cycle threshold (ct) < 40 to indicate a positive result for both INSV and TSWV (Chen et al. 2013; Boonham et al. 2002). RT-qPCR successfully amplified INSV in 90 out of 96 samples and TSWV in 8 out of 96 samples. These 8 samples tested positive for both TSWV and INSV, showing that INSV and TSWV co-infected lettuce plants. Thus overall, ∼ 95% of symptomatic plants were infected with INSV alone, and ∼ 8% were co-infected with TSWV and INSV. Amplicons of 4 samples testing positive for TSWV were sent for Sanger sequencing (Eurofins Genomics, Louisville, KY). All were identified as TSWV. One amplicon with TSWV was sequenced for INSV and double infection was confirmed. BLAST results from the NCBI nt database show 100% (138 bp) identity to TWSV (MW519211) for the 4 TWSV amplicons and 99.22% (137 bp) identity to INSV (KX790323) for the INSV amplicon. Sanger sequence data are in the GenBank (accession: OQ685940-OQ685944). Based on RT-qPCR results, all TSWV infected plants were also infected with INSV. INSV may have been introduced to Yuma by infected plants or thrips from lettuce transplants produced in California (Hasegawa et al. 2022). TSWV could have been introduced similarly. To our knowledge, this is the first report of TSWV infecting lettuce in Yuma and the first report of INSV and TSWV co-infecting lettuce. TSWV and INSV infections have remained low since their discovery in Yuma, in part due to effective cultural and chemical management by lettuce growers (Palumbo, 2022). However, an increase in disease incidence and severity in the future could have a significant negative impact on production of romaine lettuce in the region.

Analysis of genetic mutation in dilated cardiomyopathy

BackgroundCardiomyopathies (CM) are heterogeneous group of genetic disorders, specifically damage the heart muscles that leads to the progression of many clinical symptoms such as congestive heart failure (CHF). Dilated Cardiomyopathy (DCM) is the most common cause of heart failure. This study focuses on genetic variations in DCM, which is chronic, generally nonreversible heart muscles disease characterized by heart chamber dilation mostly left ventricle and decreased systolic efficiency. MethodologyThis study includes two parts, First we identified 100 genes associated with DCM using deep literature survey. Genetic variants were extracted from 1000 Genome Project (1000 GP), ExAC and Simon Genome Diversity project, for the identification of deleterious variants using in silico prediction tools followed by population genetic differentiation. In Second part, we examined 21,469 patients related to cardiac diseases. Of these, 100 DCM patients were recruited for the identification of rare pathogenic variants rs143187236 in MYOM3 and an INDEL rs36212066 in MYBPC3 using ARMS-PCR and validation through Sanger sequencing. ResultsWe identified 1658, 15,058 and 461 deleterious variants in 1000 GP, ExAC and Simon Genomes project respectively predicted by both SIFT and PolyPhen2 prediction tools. Clinvar database revealed 10, 69 and one pathogenic and likely pathogenic variants associated with DCM in 1000 GP, ExAC and Simon genome. Population genetic differentiation analysis revealed South Asian populations are closely related to European and American while highly differentiated to East Asian and African populations. ARMS-PCR and Sanger sequencing data of enlisted DCM patients showed that rs143187236 in MYOM3 was found in heterozygous state in three patients while rs36212066 in MYBPC3 was not detected in studied patients. ConclusionBased on bioinformatics analysis and investigation of deleterious variants through Sanger sequencing, it is deduced that large sample size is necessary for further detection of variants that will aid in the early diagnosis of DCM.

Evolutionary history, traits, and weediness inDigitaria(Poaceae: Panicoideae)

AbstractDigitaria is a large pantropical genus, which includes a number of economically problematic agricultural weeds. Difficulties in species identification and the circumscription of the genus have previously hindered progress in understanding its evolution and developing a stable classification. We investigate the evolutionary history of Digitaria by combining two approaches: (i) phylogenetic analyses of the currently available Sanger sequence data with previously unpublished sequences mainly from South East Asia, and (ii) newly generated genomic data from low-copy nuclear genes obtained using the Angiosperms353 probe kit, with new data mainly from Madagascar. Both methods gave congruent results. The combination of these results allowed us to confirm a new clade structure within Digitaria lineages, in partial agreement with the previous morphology-based classifications. We confirm that Digitaria as traditionally circumscribed is not monophyletic, comprising Anthephora, Baptorhachis, Chaetopoa, and Chlorocalymma. We present an inference of key character evolution with morphological support for clades with spikelets arranged in binate and ternate groups, with variable trichome morphology, spikelet length, and culm length. Geographical clade structure is also observed. No phylogenetic signal of weedy function was observed.

Sequencing Analysis of MUC6 and MUC16 Gene Fragments in Patients with Oropharyngeal Squamous Cell Carcinoma Reveals Novel Mutations: A Preliminary Study.

The growing incidence of oropharyngeal squamous cell carcinoma (OPSCC) calls for better understanding of the mutational landscape of such cases. Mucins (MUCs) are multifunctional glycoproteins expressed by the epithelial cells and may be associated with the epithelial tumour invasion and progression. The present study aimed at the analysis of the sequence of selected MUC6 and MUC16 gene fragments in the tumour, as well as the margin, samples obtained from 18 OPSCC patients. Possible associations between the detected mutations and the clinicopathological and demographic characteristics of the study group were analysed. Sanger sequencing and bioinformatic data analysis of the selected MUC6 and MUC16 cDNA fragments were performed. Our study found 13 and 3 mutations in MUC6 and MUC16, respectively. In particular, one novelty variant found that the MUC6 gene (chr11:1018257 A>T) was the most frequent across our cohort, in both the tumour and the margin samples, and was then classified as a high impact, stop-gain mutation. The current study found novel mutations in MUC6 and MUC16 providing new insight into the genetic alternation in mucin genes among the OPSCC patients. Further studies, including larger cohorts, are recommended to recognise the pattern in which the mutations affect oropharyngeal carcinogenesis.

Sequence Variants in the WNT10B Underlying Non-Syndromic Split-Hand/Foot Malformation

Introduction: Split hand and foot malformation (SHFM) or ectrodactyly is a rare limb deformity characterized by median cleft of the hand and foot with impaired or missing central rays. It can occur as an isolated anomaly or in association with abnormalities of other body parts. Methods: After delineating the clinical features of two families (A–B), with non-syndromic SHFM, exome and Sanger sequencing were employed to search for the disease-causing variants. Results: Analysis of exome and Sanger sequencing data revealed two causative variants in the WNT10B gene in affected members of the two families. This included a novel missense change [c.338G>C; p.(Gly113Ala)] in family A and a previously reported frameshift variant [c.884-896delTCCAGCCCCGTCT; p.(Phe295Cysfs*87)] in family B. Conclusion: Our findings add a novel variant in WNT10B gene as the underlying cause of SHFM. The finding adds to the growing body of knowledge about the genetic basis of developmental disorders and provides valuable insights into the molecular mechanisms that regulate limb development.

Interlaboratory study to assess the practical utility of OSAC proposed standard 2021-S-0006: Standard for the use of GenBank for taxonomic assignment of wildlife

Wildlife forensic laboratories are frequently asked to confirm the taxonomic identity of material submitted as evidence. When taxonomic assignment based on morphology is not possible, DNA-based approaches are typically employed and practitioners often rely on comparing unknown evidence sequences to public sequence databases such as GenBank. To address the use of GenBank in casework, the Organization of Scientific Area Committees for Forensic Science (OSAC) Wildlife Forensic Biology Subcommittee developed 2021-S-0006 Standard for the Use of GenBank for Taxonomic Assignment of Wildlife (currently an OSAC proposed standard). An interlaboratory study was conducted to assess the practical utility of 2021-S-0006. Sanger sequence data for taxonomically informative mitochondrial loci from forty diverse wildlife species were generated. Eleven wildlife forensic laboratories across the globe were recruited, and were sent ten unknowns to analyze using 2021-S-0006. Taxonomic assignments generated by each laboratory were compared a) to the true identity to determine correctness, and b) across laboratories to assess congruence. When following 2021-S-0006, laboratories that reported to the species level were correct in 98.3% of cases and 100% congruence was observed among laboratories. Additionally, criteria selected by laboratories for interpreting GenBank accession(s) (section 4.3 of 2021-S-0006) were independently evaluated, given they form the backbone for determining the appropriate level for taxonomic reporting. Criteria were correctly selected in 84% of cases and laboratories reported to an appropriate taxonomic level in the majority of cases. Overall, this interlaboratory study demonstrated that the OSAC proposed standard 2021-S-0006 provides an excellent framework for using GenBank for taxonomic assignment in wildlife forensic casework.