Rye: genetic ancestry inference at biobank scale.

Biobank projects are generating genomic data for many thousands of individuals. Computational methods are needed to handle these massive data sets, including genetic ancestry (GA) inference tools. Current methods for GA inference do not scale to biobank-size genomic datasets. We present Rye-a new algorithm for GA inference at biobank scale. We compared the accuracy and runtime performance of Rye to the widely used RFMix, ADMIXTURE and iAdmix programs and applied it to a dataset of 488221 genome-wide variant samples from the UK Biobank. Rye infers GA based on principal component analysis of genomic variant samples from ancestral reference populations and query individuals. The algorithm's accuracy is powered by Metropolis-Hastings optimization and its speed is provided by non-negative least squares regression. Rye produces highly accurate GA estimates for three-way admixed populations-African, European and Native American-compared to RFMix and ADMIXTURE (${R}^2 = \ 0.998 - 1.00$), and shows 50× runtime improvement compared to ADMIXTURE on the UK Biobank dataset. Rye analysis of UK Biobank samples demonstrates how it can be used to infer GA at both continental and subcontinental levels. We discuss user consideration and options for the use of Rye; the program and its documentation are distributed on the GitHub repository: https://github.com/healthdisparities/rye.

Multistate Outbreak of SARS-CoV-2 Infections, Including Vaccine Breakthrough Infections, Associated with Large Public Gatherings, United States.

During July 2021, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) B.1.617.2 variant infections, including vaccine breakthrough infections, occurred after large public gatherings in Provincetown, Massachusetts, USA, prompting a multistate investigation. Public health departments identified primary and secondary cases by using coronavirus disease surveillance data, case investigations, and contact tracing. A primary case was defined as SARS-CoV-2 detected <14 days after travel to or residence in Provincetown during July 3-17. A secondary case was defined as SARS-CoV-2 detected <14 days after close contact with a person who had a primary case but without travel to or residence in Provincetown during July 3-August 10. We identified 1,098 primary cases and 30 secondary cases associated with 26 primary cases among fully and non-fully vaccinated persons. Large gatherings can have widespread effects on SARS-CoV-2 transmission, and fully vaccinated persons should take precautions, such as masking, to prevent SARS-CoV-2 transmission, particularly during substantial or high transmission.

Native American admixture recapitulates population-specific migration and settlement of the continental United States.

European and African descendants settled the continental US during the 17th-19th centuries, coming into contact with established Native American populations. The resulting admixture among these groups yielded a significant reservoir of Native American ancestry in the modern US population. We analyzed the patterns of Native American admixture seen for the three largest genetic ancestry groups in the US population: African descendants, Western European descendants, and Spanish descendants. The three groups show distinct Native American ancestry profiles, which are indicative of their historical patterns of migration and settlement across the country. Native American ancestry in the modern African descendant population does not coincide with local geography, instead forming a single group with origins in the southeastern US, consistent with the Great Migration of the early 20th century. Western European descendants show Native American ancestry that tracks their geographic origins across the US, indicative of ongoing contact during westward expansion, and Native American ancestry can resolve Spanish descendant individuals into distinct local groups formed by more recent migration from Mexico and Puerto Rico. We found an anomalous pattern of Native American ancestry from the US southwest, which most likely corresponds to the Nuevomexicano descendants of early Spanish settlers to the region. We addressed a number of controversies surrounding this population, including the extent of Sephardic Jewish ancestry. Nuevomexicanos are less admixed than nearby Mexican-American individuals, with more European and less Native American and African ancestry, and while they do show demonstrable Sephardic Jewish ancestry, the fraction is no greater than seen for other New World Spanish descendant populations.

Population structure and pharmacogenomic risk stratification in the United States.

BACKGROUND: Pharmacogenomic (PGx) variants mediate how individuals respond to medication, and response differences among racial/ethnic groups have been attributed to patterns of PGx diversity. We hypothesized that genetic ancestry (GA) would provide higher resolution for stratifying PGx risk, since it serves as a more reliable surrogate for genetic diversity than self-identified race/ethnicity (SIRE), which includes a substantial social component. We analyzed a cohort of 8628 individuals from the United States (US), for whom we had both SIRE information and whole genome genotypes, with a focus on the three largest SIRE groups in the US: White, Black (African-American), and Hispanic (Latino). Our approach to the question of PGx risk stratification entailed the integration of two distinct methodologies: population genetics and evidence-based medicine. This integrated approach allowed us to consider the clinical implications for the observed patterns of PGx variation found within and between population groups. RESULTS: Whole genome genotypes were used to characterize individuals' continental ancestry fractions-European, African, and Native American-and individuals were grouped according to their GA profiles. SIRE and GA groups were found to be highly concordant. Continental ancestry predicts individuals' SIRE with > 96% accuracy, and accordingly, GA provides only a marginal increase in resolution for PGx risk stratification. In light of the concordance between SIRE and GA, taken together with the fact that information on SIRE is readily available to clinicians, we evaluated PGx variation between SIRE groups to explore the potential clinical utility of race and ethnicity. PGx variants are highly diverged compared to the genomic background; 82 variants show significant frequency differences among SIRE groups, and genome-wide patterns of PGx variation are almost entirely concordant with SIRE. The vast majority of PGx variation is found within rather than between groups, a well-established fact for almost all genetic variants, which is often taken to argue against the clinical utility of population stratification. Nevertheless, analysis of highly differentiated PGx variants illustrates how SIRE partitions PGx variation based on groups' characteristic ancestry patterns. These cases underscore the extent to which SIRE carries clinically valuable information for stratifying PGx risk among populations, albeit with less utility for predicting individual-level PGx alleles (genotypes), supporting the concept of population pharmacogenomics. CONCLUSIONS: Perhaps most interestingly, we show that individuals who identify as Black or Hispanic stand to gain far more from the consideration of race/ethnicity in treatment decisions than individuals from the majority White population.

Rapid Nanopore Whole-Genome Sequencing for Anthrax Emergency Preparedness.

Human anthrax cases necessitate rapid response. We completed Bacillus anthracis nanopore whole-genome sequencing in our high-containment laboratory from a human anthrax isolate hours after receipt. The de novo assembled genome showed no evidence of known antimicrobial resistance genes or introduced plasmid(s). Same-day genomic characterization enhances public health emergency response.

Ancestry effects on type 2 diabetes genetic risk inference in Hispanic/Latino populations.

BACKGROUND: Hispanic/Latino (HL) populations bear a disproportionately high burden of type 2 diabetes (T2D). The ability to predict T2D genetic risk using polygenic risk scores (PRS) offers great promise for improved screening and prevention. However, there are a number of complications related to the accurate inference of genetic risk across HL populations with distinct ancestry profiles. We investigated how ancestry affects the inference of T2D genetic risk using PRS in diverse HL populations from Colombia and the United States (US). In Colombia, we compared T2D genetic risk for the Mestizo population of Antioquia to the Afro-Colombian population of Chocó, and in the US, we compared European-American versus Mexican-American populations. METHODS: Whole genome sequences and genotypes from the 1000 Genomes Project and the ChocoGen Research Project were used for genetic ancestry inference and for T2D polygenic risk score (PRS) calculation. Continental ancestry fractions for HL genomes were inferred via comparison with African, European, and Native American reference genomes, and PRS were calculated using T2D risk variants taken from multiple genome-wide association studies (GWAS) conducted on cohorts with diverse ancestries. A correction for ancestry bias in T2D risk inference based on the frequencies of ancestral versus derived alleles was developed and applied to PRS calculations in the HL populations studied here. RESULTS: T2D genetic risk in Colombian and US HL populations is positively correlated with African and Native American ancestry and negatively correlated with European ancestry. The Afro-Colombian population of Chocó has higher predicted T2D risk than Antioquia, and the Mexican-American population has higher predicted risk than the European-American population. The inferred relative risk of T2D is robust to differences in the ancestry of the GWAS cohorts used for variant discovery. For trans-ethnic GWAS, population-specific variants and variants with same direction effects across populations yield consistent results. Nevertheless, the control for bias in T2D risk prediction confirms that explicit consideration of genetic ancestry can yield more reliable cross-population genetic risk inferences. CONCLUSIONS: T2D associations that replicate across populations provide for more reliable risk inference, and modeling population-specific frequencies of ancestral and derived risk alleles can help control for biases in PRS estimation.

GlobAl Distribution of GEnetic Traits (GADGET) web server: polygenic trait scores worldwide.

Human populations from around the world show striking phenotypic variation across a wide variety of traits. Genome-wide association studies (GWAS) are used to uncover genetic variants that influence the expression of heritable human traits; accordingly, population-specific distributions of GWAS-implicated variants may shed light on the genetic basis of human phenotypic diversity. With this in mind, we developed the GlobAl Distribution of GEnetic Traits web server (GADGET http://gadget.biosci.gatech.edu). The GADGET web server provides users with a dynamic visual platform for exploring the relationship between worldwide genetic diversity and the genetic architecture underlying numerous human phenotypes. GADGET integrates trait-implicated single nucleotide polymorphisms (SNPs) from GWAS, with population genetic data from the 1000 Genomes Project, to calculate genome-wide polygenic trait scores (PTS) for 818 phenotypes in 2504 individual genomes. Population-specific distributions of PTS are shown for 26 human populations across 5 continental population groups, with traits ordered based on the extent of variation observed among populations. Users of GADGET can also upload custom trait SNP sets to visualize global PTS distributions for their own traits of interest.

A Nationwide Screen of Carbapenem-Resistant Klebsiella pneumoniae Reveals an Isolate with Enhanced Virulence and Clinically Undetected Colistin Heteroresistance.

The convergence of hypervirulence and multidrug resistance in Klebsiella pneumoniae is a significant concern. Here, we report the first screen for hypermucoviscosity, a trait associated with increased virulence, using a U.S. surveillance collection of carbapenem-resistant (CR) K. pneumoniae isolates. We identified one hypermucoviscous isolate, which carried a gene encoding the KPC-3 carbapenemase, among numerous resistance genes. The strain further exhibited colistin heteroresistance undetected by diagnostics. This convergence of diverse resistance mechanisms and increased virulence underscores the need for enhanced K. pneumoniae surveillance.

Genetic ancestry, admixture and health determinants in Latin America.

BACKGROUND: Modern Latin American populations were formed via genetic admixture among ancestral source populations from Africa, the Americas and Europe. We are interested in studying how combinations of genetic ancestry in admixed Latin American populations may impact genomic determinants of health and disease. For this study, we characterized the impact of ancestry and admixture on genetic variants that underlie health- and disease-related phenotypes in population genomic samples from Colombia, Mexico, Peru, and Puerto Rico. RESULTS: We analyzed a total of 347 admixed Latin American genomes along with 1102 putative ancestral source genomes from Africans, Europeans, and Native Americans. We characterized the genetic ancestry, relatedness, and admixture patterns for each of the admixed Latin American genomes, finding a spectrum of ancestry proportions within and between populations. We then identified single nucleotide polymorphisms (SNPs) with anomalous ancestry-enrichment patterns, i.e. SNPs that exist in any given Latin American population at a higher frequency than expected based on the population's genetic ancestry profile. For this set of ancestry-enriched SNPs, we inspected their phenotypic impact on disease, metabolism, and the immune system. All four of the Latin American populations show ancestry-enrichment for a number of shared pathways, yielding evidence of similar selection pressures on these populations during their evolution. For example, all four populations show ancestry-enriched SNPs in multiple genes from immune system pathways, such as the cytokine receptor interaction, T cell receptor signaling, and antigen presentation pathways. We also found SNPs with excess African or European ancestry that are associated with ancestry-specific gene expression patterns and play crucial roles in the immune system and infectious disease responses. Genes from both the innate and adaptive immune system were found to be regulated by ancestry-enriched SNPs with population-specific regulatory effects. CONCLUSIONS: Ancestry-enriched SNPs in Latin American populations have a substantial effect on health- and disease-related phenotypes. The concordant impact observed for same phenotypes across populations points to a process of adaptive introgression, whereby ancestry-enriched SNPs with specific functional utility appear to have been retained in modern populations by virtue of their effects on health and fitness.

Genome-wide map of Apn1 binding sites under oxidative stress in Saccharomyces cerevisiae.

The DNA is cells is continuously exposed to reactive oxygen species resulting in toxic and mutagenic DNA damage. Although the repair of oxidative DNA damage occurs primarily through the base excision repair (BER) pathway, the nucleotide excision repair (NER) pathway processes some of the same lesions. In addition, damage tolerance mechanisms, such as recombination and translesion synthesis, enable cells to tolerate oxidative DNA damage, especially when BER and NER capacities are exceeded. Thus, disruption of BER alone or disruption of BER and NER in Saccharomyces cerevisiae leads to increased mutations as well as large-scale genomic rearrangements. Previous studies demonstrated that a particular region of chromosome II is susceptible to chronic oxidative stress-induced chromosomal rearrangements, suggesting the existence of DNA damage and/or DNA repair hotspots. Here we investigated the relationship between oxidative damage and genomic instability utilizing chromatin immunoprecipitation combined with DNA microarray technology to profile DNA repair sites along yeast chromosomes under different oxidative stress conditions. We targeted the major yeast AP endonuclease Apn1 as a representative BER protein. Our results indicate that Apn1 target sequences are enriched for cytosine and guanine nucleotides. We predict that BER protects these sites in the genome because guanines and cytosines are thought to be especially susceptible to oxidative attack, thereby preventing large-scale genome destabilization from chronic accumulation of DNA damage. Information from our studies should provide insight into how regional deployment of oxidative DNA damage management systems along chromosomes protects against large-scale rearrangements. Copyright © 2017 John Wiley & Sons, Ltd.

Compound cis-regulatory elements with both boundary and enhancer sequences in the human genome.

MOTIVATION: It has been suggested that presumably distinct classes of genomic regulatory elements may actually share common sets of features and mechanisms. However, there has been no genome-wide assessment of the prevalence of this phenomenon. RESULTS: To evaluate this possibility, we performed a bioinformatic screen for the existence of compound regulatory elements in the human genome. We identified numerous such colocated boundary and enhancer elements from human CD4(+) T cells. We report evidence that such compound regulatory elements possess unique chromatin features and facilitate cell type-specific functions related to inflammation and immune response in CD4(+) T cells.

Epigenetic regulation of human cis-natural antisense transcripts.

Mammalian genomes encode numerous cis-natural antisense transcripts (cis-NATs). The extent to which these cis-NATs are actively regulated and ultimately functionally relevant, as opposed to transcriptional noise, remains a matter of debate. To address this issue, we analyzed the chromatin environment and RNA Pol II binding properties of human cis-NAT promoters genome-wide. Cap analysis of gene expression data were used to identify thousands of cis-NAT promoters, and profiles of nine histone modifications and RNA Pol II binding for these promoters in ENCODE cell types were analyzed using chromatin immunoprecipitation followed by sequencing (ChIP-seq) data. Active cis-NAT promoters are enriched with activating histone modifications and occupied by RNA Pol II, whereas weak cis-NAT promoters are depleted for both activating modifications and RNA Pol II. The enrichment levels of activating histone modifications and RNA Pol II binding show peaks centered around cis-NAT transcriptional start sites, and the levels of activating histone modifications at cis-NAT promoters are positively correlated with cis-NAT expression levels. Cis-NAT promoters also show highly tissue-specific patterns of expression. These results suggest that human cis-NATs are actively transcribed by the RNA Pol II and that their expression is epigenetically regulated, prerequisites for a functional potential for many of these non-coding RNAs.

Multiplex Real-time PCR Assay for the Detection of all Chlamydia Species and Simultaneous Differentiation of C. psittaci and C. pneumoniae in Human Clinical Specimens.

We developed and assessed the performance of a new multiplex real-time PCR assay for the detection of all Chlamydia species and simultaneous differentiation of Chlamydia psittaci and Chlamydia pneumoniae-two important human respiratory pathogens-in human clinical specimens. Next-generation sequencing was used to identify unique targets to design real-time PCR assays targeting all Chlamydia species, C. psittaci, and C. pneumoniae. To validate the assay, we used a panel of 49 culture isolates comprising seven C. psittaci genotypes, eight C. pneumoniae isolates, seven other Chlamydia species, and 22 near-neighbor bacterial and viral isolates, along with 22 specimens from external quality assessment (EQA) panels and 34 nasopharyngeal and oropharyngeal swabs and cerebrospinal fluid, stool, and sputum specimens previously identified as positive or negative for C. psittaci or C. pneumoniae. The assays were 100% specific, with limits of detection of 7.64- 9.02 fg/µL. The assay results matched with historical assay results for all specimens, except for one owing to the increased sensitivity of the new C. psittaci assay; the results of the EQA specimens were 100% accurate. This assay may improve the timely and accurate clinical diagnosis of Chlamydia infections and provide a greater understanding of the burden of disease caused by these agents.

Ethnic disparities in mortality and group-specific risk factors in the UK Biobank.

Despite a substantial overall decrease in mortality, disparities among ethnic minorities in developed countries persist. This study investigated mortality disparities and their associated risk factors for the three largest ethnic groups in the United Kingdom: Asian, Black, and White. Study participants were sampled from the UK Biobank (UKB), a prospective cohort enrolled between 2006 and 2010. Genetics, biological samples, and health information and outcomes data of UKB participants were downloaded and data-fields were prioritized based on participants with death registry records. Kaplan-Meier method was used to evaluate survival differences among ethnic groups; survival random forest feature selection followed by Cox proportional-hazard modeling was used to identify and estimate the effects of shared and ethnic group-specific mortality risk factors. The White ethnic group showed significantly worse survival probability than the Asian and Black groups. In all three ethnic groups, endoscopy and colonoscopy procedures showed significant protective effects on overall mortality. Asian and Black women show lower relative risk of mortality than men, whereas no significant effect of sex was seen for the White group. The strongest ethnic group-specific mortality associations were ischemic heart disease for Asians, COVID-19 for Blacks, and cancers of respiratory/intrathoracic organs for Whites. Mental health-related diagnoses, including substance abuse, anxiety, and depression, were a major risk factor for overall mortality in the Asian group. The effect of mental health on Asian mortality, particularly for digestive cancers, was exacerbated by an observed hesitance to answer mental health questions, possibly related to cultural stigma. C-reactive protein (CRP) serum levels were associated with both overall and cause-specific mortality due to COVID-19 and digestive cancers in the Black group, where elevated CRP has previously been linked to psychosocial stress due to discrimination. Our results point to mortality risk factors that are group-specific and modifiable, supporting targeted interventions towards greater health equity.

Genomic analysis of Chlamydia psittaci from a multistate zoonotic outbreak in two chicken processing plants.

Four Chlamydia psittaci isolates were recovered from clinical specimens from ill workers during a multistate outbreak at two chicken processing plants. Whole genome sequencing analyses revealed high similarity to C. psittaci genotype D. The isolates differed from each other by only two single nucleotide polymorphisms, indicating a common source.

Investigation of multidrug-resistant plasmids from carbapenemase-producing Klebsiella pneumoniae clinical isolates from Pakistan.

Objectives: The study aim was to investigate multidrug-resistant (MDR) plasmids from a collection of 10 carbapenemase-producing Klebsiella pneumoniae clinical isolates identified within the same healthcare institution in Pakistan. Full characterization of the MDR plasmids including structure, typing characteristics, and AMR content as well as determination of their plasmid-based antimicrobial susceptibility profiles were carried out. Methods: Plasmids were isolated from 10 clinical isolates of Klebsiella pneumoniae, and from a corresponding set of Escherichia coli transconjugants, then sequenced using Nanopore/Illumina technology to generate plasmid hybrid assemblies. Full characterization of MDR plasmids, including determination of next generation sequencing (NGS)-based AMR profiles, plasmid incompatibility groups, and types, was carried out. The structure of MDR plasmids was analyzed using the Galileo AMR platform. For E. coli transconjugants, the NGS-based AMR profiles were compared to NGS-predicted AMR phenotypes and conventional broth microdilution (BMD) antimicrobial susceptibility testing (AST) results. Results: All carbapenemase-producing K. pneumoniae isolates (carrying either blaNDM-1, or/and blaOXA-48) carried multiple AMR plasmids encoding 34 antimicrobial resistance genes (ARGs) conferring resistance to antimicrobials from 6 different classes. The plasmid incompatibility groups and types identified were: IncC (types 1 and 3), IncFIA (type 26) IncFIB, IncFII (types K1, K2, K7, and K9), IncHI1B, and IncL. None of the blaNDM-1 and blaESBL-plasmids identified in this study were previously described. Most blaNDM-1-plasmids shared identical AMR regions suggesting potential genetic material/plasmid exchange between K. pneumoniae isolates of this collection. The majority of NGS-based AMR profiles from the E. coli transconjugants correlated well with both NGS-based predicted and conventional AST results. Conclusion: This study highlights the complexity and diversity of the plasmid-based genetic background of carbapenemase-producing clinical isolates from Pakistan. This study emphasizes the need for characterization of MDR plasmids to determine their complete molecular background and monitor AMR through plasmid transmission between multi-resistant bacterial pathogens.

Effects of genetic ancestry and socioeconomic deprivation on ethnic differences in serum creatinine.

The inclusion of ethnicity in equations for estimating the glomerular filtration rate (eGFR) from serum creatinine levels has been challenged since ethnicity is socially defined and therefore a poor proxy for biological differences. We hypothesized that genetic ancestry (GA) would be more strongly associated with creatinine levels among healthy individuals than self-identified ethnicity. We studied a diverse cohort of 35,590 participants characterized as part of the UK Biobank, grouped by self-reported ethnicity: Black, East Asian, Mixed, Other, South Asian, and White. We used multivariable modeling to test for associations between ethnicity, GA, socioeconomic deprivation, and serum creatinine levels, including covariates for age, sex, height, and body mass index. Model fit comparisons and relative importance analysis were used to compare the effects of ethnicity and GA on creatinine levels. Black ethnicity shows a positive effect on participant serum creatinine levels (ß = 9.36 ± 0.38), whereas East Asian (ß = -1.80 ± 0.66) and South Asian (ß = -0.28 ± 0.36) ethnicity show negative effects on creatinine. Male sex (ß = 17.69 ± 0.34) and height (ß = 0.13 ± 0.02) also show high positive associations with creatinine levels, while socioeconomic deprivation (ß = -0.04 ± 0.04) shows no significant association. African ancestry has the highest association (ß = 13.81 ± 0.52) with creatinine levels. Overall, GA (9.06%) explains significantly more of the variation in creatinine levels than ethnicity (4.96%), with African ancestry (6.36%) alone explaining more of the variation than ethnicity. We found that GA explains more of the variation in serum creatinine levels than socioeconomic deprivation, suggesting the possibility that ethnic differences in creatinine are shaped by genetic rather than social factors.

Genome sequences for five strains of the emerging pathogen Haemophilus haemolyticus.

We report the first whole-genome sequences for five strains, two carried and three pathogenic, of the emerging pathogen Haemophilus haemolyticus. Preliminary analyses indicate that these genome sequences encode markers that distinguish H. haemolyticus from its closest Haemophilus relatives and provide clues to the identity of its virulence factors.

Human cis natural antisense transcripts initiated by transposable elements.

The capacity of human transposable elements (TEs) to promote cis natural antisense transcripts (cis-NATs) is revealed by the discovery of 48,718 human gene antisense transcriptional start sites (TSSs) within TE sequences. TSSs that yield cis-NATs are overrepresented among TE sequences, and TE-initiated cis-NATs are more abundant close to the 3' ends of genes. The TE sequences that promote antisense transcription within human genes are relatively ancient, suggesting that selection has acted to conserve their function.

A computational genomics pipeline for prokaryotic sequencing projects.

MOTIVATION: New sequencing technologies have accelerated research on prokaryotic genomes and have made genome sequencing operations outside major genome sequencing centers routine. However, no off-the-shelf solution exists for the combined assembly, gene prediction, genome annotation and data presentation necessary to interpret sequencing data. The resulting requirement to invest significant resources into custom informatics support for genome sequencing projects remains a major impediment to the accessibility of high-throughput sequence data. RESULTS: We present a self-contained, automated high-throughput open source genome sequencing and computational genomics pipeline suitable for prokaryotic sequencing projects. The pipeline has been used at the Georgia Institute of Technology and the Centers for Disease Control and Prevention for the analysis of Neisseria meningitidis and Bordetella bronchiseptica genomes. The pipeline is capable of enhanced or manually assisted reference-based assembly using multiple assemblers and modes; gene predictor combining; and functional annotation of genes and gene products. Because every component of the pipeline is executed on a local machine with no need to access resources over the Internet, the pipeline is suitable for projects of a sensitive nature. Annotation of virulence-related features makes the pipeline particularly useful for projects working with pathogenic prokaryotes. AVAILABILITY AND IMPLEMENTATION: The pipeline is licensed under the open-source GNU General Public License and available at the Georgia Tech Neisseria Base (http://nbase.biology.gatech.edu/). The pipeline is implemented with a combination of Perl, Bourne Shell and MySQL and is compatible with Linux and other Unix systems.