Evaluation of SARS-CoV-2 ORF7a Deletions from COVID-19-Positive Individuals and Its Impact on Virus Spread in Cell Culture.

The spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), causing the COVID-19 outbreak, posed a primary concern of public health worldwide. The most common changes in SARS-CoV-2 are single nucleotide substitutions, also reported insertions and deletions. This work investigates the presence of SARS-CoV-2 ORF7a deletions identified in COVID-19-positive individuals. Sequencing of SARS-CoV-2 complete genomes showed three different ORF7a size deletions (190-nt, 339-nt and 365-nt). Deletions were confirmed through Sanger sequencing. The ORF7a∆190 was detected in a group of five relatives with mild symptoms of COVID-19, and the ORF7a∆339 and ORF7a∆365 in a couple of co-workers. These deletions did not affect subgenomic RNAs (sgRNA) production downstream of ORF7a. Still, fragments associated with sgRNA of genes upstream of ORF7a showed a decrease in size when corresponding to samples with deletions. In silico analysis suggests that the deletions impair protein proper function; however, isolated viruses with partial deletion of ORF7a can replicate in culture cells similarly to wild-type viruses at 24 hpi, but with less infectious particles after 48 hpi. These findings on deleted ORF7a accessory protein gene, contribute to understanding SARS-CoV-2 phenotypes such as replication, immune evasion and evolutionary fitness as well insights into the role of SARS-CoV-2_ORF7a in the mechanism of virus-host interactions.

Genetic diversity of the melanocortin-1 receptor in an admixed population of Rio de Janeiro: Structural and functional impacts of Cys35Tyr variant.

The melanocortin-1 receptor (MC1R) is one of the key proteins involved in the regulation of melanin production and several polymorphisms have been associated with different phenotypes of skin and hair color in human and nonhuman species. Most of the knowledge is centered on more homogeneous populations and studies involving an admixed group of people should be encouraged due to the great importance of understanding the human color variation. This work evaluates the MC1R diversity and the possible impacts of MC1R variants in an admixed sample population of Rio de Janeiro, Brazil, which is a product of Native American, African, and European miscegenation. Sequencing of complete coding region and part of the 3´UTR of MC1R gene identified 31 variants including one insertion and three novel synonymous substitutions in sample population grouped according to skin, hair and eye pigmentation levels. In nonmetric multidimensional scaling analysis (NMDS), three main clusters were identified, in which the Brazilian dark skin group remained in the African cluster whereas the intermediate and the light skin color phenotype in the European one. None gathered with Asians since their immigration to Brazil was a recent event. In silico analyses demonstrated that Cys35Tyr, Ile155Thr and Pro256Ser, found in our population, have a negative effect on receptor function probably due to changes on the receptor structure. Notably, Cys35Tyr mutation could potentially impair agonist binding. Altogether, this work contributes to the understanding of the genetic background of color variation on an admixed population and gives insights into the damaging effects of MC1R variants.

Frequency distribution of HCV resistance-associated variants in infected patients treated with direct-acting antivirals.

BACKGROUND: Hepatitis C virus (HCV) is a global public health problem. Second-generation direct-acting antivirals targeting non-structural regions on the viral genome are the cornerstone for treatment of chronic infection. However, resistance-associated variants (RAVs) have been reported to be associated with therapeutic failure. The aim of this study was to assess the frequency of variants, including RAVs, in the NS3, NS5A and NS5B regions at baseline in Brazilian patients with chronic hepatitis C with HCV genotypes 1a, 1b and 3a. METHODS: Serum samples from 13 patients were used to obtain viral RNA. Massively parallel sequencing was performed using genotype-specific amplicons and a panel of Ampliseq technology for all genotypes. RESULTS: Several non-synonymous substitutions were detected at baseline for 11 responders and pre-/post-treatment for two non-responders. HCV genotype 3a was found to have significantly more non-synonymous substitutions than HCV genotype 1 in the NS3 and NS5A regions. Analyses were conducted using quantitative and qualitative inter- and intrapatient comparisons. Variants that confer resistance to the treatment used by the patients were found in both responders and non-responders. CONCLUSIONS: A wide frequency distribution of RAVs was found at baseline, and this did not interfere with the achievement of a sustained response. Evaluation of the presence of RAVs requires additional study in order to determine clinical relevance.

Evaluation of Promega PowerSeq™ Auto/Y systems prototype on an admixed sample of Rio de Janeiro, Brazil: Population data, sensitivity, stutter and mixture studies.

Forensic DNA typing typically relies on the length-based (LB) separation of PCR products containing short tandem repeat loci (STRs). Massively parallel sequencing (MPS) elucidates an additional level of STR motif and flanking region variation. Also, MPS enables simultaneous analysis of different marker-types - autosomal STRs, SNPs for lineage and identification purposes, reducing both the amount of sample used and the turn-around-time of analysis. Therefore, MPS methodologies are being considered as an additional tool in forensic genetic casework. The PowerSeq™ Auto/Y System (Promega Corp), a multiplex forensic kit for MPS, enables analysis of the 22 autosomal STR markers (plus Amelogenin) from the PowerPlex® Fusion 6C kit and 23 Y-STR markers from the PowerPlex® Y23 kit. Population data were generated from 140 individuals from an admixed sample from Rio de Janeiro, Brazil. All samples were processed according to the manufacturers' recommended protocols. Raw data (FastQ) were generated for each indexed sample and analyzed using STRait Razor v2s and PowerSeqv2.config file. The subsequent population data showed the largest increase in expected heterozygosity (23%), from LB to sequence-based (SB) analyses at the D5S818 locus. Unreported allele was found at the D21S11 locus. The random match probability across all loci decreased from 5.9 × 10-28 to 7.6 × 10-33. Sensitivity studies using 1, 0.25, 0.062 and 0.016 ng of DNA input were analyzed in triplicate. Full Y-STR profiles were detected in all samples, and no autosomal allele drop-out was observed with 62 pg of input DNA. For mixture studies, 1 ng of genomic DNA from a male and female sample at 1:1, 1:4, 1:9, 1:19 and 1:49 proportions were analyzed in triplicate. Clearly resolvable alleles (i.e., no stacking or shared alleles) were obtained at a 1:19 male to female contributor ratio. The minus one stutter (-1) increased with the longest uninterrupted stretch (LUS) allele size reads and according to simple or compound/complex repeats. The haplotype-specific stutter rates add more information for mixed samples interpretation. These data support the use of the PowerSeqTM Auto/Y systems prototype kit (22 autosomal STR loci, 23 Y-STR loci and Amelogenin) for forensic genetics applications.

STRait Razor Online: An enhanced user interface to facilitate interpretation of MPS data.

Since 2013, STRait Razor has enabled analysis of massively parallel sequencing (MPS) data from various marker systems such as short tandem repeats, single nucleotide polymorphisms, insertion/deletions, and mitochondrial DNA. In this paper, STRait Razor Online (SRO), available at https://www.unthsc.edu/straitrazor, is introduced as an interactive, Shiny-based user interface for primary analysis of MPS data and secondary analysis of STRait Razor haplotype pileups. This software can be accessed from any common browser via desktop, tablet, or smartphone device. SRO is available also as a standalone application and open-source R script available at https://github.com/ExpectationsManaged/STRaitRazorOnline. The local application is capable of batch processing of both fastq files and primary analysis output. Processed batches generate individual report folders and summary reports at the locus- and haplotype-level in a matter of minutes. For example, the processing of data from ∼700 samples generated with the ForenSeq Signature Preparation Kit from allsequences.txt to a final table can be performed in ∼40 min whereas the Excel-based workbooks can take 35-60 h to compile a subset of the tables generated by SRO. To facilitate analysis of single-source, reference samples, a preliminary triaging system was implemented that calls potential alleles and flags loci suspected of severe heterozygote imbalance. When compared to published, manually curated data sets, 98.72 % of software-assigned allele calls without manual interpretation were consistent with curated data sets, 0.99 % loci were presented to the user for interpretation due to heterozygote imbalance, and the remaining 0.29 % of loci were inconsistent due to the analytical thresholds used across the studies.

Evaluation of 16S rRNA Hypervariable Regions for Bioweapon Species Detection by Massively Parallel Sequencing.

Molecular detection and classification of the bacterial groups in a sample are relevant in several areas, including medical research and forensics. Sanger sequencing of the 16S rRNA gene is considered the gold standard for microbial phylogenetic analysis. However, the development of massively parallel sequencing (MPS) offers enhanced sensitivity and specificity for microbiological analyses. In addition, 16S rRNA target amplification followed by MPS facilitates the combined use of multiple markers/regions, better discrimination of sample background, and higher sample throughput. We designed a novel set of 16S rRNA gene primers for detection of bacterial species associated with clinical, bioweapon, and biohazards microorganisms via alignment of 364 sequences representing 19 bacterial species and strains relevant to medical and forensics applications. In silico results indicated that the hypervariable regions (V1V2), (V4V5), and (V6V7V8) support the resolution of a selected group of bacteria. Interspecies and intraspecies comparisons showed 74.23%-85.51% and 94.48%-99.98% sequencing variation among species and strains, respectively. Sequence reads from a simulated scenario of bacterial species mapped to each of the three hypervariable regions of the respective species with different affinities. The minimum limit of detection was achieved using two different MPS platforms. This protocol can be used to detect or monitor as low as 2,000 genome equivalents of bacterial species associated with clinical, bioweapon, and biohazard microorganisms and potentially can distinguish natural outbreaks of pathogenic microorganisms from those occurring by intentional release.

A novel phylogenetic approach for de novo discovery of putative nuclear mitochondrial (pNumt) haplotypes.

Current approaches for parsing true variation (i.e. signal) from noise, broadly involve estimating a baseline value of the latter, below which all sequence data are ignored. In an effort to deliver a more objective criterion for setting such thresholds, a novel approach based on phylogenetic principles is presented here., Our method deconstructs a special category of noise from true mitochondrial genome data, namely nuclear insertions of mitochondrial DNA (Numts). This bioinformatic approach leverages the relationship of massively parallel sequence reads and is capable of discovering putative Numts (pNumts) in absence of a reference genome. The new method was tested on a whole mitochondrial genome dataset (n = 41 individuals from an admixed population sample from Rio de Janeiro) and led to the discovery of 451 pNumt variants. Comparison of these pNumts haplotypes against an existing Numt database revealed 147 exact matches to previously discovered Numts, while 122 haplotypes differed only by a single base pair and none matched exclusively to the mitochondrial genome. In general, these sequences were considerably more divergent from the mitochondrial genome than from those of the Numt database, supporting that the novel pNumts were probably hitherto uncatalogued variants. Unlike previous techniques, our method appears to be able to detect both polymorphic and fixed Numt sequences. It was also found that the region containing the D-Loop and associated Promoters (DLP) in the human mitochondrial genome, which harbors markers of forensic genetics importance, is the origin of several Numts. Though currently designed for the mitochondrial genome, our novel approach has the potential to be expanded to other scenarios that might require construing signal from noise, including the deconvolution of mixtures, thus significantly improving how analytical thresholds may be established.

Evaluation of mitogenome sequence concordance, heteroplasmy detection, and haplogrouping in a worldwide lineage study using the Precision ID mtDNA Whole Genome Panel.

The emergence of Massively Parallel Sequencing technologies enabled the analysis of full mitochondrial (mt)DNA sequences from forensically relevant samples that have, so far, only been typed in the control region or its hypervariable segments. In this study, we evaluated the performance of a commercially available multiplex-PCR-based assay, the Precision ID mtDNA Whole Genome Panel (Thermo Fisher Scientific), for the amplification and sequencing of the entire mitochondrial genome (mitogenome) from even degraded forensic specimens. For this purpose, more than 500 samples from 24 different populations were selected to cover the vast majority of established superhaplogroups. These are known to harbor different signature sequence motifs corresponding to their phylogenetic background that could have an effect on primer binding and, thus, could limit a broad application of this molecular genetic tool. The selected samples derived from various forensically relevant tissue sources and were DNA extracted using different methods. We evaluated sequence concordance and heteroplasmy detection and compared the findings to conventional Sanger sequencing as well as an orthogonal MPS platform. We discuss advantages and limitations of this approach with respect to forensic genetic workflow and analytical requirements.

17 Y-STR haplotype diversity in São Paulo state (southeast of Brazil).

A sample of 158 Brazilian males from São Paulo (SP), Brazilian southeast, was typed for 17 Y-STR loci (DYS19, DYS389I, DYS389II, DYS390, DYS391, DYS392, DYS393, DYS437, DYS438, DYS439, DYS448, DYS456, DYS458, DYS635, YGATA_H4.1, and DYS385ab). A total of 158 haplotypes were identified, of which all were unique. The haplotype diversity and discrimination capacity were calculated in 1.0 and the genetic diversity was 67.4%. Pairwise haplotype distances showed that the São Paulo population is not significantly different from Rio de Janeiro and Portugal, but is different from African and Native American.

Evaluation of the precision ID mtDNA whole genome panel on two massively parallel sequencing systems.

Sequencing whole mitochondrial genomes by capillary electrophoresis is a costly and time/labor-intensive endeavor. Many of the previous Sanger sequencing-based approaches generated amplicons that were several kilobases in length; lengths that are likely not amenable for most forensic applications. However, with the advent of massively parallel sequencing (MPS) short-amplicon multiplexes covering the entire mitochondrial genome can be sequenced relatively easily and rapidly. Recently, the Precision ID mtDNA Whole Genome Panel (Thermo Fisher Scientific by Applied Biosystems™) has been introduced. This panel is composed of 162 amplicons (in two multiplexes) that are considerably smaller in length (∼163bp) and thus are more amenable to analyzing challenged samples. This panel was evaluated on both the Ion S5™ System (Thermo Fisher Scientific) and the MiSeq™ FGx Desktop Sequencer (Illumina). A script was developed to extract phased haplotypes associated with these amplicons. Levels of read-depth were compared across sequencing pools and between sequencing technologies and haplotype concordances were assessed. Given modest thresholds on read depth, the haplotypes identified by either technology were consistent. Nuclear mitochondrial sequences (Numts) were also inferred, and the effect of different mapping strategies commonly used to filter out Numts were contrasted. Some Numts are co-amplified with this amplification kit, and while the choice of reference sequence can mitigate some of these effects, some data from the mitochondrial genome were lost in the process in this study. This study demonstrates that the Ion and MiSeq platforms provide consistent haplotype estimation of the whole mitochondrial genome, thus providing further support for the reliability and validity of the Precision ID mtDNA Whole Genome Panel.

Planktonic microbial profiling in water samples from a Brazilian Amazonian reservoir.

Our comprehension of the dynamics and diversity of freshwater planktonic bacterial communities is far from complete concerning the Brazilian Amazonian region. Therefore, reference studies are urgently needed. We mapped bacterial communities present in the planktonic communities of a freshwater artificial reservoir located in the western Amazonian basin. Two samples were obtained from rainy and dry seasons, the periods during which water quality and plankton diversity undergo the most significant changes. Hypervariable 16S rRNA and shotgun sequencing were performed to describe the first reference of a microbial community in an Amazonian lentic system. Microbial composition consisted mainly of Betaproteobacteria, Cyanobacteria, Alphaproteobacteria, and Actinobacteria in the dry period. The bacteria distribution in the rainy period was notably absent of Cyanobacteria. Microcystis was observed in the dry period in which the gene cluster for cyanotoxins was found. Iron acquisition gene group was higher in the sample from the rainy season. This work mapped the first inventory of the planktonic microbial community of a large water reservoir in the Amazon, providing a reference for future functional studies and determining other communities and how they interact.

Green Tobacco Sickness among Brazilian farm workers and genetic polymorphisms.

OBJECTIVE: Green Tobacco Sickness (GTS) is an occupational illness caused by dermal absorption of nicotine from tobacco leaves. It affects thousands of farm workers worldwide. Brazil is the second tobacco producer in the world; despite this, there are few studies on GTS among Brazilian harvesters. This study aimed to determine the prevalence of GTS among a population of tobacco workers from a producing area in northeastern Brazil and investigate whether the occurrence of the disease was influenced by factors such age, gender and smoking status. In addition, it was investigated if there was association between the onset of GTS and genetic polymorphisms in genes that encode some detoxification enzymes. A semi-structured questionnaire was used to collect demographic, behavioral and occupational data from the referred workers. Polymorphisms were tested through the Polymerase Chain Reaction technique. RESULTS: The total prevalence of GTS found was 56.9%, with a significant difference between genders (71.7% for women and 35.3% for men, p < 0.0001). No association was identified between the investigated polymorphisms and GTS. This study confirms the occurrence of GTS among tobacco harvesters in Brazil with high prevalence. The investigation suggests the need to take preventive measures to protect tobacco workers against this disease.

Genetic data for 26 autosomal STR markers from Brazilian population.

The allelic frequency distributions and statistical forensic parameters of 26 mini short tandem repeat (mini-STR) loci in a sample of 1575 unrelated individuals from five different Brazilian regions were obtained. All the analyzed loci showed great diversity and were highly informative. The results were compared with those of the US Caucasian, African American, and Hispanic population studies. This study aimed to contribute to forensic analysis for human identification and inference of the evidential value in familial bond tests.

Full-gene haplotypes refine CYP2D6 metabolizer phenotype inferences.

CYP2D6 is a critical pharmacogenetic target, and polymorphisms in the gene region are commonly used to infer enzyme activity score and predict resulting metabolizer phenotype: poor, intermediate, extensive/normal, or ultrarapid which can be useful in determining cause and/or manner of death in some autopsies. Current genotyping approaches are incapable of identifying novel and/or rare variants, so CYP2D6 star allele definitions are limited to polymorphisms known a priori. While useful for most predictions, recent studies using massively parallel sequencing data have identified additional polymorphisms in CYP2D6 that are predicted to alter enzyme function but are not considered in current star allele nomenclature. The 1000 Genomes Project data were used to produce full-gene haplotypes, describe their distribution in super-populations, and predict enzyme activity scores. Full-gene haplotypes generated lower activity scores than current approaches due to inclusion of additional damaging polymorphisms in the star allele. These findings are critical for clinical implementation of metabolizer phenotype prediction because a fraction of the population may be incorrectly considered normal metabolizers but actually may be poor or intermediate metabolizers.

Are your results valid? Cellular authentication a need from the past, an emergency on the present.

The cultures of immortalized cells have been established in the 50s and become popular as a biological model for in vitro assays. The success and popularization brought side effects. Still, in the 60 years emerge the first cases of misidentification/contamination of cell line. Because of that, the scientific community has been oriented to authenticate their lines before performing assays. The use of cells with incorrect identification or contamination has been identified as responsible for an increasing number of unmatched results and a waste of resources. For this reason, we implemented the Cell Line Authentication Service at Brazilian Metrology Institute (Inmetro), open to Brazilian scientific community and society in general. From 2012 to 2014 were conducted 111 cell line authentication test, of which 13.8% had some problem. Here are the description and discussion of these data and simple guidelines to minimize the risk of contamination and misidentification, and invite the scientific community to maintain an alert system to avoid spending unnecessary resources and produce unreliable data.

Hybridization Capture-Based Next-Generation Sequencing to Evaluate Coding Sequence and Deep Intronic Mutations in the NF1 Gene.

Neurofibromatosis 1 (NF1) is one of the most common genetic disorders and is caused by mutations in the NF1 gene. NF1 gene mutational analysis presents a considerable challenge because of its large size, existence of highly homologous pseudogenes located throughout the human genome, absence of mutational hotspots, and diversity of mutations types, including deep intronic splicing mutations. We aimed to evaluate the use of hybridization capture-based next-generation sequencing to screen coding and noncoding NF1 regions. Hybridization capture-based next-generation sequencing, with genomic DNA as starting material, was used to sequence the whole NF1 gene (exons and introns) from 11 unrelated individuals and 1 relative, who all had NF1. All of them met the NF1 clinical diagnostic criteria. We showed a mutation detection rate of 91% (10 out of 11). We identified eight recurrent and two novel mutations, which were all confirmed by Sanger methodology. In the Sanger sequencing confirmation, we also included another three relatives with NF1. Splicing alterations accounted for 50% of the mutations. One of them was caused by a deep intronic mutation (c.1260 + 1604A > G). Frameshift truncation and missense mutations corresponded to 30% and 20% of the pathogenic variants, respectively. In conclusion, we show the use of a simple and fast approach to screen, at once, the entire NF1 gene (exons and introns) for different types of pathogenic variations, including the deep intronic splicing mutations.

Common NOD2/CARD15 and TLR4 Polymorphisms Are Associated with Crohn's Disease Phenotypes in Southeastern Brazilians.

AIM: To investigate whether variants in NOD2/CARD15 and TLR4 are associated with CD and ulcerative colitis (UC) in a genetically admixed population of Rio de Janeiro, where IBD has continued to rise. METHODS: We recruited 67 consecutive patients with CD, 61 patients with UC, and 86 healthy and ethnically matched individuals as controls. DNA was extracted from buccal brush samples and genotyped by PCR with restriction enzymes for G908R and L1007finsC NOD2/CARD15 single-nucleotide polymorphisms (SNPs) and for T399I and D299G TLR4 SNPs. Clinical data were registered for subsequent analysis with multivariate models. RESULTS: NOD2/CARD15 G908R and L1007finsC SNPs were found in one and three patients, respectively, with CD. NOD2/CARD15 G908R and L1007finsC SNPs were not found in any patients with UC, but were found in three and three controls, respectively. With regard to the TLR4 gene, no significant difference was detected among the groups. Overall, none of the SNPs investigated determined a differential risk for a specific diagnosis. Genotype-phenotype associations were found in only CD, where L1007finsC was associated with colonic localization; however, TLR4 T399I SNP was associated with male gender, and D299G SNP was associated with colonic involvement, chronic corticosteroid use, and the need for anti-TNF-alpha therapy. CONCLUSION: Variants of NOD2/CARD15 and TLR4 do not confer susceptibility to IBD, but appear to determine CD phenotypes in this southeastern Brazilian population.

Selection of highly informative SNP markers for population affiliation of major US populations.

Ancestry informative markers (AIMs) can be used to detect and adjust for population stratification and predict the ancestry of the source of an evidence sample. Autosomal single nucleotide polymorphisms (SNPs) are the best candidates for AIMs. It is essential to identify the most informative AIM SNPs across relevant populations. Several informativeness measures for ancestry estimation have been used for AIMs selection: absolute allele frequency differences (δ), F statistics (F ST), and informativeness for assignment measure (In). However, their efficacy has not been compared objectively, particularly for determining affiliations of major US populations. In this study, these three measures were directly compared for AIMs selection among four major US populations, i.e., African American, Caucasian, East Asian, and Hispanic American. The results showed that the F ST panel performed slightly better for population resolution based on principal component analysis (PCA) clustering than did the δ panel and both performed better than the In panel. Therefore, the 23 AIMs selected by the F ST measure were used to characterize the four major American populations. Genotype data of nine sample populations were used to evaluate the efficiency of the 23-AIMs panel. The results indicated that individuals could be correctly assigned to the major population categories. Our AIMs panel could contribute to the candidate pool of AIMs for potential forensic identification purposes.