A band analysis

This tutorial provides simple instructions on performing a new analysis with MyoProfiler. Clicking on any of the images on this page will open a larger version in a new browser window. The preparation shown here was probed for myosin heavy chain and an associated protein.

Getting started

Using the MyoProfiler through MATLAB:
- Launch MATLAB and navigate to the Apps tab on the top menu. Find the MyoProfiler under My Apps and start the application by clicking it. The instructions on how to locate the Apps tab can be found here.
Using the MyoProfiler as a stand-alone application:
- Locate your MyoProfiler.exe shortcut and start the application by double-clicking it.

After a few seconds, you should see the program window, given below.

User panels

The interface is divided into multiple panels:

Controls: The user controls are accessed in this panel:
- Channel Colormap: Users can either use the emission wavelengths of each channel or a generic colormap for pseudo coloring and trace colors.
- ROI: The regions of interest (ROI) is selected through this panel. By default, the ROIs are derived from splines, but users can expand them into quadrilaterals to average around the spline.
- Profile: The software uses a peak finding algorithm to identify the crescents and troughs in the profile. Users can change the default values through this panel. More detail on these options can be found below.
Image Display: Loaded images are displayed in this panel. Channel tabs and an additional tab for the merged image is automatically generated.
Profiler Panel: Extracted profiles are shown in this panel. Since the profiles are two-dimensional, there are additional displays for “views” from horizontal and vertical image axes.
Sarcomere Panel: Analyzed sarcomere profiles are shown in this panel. Sarcomere length and full half max width are calculated and tabulated.

Loading images

Users start their analysis by loading their images to the software environment. They can either use the commercial image formats, such as Nikon originated “.ND2” files, or standard image formats, for instance “.TIF” or “.PNG” files.

The MyoProfiler is built to support Nikon images. The image support will be extended in the future and the tutorial will be updated as needed.

Loading Nikon .ND2 images

The Load Image option is located under the Menu on the top left corner, shown in the red rectangle. Click on the Menu and access to the dropdown menu.

As you hover over the Load Image option, additional options appear. Click on the Nikon option, highlighted in the red rectangle.

The Nikon option opens the below file dialog. The dialog automatically looks for the .ND2 files. Navigate to your folder with the image files and click Open.

.ND2 files packages data and information from utilized channels under a single file. MyoProfiler directly reads into the files and extract taken images from individual channels, pixels to length scale (um) calibration, and emission wavelengths.

Loaded images appear under the Image Display. The channel tabs are automatically generated. Below image shows the image from the Channel 1, shown in red rectangle.

Also, note that the pixels to length scale (micron) calibration is loaded from the loaded file.

Below image shows the image from the Channel 2, shown in red rectangle.

Although the software uses these monochromatic images for analysis, a merged image with pseudo coloring is provided for users to review and qualitative analysis. The pseudo coloring is based on the emitted wavelengths, users can change the coloring using the Channel Colormap dropdown, shown in red rectangle.

You can hover over the images to reveal the toolbar to Zoom In, Zoom Out, and Reset View. The toolbar is shown in the red rectangle.

A zoomed in view is useful for qualitative localization and judiciously select ROIs.

Loading standard format images

The Load Image option is located under the Menu on the top left corner, shown in the red rectangle. Click on the Menu and access to the dropdown menu.

As you hover over the Load Image option, additional options appear. Click on the Standard Formats option, highlighted in the red rectangle.

Once the Standard Formats option is clicked, the software opens a dialog box, shown below in red rectangle, for users to enter number of channels in their analysis. Although optional, users can provide the emission wavelengths.

Click OK after the fields are filled out.

In this tutorial the user provided that there are 2 images with 670 and 490 nm emission wavelengths, respectively for Channel 1 and 2. Please note that these image are exported from Nikon example.

MyoProfiler asks for users to load the image for Channel 1 first. Navigate to your folder with the image files and click Open.

Then, repeat it for the same for Channel 2.

The channel tabs are automatically generated. Below image shows the image from the Channel 1, shown in red rectangle.

Below image shows the image from the Channel 2, shown in red rectangle.

The pseudo colored merged image appears as red and cyan, instead of green and orange, with respect to the provided wavelengths.

Please note that in the Standard Formats case, users are expected to provide the calibration constant from pixels to microns.

Although the rest of the tutorial uses the loaded Nikon image, users can get the same results using the standard format images.

Region of interests

Now that the image is loaded, users can define their ROI that they want to analyze. Click the Select Points button under the ROI panel shown in the red rectangle.

Once clicked the mouse cursor changes to a crosshair for users to pick points. Users can pick as many points as they want as long as the total number of points are greater than 2. These points are connected with a thin blue line on the image axes. Once you are done with the selection, use the right click of your mouse to finalize the process.

Once the ROI is determined, MyoProfiler defines a spline using the selected points as the controls. The generated spline appears as magenta and in this case it resembles a line as there are only 2 control points. Since all the channels are simultaneously analyzed, software copies the ROI and the splines onto the other channels. Users can define the ROIs in any channel they prefer.

Here is the ROI and spline in the Channel 1 axis.

Here is the ROI and spline in the Channel 2 axis.

The ROI and spline can be also found in the Merged image axis. Once again, this tab is used for qualitative review.

Profile extraction

Once the ROIs are defined, the software initiates the analysis. The pixel intensity is extracted along the crosssection of the generated spline. Since the spline is defined in the 2 dimensional space of the horizontal and vertical image axes. Therefore, the extracted profiles are shown in a 3 dimensional axes. The trace colors follow the Channel Colormap scheme. Channel 1 results are shown in green and the Channel 2 results are shown in orange throughout this tutorial.

The extracted profiles are formed by patterns of troughs and crescents and the software uses a peak finding algorithm to identify them. The software identifies all the troughs in the profile and then judiciously divide them into possible Z and M line points with respect to their widths. Then, the software performs a t-test to see whether the groups have different widths. If the test is successful, sarcomere profiles can be extracted. Z lines locations are marked on respective x and y axes. The Z-lines are shown with dashed lines on the horizontal and vertical image axes figures.

Sarcomere analysis

One sarcomere spans between 2 adjacent Z disks. It can be seen from the Profiler Panel that there are two full sarcomeres in the extracted profile. The sarcomere length is defined as the distance between 2 Z-lines in 2 dimensional space.

The software measures the sarcomere length through an arc length calculation. Once the sarcomere length is defined, each sarcomere profile is extracted and the M-lines (magenta diamonds) are identified. Then, the sarcomere profiles are mapped onto a new spatial axis centered on the M-line with respect to their sarcomere length. The A band peaks are marked with red diamonds. The intensity value at each peak is used to find the location at which the intensity is at 50% of the peak value. This is done for both right and left hand side of the M-line. The full width half maximum (FWHM) is calculated as the distance between the left and right side 50% values. All the calculated sarcomere lengths and FWHM values are tabulated in the Sarcomere Panel. The average sarcomere profiles are shown in the Average Sarcomere Intensity plot. Each average profile is shown with the standard deviation envelope.

Profiler options

Channel colormap

As mentioned earlier, the color scheme follows the emission wavelengths of interest. If the color scheme is not feasible for visualization, users can switch to the Parula colormap. This option is available under the Channel Colormap drop down shown in the red rectangle.

Parula colormap is applied all the traces as well as the pseudo coloring of the merged image.

ROI update

Each ROI is attached to a callback at the backend to update the calculations once the position is changed. This feature is executed once users are done with repositioning. The blue circles are the anchors of the ROI and can be moved anywhere around the image. Any available anchor point can be used for this purpose. The ROIs and generated splines are updated in all the channels. Users can perform this from any image axis they want. The below video provides and example.

ROI width (expansion)

In the above video, the images are analyzed using along a spline. This method is similar to “line scan” in imaging. Users can expand their chosen ROI and transform it to a quadrilateral. As it can be seen from the images in this tutorial, the extracted profiles are somewhat noisy. In a quadrilateral ROI case, the profiles along the width of the ROI is averaged. The software rotates the image in the backend until the line segment is perpendicular to the vertical image axis, then records the intensity values centered around the spline. This option can be accessed in the ROI panel. Please keep in mind that the width values are entered in pixels. Once the value is changed, a progress window appears to relay progress for each images. This process is robust but please keep in mind that the length of this process depends on the length of the ROI. The below video shows an example in real time on a standard laptop.

As a result of the averaging, sarcomere profiles are smoothened and the initial values are slightly changed. The analyzed area appears as a transparent patch, shown in red rectangle, on the all the axes following the calculations.

Primary prominence

Since the extracted profiles are cyclic, it is important to distinguish the Z lines from the M lines. The software first identifies all the dips in the profile. Since some of these dips are relatively smaller, the software uses an option called prominence to go around this. The default prominence value is 0.05, which means that MyoProfiler is looking for values that are at least 5% smaller of the sampled space measured from a reference point. Users can change this value through the Profiler panel, shown in red rectangle.

The primary prominence value is increased to 0.5 and none of the troughs are labeled as Z lines. Instead, the software returned a warning message saying that the none of the profiles have a distinguishable pattern. This option might be useful for troubleshooting purposes, but it resulted in with an unwanted results in this case.

Secondary prominence

Similar to the primary prominence, secondary prominence helps to distinguish the peaks in the sarcomere profiles. This value is relatively higher than the primary prominence since they are close in intensity and slightly separated by M-lines with similar intensities.

In this example, increase in the value resulted in missing an A-band peak which is quite close to the M-line pixel intensity.

Pixel to length scale calibration

The reported sarcomere length and FWHM are in microns. Users can change the calibration constant from the Calibration field in the Profile Panel, shown in red rectangle.

Please not that the only sarcomere length and FWHM values are changed.

Export results

Once users are done with their analysis, they can export the analysis into their host computers. A summary Excel sheet and an analysis file with .myoprof extension are saved. The exported analysis files can be loaded back into the software to revisit the analysis.

Export Analysis option is accessed through the top menu, shown in the red rectangle.

Once it is clicked, it opens a file dialog. Navigate to the folder that you want to use.

The Excel file and the .myoprof file share the same name. Enter the name and hit Save. You can find both of the files under the designated folder.

The exported Excel file has multiple sheets. The first sheet is the Analysis Summary sheet. It consists of the name of the analyzed image, used calibration, and number of channels. The summary stats for sarcomere length and FWHM are given with mean, standard deviation, and standard error of the mean.

The second sheet is the Sarcomere Summary. It is the exported version of the sarcomere table in the Sarcomere Panel. You can find the individual values from each sarcomere here.

The following sheets are available for all the analyzed channels. Summary profiles sheet stores the mean, standard deviation, and standard error of the mean sarcomere intensity profiles.

The Sarcomere Profiles has the individual intensity profiles of each sarcomere.